ABAQUS与你我的约定

weixin

嗯，没错，ABAQUS和我们有个约定，约定了什么呢？比如自由度、单位制、应力应变等符号的含义。其实对应于ABAQUS帮助文档Abaqus Analysis User's Guide 1.2.2 Conventions，如果没有遵守这些约定，冥冥之中我们就种下错误的小种子！约定的内容总共包含：自由度、坐标系统、单位、时间尺度、曲面方向、应力与应变、旋转、自由度。

      自由度      

除了轴对称单元、流体单元和电磁单元以外的其他单元自由度进行如下的约定：

其中，x、y、z默认情况下是分别与系统的整体坐标系X、Y、Z相一致的，但如果使用*Transform对结点进行局部坐标系转化的话，它们将与局部坐标系中的相关坐标轴一致。Abaqus/Standard中对于壳或梁单元最多可以20个温度自由度（自由度从11到30）。

01、轴对称单元

对轴对称单元的平动与旋转自由度如下规定：

用*transform进行结点坐标系转换的自由度改变同上。

02、可用的自由度

上述所列自由度并不是同时都能用在某一单元结点上的，不同的分析、不同的单元将拥有适合其分析的自由度，而其他自由度在某些情况下是失效的。

      坐标系统      

同一般规定一样，ABAQUS采用的基本坐标系（系统整体坐标系）是直角坐标系，方向遵循右手法则。为便于各种分析，用户可以自行定义局部坐标系以便于建模、定义材料、定义载荷以及变量输出等。

  · 建模中，常在node/ngen中加入*system；

  · 材料中，常以*orientation进行定义局部坐标系（尤其对于壳、梁单元）；

  · 载荷中，用*transform可以定义局部坐标系下的载荷；

  · 输出中，先前定义的*transform（用于结点变量）与*orientation可以发挥作用。

         单位         

实际上，ABAQUS没有单位的概念，它仅是通过有限元方法对矩阵进行数学运算得到结果，理论上没有什么物理意义，但各种变量从人为地角度赋予物理意义以及物理定理的数学表示，从而使得ABAQUS求解出“有意义的”的结果来，由此看来，结果是否有效，人对各种数据变量的主观把握是重要的。单位的一致性可以保证结果运算不会产生与之相关的问题。一般，ABAQUS建议用一套认可的单位制进行单位定义，比方说，国际单位制。当然，如果你愿意通过一系列的转化（转化因子可能复杂）的话，可以不遵守单位一致性的约定。

01、旋转及角度表示

ABAQUS中，旋转自由度 (4、5、6) 以弧度来表示，而其他角度相关的一般都以真实角度表示（如相角、*ncopy,shift中的旋转角度），其实，便于好记的说法就是：与长度或三角运算相关的用弧度表示，与旋转相关的用角度表示。

02、国际单位制 (SI)

国际单位制是我们最常用的一套符合一致性要求的单位系统 (SI)。总共有七个基本单位，这里主要用到其中的五个：长度 (m)，质量 (kg)、时间 (s)、温度 (K)、电流 (A)，其他相关单位均在此单位基础上组合或推导出来。比如力，其单位为牛顿 (N=kg·m/s2)；焦耳 (J=N·m=kg·m2/s2)。

然而，有时国际标准单位在分析中使用并不方便。比如说，结构分析中常用到的杨氏模量，其单位常用为MPa(=N/mm2)，此时，为单位一致，与之相对应的其他基本单位应为（吨、毫米、秒）。

03、美（英）式单位

对我们而言，总是不太习惯使用美式或英式单位制，这是因为其单位命名规则不像国际单位制中的表示的那样清楚。比方说，1磅力等于1磅质量 (lbm) 乘以重力加速度。如果以磅力、英尺 (ft) 和秒作为基本单位的话，则质量lbm 可以表示为：ft/sec2。但在一般情况下，密度的单位为lbm/in3，所以若要以上述基本单位表示的话，就必须转化为：sec2/ft4。是不是感觉非常不方便？

而且，从手册中也容易让我们在lb 到底代表的是lbm 还是lbf 而感到模糊。我们必须通过查看其由哪些单位组成推导出才能确认是什么。

另外，难理解的还有两个单位，一个是斯（勒格）-slug(＝lbf sec2/ft)，其被定义为在受到一磅重的力作用时产生每秒钟一英尺的加速度的质量单位。另一个是磅达-poundal，被定义为使质量1磅的质点发生1尺/秒加速度的力。以下两个有用的转化为slug=g·lbm 和lbf=g·poundals，这里g 为重力加速度值。

04、ABAQUS单位符号表示

时间尺度      

