声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

声振论坛 展示 振动理论 非线性振动 查看内容

这是一个非线性的世界(Nonlinear World)

2016-8-17 14:04| 发布者: tangwx| 查看: 1790| 评论: 0|原作者: 不详|来自: 有限元分析理论与方法微信公众账号

摘要:   线形(linear)和非线性(nonlinear)分析是有限元分析中两大“双塔”,线形问题总是我们所擅长或者乐意研究的,因为其够简单,分析成本低。但本质上说,这个世界主要还是非线性居多,线形总隐含着小变形、小应变、 ...
5.webp.jpg
  线形(linear)和非线性(nonlinear)分析是有限元分析中两大“双塔”,线形问题总是我们所擅长或者乐意研究的,因为其够简单,分析成本低。但本质上说,这个世界主要还是非线性居多,线形总隐含着小变形、小应变、小的旋转、小的温度变化等等前提条件,但即使是“小”,也不会是数学上严格意义的真正线性,不会是一条完美的直线。不过你也不必惊慌,虽然不是完美的线形,但我们并不需要那么完美同样可以在“小”的前提下得到较为精确的结果,所以很多工程师都乐意把问题简化为线形分析。
  一旦我们进入非线性的世界,那么我们将踏上“未知”的旅程。下图展现了一个典型的非线性的历程。横坐标为位移,纵坐标为载荷,在非线性计算中,载荷的施加也是离散叠加的,如图中所示我们先施加10%的载荷(实际情况中没这么大,往往很小如0.1%),第一个线形逼近(有限元分析是数值算法,是多项式的线形逼近过程)没有收敛(没达到指定的收敛阈值),程序要继续寻找收敛点,经过几次逼近后才能找到收敛点,这才完成第一次的载荷收敛,下一步再增加10%的载荷,这个逼近过程也是类似的,直到完成100%的载荷加载。

  这个过程是非常耗费成本的,因为每一步都需要对刚度矩阵(stiffness matrix)进行迭代,而一般线形分析只需要一个刚度矩阵即可,也就是说,如这个例子迭代十次的工作量(理想化的时间,成本)相当于线形分析的10倍,实际情况往往更多。

6.webp.jpg

  有些时候,我们不得不接受非线性,这一节我们主要讨论常见的几种非线性:
  几何非线性(geometic nonlinaerity)
  材料非线性(material nonlinearity)
  接触非线性(contact nonlinearity)

  1.几何非线性(geometic nonlinaerity)
  几何非线性有时候又被叫做大变形分析(large displacement,但有些人并不同意这样叫,因为多大才算大变形呢并没有一个定量的约定),因为几何非线性往往涉及到大的位移,大的应变等等。

  下面我们通过一个理论实例来说明几何非线性。如下图所示,一个杆件一端安装在铰链上,铰链上安装有扭簧(torsional spring)抵抗变形,杆件另一端施加有向上的一定大小的力F。我们来看看杆件旋转角度Θ和力F的关系。

7.webp.jpg

8.png
  此时F和Θ呈线形关系。

  如果我们一开始就把问题简化为线形,那么如下图所示,线形结果和正确的非线性就有非常大的偏差,从图中可以看出,当旋转角度大于30°时,两种分析结果就分道扬镳了。

9.webp.jpg

  几何非线性在工程中非常常见。比如屈曲过程(buckling)就伴随着几何非线性,如下图所示。

10.webp.jpg

  又比如一根梁两端约束的不同,受力情况也不同。如下图所示,一种情况是简支梁一端约束,另一端可以轴向活动;另一种情况是两段都受到约束。当施加垂直方向的载荷时,如果运用线形分析,那么结果没太大差比;当考虑非线性(几何非线性)时,结果将非常不同。

11.webp.jpg

  当考虑几何非线性时,简支梁的位移情况:

12.webp.jpg

  2.材料非线性(material nonlinearity)
  材料非线性应该是我们最熟悉的,在材料力学中主要讲的就是材料的非线性。当材料通过线形区域的屈服点(yield point)时,非线性的塑性(plastic)变形开始接管材料的变形形式,此时应力-应变将不再符合线形关系,

13.webp.jpg

  非线性分析中,步长大小的设定对非线性分析有很大的影响,过大的步长、高度非线性的事件或者接触条件的巨大改变,以及涉及到碰撞、断裂和高度屈曲的情况,求解器就需要长时间寻找到收敛的平衡点,有时候甚至找不到。但为了简化,对材料非线性来说,有些时候,在进入塑形变形后,我们同样可以进行线形简化。如下图所示,塑形阶段应力-应变的变化还是近似符合线形的,我们仍然可以通过线形拟合来进行替代,或者我们可以通过输入来自excel的数据来进行更精确的拟合。

14.webp.jpg

  同样我们可以用上面几何非线性的例子来说明材料非线性,在上面的计算中,我们是假定扭簧在旋转过程中材料属性是不会有变化的,但实际上随着纽簧的旋转会不断硬化,材料属性多少会有些变化(但很小)。为了说明材料非线性,我们假定扭簧系数Kt和角度成函数关系:

  kt=k0+k1Θ
  此时有扭簧的扭矩为:
15.png
  外部弯矩不变,此时可得:
16.png
  可以看出,此时公式既包含材料的非线性,也包含几何非线性。
  同样的,如果旋转角度很小,此时有:
17.png
  此时公式只包含材料非线性
  下图综合绘制了几何+材料非线性、材料非线性、几何非线性和线性4种不同条件下力和角度的关系。

18.webp.jpg

  3.接触非线性(contact nonlinearity)
  接触非线性顾名思义就是当两个对象发生接触时的现象,接触非线性和几何非线性类似,因为也是关于位移(变形)的函数,所以有些介绍中也把其归为几何非线性之中。

  接触非线性的计算公式比较多,不同的接触类型适用于不同的公式,但在有限元分析时,一般都会引入“惩罚因子”(Penality Factor)这个物理量,以后我们会详细讲到。

  下面列举了一些常见的非线性接触的例子:
19.webp.jpg
  卡扣的安装过程

20.webp.jpg
  球状体装入夹持件的过程

21.gif
  轮胎和地面接触的过程

22.gif
  球状体和平面接触

23.webp.jpg
  尖锐物与面接触的过程

  非线性分析(nonlinear analysis)是有限元分析中非常重要的一个大类,我们以后会通过实例重点关注和讲解。




最新评论

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 07:29 , Processed in 0.046251 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

返回顶部