001 结构极限强度分析的三种有限元解法研究_彭大炜

第25卷第2期2010年4月中国海洋平台CHI NA O FFSH OR E PL AT FO RM V ol.25N o.2A pr.,2010收稿日期:2009-11-03; 修改稿收到日期:2009-12-25作者简介:彭大炜(1984-),男,硕士研究生,主要从事船体结构极限强度分析研究。

文章编号:1001-4500(2010)02-0001-05结构极限强度分析的三种有限元解法研究彭大炜, 张世联(上海交通大学,上海200030)摘 要:介绍了弧长法、阻尼因子法和准静态法等3种计算结构极限强度的非线性有限元解法及其求解思路。

通过Reckling N o.23模型数值计算,对上述3种解法各自的计算过程、求解特点和要点进行了比较和归纳。

为更准确、高效地运用非线性有限元程序计算分析结构极限强度提出了合理的建议。

关键词:弧长法;阻尼因子法;准静态法;极限承载能力中图分类号:U 661.41 文献标识码:AStudy on Three FEA Methods of Structure Ultimate Strength AnalysisPENG Da -w ei, ZH ANG Sh-i lian(Shang hai Jiao T ong Univer sity,Shanghai 200030,China)Abstract:Three alg orithm s,Arc -Length metho d,Dam ping Factor metho d and Quas-iStatic method in no nlinear FEA method of ultim ate strength analy sis are presented in this pa -per.The three algorithms have differ ent analysis processes,features and key points,w hichare discussed and com pared w ith one another through a num eral calculation by Reckling N o23test m odel.The result of the comparisons gives a r easonable sugg estio n for better use of non -linear FEA m ethod of ultimate streng th analysis.Key words:Arc -Leng th method,Dam ping Factor m ethod,Quas-i Static method,ult-imate streng th图1 非线性有限元分析流程0 前言近年来,结构设计的思路已经逐步从传统的许用应力设计向极限状态设计发展[1]。

研究结构的极限强度水平,确定其与设计载荷水平之间的确切裕度,已成为结构理性设计的核心和基础。

目前,船体结构极限承载能力研究的方法主要有[2]:实船事故调查和模型试验法、直接法(解析法)、以Smith 法为代表的逐步崩溃法和非线性有限元法。

随着计算机技术的不断发展,非线性有限元法现已成为计算和评估结构极限承载能力最理想的方法,其分析问题的主要流程如图1所示。

在非线性有限元法中,有3种不同的分析结构极限强度的解法。

它们分别是弧长法、阻尼因子法和准静态法。

本文基于ABAQU S 大型通用非线性有限元程序,以Reckling -N o.23试验模型为对象,分别采用上述3种解法进行极限强度的数值分析,并研究不同解法的求解过程,对不同解法的应用特点进行了归纳和总结,且还针对如何更准确、高效地运用非线性有限元法计算分析结构极限强度的问题提出合理的建议。

1 理论介绍从求解本质而言,弧长法和阻尼因子法是以结构非线性静态平衡方程式(1)的求解(牛顿拉-普森迭代)为基础。

而准静态方法是以结构非线性运动方程式(2)的显式求解(中心差分法)为基础。

{P }-{I }=0(1)[M]{&u }={P }-{I }(2)式中:{P}为载荷列阵;{I }为内力列阵;[M]为质量矩阵;{&u }为加速度列阵。

以下分别简要介绍这3种解法的基本思路。

1.1 弧长法弧长法的基本思路是通过设置一个参数(弧长l )来控制平衡方程的增量迭代和收敛。

将(1)式写成增量格式:[K T ]{$u}={$P}-{R}(3)式中:[K T ]为切线刚度矩阵;{$u}为位移增量;{$P}为载荷增量;{R}为残差力。

在弧长法中,第i 步迭代的载荷增量{$P}i ,由载荷增量因子$K i 和参考载荷{P ref }来控制,即{$P}i =$K i {P ref }(4)将(4)式代入(3)式即得到弧长法的第i 步迭代的增量形格式:[K T ]{$u}i =$K i {P ref }-{R}i(5)图2 弧长法示意图弧长法求解时,是以前一步增量计算得到的平衡点为圆心,以弧长增量$l i 为半径,通过牛顿拉-普森迭代搜索下一个平衡点,如图2所示。

每一步的弧长增量$l i 、载荷增量因子$K i 和位移增量{$u}i 通过下面的约束方程来控制:|$l i |=|{$u i }|2+|$K i {P ref }|2(6)通过迭代,直至残差力在容差{R}i 范围内。

当第i 步迭代完成时,有载荷列阵{P }i =E i {$P}i =E i $K i {P ref }(7)位移列阵{u}i =E i {$u}i(8)弧长l i =E i $li(9)由于在弧长增量$l i 中同时包含了载荷增量{$P i }和位移增量{$u i }的信息,运用弧长法能够全程跟踪结构在/加载0($K i >0)和/卸载0($K i <0)过程中的载荷-位移平衡路径。