ABAQUS共有两种时间计法，一种是分析步时间 (step time)，另一种是总分析时间 (total time)。除了线性扰动分析（它不考虑时间，频域分析），分析步时间是从每一分析步开始计算，而总分析时间则是从第一个分析步开始计算的所有分析步的时间积累（包括*restart步）。

      曲面方向      

完全的空间曲面定义需要有局部面方向的定义，这样在诸如基于单元接触曲面的切向滑移方向或是壳单元的应力应变方向均能得到一致的表示。为此，ABAQUS作如下规定：

  · 默认的局部1方向是整体x 轴向曲面的投影，如果该x 轴垂直于曲面（与曲面法线夹角小于0.1度），其局部1方向是整体z轴向曲面的投影。

  · 局部2方向则是局部1方向依右手法则形成，故而局部1、2和曲面的正法线的构形局部坐标轴方向（依右手法则），如下图。

  · 曲面的正法线方向是通过构成该单元的结点以右手法则旋转而成，局部曲面方向可以通过*orientation定义。

  · 当考虑垫片单元或与*section print和*section file相关的局部坐标系统定义时，曲面的局部1-、2-方向变成局部2-、3-方向。

  · 对于线型单元如梁beam、管pipe、杆truss的空间方向，其默认的局部1-方向和2-方向该单元的切向与横向（其具体方向与结点定义顺序有关），当然也可以通过*orientation进行修改。

01、局部旋转方向

对于几何线性分析而言，以默认的材料方向（初始参考中定义）就可以将应力应变表示出来。

对于几何非线性分析，在ABAQUS/Standard中的小应变壳单元 (S4R5, S8R, S8R5,S8RT, S9R5, STRI3, 和STRI65) 使用总体拉格朗日应变算法，应力应变可以相对于参考构型的材料方向绝定。垫片单元是小应变小位移单元，默认情况下其应力应变值也是以初始参考构型定义的行为方向输出。

对于有限膜应变单元（所有的膜单元以及S3/S3R, S4, S4R, SAX, 和 SAXA单元）和在ABAQUS/Explicit中的小应变单元，其材料方向是随着曲面的平均刚性旋转运动，而变形成当前构型的材料方向。此时，这些单元的应力应变则是根据当前的参考构型中的材料方向给出的（更详细地说明可以参考ABAQUS相关手册）。用户可以决定与*section print和*section file相关的局部坐标系统是固定不动还是随着曲面的平均刚性运动而旋转。

     应力与应变     

在定义材料特性时，在ABAQUS中应力应变的各个分量定义如下顺序：

比如，一个完全各向异性、线弹性矩阵为：

1-、2-和3-方向的确定依靠所选择的单元类型，对于实体单元就是整体坐标系的空间方向。对于壳单元和膜单元，1-、2-方向默认上是壳或膜曲面的局部方向，可以通过*orientation对三个方向进行适合自己分析的修改。

对于实体单元的几何非线性分析，默认的方向并不会随着材料的旋转而变化。然而，如果通过*orientation定义的方向则会旋转。

ABAQUS/Explicit在分析中，在内部以不同的顺序完成对各应力的存储。对于几何非线性分析，无论是否使用*orientation，这些被存储的变量总是随着材料的方向旋转的，这一点在子程序VUMAT被使用时尤显重要。

01、各向异性材料行为

在连续体单元中定义各向异性材料行为时，用*orientation所定义的材料方向是很重要且必要的。

02、零值应力

在分析中，如果应力值一直为零，在矩阵存储时将被忽略。例如，对于平面应力分析，ABAQUS仅存储两个面内的正应力和一个面内的剪应力。

03、剪应变

ABAQUS输出工程剪应变：

04、应力应变测量

在ABAQUS中，应力使用的是柯西应力或真实应力值，即每单位面积（变形后的面积）上的力大小。详细地请看手册"Stress measures" Section 1.5.2 of theABAQUS Theory Manual。

对于几何非线性分析，存在多种不同的应变计量方法。不像真应力，它没有很清楚的真应变的计算方法。对同一种物理变形，在大应变分析中不同的应变测量方法会给出不同的应变值，当然其值所能反应的实际性也就不同。如何选择最好的应变测量方法，这依靠分析类型、材料行为以及在一定程度上也依靠个人喜好！详见"Strain measures" Section 1.4.2 of the ABAQUS Theory Manual。

默认情况下，在ABAQUS/Standard中应变以“综合”总应变 (E) 输出，对于大应变的壳、膜和实体单元，还有两种其他的总应变计法可以输出：自然应变 (LE) 和公称应变 (NE)。

