有限元分析中的应力.pdfVIP

你真的了解有限元分析中的“应力”吗 原创 2016-07-09Feaforall ?????? 虽然在有限元分析中我们常常会用到软件后处理程序得出的应力值 （stress），但其实应力有很多值得我们研究的地方。 ?????? 如果我们把作用于物体的力产生的各处应力汇总起来，那么应力也就像 流体分析CFD 中的速度或者压力一样形成应力场 “流过”物体，我们抓取感兴趣 的地方来进行强度的评估。然而，由于应力状态变化复杂，并不好在3D 单元中 进行可视化，所以我们更需要根据软件已有的功能来探究应力的意义。 ??????1. 几乎所有的有限元分析结果中，默认的应力结果是冯米斯应力（Von Mises ），冯米斯应力是一个标量结果，并没有力的方向性指示。学过材料力学 的应该知道还有一种应力是主应力（principle stress），主应力是矢量，某些情 况下也是非常有用的，那么他们之间有什么区别？ ?????? 2.物理内部的受力在不同部位都不一样，我们怎样尽可能多的去研究内 部力场的不同特性并且通过软件可视化出来呢？ 下面我们将探究上面的两个问题。 什么是应力？ ??????? 首先我们先说说什么是应力。众所周知，应力（stress）是单位面积上 作用的力（forces）。我们并不好感知或者测量应力，但力（force）是实实在在 的，我们可以很好的感知和测量。物质总是由原子构成的，从原子的维度看，原 子之间相吸或者相斥。物体在没有受力的状态下，原子处于自然状态，所有的力 互相平衡，如果物体受到外部力的作用，原子就会偏离平衡位置去寻找新的平衡 位置来平衡外部力。如下图所示，相同长度L 上分别有两排5 对的原子和两排6 对的原子，如果假设原子之间的吸引力相同，那么单位长度上6 对原子的应力要 比5 对的大，扩展到宏观的3D 情形同样适用。 力和应力单元 ?????? ?微积分学科的发展可以使我们通过数学运用无限（无限大或者无限小） 的原理来处理很多实际问题，宏观物体的受力是微观单元的叠加。在材料力学中， 我们把一个无限小的立方体（cube ）单元来描述某一点的受力情况。为什么无 限小呢？因为由于无限小，小到物体内部力是均匀的，没有应力变化，只有一种 应力状态。如下图所示，六个面分别受到法向力平衡。 ????? ?上图是垂直于截面的法向力（normal force）情形，那么自然还有一种 剪切力（shear force）。如下图所示，如果X 方向截面受到剪切力Fxy （下标x 代表x 方向截面，y 代表受力朝向y 方向），那么为了使单元平衡，就会产生其 余三个力Fyx ，Fxy ，Fyx （如果想当然觉得只有Fxy 产生，那么立方体将受到 弯矩无法平衡） ?????? 如果将法向力和剪切力汇总到一个立方体中（为了便于图形呈现，其它 三个面的受力状态这里没表示出来）： ?????? 有限元模型中，每个单元受到的力，包括法向力、剪应力的合力都是和 外力通过节点传递到该单元的力平衡的，这种微观的平衡是力学平衡的微观表现。 有力（force）就有应力（stress），相应的应力单元如下图所示： ?????? 下面我们通过一个实例来研究物体在受力状态下的力的多种观察视角。 ?????? 如下图所示，一个斜十字交叉的简易桁架模型左端两个孔完全约束，右 边两个吊耳孔分别受到向下和左/右方向的力（大小任意）。 ?????? 整体的Von Mises 应力状态如下图所示，一般情况下软件都默认得到这 个应力云图，我们看看最大的受力区域和值就可以了，但今天我们不关注这个， 我们更关注力在不同区域，不同方向的不同形式，von mises 是标量，没有方向， 得出的数值也没有正负，得不到这些细节信息。 首先我们来看看X 方向的受力情况： ?????? 从上图中我们可以看出来，上面部分主要受到拉伸力（数值是正值）， 下面主要受压缩力（数值是负值），为了证明我们的观察，我们将上面受到拉伸 的区域C 和压缩力的区域E 局部放大得到如下的结果。 放大区域C ： ?????? 从上图可以看出在图中虚线方向上，力的变化都是正值，还可以看出这 种变化是线性的。 再看区域E ： ?????? 从上图可以看出在图中虚线方向上，力的变化都是负值，同样是线性变 化的。 ?????? 此外，从X 方向的应力分布云图中还可以看出，区域A 似乎是拉伸最 大点，区域B 似乎是压缩最大点，但这只是X 方向的情形并不能告诉我们全部 的信息。 我们再看看Y 方向的应力分布： ?????? 从Y 方向的应力分布来看，最感兴趣的应该是区域D 。绘制区域D