但在实际运用时,存在跟踪失败[3]和收敛性问题。

目前弧长法还在不断地研究改进。

1.2 阻尼因子法阻尼因子法仍采用传统的载荷步长增量来进行非线性迭代求解。

其思路是通过在平衡方程式(1)中引入阻尼力项来求解结构不稳定崩溃的问题,求解的平衡方程为{P}-{I }-{F v }=0(10)其中{F v }为阻尼力列阵,取决于求解时结构变形的速度,由广义节点速度{v }来表征:#2#中国海洋平台 第25卷 第2期{F v }=c[M *]{v }(11)式中:c 为阻尼系数[4];M *为人工质量矩阵;广义节点速度{v}={$u}/$t 。

加载的初始阶段,结构尚处于稳定状态,此时广义节点速度{v}很小,故阻尼力项{F v }对平衡方程式(10)几乎没有影响,可以忽略。

随着载荷的不断增加,结构趋向于不稳定,当外载{P}已经不能完全由结构内力{I }来平衡时,结构达到极限状态,相当部分的应变能将释放转化为动能,广义节点速率{v}迅速增大。

此时,阻尼力项{F v }通过做功消耗释放的应变能,在平衡方程式(10)中起到维持求解系统的/平衡0作用。

由于采用传统的载荷步长来进行非线性迭代,阻尼因子法无法继续有效跟踪结构在/卸载0过程中的载荷-位移路径,取而代之的是一条几乎水平(斜率为0)的直线。

通过考察阻尼力项{F v }为维持系统/平衡0所做的功所占结构应变能百分比的历史变化曲线,就能够确定结构的极限承载能力。

相比弧长法、阻尼因子法的数值收敛性要好一些,能够解决的问题也更广泛。

1.3 准静态法与弧长法和阻尼因子法的静态求解方法不同,准静态法从本质上讲是一个结构动态求解的过程。

在求解时,对式(2)采用中心差分法进行显示的时间积分,由一个增量步的动力学条件计算下一个增量步的动力学条件,直至求解时间结束。

准静态法的基本思路就是用慢速加载的动态分析来模拟静态问题,所以,求解关键是要设置合适加载速率。

加载速率过快会导致求解结果的局部性(剧烈的结构局部变形),使计算结果偏离/准静态0的要求;而过慢的加载速率意味着较长的加载时间,从而使计算时间大幅度增加。

故分析时,通常会取从快到慢多个加载速率进行分析比较,以选定一个合适的加载速率。

判断加载速率是否合适的一个重要标准就是分析过程中结构模型的动能与其应变能之比,一般准静态的要求是小于5%。

由于中心差分法是条件稳定的算法[5],在分析时,时间步长$t 必须小于稳定性限制$t stable 才能保证求解的稳定性:$t s table =2/X max \$t e(12)式中:X max 为结构最高阶固有频率;$t e 为结构模型中最小尺寸的单元的稳定时间步长。

$t e 与单元的特征尺度L e 、弹性模量E 和材料密度Q 有关,是稳定性限制$t stable 的上界[6]:$t e =L e /E/Q(13) 在计算中,由于结构最高阶频率X max 不易求得,故时间步长$t 就取为$t e 。

相比弧长法和阻尼因子法,准静态方法最大的优势在于采用中心差分法进行显式时间积分不存在收敛性的问题。

因此,准静态法能够很好地求解更复杂结构崩溃问题,如结构的自接触和材料的失效问题。

在求解极限状态问题时,时间步长$t 往往较小,准静态/缓慢0加载可能导致求解的机时很长,这时可采用质量放大等方法进行调整。

2 数值计算分析2.1 Reckling N o.23模型试验Reckling (1979)[7]采用箱型剖面模拟实船剖面进行了系列总纵极限承载能力试验,以研究船体梁在极限状态下各种崩溃模式和剖面应力分布。

本文选取加强筋较多的Reckling No.23号模型作为分析对象,模图3 Reckling 有限元模型一阶屈曲模态型跨长l =500mm 。

表1列出了模型截面尺寸和材料属性。

表1 R eckling N o.23模型截面尺寸和材料属性加筋板尺寸/mm 屈服应力/(N/mm 2)弹性模量/(N /mm 2)压缩翼板(甲板)600@2.5246210000拉伸翼板(底板)600@2.5246210000舷侧板400@2.5246210000翼板加强筋30@20@2.5246210000舷侧加强筋30@2.5246210000#3#第2期 彭大炜等 结构极限强度分析的三种有限元解法研究在进行极限强度分析时,可将如图3所示模型的第一阶线性屈曲模态的变形形式作为结构的初始缺陷图4 载荷-位移平衡路径曲线(弧长法)引入模型网格。

2.2 弧长法解采用弧长法时,在模型两端设置大小相同、方向相反的参考载荷(中垂弯矩),对应各弧长增量步的载荷与位移的表达式如式(7)~式(9)所示。

得到的端面弯矩-转角曲线如图4所示:考察曲线峰值点对应的弯矩,即为采用弧长法得到的极限载荷。