本帖最后由 wdhd 于 2016-3-16 10:44 编辑
==6==应力应变的规定
在定义材料特性时,在ABAQUS中应力应变的各个分量定义如下顺序:
:1-方向的正应力
:2方向的正应力
:3方向的正应力
:1-2面的剪应力
:1-3方向的剪应力
:2-3方向的剪应力
比如,一个完全各向异性、线弹性矩阵为:
1-、2-和3-方向的确定依靠所选择的单元类型。对于实体单元就是整体坐标系的空间方向。对于壳单元和膜单元,1-、2-方向默认上是壳或膜曲面的局部方向。可以通过*orientation对三个方向进行适合自己分析的修改。
对于实体单元的几何非线性分析,默认的方向并不会随着材料的旋转而变化。然而,如果通过*orientation定义的方向则会旋转。
ABAQUS/Explicit在分析中,在内部以不同的顺序完成对各应力的存储。对于几何非线性分析,无论是否使用*orientation,这些被存储的变量总是随着材料的方向旋转的,这一点在子程序VUMAT被使用时尤显重要。
● 各向异性材料行为
在连续体单元中定义各向异性材料行为时,用*orientation所定义的材料方向是很重要且必要的。
● 零值应力
在分析中,如果应力值一直为零,在矩阵存储时将被忽略。例如,对于平面应力分析,ABAQUS仅存储两个面内的正应力和一个面内的剪应力。
● 剪应变
ABAQUS输出工程剪应变 :
● 应力应变测量
在ABAQUS中,应力使用的是柯西应力(?)或真实应力值,即每单位当前面积上的应大小。详细地请看手册。“Stress measures,” Section 1.5.2 of the ABAQUS Theory Manual。
对于几何非线性分析,存在多种不同的应变计量方法。不像‘真实’应力,它没有很清楚的真实应变的计法。对同一种物理变形,在大应变分析中不同的应变测量方法会给出不同的应变值,当然其值所能反应的实际性也就不同。如何选择最好的应变测量方法,这依靠分析类型、材料行为以及在一定程度上也依靠个人喜好!详见“Strain measures,” Section 1.4.2 of the ABAQUS Theory Manual
默认情况下,在ABAQUS/Standard中应变以“综合”总应变(E)输出,对于大应变的壳、膜和实体单元,还有两种其他的总应变计法可以输出:自然应变(LE)和公称应变(NE)。
在ABAQUS/Explicit中,自然应变(LE)是默认的应变输出,也可以要求公称应变输出。而“综合”应变在ABAQUS/Explicit中是无法得到的。
● 总(综合)应变
默认情况下,ABAQUS/Standard输出到.dat文件和.fil文件的应变为综合应变,这是对于所有可以在材料体下将应变率数值积分获得有限应变的单元都适用的:
其中 和 分别是分析中第n+1次和第n次的总应变, 为旋转张量的增量, 为从第n次增量到第n+1次增量的总应变增量。对于使 用正转坐标系统(右手法则)的单元(使用*orientation有限应变壳单元、膜单元和实体单元)上式可以简化为:
应变增量可以通过对变形率D在整个时间增量上积分得到:
这种应变计法对于弹塑性(粘弹性材料)或弹性-蠕变材料都是适合的,这是因为塑性应变和蠕变应变都是通过相同的积分方式得到的。在这样的材料里,弹性应变是很小的(因为屈服应力相对于弹模来说是较小的),此时总(综合)应变可以直接地与塑性应变和蠕变应变相对比!
如果应变的主方向关于材料方向的变化而旋转变化的,那么最后所得到的应变是不能和总变形相联系的,此时无论采用的是何种坐标系统。如果主应变方向保持固定,那么应变就是变形率的积分:
这与稍后将要讨论的自然应变等效。
● 格林应变
在ABAQUS/Standard中,对于小应变壳和梁单元,默认的应变计法E为格林应变:
此时,F为变形适量梯度而I为特征张量。这种应变计法适合于在小应变、大旋转分析中使用这些单元。分量 代表沿原定义方向的应变,不能在有弹塑性或超弹性材料行为的有限应变分析中使用小应变壳和梁单元,因为可能会导致不正确的分析结果!
●公称应变
公称应变NE被定义为:
其中, 为左拉伸张量, 为主拉伸, 为在当前材料参数下的主拉伸方向。因此,主公称应变也就是在主材料参数方向下长度变化对原长度的比值,从而直观地解释了变形。
●自然应变
自然应变LE被定义为:
(式中变量同上)。
● 应力不变量
在ABAQUS中许多结构模型都是根据应力不变量进行阐述的,这些应力不变量为:
等效压应力:
Mises等效应力:
第三偏应力不变量:
这中间,S为偏应力,被定义成:
==七==有限旋转
空间中,对于有限旋转作如下规定:定义 、 、 为关于整体X、Y和Z轴的旋转量(也就是指在一个节点的4、5、6自由度)。
定义 ,其中 ,方向P为旋转同方向, 为右手法则得出的旋转角度(弧度表示),见下图。 值不是唯一确定的。在大旋转问题中,如果总旋转量超过 ,可以加或减任意倍的 (以得到在 内的值),而这可能会导致旋转分量的不连续输出。
这个规则提供了在大多数分析中对于运动边条和弯矩的简便输入以及输出的简单解释。由ABAQUS产生的对旋转的输出代表的是关于一固定轴从原参考构型到当前构型的旋转。此输出不同于对结点的旋转历史输出。而且,这个规则成为了小旋转分析甚至是施加到初始有限旋转上的小旋转情况(这种情况可在关于一个预变形状态下小振动研究考虑)最通常的规则。
●组合旋转
因为有限旋转的方向性,它不是累积的,它们被指定的方式与其它种边条被指定的方式有所不同:在整个step上指定的旋转增量必须是从分析步开始的构型到该步末所需要的构型旋转结点所需要的旋转量。对于总旋转量,在分析中它将结点旋转到最终的位置,但如果它是在其他的几个初始参考构型中施加的,那么在一个分析步中就使结点旋转到总量来说是不够的(也可以说是不正确的)。如果对于某结点施加的一个旋转增量 是在该分析步开始(也就是上一分析步结束时)时使该结点从其旋转边界条件旋转并最终在该步结束时到其最终位置的话,那么在该步末的模型此结点的边界条件诸如旋转向量就应该是 。如果旋转向量的方向保持不变,那么定义旋转边条和总旋转向量的方法是相同的。
下面以梁的旋转为例来说明如何指定组合有限旋转以及说明有限旋转输出:
梁最初位一x轴上,现在想进行如相复合旋转:最行使梁绕z轴旋转旋转60度,随后梁自旋(以自己梁轴)90度,最后使梁绕在x-y平面内垂直于梁体的轴旋转90度。
该分析的组合旋转分为三步,分别施加旋转微量 、 和 ,其中:
在这个例子中, 、 和 。这里 代表每一个绕旋转轴有限旋转的大小。上述的三个旋转向量分别在第一个步开始的构型定义时施加。用*Boundary,type=velocity可以很简便地进行定义。为方便起见,在ABAQUS中对于velocity类型的边条其默认的amplitude值为常量。
如下为一个典型的ABAQUs分析步,将结点1固定在原点而将旋转加于结点2上:
- STEP, NLGEOM
- Step 1: Rotate 60 degrees about the z-axis
- *STATIC
- *BOUNDARY, TYPE=VELOCITY
- 2, 4, 5
- 2, 6, 6, 1.047198
- *END STEP
- **
- *STEP, NLGEOM
- Step 2: Rotate 90 degrees about the beam axis
- *STATIC
- *BOUNDARY, TYPE=VELOCITY
- 2, 4, 4, 0.785398
- 2, 5, 5, 1.36035
- 2, 6, 6
- *END STEP
- **
- *STEP, NLGEOM
- Step 3: Rotate beam onto z-axis
- *STATIC
- *BOUNDARY, TYPE=VELOCITY
- 2, 4, 4, 1.36035
- 2, 5, 5, -0.785398
- 2, 6, 6
- *END STEP
复制代码
推荐使用上述方示施加有限旋转边界条件(*boundary,type-velocity为默认的amplitude定义类型)。然而,如果旋转边条是作为位移类型边条加载的话,上面的输入格式需要进行改变,在ABAQUS/Standard中,对于某一分析步的边条是在分析步内定义整体或该步最终的边界状态。这种情况下,所有以前分析步中所指定的边条必须要被添加到增量的旋转向量分量上。ABAQUS/Standard的分析步定义将变为:
- *STEP, NLGEOM
- Step 1: Rotate 60 degrees about the z-axis
- *STATIC
- *BOUNDARY
- 2, 4, 5
- 2, 6, 6, 1.047198
- *END STEP
- **
- *STEP, NLGEOM
- Step 2: Rotate 90 degrees about the beam axis
- *STATIC
- *BOUNDARY
- 2, 4, 4, 0.785398
- 2, 5, 5, 1.36035
- 2, 6, 6, 1.047198
- *END STEP
- **
- *STEP, NLGEOM
- Step 3: Rotate beam onto z-axis
- *STATIC
- *BOUNDARY
- 2, 4, 4, 2.145748
- 2, 5, 5, 0.574952
- 2, 6, 6, 1.047198
- *END STEP
复制代码
可以看出,第二、三步的边条是分别加上了前几步已定义的旋转边界条件!
|