直接积分解法(通用7篇)
直接积分解法 篇1
0 引言
在板材冲压成形中, 圆锥形件是典型的曲面类零件, 在成形过程中既有径 (经) 向受拉、周向受压的拉深成形区域, 一般又存在双向拉应力的胀形区域, 在凹模圆角和凸模圆角区域还会产生弯曲变形。变形过程中, 可能出现的两种失稳形式是破裂和起皱, 其中, 起皱失稳可能出现在凸缘区和悬空的锥壁区。圆筒形件可以认为是当圆锥形件的大端尺寸接近小端尺寸时的特例, 而其他轴对称曲面类零件又可以看作由一系列的微锥组成。此外, 圆锥形零件在板材冲压成形中也有广泛应用, 如汽车灯罩、消音器后盖及轮毂防尘罩、电振动喇叭筒和圆锥滚子轴承保持架等都是圆锥形零件[1]。因此, 对圆锥形件成形的研究是有理论意义和实际意义的。
许多学者对包括圆锥形件的轴对称成形问题进行了研究, 涉及拉深[2,3,4]、胀形[5]和内孔翻边[6]等问题。在分析轴对称成形问题时, 主要采用两种假设模型[2]:一是板坯厚度不变的平面应变假设, 二是板平面 (或曲面) 法线方向应力为零的平面应力假设。当假设板坯厚度不变时, 可使分析过程大大简化, 所以许多学者常用这种假设模型[2,3,4,7,8,9]。
由于板材成形在很多情况下, 垂直于板平面或曲面方向上无力的作用或作用力很小 (即使采用液压成形方式, 板坯法线方向的应力也很小[10]) , 此时板材成形过程更接近平面应力状态, 因此, 按平面应力模型进行分析更接近实际情况。
关于平面内的轴对称成形问题, 在平面应力和比例加载等假设条件下, 李敏华[11]给出了参数方程的直接积分求解方法, 改进了Nadai[12]、Millenson等[13]给出的需要不断迭代才能得到解析解的方法, 显著简化了求解过程。文献[8, 14]将轴对称问题的参数方程分析方法用于求解轴对称成形板平面内冲压成形问题, 给出了不考虑厚度变化时用参数方程形式表示的凸缘区应力显式表达式。但由于轴对称成形问题多数情况下都涉及凸缘区以外的成形区域, 这些区域的形状是锥面或曲面, 无法采用平面应力假设条件下的直接积分解法。这也是该方法没有得到推广应用的原因之一。
本文将求解平面内轴对称问题的分析方法推广至求解圆锥形件的成形问题, 给出了基于参数方程的积分解法。针对圆锥形零件的拉深成形, 将变形区主要划分为凸缘区、凹模圆角区、悬空锥壁区和凸模圆角区。对前三个区域都给出了基于积分解法的应力、应变解, 采用厚0.87mm的ST16板材进行了圆锥形件拉深实验, 测量了应变分布, 理论计算值与实验结果一致。采用本文方法, 还可将积分解法推广至求解一般轴对称成形问题。
1 变形协调方程及平衡方程
1.1 变形协调方程
如图1所示, 从圆盘或圆环的平板毛坯变形为圆锥形壳体, 变形后某一微锥形环内外半径分别为r和r+dr, 若对应半径是r的质点变形后的径向位移为u, 则变形前微圆环的内外缘半径分别为r-u和r-u+dr-du。
α为锥面母线与竖直方向的夹角, 又称为半锥角。由应变的定义, 可以得到经向应变ερ和周向应变εθ[11]:
消去u, 可得
对于平面内的轴对称成形, 可令α=π/2。
1.2 微分平衡方程
在圆锥形件的成形中, 每一变形质点的主轴方向为经向、纬向及法向, 对应的三个方向的应力分别表示为σρ、σθ和σz。截取一微锥壳环, 然后用过轴对称线的平面将其剖开, 得到一半微锥环, 其受力情况如图2所示。
图2中符号含义如下:O O′为对称轴线;r、α的含义与图1相同;σρ、σρ+dσρ分别为上下纬端面上作用的经向应力;σθ为作用于微锥壳体上的纬向应力;t、t+dt分别为上下纬端面处的厚度;ds为微锥壳体的母线长度;p为作用于壳体下表面的单位压力。
分别在轴线方向和水平方向列平衡方程, 有
由式 (4) 、式 (5) , 并注意到ds=dr/sinα, 得
式 (6) 与轴对称平面内成形的微分方程在形式上完全相同, 也可表示成
2 参数方程及边界条件
对于圆锥形件成形的锥壁区, 忽略板厚方向的应力, 根据等效应力的定义, 有
式中, σ为等效应力;R为板厚方向性系数。
参照平面内的轴对称成形问题的解法[11], 应力可以用参数方程的形式表示:
式中, ω为参数。
比例加载条件下, 等效应变ε的参数表达式为
由应力参数方程式 (9) , 得
将式 (9) 、式 (11) 代入式 (7) , 并利用体积不变条件dt/t=-dερ-dεθ, 得
设材料的等效应力应变关系满足σ=Bεn (B为强度系数, n为硬化指数) , 则
将式 (13) 代入式 (12) , 得
由参数方程式 (10) , 得
将式 (10) 、式 (15) 代入式 (3) , 得
在微分方程 (式 (14) 和式 (16) ) 中, 变量ε、r是参变量ω的函数, 当边界条件已知时, 方程的解是确定的。
图3所示为圆锥形件拉深至某一位置时的各变形区, AB为凸缘区, BC为凹模圆角区, CD为锥壁区, DE为凸模圆角区, EF为锥底。
当拉深位置一定时, 若采用直母线假设, 则DC为直线, 锥形件的锥角α已知。设α=α0, α0可由拉深位置、模具几何尺寸以及板厚变化情况确定。
板材平面内轴对称问题可认为是圆锥形件成形的一个特例, 因式 (14) 不含变量α, 对凸缘区的变形, 在求解时令α=π/2, 代入式 (16) 即可。
圆锥形件拉深时, 因凹模圆角对板坯的法向作用力很小, 凹模圆角区板坯的变形可以近似认为只有经向和纬向应力。如图3所示, 设凹模圆角处的板坯形状与凹模圆角轮廓形状一致, 在凹模圆角处, BC与AB、CD分别在B、C点相切。将圆弧BC分成若干个弧段, 当弧段尺寸较小时, 该弧段接近直线段, 其对应的变形区可作为一个小锥环处理。设小弧段的长度为ΔS, 该处切线与竖直方向夹角为αi, 令α=αi, 应变微分方程和求解过程与锥壁区相同。
对于圆锥形件的拉深成形问题, 采用积分解法的计算次序是凸缘区、凹模圆角区和悬空侧壁区。在凸缘变形区, 设F为压边力, μ为摩擦因数, tw为凸缘外缘对应的板坯厚度, ω0为参数初始值, ε0为应变初始值。则边界条件如下:当ρ=Rw (Rw为凸缘外缘半径) 时
当不考虑摩擦时, σe初值为0, ω0=3π/2-β。
在凹模圆角区和悬空侧壁区, 只要给出初始点的应变作为边界条件, 由式 (14) 和式 (16) , 应变与坐标位置的关系是确定的。凸缘区与凹模圆弧交接的位置对应的应变可作为凹模圆角区的初始值, 凹模圆角区与悬空侧壁区交接位置对应的应变可作为悬空侧壁区的初值。
3 直接积分解法
将式 (14) 、式 (16) 联立, 消去dr/r, 设方程为
对于给定的边界条件, 一般需要采用迭代解法才能得到微分方程式 (17) 的解。
对于等锥角的圆锥形件成形问题, 已知边界条件, 即当ω=ω0时, ε=ε0, 采用直接积分解法求任意ω对应的应变ε的计算过程和收敛性可说明如下:
(1) 将区间 (ω0, ω) N等分, 各值为ωi (i=0, 1, …, N) , 对应的应变值为εi (i=0, 1, …, N) , 每一区间的增量Δω= (ω-ω0) /N。
(2) 根据微分学的知识, 有
式中, o (Δω) 为比Δω高级的微分项。
因此, 当参变量由ω0增加为ω1时, 对应的应变ε1可以近似表示为
(3) 一般地, 当已知ωi-1 (i≥1) 对应的应变为εi-1时, 其临近点ωi对应的应变εi近似表示为
(4) 对于任意ω, 对应的应变值近似为
在步骤 (2) ~ (4) 中, 当ω和ε已知时, 由式 (17) 可得到该点的一阶导数。因此, 式 (21) 又可写为
根据积分的定义, 有
N值越大, 计算精度越高。步骤 (1) ~ (4) 的求解过程实际上是用数值方法计算积分问题, 即
在计算过程中, 每次求出ε后, 代入式 (14) 或式 (16) 可求出dr/r, 积分后可得到r。这样就可得到等效应变ε与位置r之间的关系, 进而可求出应力、应变。
4 圆锥形件拉深成形应力应变分布
4.1 圆角区和锥壁区α值的确定
4.1.1 凹模圆角区α值的确定
拉深过程中, 板坯经过凹模圆角时的变形非常复杂, 为了简单起见, 设板坯经过凹模圆角时, 板坯内层保持凹模圆角的外轮廓形状。如图4所示, 根据几何关系可建立质点位置半径r和α的关系:
4.1.2 锥壁区α值的确定
分析表明, 悬空侧壁区锥角对应变分布的影响很大。如图3所示, 设从凸缘区边缘A至凸模圆角与锥底连接点E为变形区。
在计算出凸缘区、凹模圆角区的应变分布之后, 可按体积不变条件确定锥角。实际计算时, 暂以变形前后总面积不变, 初步确定锥角α0, 再根据计算后的厚度变化按体积不变进行修正。计算表明, 采用面积不变假设与体积不变假设结果相差不大, 一般需修正2~3次就能使体积计算值收敛。
采用积分法求解, 实际是将板坯分成若干微段, 成形后的体积很容易得到, 体积可按下式进行计算:
式中, Vi为各部分的体积;rai、rbi分别为各计算区间的起止点对应的径向坐标;ri、ti分别为位置半径和相应位置的板坯厚度;αi为对应计算段的板坯所在弧或平面与轴对称线夹角。
4.2 应力、应变分布
取凹模内半径Rd=26.82mm, 凸模外半径Rp=20mm, 凹模圆角rd=13mm, 凸模圆角rp=5 mm, 原始毛坯半径R0=55 mm, 板厚0.87mm。压边力取5kN, 板坯与凹模、板坯与压边圈之间的摩擦因数都取0.06, 不考虑板坯在凹模圆角处的摩擦。材料选用ST16板材, 性能参数经拉伸实验确定:B=511.4 MPa, n=0.26, R=2.243。拉深相对位置Rw/R0分别为0.919、0.878和0.838, 采用上述计算方法得到的应变分布和应力分布分别如图5~图10所示。
由前面的理论分析可知, 应力值与材料强度系数B成比例关系, 为了简明起见, 图8~图10采用了应力与强度系数的比值。
图5~图7显示, 经向应变为正值, 除接近凸模圆角附近的很小区域外, 其他区域的纬/周向应变都为负值。从凸缘外缘到凸模圆角, 板坯厚度应变由正变为负。随着拉深过程的进行, 变形程度逐渐增大, 经向应变、周向应变绝对值最大值也逐渐增大。
图8~图10表明, 经向应力都为拉应力, 周向应力分为拉和压两个区域。经向应力最大值随拉深变形程度的增大而逐渐增大。周向应力在大部分区域都为负值 (压应力) , 但在接近凸模圆角的区域为正值 (拉应力) , 且随着变形程度的增大, 拉应力值也增大。接近凸缘区变形性质为经向受拉、周向受压的拉深, 接近凸模圆角区变形为双向受拉的胀形, 且随着变形程度的增大, 双向应力越来越接近。
5 圆锥形件拉深成形实验
圆锥形件拉深成形实验选取ST16板坯, 板坯几何尺寸、性能参数、拉深工艺参数和模具几何参数与前面理论计算相同。实验过程中先在圆形毛坯上印制网格, 再进行拉深成形实验, 测量成形后板坯的应变分布。
拉深成形时, 将板坯印制有网格的面放置到凹模一侧, 变形后这个表面即为图4所示的内层, 其与模具轮廓形状最接近。为了减小摩擦的影响, 选择了薄膜润滑。如图11所示, 三个试件拉深高度H从左至右依次为22.87mm、27.26mm和32.53mm, 对应的法兰外缘相对半径Rw/R0分别为0.919、0.878和0.838。
成形后的试件用自动应变测量系统Vialux AutoGrid进行测量。为了真实地反映应变的分布情况, 将拉深试件的整个凸缘区、凹模圆角区和悬空侧壁区所有网格节点都作为测量点, 将经向应变、周向应变测量值与相应各点所在位置半径之间的关系表示在图12~图14上。图中光滑曲线为理论值, 离散点为实验值。实验得到的经向应变、周向应变应变值都分别集中在一个条带区域内。
实验和计算结果表明, 采用直接积分法得到的理论计算值与实测的周向应变值一致, 经向应变值的理论计算结果与实测值在凹模圆角区、接近凹模圆角的凸缘, 以及接近凸模圆角的悬空侧壁区相差稍大, 但变化趋势基本一致。由于理论计算采用了薄膜理论的平衡方程式, 没有考虑弯曲的影响, 另外采用直母线假设等, 这些都可能造成了理论计算与实际存在一定的偏差。板坯在凹模圆角的弯曲变形主要是沿经向的, 因而对周向应变影响不大。
6 结论
(1) 对圆锥形零件的成形问题, 给出了基于参数解法的直接积分解法和具体求解过程, 这一方法避免了进行复杂的迭代求解。
(2) 选取0.87 mm厚的ST16板材和一定的模具尺寸及工艺条件, 对圆锥形件拉深成形, 求解得到了板坯拉深至不同位置时的应力应变分布。
(3) 进行了圆锥形件拉深成形实验, 实验测量了凸缘区、凹模圆角区和悬空侧壁区的应变分布。结果表明, 在相同变形条件下采用直接积分法计算的周向应变与实测结果较接近, 而经向应变在某些区域差别稍大, 但趋势基本一致。
浅析不定积分的解法 篇2
1 直接积分法
直接积分法是求不定积分的基本方法, 是基本途径, 也是其他积分方法的基础, 这一方法是直接利用积分法则和公式得出结果, 或将被积函数做恒等变形, 使之符合基本法与公式, 然后再利用积分法则与公式作出结果。例如求把该式分子乘得到分项后得到然后, 利用基本公式求得结果为
注:在分项积分后, 每个不定积分的结果都含有任意常数.由于任意常数的代数和仍为任意常数, 故只需在最后一个积分符号消失的同时, 加上一个积分常数就可以了.
2 第一换元积分法 (凑微分法)
求一个函数的不定积分是积分学的一个基本问题, 解决这类问题的方法多种多样, 其中有一种方法就是第一换元法, 换元法是求不定积分的基本方法。例如求∫x2exdx分析:可以用两种方法凑微分, 但用哪一种行得通?要试试看。
∫x2exdx=∫x2d (ex) =x2ex-∫ex d (x2) =x2 ex-2∫xexdx虽然还不能得到结果, 但次数降低了, 越变越简单。再进行一次分部积分得到:
有些积分, 用一次分部积分不行的话, 可进行两次、三次或更多次的分部积分法, 只要能给它凑成基本公式的形式, 分部积分法就可用。
一般地, 形如∫xk ex dx的积分, 都可以用分部积分法计算, 并且计算方法和例题类似。
计算熟练以后, 就可以省略“设”的步骤, 把所设的式子当作一个整体, 在心里面想着它是一个变数, 就可以使书写简化。
注:给分式的分子加一项再减一项分成两项, 或将被积函数乘一项再除一项是常用的求积分技巧。
本例的 (2) 有多种解法, 各解法所得结果有所不同 (请读者试着用其他方法求此不定积分) , 它们可以通过三角函数间的关系化为相同。但多数情况下, 把各解法所得结果化为相同是很复杂的, 此时, 检验积分结果是否正确, 只需对结果求导, 如果导数等于被积函数, 则结果正确, 否则结果错误。
3 第二换元积分法
第一换元积分法是选择新变量u, 令u={ (x) 进行换元。但对于某些被积函数来说, 用第一换元积分法很困难, 用第二换元积分法求不定积分最常见的是被积函数中含有下列表达式的不定积分:例如, 而用相反的方法: (1) 令x=asint, (2) 令x=atant, (3) 令x=asect进行换元, 就能比较顺利地求出结果。这是一种与凑微分法的过程刚好相反的计算不定积分的方法。
注:第二换元积分法的换元表达式中, 新变量t处于自变量的地位, 而在第一换元积分法的换元表达式中, 新变量则处于因变量的地位。此外, 在使用第二换元积分法时, 为保证x=ψ (t) 的反函数确实存在及原来的积分有意义, 通常要求x=ψ (t) 是单调函数、有连续导数且ψ' (t) ≠0。
4 分部积分法
微积分中对于那些由两个不同函数组成的被积函数, 不便于进行换元的组合分成两部分进行积分, 其原理是复合函数求导的逆用。根据组成积分函数的基本函数将积分顺序整理为口诀:“反对幂指三”。分别带指五类基本函数:反三角函数、对数函数、幂函数、指数函数、三角函数的积分次序。具体操作如:被积函数由幂函数和三角正弦组成则按口诀先积三角函数。分部积分也是求不定积分的基本方法:
例∫x 2e x dx分析:可以用两种方法凑微分, 但用哪一种行得通?要试试看。
虽然还不能得到结果, 但次数降低了, 越变越简单。再进行一次分部积分, 应该行得通。
解:
有些积分, 用一次分部积分不行的话, 可进行两次、三次或更多次的分部积分。直到能用基本公式求出或是能转化成所求式子即可。
不过, 在进行这种涉及繁复的代数计算时, 一定要注意掌握一个原则, 就是动手之前仔细观察, 根据经验判断是否存在更为简单的方法, 只有在确实找不到简单方法之后, 再开始根据这种确定的计算程式来进行计算。
一道不定积分的多种解法 篇3
解法2在解法1的基础上,在第二步时可采用整体凑微分的方法.
解法3解法1与2均是对三角函数采用降幂方法先进行化简,再凑微分,然后利用分部积分法求得结果,那也可以不化简而直接凑微分,但也只能将一个cosx凑进去,然后利用分部积分法求解,解法如下:
解法4同上,依然是优先考虑三角函数,先变形再求解.
解法5以上四个解法,都是从三角函数出发,那也可逆向思维,打破常规,考虑幂函数凑进去作为v,再采用分部积分法进行求解.
注:解法5中,虽然在第一步处理时是优先选择了将幂函数凑进去作为v,但是注意后面处理相同情况时,还是优先选择了三角函数,否则是解不出来结果的.
以上,我们给出了一道不定积分的五种计算方法,当然这方法都是分部积分法的应用,虽然是同一个应用,但是我们处理方式不同,方法过程就有所不同,今后,在做不定积分题目时,在掌握基本方法和原理的基础上,我们要学会发散思维,灵活处理.
摘要:不定积分是高等数学的重要内容之一,它的求解方法有很多,熟练掌握不定积分的求解方法是后期学习的坚实基础,本文主要是基于分部积分法就一道不定积分的求解方法就行探讨.
一个不定积分问题的多种解法 篇4
求函数不定积分有公式法、代数变形方法、换元法、分部积分等方法, 灵活性强, 知识的综合应用较多, 现给出一道不定积分题的多种解法, 以体现不定积分各种解法的优劣及各种思维方法的应用.
例 求不定积分∫
分析 该被积函数显然是有理真分式, 故首选按有理函数的积分计算:
解法一 设
则
∫
其次考虑换元积分法, 将幂次数降下来解决问题, 及三角换元法.
解法二 令x2=u, 则
∫
∫
其中, u=x2, 则原式
解法三 设x=tant, dx=sec2tdt,
=∫tantdt-∫sintcostdt
本方法虽然过程较复杂但体现出逆向思维, 即三角函数的有理式可化为有理函数积分解决, 而有理式也同样地可化为三角函数的积分问题解决.
解法四
∫
解法五
∫
利用拆项凑微分和先凑微分再拆项能培养学生思维的创新.用第一换元法与分部积分综合应用, 也是既快捷又理想求不定积分的方法.
解法六
∫
解法七
一道微积分题的新解法 篇5
在数学分析中, 极限占有非常重要的地位, 所以学好极限是非常重要的, 本文对钱吉林编著的《数学分析题解精粹》中一道考研题力图以具用高中水平就能进行求解, 然后对其加以推广, 以达到举一反三的目的。它不仅对考研朋友们有一定的帮助, 而且也可以为高三的学生、中学教师以及从事《数学分析》或《高等数学》的高等教育研究者提供参考。原题如下:求求吉林用下面方法1求解, 但没有证明1援1, 我们在讲新内容之前, 不妨回顾以下几个高中数列中的公式:
以上几个公式中第一个可以用等差数列求和公式来直接来求, 其他可以用数学归纳法来证明, 我们这里不作证明, 下面只针对 (4) 用另一种方法加以证明。
1 主要结果
1.1 证明
证 因为 (n+1) 5-n5=5n4+10n3+10n2+5n+1
所以
又因为
所以
14+24+34+44+…+n4
1.2 求 (这里的琢是正整数) n寅肄n
方法1根据2援1的证明可以知道
当α<5时, 。
方法2我们知道闭区间[0, 1]上函数y=x2的定积分为, 这是通过利用定积分的定义通过极限的思想来求的得, 而极限的思想我们在高中立体几何求球的表面积就开始接触, 并不是陌生的东西, 因此受此启发我们用定积分的定义解本题, 事实上,
下面对α分三种情况进行讨论:
当α>5时, 原式=0;
当α=5时, ;
当α<5时, 原式=∞。
从上面两种方法可以看出, 这两种方法的极大地减少了用极限定义带来的繁琐。一般地, 1n+2n+3n+4n+…+nn (这里的n是正整数) 等于多少我们并不知道, 因此方法1有局限性, 这时我们可以用方法2就是对这一类问题进行推广的解法。请看下面例题:
1.3 求 (这里的α、β是正整数)
解
下面对α分三种情况进行讨论:
当α>β+1时, 原式=0;
当α=β+1时, ;
当α<β+1时, 原式=∞。
参考文献
[1]Polya G (李心灿等译) .教学与猜想[M].北京:科学出版社, 1984.
[2]徐利治.数学方法论选讲[M].武汉:华中工学院出版社, 1983:177-182.
[3]裴礼文.数学分析中的典型问题与方法[M].北京:高等教育出版社, 1993.
一道不定积分题的若干种解法 篇6
在学习不定积分时, 很多同学不能很好地理解不定积分的定义, 总是喜欢求出被积函数的一个原函数后不在后面加上任意常数C, 这是对不定积分概念的本质不理解造成的. 不定积分实际上就是其原函数族. 又因为原函数是不唯一的, 从而造成了其结果形式的不统一, 同时也给学生判断不定积分结果的正确性带来一定的困难. 因为即便结果不一样也可能都是正确的, 因此要验证其结果的正确性, 必须要通过对结果求导看其是不是被积函数来判断. 对于计算不定积分来说, 求原函数就是关键, 不同的解法可以得到不同的结果, 然而都是正确的. 下面就通过给出一道不定积分的几种解法, 帮助学生更好地理解不定积分的概念和掌握计算不定积分的方法.
二、有关的几种解法
1. 第一换元法
2. 第二换元法
3. 无理式代换
参考文献
[1]同济大学应用数学系.高等数学:上册[M].5版.北京:高等教育出版社, 2002.
例谈求广义积分的三种不常见解法 篇7
关键词:留数定理,Fourier积分,Fourier变换,Frullani公式
在科学研究中, 数学是重要的基本工具, 已经广泛地渗透到各个学科, 是培养各类科技人才的必要基础。《高等数学》这门课在内容上极其丰富, 方法上灵活多样。在多年的教学实践中, 作者发现, 许多学生在学习求广义积分时感到非常困难, 分析其原因在于广义积分的概念不仅抽象, 而且计算方法灵活, 不易掌握。
因此, 熟练准确地计算各种广义积分, 对于进一步学好《高等数学》的其它相关内容是十分必要的。本文提出了三种不常见的求广义积分的解题方法, 目的是开阔学生的解题思路, 从而提高解题能力。
一、利用复变函数的留数定理计算广义积分
二、利用Fourier积分和Fourier变换求广义积分值
三、利用Frullani公式计算广义积分
参考文献
[1]丁家泰.微积分解题方法[M].北京师范大学出版社, 1981.