首页/文章/ 详情

屈曲分析是怎么回事-非线性屈曲分析(一)

2月前浏览1751

非线性屈曲分析

   

   



      前面我们已经简要说明了屈曲的概念以及欧拉杆公式的推导,并通过仿真校核了线性屈曲理论计算的结果↑。今天我们来说一下非线性屈曲。

01    
非线性屈曲简介    

 

      前面我们采用求特征值的方法求解结构线性屈曲的临界载荷,但是该方法只适用于弹性屈曲问题。如果在变形过程中,结构材料超过了材料的弹性范围,到达塑性阶段,求解结构屈曲的临界载荷则需要使用非线性求解方法。同时在线性屈曲计算时,我们只能看到结构发生屈曲时的临界载荷以及发生屈曲时结构形状,也就是模态振型,从而看不到屈曲变形的详细历程。

图1.屈曲历程:载荷位移曲线

      图1中的A点为结构发生屈曲的临界点,在发生屈曲后,结构跳到另一个平衡位置C,在这两点之间,结构的承载能力急剧下降,并且力和位移曲线开始偏离线性。从图1可以看到,非线性屈曲分析包含了屈曲发生前后的完整过程。其中前屈曲部分主要是线性屈曲,而后屈曲分析主要是计算结构发生屈曲后的变形和新的承载能力。

02    
非线性屈曲的计算方法    

 

      在有限元分析计算时,不能采用逐步增加力的方法。这是因为当发生屈曲行为时,载荷无法确定对应的位移,也就是载荷与位移无法一一对应。在这种情况下,若使用位移控制,就不会出现这样的问题,即位移与载荷是一一对应的。具体情况见图2。

图2

      对于瞬间翻转型屈曲问题,其位移和载荷曲线如图3所示,此时每个力可能对应多个位移,每个位移也可能对应多个载荷,无论是位移还是力加载都无法分析收敛,这时候只能采用弧长法进行分析。

图3.瞬间翻转型屈曲

        下面我们通过abaqus来模拟几个非线性屈曲的实例以作说明。

03    
压杆结构后屈曲分析
   

 

      在非线性屈曲分析中,首先需要进行线性屈曲分析,根据线性屈曲分析的模态输出文件并结合缺陷的引入进行后屈曲非线性分析。本次选用的结构还是上篇文章的结构尺寸:截面为10mmX10mm,长度为400mm,材料选为Q235材料,其密度、弹性、塑性参数见图4:

A1线性屈曲分析-材料设置

图4.材料属性

A2线性屈曲分析-约束设置:在杆件两端建立两个参考点,在两个参考及附近的面间分别建立耦合约束:

图5.耦合约束

A3线性屈曲分析-载荷步设置:创建载荷步,并在Procedure type中选择Linear perturbation-Buckle,模数设置3即可,最大迭代次数设为300,其他默认即可。

图6.载荷步设置

A4线性屈曲分析-边界条件:在右侧参考点处加Z向-1N的压力,并在左侧参考点处建立固定约束,具体见图7:


图7.边界设置

A5线性屈曲分析-网格划分:对结构进行网格划分时,网格大小设置为1mm,并且为六面体网格:

图8.网格划分

A6线性屈曲分析-输出文件设置:右键所建立的模型,然后选择Edit-Keywords设置关键字,输出文件,该文件包含线性屈曲模态的各节点位移,在非线性分析时,以此作为材料缺陷进行后屈曲分析。


图9.模态节点位移输出文件设置

      提交计算后,在下图中我们可以看到最小的临界载荷为2667.8N(与上篇文章稍有不同是因为两个分析其边界条件稍有不同)。

图10.一阶模态及载荷因子

      首先我们来说下为什么要引入缺陷:如果一个结构是完全对称的,包括载荷、边界、材料等,则在数值模拟分析上不会发生非线性屈曲的行为,无法获得屈曲分析中比较合理的载荷-位移曲线。因此,为了解决这个问题,必须引入小的几何缺陷来触发屈曲变化行为。

      在前面已经说明了由于压杆稳定中的某些情况其位移-载荷是一一对应的,因此可以采用静力通用算法进行非线性分析,或者采用弧长法分析。因此本文分别应用这两种方法分析计算。

      在计算前,可以对前面线性分析模型进行复 制,并且重命名,这样就可以减小一定的工作量:

图11.创建新的非线性屈曲分析模型

B1非线性屈曲分析-载荷步设置:首先建立一个新的静力通用载荷步,并开启大变形,在步长设置方面,初始步长设置为0.01,最大步长设置为0.05(设置小一些,便于计算的收敛)。

图12.非线性屈曲分析载荷步设置-静力通用

B2非线性屈曲分析-边界设置:在结构右侧的参考点处,去掉原有的载荷,并施加位移约束:Z向压缩100mm:

图13.非线性屈曲分析边界设置-静力通用

B3非线性屈曲分析-缺陷因子设置:同样的方法,首先点击Discard All Edit去掉原有添加的关键字,然后在step后面添加以下命令:

*imperfection,file =liner_buckling_L,step=1

1,0.4

2,0.4

3,0.4

