《ABAQUS个人使用笔记》由会员分享，可在线阅读，更多相关《ABAQUS个人使用笔记（7页珍藏版）》请在金锄头文库上搜索。
1、ABAQUS 学习笔记一AQUS .inp 编码介绍（一） ABAQUS 头信息文件段（ 14）1.*PREPRINT 输出求解过程所要求的信息（在 dat 文件中）ie:*PREPRINT, ECHO=YES, HISTORY=YES, MODEL=YES2.*HEADING 标题输出文件（出现在 POST/VIEW 窗口中，且出现在结果输出文件中）ie:*HEADINGSTRESS ANALYSIS FOR A PLATE WITH A HOLE3.*RESTART 要求 abaqus/standard 输出其 POST/view 模块所需要的.res 文件。其中的FREQ？控制结果在每次
2、迭代（或载荷步）输出的次数。ie:*RESTART, WRITE, FREQ=14.*FILE FORMAT 要求 abaqus/standard 输出到.fil 中的某些信息。它也用于 post。对于在后处理中得到 x-y 形式的诸如应力时间、应力应变图有用！ie: *FILE FORMAT, ZERO INCREMENT（二） ABAQUS 网格生成段定义结点、单元，常用的命令有：结点定义（*NODE，*NGEN） ，单元定义（*ELEMENT，*ELGEN 等） 。1.*NODE 定义结点，其格式为：*NODE 结点号，x 轴坐标，y 轴坐标， （z 轴坐标）2.*NGEN 在已有结点的
3、基础上进行多个结点的生成，一般是在两结点间以某种方式（直线、圆）产生一定分布规律的结点。如：*NGEN, LINE=C, NSET=HOLE，119, 1919, 100, 101 在两结点（结点号为 119，1919）间以圆弧形式生成多个结点，100 为任意相邻结点的单元号增量，101 为圆弧形成时圆心位置的结点（对于直线形式生成没有此结点） 。所有这些生成的结点（包括 119，1919）被命名成 HOLE 的集合（这样做的目的是以后的命令中使用到它，比如说对这些结点施加同等条件的边界条件或载荷等，HOLE 就是这些结点的代称） 。*NGEN 使用的前提就是必须存在已有结点。*NGEN, N
4、SET=OUTER131, 1031, 100 以线形式形成结点，结点号增量 100，结点集合名为 OUTER。*NGEN, NSET=OUTER1031, 1931, 100 同上生成结点，可以同上结点集合名，这样 OUTER 就包括这两次生成的所有结点3.*NFILL 在如上生成的结点集（实际上，代表两条几何意义上的边界线）之间按一定规律（BIAS？）填充结点。这样所有生成的结点构成一定形状的实体（面） 。如：*NFILL, NSET=PLATE, BIAS=0.8HOLE, OUTER, 12, 1 以 HOLE 为第一条边界，OUTER 为第二条边界（终止边） ，以从疏到密的规律（BI
5、AS 小于 1）分布，其生成结点数在两内外对应结点间为12，1 为每组结点号的增量。所有这些结点被置于 PLATE 的集合中。下面以上面生成的结点来生成单元：4.*ELEMENT 定义单元所使用的类型（TYPE=？） ，然后另行定义通过联结结点形成单元，其结点数目依靠单元类型而变。*ELEMENT, TYPE=CPS4 /采用四单元的平面应力单元19, 119, 120, 220, 219 /定义顺序：单元号，以逆时针方向形成单元的各结点号（三）ABAQUS 单元注意：分析前要选择合适的元素，这时要考虑的问题就是：使用什么样类型的单元？有限元的基本思路就是将实际中的连续体离散化，实际结果是将众
6、多离散分析结果的集合，这似乎有点像积分的概念。选择元素种类最重要考虑的是分析必要的现象，满足必要的准确度基础上去掉不必要的细节与准确度。是选择 1-D, 2-D or 3-D 单元、用于何种分析的单元、是否高阶单元等。（四）ABAQUS 材料ABAQUS 本身提供了丰富的材料库供分析使用，并已能满足常用的分析。但对于新型本构关系的材料 abaqus 本身是无法体现的，UMAT 则为这个问题提供了解决。自己编程将材料的应力应变本构表示出来，ABAQUS 调用完成分析。ABAQUS 的材料行为模式主要分为弹性材料： Linear elasticity （线弹性） No compression or
7、 tension elasticity （无压缩或位伸弹性材料，即单力性材料） Plane stress orthotropic failure （平面应力单元） Porous elasticity （多孔弹性） Hypoelasticity （亚弹性） Hyperelasticity （超弹性） Foam elasticity （泡沫单元） Viscoelasticity （粘弹性）非弹性材料 Classical metal plasticity （塑性） Metals subjected to cyclic loading （受周期荷载金属单元） Rate-dependent yield（
8、率相关屈服单元） Creep and Swelling （蠕变） Anisotropic yield and creep （各向异性） Porous metal plasticity (多孔塑性) Deformation plasticity （塑变单元） Granular materials or polymers （粒状材料或复合材料） Clay plasticity （粘土塑性） Crushable foam plasticity （可压泡沫塑性） Jointed material （？） Concrete （混凝土）（五）ABAQUS 求解对于一个 inp 文件，不进入 CAE 时，需
9、要这样做：1. 检查 inp 文件的正确性（当然主要是指 keyword 的使用） ，自己能做检查最好，否则可以通过：ABAQUS datacheck job=yourjobname2. 检查确认修正后进行计算：通过：ABAQUS job=yourjobname3. 检验分析结果的合理性：不只是会算，更要会对分析结果进行确认。首先要对整个分析及分析的并键之处成竹在心。然后可以通过以下途径作结果确认：自已能够得到的解析解实验数据其它数值解别人的求解结果（当然你得信任他）直觉与经验4. 如果迭代无法收敛：需要通过.msg,.sta 文件查看出错信息并做出判断（在 CAE 中submit 分析时可以
10、通过 monitor 查看） ，判断依据为：结构约束是否足够或过多材料数据是否正确单元是否适合此分析网格有没有过扭曲、奇异接触单元是否足够步长是否过大二有限元理论（一）关于应力应变金属的工程应力（未变形单位面积上的力）称为名义应力，与之相对应的为名义应变（每单位未变形长度的伸长） 。 -名义应力 -名义应变0/FA0/Ll在只考虑 的情况下，拉伸和压缩应变是相同的，即：ld，其中 l 是当前长度， 是原始长度， 为真实应变或对数00ln()l 0l应变。与真实应变对应的真实应力： ，F 为材料受力，A 是当前面积。 在 ABAQUS 中必须用真实应力和真实应变定义塑性.ABAQUS 需要这些值
11、并对应地在输入文件中解释这些数据。然而，大多数实验数据常常是用名义应力和名义应变值给出的。这时，必须应用公式将塑性材料的名义应力（变）转为真实应力（变） 。考虑塑性变形的不可压缩性，真实应力与名义应力间的关系为：，0lA当前面积与原始面积的关系为：0l将 A 的定义代入到真实应力的定义式中，得到：00()nomFllA其中 也可以写为 。0l1nom这样就给出了真实应力和名义应力、名义应变之间的关系：(1)nomno真实应变和名义应变间的关系很少用到，名义应变推导如下： 01nomll上式各加 1，然后求自然对数，就得到了二者的关系：l()nomABAQUS 中的*PLASTIC 选项定义了大
12、部分金属的后屈服特性。ABAQUS 用连接给定数据点的一系列直线来逼近材料光滑的应力应变曲线。可以用任意多的数据点来逼近实际的材料性质;所以，有可能非常逼真地模拟材料的真实性质。在*PLASTIC 选项中的数据将材料的真实屈服应力定义为真实塑性应变的函数。选项的第一个数据定义材料的初始屈服应力，因此，塑性应变值应该为零。在用来定义塑性性能的材料实验数据中，提供的应变不仅包含材料的塑性应变，而是包括材料的总体应变。所以必须将总体应变分解为弹性和塑性应变分量。弹性应变等于真实应力与杨氏模量的比值，从总体应变中减去弹性应变，就得到了塑性应变，其关系为：/plteltE其中 是真实塑性应变， 是总体真
13、实应变， 是真实弹性应变。plt el总体应变分解为弹性与塑性应变分量实验数据转换为 ABAQUS 输入数据的示例下图中的应力应变曲线可以作为一个例子，用来示范如何将定义材料塑性特性的实验特性的实验数据转换为 ABAQUS 适用的输入格式。名义应力应变曲线上的 6 个点将成为*PLASTIC 选项中的数据。第一步是用公式将名义应力和名义应变转化为真实应力和应变。一旦得到这些值，就可以用公式不确定与屈服应力相关联的塑性应变。下面给出转换后的数据。在小应变时，真实应变和名义应变间的差别很小，而在大应变时，二者间的就会有明显的差别;因此，如果模拟的应变比较大，就一定要向 abaqus 提供正确的应力
14、应变数据。定义这种材料的输入数据格式在图中给出。（二）. 对于受力的大小,受力的方式,还有本构方程参数的选择对于模型是否收敛影响很大.泊松比的影响：材料的泊松比的大小对于网格的扰动影响很大，在 foam 中，由于其泊松比是 0，所以它对于单元的扰动不是很大。所以在考虑到经常出现单元节点被翻转过来的现象，可以调整泊松比的大小。REMESH：对于 creep 的，特别是材料呈现非线性的状态下，变形很大，就有必要对其进行重新划分网格，用 map solution 来对其旧网格进行映射。这就要决定何时进行重新划分网格，这个就要看应变的增长幅度了，通过观察网格外形的变化曲线来决定是否要进行重新划分区域。
15、接触表面的 remesh 时，网格类型，单元数目等必须和原有的 mesh 保持一致，这个对于 contact 的计算十分重要。但是对于刚体表面的 remesh 没有这个必要的，单元数目可以减少，网格可以粗化，但是对于非刚体，一般将网格进行细化。对于 NIGEOM（非线性）：the load must be applied gradually. We apply the load gradually by dividing the step into increments。Omit this parameter or set NLGEOM=NO to perform a geometrically linear analysis during the current step. Include this parameter or set NLGEOM=YES to indicate that geometric nonlinearity should be accounted for during the step (stress analysis and fully coupled thermal-stress analysis only). Once the NLGEOM option has been switched on, it will be a
