《ABAQUS规则.doc》由会员分享，可在线阅读，更多相关《ABAQUS规则.doc（9页珍藏版）》请在三一文库上搜索。
1、ABAQUS勺定相关 每种软件在顺利运行中都有自己的一套在诸如单位、 符号、 变量值表示等方面的约 定用法， 如果想用此种软件进行适合自己的分析， 自己进行主观操作之外， 对它的这种勺定 我们也要提起注意，否则很容易产生我们觉察不到的问题。（参考 abaqus analysis manual中 1.2.2 Conventions ）1.自由度2.坐标系统3.单位4.时间尺度5.曲面方向6.应力与应变7.旋转=1 =自由度Abaqus中对单位的认定与其他软件（如 ANSY$稍微有点不同就在于默认情况下abaqus是以 1、2、3 等数字来表示各种自由度的标符的，在手写inp 中，只能以它们表示自
2、由度。除了轴对称单元（ .ax. ）以外，其它单元对自由度进行如下勺定：1. x 方向（平动自由度）2. y 方向。3. z 方向。4. 绕 x 轴旋转的旋转自由度（以弧度表示）5.6.绕 y 。 绕 z 。7. 翘曲（对于开口截面梁单元）8. 孔隙压力（或静水压）9. 电势11温度（或质量扩散分析中的归一化浓度）12. 第二温度（对于壳、梁）13. 第三温度。14. 其他其中，x、y、z默认情况下是分别与系统的整体坐标系X、Y、Z相一致的，但如果使用*Transform对结点进行局部坐标系转化的话，那么它们将与局部坐标系中的相关坐标轴一致。 对轴对称单元的平动与旋转自由度如下规定：1. r
3、方向（径向）位移2. z 方向（轴向）位移5. 绕 z 轴旋转（用于带扭曲的轴对称单元） ，以弧度表示6. r-z 平面的旋转（用于轴对称壳单元） ，以弧度表示用*tra nsform 进行结点坐标系转换的自由度改变同上。 可用的自由度 上述所列自由度并不是同时都能用在某一单元结点上的，不同的分析，不同的单元自会 有适合其分析的自由度，而其他则在此是失效的。 ABAQUS/Standard中的内部变量除了上述所列的自由度外，ABAQUS/Standard对某些单元还内定了内部自由度变量（如用于施加约束的拉格朗日乘子），一般情况下，使用 ABAQUS分析并不需要去了解这 些变量，但在进行分析过程
4、中，当迭代中对非线性约束的满足进行检验时常用到这些内 部变量，这从msg文件中的错误警告信息中可以看到。内变量与内部结点相关，内部结 点在ABAQUS为系统分析所用，以负的结点号出现以别于我们所定义的结点。=2=坐标系同一般规定一样，ABAQU睬用的基本坐标系（系统整体坐标系）是直角坐标系，方向遵 循右手法则。为便于各种分析，用户可以自行定义局部坐标系以便于建模、定义材料、定义载荷以及变量输出等。建模中，常在 node/ngen中加入*system材料中，常以*orientation进行定义局部坐标系（尤其对于壳、梁单元）载荷中，用*transform 可以定义局部坐标系下的载荷输出中，先前定
5、义的*tran sform （用于结点变量）与*orie ntation可以发挥作用。=3=单位实际上，ABAQU毀有单位的概念，它仅是通过有限元方法对矩阵进行数学运算得到结果， 理论上没有什么物理意义，但各种变量从人为地角度赋与物理意义以及物理定理的数学表 示，从而发动 ABAQU进行求解出有意义的的结果来，由此看来，结果是否有效，人对 各种数据变量的主观把握是重要的。单位的一致性可以保证结果运算不会产生与之相关的问题。一般，ABAQUS建议用一套认可的单位制进行单位定义，比方说，国际单位制。当然， 如果你原意通过一系列的转化（转化因子可能复杂）的话，可以不遵守单位一致性的约定。（）旋转及角
6、度表示ABAQUS，旋转自由度（4、5、6）以弧度来表示，而其他角度相关的一般都以真实角 度表示（如相角、*n copy,shift中的旋转角度），其实，便于好记的说法就是，与长度或三角运算相关的用弧度，与旋转相关的用角度表示。（个人理解）国际单位制国际单位制是我们最常用的一套符合一致性要求的单位系统（SI）。其基本的单位共五个：长度（m）,质量（kg ）、时间（s）、温度（K）、电流（A）。其他相关单位均在此单 位基础上组合或推导出来。如力，其单位为牛顿（N）：1 Newton = 1 kg in/s2,类似地，库伦在国际单位制中为：1 Coulomb = 1 A s.焦耳（J）:1 Jou
7、le = 1 N m = I A olt s = 1 kg nr/sJ.电势在国际单位制中为伏特（ V）,其单位标准为：1 Joule = 1 Vult C = 1 Wilt A s,然而，有时国际标准单位在分析中使用并不方便。比如说，结构分析中常用到的杨氏模 量，其单位常用为 MPa（= N/mmA2）,此时，为单位一致，与之相对应的其他基本单位应 为（1吨=1000KG千米、秒）。美（英）式单位对我们而言，总是不太习惯使用美式或英式单位制，这是因为其单位命名规则不像国际单位制中的表示的那样清楚。比方说，1磅力等于1磅质量（Ibm）乘以重力加速度（单位为）。如果以磅力、英尺（ft）和秒作为基
8、本单位的话，则质量Ibm可以表示为：述9它。但在一般情况下，密度的单位为，所以若要以上述基本单位表1 bm/in3 = Ibf示的话，就必须转化为：（不方便？）而且，从手册中也容易让我们在lb倒底代表的是lbm还是Ibf而感到模糊。我们必须通过查看其由哪些单位组成推导出才能确认是什么。另外，难理解的还有两个单位，一个是斯（勒格）slug （= Ibf secA2/ft ）,其被定义为在受到一磅重的力作用时产生每秒钟一英尺的加速度的质量单位。另一个是磅达一poundal （是这么称呼的吧，？），被定义为使质量1磅的质点发生1尺/秒加速度的力。 以下两个有用的转化为：訂，和1.门！。这里g为重力加
9、速度值（单位为：）。 ABAQUS!位符号表示单位符号SI单位len gthLmetermassMkilogramtimeTsecondtemperature0degree Celsiuselectric curre ntAAmpereforceFNewt onen ergyJJouleelectric chargeCCoulombelectric pote ntialPIvoltmass concen trati onPParts per milli on=4=时间ABAQUS共有两种时间计法，一种是步时间（ step time ）另一种是总分析时间（total time ）。除了线性扰动分
10、析（它不考虑时间），步时间是从每一分析步开始计算，而总分析时间则是从第一个 step开始计算起的所有 step的时间积累（包括*restart 步）。=5=关于空间曲面的局部方向完全的空间曲面定义需要有局部面方向的定义以完成定义诸如基于单元接触曲面的切向 滑移方向或是壳单元的应力应变方向。对此类方向ABAQUS作如下规定：默认的局部1方向是整体x轴向曲面的投影，如果该 x轴垂直于曲面（与曲面法线夹角 小于0.1度），其局部1方向是整体z轴向曲面的投影。而局部 2方向则是局部1方向依右 手法则形成，故而局部1、2和曲面的正法线的构形局部坐标轴方向（依右手法则）。如下图。曲面的正法线方向是通过构成
11、该单元的结点以右手法则旋转而成。局部曲面方向可以通过*orie ntatio n定义。surface normal当考虑垫片单元或与*section print和*section file 相关的局部坐标系统定义时， 曲面的局部1-、2-方向变成局部2-、3 方向。对于线型单元如梁 beam管pipe、杆truss的空间方向，其默认的局部 1-方向和 2 方向该单 元的切 向与 横向（其具 体方向 与结点定义 顺序 有关），当然 也可 以通过 *orie ntation 进行修改。局部旋转方向对于几何线性分析而言，以默认的材料方向（初始参考中定义）就可以将应力应变表示 出来。对于几何非线性分析
12、，在ABAQUS/Standard中的小应变壳单元 （S4R5, S8R, S8R5, S8RT,S9R5, STRI3,和STRI65 ）使用总体拉格朗日应变算法，应力应变可以相对于参考构型 的材料方向改定。垫片单元是小应变小位移单元，默认情况下其应力应变值也是以初始 参考构型定义的行为方向输出。对于有限膜应变单元（所有的膜单元以及S3/S3R, S4, S4R, SAX,和SAXA单元）和在ABAQUS/Explicit中的小应变单元，其材料方向是随着曲面的平均刚性旋转运动而变以 形成当前构型的材料方向。此时这些单元的应力应变则是根据当前的参考构型中的材料方向给出的。（更详细地说明可以参考
13、ABAQUS相关手册）。用户可以决定与*sectionprint和*section file相关的局部坐标系统是固定不动还是随着曲面的平均刚性运动而旋转。=6=应力应变的规定在定义材料特性时，在ABAQUS应力应变的各个分量定义如下顺序：:1-方向的正应力:2方向的正应力:3方向的正应力:1-2面的剪应力:1-3方向的剪应力:2-3方向的剪应力比如，一个完全各向异性、线弹性矩阵为:/11OHLJD 122Dii囲5订卫DinsQiL)2-|2”333312D3323SX 111111 1212Dins*12TnJ?131313血3/D沁u -丿1-、2-和3-方向的确定依靠所选择的单元类型。对
14、于实体单元就是整体坐标系的空间 方向。对于壳单元和膜单元，1-、2-方向默认上是壳或膜曲面的局部方向。可以通过*orie ntation 对三个方向进行适合自己分析的修改。对于实体单元的几何非线性分析，默认的方向并不会随着材料的旋转而变化。然而，如果通过*orientation定义的方向则会旋转。ABAQUS/Explicit在分析中，在内部以不同的顺序完成对各应力的存储。对于几何非线性分析，无论是否使用*orientation ，这些被存储的变量总是随着材料的方向旋转的，这一点在子程序 VUMAT被使用时尤显重要。各向异性材料行为在连续体单元中定义各向异性材料行为时，用*orie ntati
15、o n所定义的材料方向是很重要且必要的。零值应力在分析中，如果应力值一直为零，在矩阵存储时将被忽略。例如，对于平面应力分析， ABAQUS仅存储两个面内的正应力和一个面内的剪应力。剪应变ABAQU輸出工程剪应变:* ij = ij +应力应变测量在ABAQU中，应力使用的是柯西应力（？）或真实应力值，即每单位当前面积上的应大小。详细地请看手册。？ Stress measures, ? Section 1.5.2 of the ABAQUSTheory Manual。对于几何非线性分析，存在多种不同的应变计量方法。不像真实应力，它没有 很清楚的真实应变的计法。对同一种物理变形，在大应变分析中不同
16、的应变测量方法会 给出不同的应变值，当然其值所能反应的实际性也就不同。如何选择最好的应变测量方 法，这依靠分析类型、材料行为以及在一定程度上也依靠个人喜好！详见？Strainmeasures, ? Section 1.4.2 of the ABAQUS Theory Manual默认情况下，在 ABAQUS/Standard中应变以？综合？总应变（ 日输出，对于大应变 的壳、膜和实体单元，还有两种其他的总应变计法可以输出：自然应变（LE）和公称应变（NE）。在ABAQUS/Explicit中，自然应变（LE）是默认的应变输出，也可以要求公称应变 输出。而？综合？应变在ABAQUS/Explic
17、it中是无法得到的。总（综合）应变默认情况下，ABAQUS/Standard输出到.dat文件和.fil文件的应变为综合应变，这是对于所有可以在材料体下将应变率数值积分获得有限应变的单元都适用的：川川 ，其中 和分别是分析中第 n+1次和第n次的总应变，_1：为旋转张量的增量, 用正转坐标系统（右手法则） 单元）上式可以简化为： 曰+| = eft + 上.应变增量可以通过对变形率D在整个时间增量上积分得到:-为从第n次增量到第n+1次增量的总应变增量。对于使 的单元（使用 *orie ntati on有限应变壳单元、膜单元和实体f严+1= / DJ（粘弹性材料）或弹性-蠕变材料都是适合的，这
18、是因为塑性应变在这样的材料里，弹性应变是很小的（因为屈 ，此时总（综合）应变可以直接地与塑性应变和蠕变应变这种应变计法对于弹塑性 和蠕变应变都是通过相同的积分方式得到的。服应力相对于弹模来说是较小的） 相对比！如果应变的主方向关于材料方向的变化而旋转变化的，那么最后所得到的应变是不能和总变形相联系的，此时无论采用的是何种坐标系统。如果主应变方向保持固定，那么应变就是变形率的积分：Ddt.J0这与稍后将要讨论的自然应变等效。格林应变在ABAQUS/Standard中，对于小应变壳和梁单元，默认的应变计法E为格林应变:产二（F匚 F-I）*2此时，F为变形适量梯度而I为特征张量。这种应变计法适合于
19、在小应变、大旋转分析 中使用这些单元。分量代表沿原定义方向的应变，不能在有弹塑性或超弹性材料行为的有限应变分析中使用小应变壳和梁单元，因为可能会导致不正确的分析结果！公称应变公称应变NE被定义为：3 v = V - I = Z - Hmm，其中，I 1为左拉伸张量,为主拉伸，为在当前材料参数下的主拉伸方从而直观地解向。因此，主公称应变也就是在主材料参数方向下长度变化对原长度的比值, 释了变形。自然应变自然应变LE被定义为：el = In V =hi A； n, n/ .-（式中变量同上）。应力不变量在ABAQU叩许多结构模型都是根据应力不变量进行阐述的，这些应力不变量为:1等效压应力：p -
20、trace ta):oMises等效应力:r = （-S S : S）:第三偏应力不变量：_这中间，S为偏应力，被定义成：、T 丨=七=有限旋转 空间中，对于有限旋转作如下规定： 定义、 、 为关于整体X、Y和Z轴的旋转量（也就是指在一个节点的 4、5、6自由度）。定义八、，其中，方向P为旋转同方向，为右手法则得出的旋转角度（弧度表示），见下图。，值不是唯一确定的。在大旋转问题中，如果总旋转量超过 一，可以加或减任意倍的一（以得到在一内的值）， 而这可能会导致旋转分量的不连续输出。这个规则提供了在大多数分析中对于运动边条和弯矩的简便输入以及输出的简单解释。由 ABAQU萨生的对旋转的输出代表的
21、是关于一固定轴从原参考构型到当前构型的旋转。此输 出不同于对结点的旋转历史输出。而且，这个规则成为了小旋转分析甚至是施加到初始有限旋转上的小旋转情况（这种情况可在关于一个预变形状态下小振动研究考虑）最通常的规则。组合旋转因为有限旋转的方向性，它不是累积的，它们被指定的方式与其它种边条被指定的方式有所不同：在整个step上指定的旋转增量必须是从分析步开始的构型到该步末所需要的构型 旋转结点所需要的旋转量。对于总旋转量，在分析中它将结点旋转到最终的位臵，但如果它是在其他的几个初始参考构型中施加的，那么在一个分析步中就使结点旋转到总量来说是不够的（也可以说是不正确的）。如果对于某结点施加的一个旋转增
22、量2是 在该分析步开始（也就是上一分析步结束时）时使该结点从其旋转边界条件旋转并最终在该步结束时到其最终位臵的话， 那么在该步末的模型此结点的边界条件诸如旋转向量就应该是1- - - 一 。如果旋转向量的方向保持不变，那么定义旋转边条和总旋转向量的方法是相同的。下面以梁的旋转为例来说明如何指定组合有限旋转以及说明有限旋转输出：梁最初位一 x轴上，现在想进行如相复合旋转：最行使梁绕z轴旋转旋转60度，随后梁自旋（以自己梁轴）90度，最后使梁绕在 x-y平面内垂直于梁体的轴旋转90度。该分析的组合旋转分为三步，分别施加旋转微量*、和，其中：01=斷6=0 01.0471987 *仇=6*用卩1务
23、+切1 门片）=（J853981JG0350 （）1 atid0靠=t3（sinr coe勺）= 1.360350 0*785398。丁 在这个例子中，。:一匸強、和，。这里代表每一个绕旋转轴有限旋转的大小。上述的三个旋转向量分别在第一个步开始的构型定义时施加。用*Boundary,type=velocity可以很简便地进行定义。为方便起见，在ABAQU中对于velocity类型的边条其默认的 amplitude值为常量。如下为一个典型的 ABAQUs分析步，将结点1固定在原点而将旋转加于结点2 上:STEP, NLGEOMStep 1: Rotate 60 degrees about the
24、 z-axis*STATIC*BOUNDARY, TYPE=VELOCITY2, 4, 52, 6, 6, 1.047198*END STEP*STEP, NLGEOMStep 2: Rotate 90 degrees about the beam axis*STATIC*BOUNDARY, TYPE=VELOCITY2, 4, 4, 0.7853982, 5, 5, 1.360352, 6, 6*END STEP*STEP, NLGEOMStep 3: Rotate beam onto z-axis*STATIC*BOUNDARY, TYPE=VELOCITY2, 4, 4, 1.36035
25、2, 5, 5, -0.7853982, 6, 6*END STEP推荐使用上述方示施加有限旋转边界条件（*boundary,type-velocity为默认的amplitude 定义类型）。然而，如果旋转边条是作为位移类型边条加载的话，上面的输入格式需要进行 改变，在 ABAQUS/Standard 中，对于某一分析步的边条是在分析步内定义整体或该步最终的 边界状态。 这种情况下， 所有以前分析步中所指定的边条必须要被添加到增量的旋转向量分 量上。 ABAQUS/Standard 的分析步定义将变为：*STEP, NLGEOMStep 1: Rotate 60 degrees about t
26、he z-axis*STATIC*BOUNDARY2, 4, 52, 6, 6, 1.047198*END STEP*STEP, NLGEOMStep 2: Rotate 90 degrees about the beam axis*STATIC*BOUNDARY2, 4, 4, 0.7853982, 5, 5, 1.360352, 6, 6, 1.047198*END STEP*STEP, NLGEOMStep 3: Rotate beam onto z-axis*STATIC*BOUNDARY2, 4, 4, 2.1457482, 5, 5, 0.5749522, 6, 6, 1.047198*END STEP 可以看出，第二、三步的边条是分别加上了前几步已定义的旋转边界条件！