其中*imperfection,file =liner_buckling_L,表示读取线性屈曲生成的文件,而下面的1,0.4中的1表示第一阶模态,而0.4为缺陷因子,该因子的计算可以为压杆长度/1000/模态位移值,这个值也可以根据计算结果反复调整。

图14.非线性屈曲分析-缺陷因子设置

     提交计算后,我们可以看到压杆的变形情况如下:

图15.非线性屈曲分析-静力通用变形结果

      然后提取右侧约束点的位移及载荷,其变化关系如下:可以看出在屈曲时承载能力急剧减小,作图方法可参考:弄明白塑性材料拉伸试验仿真过程

图15.非线性屈曲分析-静力通用位移载荷关系

最后我们再采用弧长法进行计算:首先复 制静力通用的计算模型,然后将原有的静力通用载荷步,替换为静力-Riks计算,然后可以控制分析过程结构右侧最大位移为-100mm,DOF为3表示控制Z方向的位移为-100mm。并设置合理的步长:

图16.非线性屈曲分析-弧长法载荷步设置

     然后删除结构右端100mm的位移约束,并添加力载荷,载荷设置为最小的临界载荷:

图17.非线性屈曲分析-弧长法边界设置

     其他设置不变,最后检查关键词,为了减少出错最好去掉原有的关键词,添加缺陷因子关键词,然后提交计算,计算结果如下:

      从上面可以看出,对于压杆稳定的部分情况,静力通用和弧长法都能进行非线性屈曲计算。

   

END





来源:ADAMS及ANSYS等机械仿真
Abaqus非线性通用理论材料控制试验
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-08-25
最近编辑:2月前
奋斗的老张
硕士 签名征集中
获赞 27粉丝 61文章 55课程 1
点赞
收藏
作者推荐

壳单元与实体单元的仿真对比

平板结构仿真分析 01 壳单元和实体单元介绍 壳单元常用于薄壁件,如容器、板材等,每个节点有6个自由度,可以沿XYZ的移动和转动。壳单元的使用场景是结构在厚度方向上的尺寸远小于其他两个维度,因此可以使用壳单元来简化模型,提高计算效率。壳单元的数学原理是基于薄壁结构的理论,其中厚度方向的变形可以忽略。壳单元适用于考虑板、壳的弯曲、扭曲等变形行为。 对于实体单元,几乎任何一个部件都可以用实体单元来建模,但是有时实体单元数量难以控制,单元质量控制比较困难,四面体中容易出现小夹角单元,六面体单元容易出现纵横比过大的单元,而且实体单元应用不当还会出现沙漏、剪切自锁等问题。02 理论分析 本次仿真分析的平板结构尺寸为长600mm、宽200mm、高15mm,一段固定,一段加载垂直平板的10000N的载荷模拟弯曲,为了避免边界影响,基于梁截面纯弯时的表面应力计算方法,该平板上表面(受拉)中间跨度位置的应力,经过计算为400Mpa,具体计算公式如下:03 仿真计算 为了验证壳单元和实体单元的差别,分别划分了如下的网格,网格尺寸都是5mm,壳单元的厚度为15mm,实体单元设置为3层,实体单元类型设置为C3D8。图1.网格划分a.材料设置 由于工况简单,材料仅设置弹性属性即可(合金钢):图2.材料属性 然后分别对两部分结构设置壳和实体属性:图3.材料截面属性赋予b.载荷步设置 载荷步方面,设置静力通用,其他的默认,不开启大变形:图4.载荷步设置c.边界设置 针对壳单元,左侧节点约束6个自由度,对于实体单元,左侧节点约束3个自由度,具体如下:图5.固定约束设置 然后在右端加载10000N的载荷,根据节点数目均布:图6.载荷施加d.结果分析 经过分析,我们看下最大主应力和Mises应力,在图中可以看到平板上表面主应力和Mises应力明显不同:图7.最大主应力图8.Mises应力 然后提取平板结构(壳单元)上表面的最大主应力,仿真结果为398Mpa,与前面理论计算结果相近,由于实体单元采用的是一阶单元,并且模拟纯弯曲,所以应用该单元可能会发生剪切锁死,导致结果出现问题。图9.提取中间表面最大主应力 然后同样应用C3D8单元计算拉伸,我们可以看到应力结果是一致的:图10.拉伸工况壳单元和实体单元结果对比 然后我们再将实体单元换成高阶单元:C3D20R,采用与上面相同的弯曲载荷(右端Z向加10000N):图11.实体高阶单元 计算完成后,首先观察整体变形情况,壳单元右端变形为59.224mm,实体单元右端变形59.253mm,基本一致:图12.整体变形 然后我们再看Mises应力,从下图中可以看到Mises应力也是一致的(除边界区域稍有不同):图13.高阶实体单元与壳单元Mises应力结果 在看主应力时,针对壳单元,需要设置顶层或底层,abaqus默认底层:图14.底层主应力结果图15.顶层主应力结果 从图14和15可以看出,高阶实体单元与壳单元在纯弯状态下,主应力也是一致的。说明高阶实体单元相对于一阶实体单元能更好的模拟薄板结构的弯曲变形。 来源:ADAMS及ANSYS等机械仿真

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