在ABAQUS/Explicit中，自然应变 (LE) 是默认的应变输出，也可以要求公称应变输出。而“综合”应变在ABAQUS/Explicit中是无法得到的。

05、总（综合）应变

默认情况下，ABAQUS/Standard输出到.dat 文件和.fil 文件的应变为综合应变，这是对于所有可以在材料体下将应变率数值积分获得有限应变的单元都适用的。

其中，ε[n+1]和εn 分别是分析中第n+1次和第n 次的总应变，ΔR 为旋转张量的增量，Δε 为从第n 次增量到第n+1次增量的总应变增量。对于使用正转坐标系统（右手法则）的单元（使用*orientation有限应变壳单元、膜单元和实体单元）上式可以简化为：

应变增量可以通过对变形率D 在整个时间增量上积分得到：

这种应变计法对于弹塑性（粘弹性材料）或弹性－蠕变材料都是适合的，这是因为塑性应变和蠕变应变都是通过相同的积分方式得到的。在这样的材料里，弹性应变是很小的（因为屈服应力相对于弹模来说是较小的），此时总（综合）应变可以直接地与塑性应变和蠕变应变相对比。

如果应变的主方向关于材料方向的变化而旋转变化的，那么最后所得到的应变是不能和总变形相联系的。此时，无论采用的是何种坐标系统。如果主应变方向保持固定，那么应变就是变形率的积分：

这与稍后将要讨论的自然应变等效。

06、格林应变

在ABAQUS/Standard中，对于小应变壳和梁单元，默认的应变计法E 为格林应变：

此时，F 为变形适量梯度，而I 为特征张量。这种应变计法适合于在小应变、大旋转分析中的单元使用。分量εG 代表沿原定义方向的应变，不能在有弹塑性或超弹性材料行为的有限应变分析中使用小应变壳和梁单元，因为可能会导致不正确的分析结果！

07、公称应变

公称应变NE 被定义为：

其中，V=√F·FT 为左拉伸张量，λi 为主拉伸，ni 为在当前材料参数下的主拉伸方向。因此，公称应变也就是在主材料参数方向下长度变化对原长度的比值，从而直观地解释了变形。

08、自然应变（对数应变）

自然应变LE 被定义为：

式中，变量同上。

09、应力不变量

在ABAQUS中许多结构模型都是根据应力不变量进行阐述的，不变量的概念是指：数值不会随着坐标系的转换而变化，这些应力不变量为：

等效压应力：

Mises等效应力：

第三偏应力不变量：

这中间，S 为偏应力，被定义成：

旋转         

空间中，对于有限旋转作如下规定：定义Φx、Φy、Φz 为关于整体X、Y和Z轴的旋转量（也就是指在一个节点的4、5、6自由度）。

定义，px=Φx/Φ，py=Φy/Φ，pz=Φz/Φ，方向P 为旋转同方向，其中Φ 为右手法则得出的旋转角度（弧度表示）。值不是唯一确定的，在大旋转问题中，如果总旋转量超过2π，可以加或减任意倍的2π（以得到在-2π到2π内的值），而这可能会导致旋转分量的不连续输出。

这个规则提供了在大多数分析中对于运动边条和弯矩的简便输入以及输出的简单解释。由ABAQUS产生的对旋转的输出，代表的是关于一固定轴从原参考构型到当前构型的旋转。此输出不同于对结点的旋转历史输出，而且，这个规则成为了小旋转分析甚至是施加到初始有限旋转上的小旋转情况（这种情况可在关于一个预变形状态下小振动研究考虑）最通常的规则。

01、组合旋转

因为有限旋转的方向性，它不是累积的，它们被指定的方式与其它种边条被指定的方式有所不同。在整个分析步上，指定的旋转增量必须是从分析步开始的构型到该步末所需要的构型旋转结点所需要的旋转量。对于总旋转量，在分析中它将结点旋转到最终的位置，但如果它是在其他的几个初始参考构型中施加的，那么在一个分析步中就使结点旋转到总量来说是不够的（也可以说是不正确的）。如果对于某结点施加的一个旋转增量 (ΔΦx , ΔΦy , ΔΦz )是在该分析步开始（也就是上一分析步结束时）时使该结点从其旋转边界条件 (Φx , Φy , Φz )t 旋转并最终在该步结束时到其最终位置的话，那么，在该步末的模型此结点的边界条件诸如旋转向量就应该是(Φxt+ΔΦx , Φyt+ΔΦy , Φzt+ΔΦz )。如果旋转向量的方向保持不变，那么定义旋转边条和总旋转向量的方法是相同的。

来源：CAE爱联盟微信公众号（ID：hxhycae），文章整理自网络。