abaqus分析收敛的个人经验整理
说一下自己在分析收敛方面的一些经验
接触分析收敛不管怎么总还是一个很大的问题，而我们经常在一个地方卡了很长的时间，怎么也找不到解决和提高的办法。而aba_aba在abaqus常见问题汇总中给了我们模型改进的方向和一些方法。在我分析的过程当中，怎么找到模型中的影响收敛的关键问题所在也是一个很让我迷茫了很长时间。下面谈一下我个人的一些经验和看法。如有错误还望大家指出，也希望大家给出自己更多的经验分享。
abaqus的隐式求解的就是求算出一个很大的刚度矩阵的解，这个方程能否通过一次一次的迭代到最后达到一个系统默认的收敛准则标准的范围之内，就决定了这一次计算能否收敛。因此要收敛的话，系统与上一个分析步的边界条件区别越小的话，系统就越容易找到收敛解。针对这一点，我们可以得到下面的几种方法来尽可能的使系统的方程的解尽可能的接近上一步，以达到收敛。下面的方法的指导思想是：尽可能小的模型，前后两个分析步的改变尽可能的少。
1.接触分析真正加载之前，设置一个接触步让两个面接触上来，在这个步骤里面，接触面的过盈小一点好，比如0.001.接下去再把作用与两个接触体的力及接触方向的自由度放开。
2.如果系统的载荷很多的话，将系统的载荷分做多步进行加载，一次性全上可能使系统无法在规定的迭代次数内收敛。所以根据需要分开，让abaqus的内核慢慢消化去。少吃多餐在这边好像也是成立的。
3.系统有多个接触的话，也最好如载荷一样，分成几个step让他们接触上。这样的做法会让你以后在模型的修改中更有方向性。
4.模型还是不收敛的话，你可以看一下是在哪一步或者那个inc不收敛。对于第一步直接不收敛的话，如果模型是像我上面把载荷和接触分成很多步建立的话，可以把载荷加载的顺序换一下。如果你把第二个加载的载荷换到第一步以后，计算收敛了，那影响收敛的主要问题应该就是原来第一个加载或者接触影响的。这种情况下面一般算到这个加载的时候还是不会收敛。这个时候可以考虑是否有什么其他办法能够使步骤的变化与上一步变动小一点，比如第一点里面提到，或者
继续把这个载荷细分呢？
5.对于接触分析不收敛的情况，可以自己看一下模型的接触面。有时候是overclosure,这个时候在assemble里面将模型相对位置稍微移动下或者用接触里面的那个adjustonly toremoveoverclose,不过这一种方法会使你的网格扭曲变形。问题不大也是可以用的。有的时候是因为，模型中的两个接触面变成了一个点和一个面接触，而点或者面中有一个位置并不是很稳定。这个时候就会出现了dividing,有时候求解无法成功。这时候可以看一下是不是能够将模型该处稍微改一下呢？或者将该处的网格细化一下。
6.模型实在是比较大的话，可以修改solver的设定，将迭代次数改大一点。对于开始计算就不收敛的，而在迭代次数到了以后时间增量还不是很小的话，可以将initial和minimum改小一点。模型越大的话这边可以改的越小，特别是前后两个step变化比较大的情况下。但对于模型不是很大的情况下，太小的时间增量是意义不大的，问题应该从模型当中是否有错误去考虑。
7. 模型太大的话会导致求解的方程太大，不需要的不重要的接触最好从模型当中去除。这样的话对结果影响也不会很大，而且可以是计算时间大大的减少。
***NOTE:THE SOLUTION APPEARS TO BE DIVERGING. CONVERGENCE IS JUDGED UNLIKELY.
求解出现分歧，一般是状态不稳定，如你的某一个物体的自由度忘了约束或着分析到了一种极限状态之类的，解决方案见上。
***ERROR: TIME INCREMENT REQUIRED IS LESS THAN THE MINIMUM SP ECIFIED 步长比最小步长小，方案见上。
***NOTE: SEVERE CONTACT OVERCLOSURES EXIST. CONVERGENCE IS J UDGED UNLIKELY. 接触出现问题，用visulazation的tool的job diagnostics查看，针对那个接触进行修改。把距离稍微改大一点点，应该是很有帮助的。如果是计算中出现问题的话，应该是你前后两个step设置的有问题。
近段时间，遇到abaqus的收敛问题，相当令人头痛。本以为大的问题解决了，小问题应该不在话下的，实际情况不是这样的。本人现在用到的是不同单元的连接与接触问题，同时还用到no compression功能及材料的非线性行为。这下可热闹了，就是不收敛？！通过与高手交流，加深了如下认识：
1.平稳地建立接触关系；
2.施加过渡约束，消除刚体位移；
3.如何查找产生刚体位移的节点？通过visualization 模块：tool →job diagnostics → warning,选中highlight selection in viewport,可以显示出现了numerical singularity的节点，一般来说，将此处的节点施加较小的位移约束，是可以解决刚体位移的困难的。
4.如何施加上述发生刚体位移的节点呢？在interaction模块,点击SPECIAL
→ SPRING/DASHOTS → CREATE，设置弹簧类型为CONECT POINTS TO GROND，选择一个节点，将degree of freedom设为出现了刚体位移的自由度，将spring stiffness的值设为一个小值（太小不足以约束刚体位移，太大会影响变形）
5.另外一个方法是通过visualization模块中的时程曲线，查找不收敛处的情况。也有利于找到问题所在。
接触问题技巧整理 1、塑性材料和接触面上都不能用C3D20R和C3D20单元，这可能是你收敛问题的主要原因。如果需要得到应力，可以使用C3D8I (在所关心的部位要让单元角度尽量接近90度)，如果只关心应变和位移，可以使用C3D8R, 几何形状复杂时，可以使用C3D10M. 2、接触对中的slave surface应该是材料较软，网格较细的面。 3、接触面之间有微小的距离，定义接触时要设定“Adjust=位置误差限度”，此误差限度要大于接触面之间的距离，否则ABAQUS会认为两个面没有接触：*Contact Pair, interaction="SOIL PILE SIDE CONTACT", small sliding, adjust=0.02 4、定义tie时也应该设定类似的position tolerance: *Tie, name=ShaftBottom, adjust=yes, position tolerance=0.1 5、msg文件中出现zero pivot说明ABAQUS无法自动解决过约束问题，例如在桩底部的最外一圈节点上即定义了tie，又定义了contact, 出现过约束。解决方法是在选择tie或contact 的slave surface时，将类型设为node region, 然后选择区域时不要包含这一圈节点（我附上的文件中没有做这样的修改）。 6、接触定义在哪个分析步取决于你模型的实际物理背景，如果从一开始两个面就是相接触的，就定义在initial或你的第一个分析步中；如果是后来才开始接触的，就定义在后面的分析步中。边界条件也是这样。 7、我在前面上传的文件里用*CONTROL设了允许的迭代次数18，意思是18次迭代不收敛时，才减小时间增量步（ABAQUS默认的值是12）。一般情况下不必设置此参数，如果在msg 文件中看到opening和closure的数目不断减小（即迭代的趋势是收敛的），但12次迭代仍不足以完全达到收敛，就可以用*CONTROL来增大允许的迭代次数。 8、桩头掉在了地表下，说明接触定义得不正确。可能接触面的距离还是大于*contact pair 中的adjust=0.02, 可改为adjust=0.2 9、原则上，90度的圆弧应该划分10个单元，适当少一些可能也行。 *contact pair中的adjust=0.005，还是太小，在后处理时可以看到，接触面之间的距离大于0.005。把adjust设置大一些没关系，比如adjust=0.1。 10、网格不好也可能产生过约束问题，不要只考虑边界条件啊！ 11、NUMERICAL SINGULARITY WHEN PROCESSING NODE 15294 D.O.F. 2 RATIO = 2.48305E+11", 说明NODE 15294 所在的实体在方向2上出现无限大的刚体位移。 可以在此实体上的任意一点和地面之间定义一个很软的弹簧，以消除刚体位移。方法是：interaction模块，菜单special / springs-dashpots / create, 选connect points to ground, 选节点，Degree of freedom 设为出现了刚体位移的自由度，spring stiffness为一个较小的值（太小则不足以消除刚体位移，太大则会影响变形）。 如果多个方向上出现了刚体位移，就要分别在相应的方向上各定义一个弹簧。spring 所在的节点在弹簧方向的位移乘以spring stiffness，就是弹簧所分担的载荷，它应该远远小于在此方向上的外载荷。如果模型位移很小，我常常把spring stiffness设为1.
整个计算过程包括2个分析步，第1步做屈曲分析，笫2步做极限强度分析。 第1步：屈曲分析 载荷步定义如下： Step 1-Initial Step 2- Buckle
? Re Mbs M^nce C^wvoini live 2oc*$ *l^*?4 tjdp V ：i.Jsa&# 录 +r A AJIu fffiC? fe3 Ha ? ABAQUS学习笔记 一．AQUS－.inp编码介绍 （一）．ABAQUS头信息文件段（1－4） 1.*PREPRINT 输出求解过程所要求的信息（在dat文件中） ie:*PREPRINT, ECHO=YES, HISTORY=YES, MODEL=YES 2.*HEADING 标题输出文件（出现在POST/VIEW窗口中，且出现在结果输出文件中） ie:*HEADING STRESS ANALYSIS FOR A PLATE WITH A HOLE 3.*RESTART 要求abaqus/standard输出其POST/view模块所需要的.res文件。其中的FREQ ＝？控制结果在每次迭代（或载荷步）输出的次数。 ie:*RESTART, WRITE, FREQ=1 4.*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 在已有结点的基础上进行多个结点的生成，一般是在两结点间以某种方式（直线、圆）产生一定分布规律的结点。 如：*NGEN, LINE=C, NSET=HOLE， 119, 1919, 100, 101 在两结点（结点号为119，1919）间以圆弧形式生成多个结点，100为任意相邻结点的单元号增量，101为圆弧形成时圆心位置的结点（对于直线形式生成没有此结点）。所有这些生成的结点（包括119，1919）被命名成HOLE的集合（这样做的目的是以后的命令中使用到它，比如说对这些结点施加同等条件的边界条件或载荷等，HOLE就是这些结点的代称）。*NGEN使用的前提就是必须存在已有结点。 *NGEN, NSET=OUTER 131, 1031, 100 以线形式形成结点，结点号增量100，结点集合名为OUTER。 *NGEN, NSET=OUTER 1031, 1931, 100 同上生成结点，可以同上结点集合名，这样OUTER就包括这两次生成的所有结点 3.*NFILL 在如上生成的结点集（实际上，代表两条几何意义上的边界线）之间按一定规律（BIAS＝？）填充结点。这样所有生成的结点构成一定形状的实体（面）。 如：*NFILL, NSET=PLATE, BIAS=0.8 HOLE, OUTER, 12, 1 以HOLE为第一条边界，OUTER为第二条边界（终止边），以从疏到密的规律（BIAS小于1）分布，其生成结点数在两内外对应结点间为12，1为每组结点号的增量。所有这些结点被置于PLATE的集合中。 下面以上面生成的结点来生成单元： 4.*ELEMENT 1. Abaqus/Explicit 中的接触形式 双击Interactions，出现接触形式定义。分为通用接触（General contact）、面面接触（Surface-to-Surface contact）和自接触（Self-contact）。 1. 通用接触General contact 通用接触用于为多组件，并具有复杂拓扑关系的模型建模。 General contact algorithm ?The contact domain spans multiple bodies (both rigid and deformable) ?Default domain is defined automatically via an all-inclusive element-based surface ?The method is geared toward models with multiple components and complex topology。 ?Greater ease in defining contact model 2. Surface-to-Surface contact Contact pair algorithm ?Requires user-specified pairing of individual surfaces ?Often results in more efficient analyses since contact surfaces are limited in scope 3. 自接触（Self-contact） 自接触应用于当部件发生变形时，可能导致自己的某两个或多个面发生接触的情况。如弹簧的压缩变形，橡胶条的压缩。 1.应用背景概述 随着科学技术的发展，汽车已经成为人们生活中必不可少的交通工具。但当今由于交通事故造成的损失日益剧增，研究汽车的碰撞安全性能，提高其耐撞性成为各国汽车行业研究的重要课题。目前国内外许多著名大学、研究机构以及汽车生产厂商都在大力研究节省成本的汽车安全检测方法，而汽车碰撞理论以及模拟技术随之迅速发展，其中运用有限元方法来研究车辆碰撞模拟得到了相当的重视。而本案例就是取材于汽车碰撞模拟分析中的一个小案例―――保险杠撞击刚性墙。 2.问题描述 该案例选取的几何模型是通过导入已有的*.IGS文件来生成的（已经通过Solidworks软件建好模型的），共包括刚性墙（PART-wall）、保险杠（PART-bumper）、平板（PART-plane）以及横梁（PART-rail）四个部件，该分析案例的关注要点就是主要吸能部件（保险杠）的变形模拟，即发生车体碰撞时其是否能够对车体有足够的保护能力？这里根据具体车体模型建立了保险杠撞击刚性墙的有限元分析模型，为了节省计算资源和时间成本这里也对保险杠的对称模型进行了简化，详细的撞击模型请参照图1所示，撞击时保险杠分析模型以2000mm/s的速度撞击刚性墙，其中分析模型中的保险杠与平板之间、平板与横梁之间不定义接触，采用焊接进行连接，对于保险杠和刚性墙之间的接触采用接触对算法来定义。 1.横梁（rail） 2.平板（plane） 3.保险杠（bumper） 4.刚性墙（wall） 图2.1 碰撞模型的SolidWorks图 为了使模拟结果尽可能真实，通过查阅相关资料，定义了在碰撞过程中相关的数据以及各部件的材料属性。其中，刚性墙的材料密度为7.83×10-9，弹性模量为2.07×105，泊松比为0.28；保险杠、平板以及横梁的材料密度为7.83×10-9，弹性模量为2.07×105，泊松比为0.28，塑形应力-应变数据如表2.1所示。 表2.1 应力-应变数据表 应力210 300 314 325 390 438 505 527 应变0.0000 0.0309 0.0409 0.0500 0.1510 0.3010 0.7010 0.9010 注：本例中的单位制为：ton，mm，s。 3.案例详细求解过程 本案例使用软件为版本为abaqus6.12，各详细截图及分析以该版本为准。3.1 创建部件 （1）启动ABAQUS/CAE，创建一个新的模型数据库，重命名为The crash simulation，保存模型为The crash simulation.cae。 （2）通过导入已有的*.IGS文件来创建各个部件，在主菜单中执行【File】→【Import】→【Part】命令，选择刚刚创建保存的的bumper_asm.igs文件，弹出【Create Part From IGS File】对话框如图3.1所示，根据图3.1所示设定【Repair Options】的相关选项，其它参数默认，单击【Ok】按钮，可以看到在模型树中显示了导入的部件bumper_asm。 图3.1 Create Part From IGS File对话框 部件模态综合法 随着科学和生产的发展，特别是航空、航天事业的发展，越来越多的大型复杂结构被采用，这使得建模和求解都比较困难。一方面，一个复杂结构势必引入较多的自由度，形成高维的动力学方程，使一般的计算机在内存和求解速度方面都难以胜任，更何况一般的工程问题主要关心的是较低阶的模态。仅为了获取少数的几个模态，必须为求解高维方程付出巨大的代价也是不合适的。另一方面，正是由于结构的庞大和复杂，一个完整的结构往往不是在同一地区生产完成的，可能一个结构的各个主要零部件不得不由不同的地区、不同的厂家生产。而且由于试验条件的限制只能进行部件的模态实验，而无法对整体结构进行模态实验。针对这些主要的问题，为了获得大型、复杂结构的整体模态参数，于是发展了部件模态综合法。 部件模态综合法又叫子结构耦合法。它的基本思想是按工程观点或结构的几何轮廓，并遵循某些原则要求，把完整的结构进行人为抽象肢解成若干个子结构（或部件）；首先对子结构（或部件）进行模态分析，然后经由各种方案，把它们的主要模态信息（常为低阶主模态信息）予以保留，并借以综合完整结构的主要模态特征。它的主要有点是，可以通过求解若干小尺寸结构的特征问题来代替直接求解大型特征值问题。同时对各个子结构可分别使用各种适宜的数学模型和计算程序，也可以借助试验的方法来获得他们的主要模态信息。 对于自由振动方程在数学上讲就是固有（特征）值方程。特征值方程的解不仅给出了特征值，即结构的自振频率和特征矢量——振兴或模态，而且还能使结构在动力载荷作用下的运动方程解耦，即所谓的振型分解法或叫振型叠加法。因此，特征值问题的求解技术，对于解决结构振动问题来说吧，是非常重要的。 考虑阻尼的振型叠加法 振型叠加法的定义：将结构各阶振型作为广义坐标系，求出对应于各阶振动的结构内力和位移，经叠加后确定结构总响应的方法。 振型叠加法的使用条件： ?（1）系统应该是线性的：线性材料特性，无接触条件，无非线性几何效应。 ?（2）响应应该只受较少的频率支配。当响应中各频率成分增加时，例如撞击和冲击问题，振型叠加技术的有效性将大大降低。 ?（3）载荷的主要频率应在所提取的频率范围内，以确保对载荷的描述足够精确。 ?（4）由于任何突然加载所产生的初始加速度应该能用特征模态精确描述。 ?（5）系统的阻尼不能过大。 看本例题之前，请务必先找着文献[1]中P75——P101中提供的例题完全照做一遍，以熟悉基本的操作流程。 下面是本例题的操作过程，模拟一片砌体墙片的滞回实验。 第一步：模型部件的建立 进入ABAQUS（中文版），在左方菜单中，选择“部件”，鼠标右键点击一下，选择“创建”，进入模型的部件创建菜单。模型中往往有一个或者多个部件构成。如图1，设定部件名称，其他选项如图1所示。 图1 然后中间的主操作界面出现，如图2。 此平面默认的在平行于计算机显示器的方向为XY平面，我们如图2，按照现实中墙体立面的尺寸画出墙体立面框，然后点击“完成”，弹出的菜单如图3. 图2 图3 此处的“深度”一项就是设定墙体部件的厚度，输入0.24（墙体厚度0.24m），点击确定。于是得到了墙体部件的基本视图如图4所示。 同理，依样设定加载梁的尺寸，得到加载梁部件。 这样，第一步部件尺寸设定就完成了。 图4 第二步：部件使用材料的设定 加载梁使用c50混凝土，砌体使用与实验相对应的材料参数。 由于模拟是针对砌体，所以不考虑加载梁的塑性，因此加载梁只设定密度和弹性。 而砌体则以混凝土塑性损伤本构模型来模拟，要设定密度、弹性、混凝土损伤塑性。 以上内容中混凝土材料参数的设定参见资料[2]，砌体材料参数的设定参见资料[3]。 第三步：将材料属性赋予模型 设定了材料参数后，还要对将材料参数“赋予”模型。其操作菜单如图5 图5 首先建立两个界面SECTION-1和SECTION-2如图5左边红框所示，将两种材料（C50混凝 土与砌体材料）“注入”SECTION-1和SECTION-2中，然后点击右边红框中的图标，选择截面所要“赋予”的对象，即可完成材料参数赋予模型的操作。 第四步：安装配件 在第一步谁定好了部件后，实际上部件就像积木玩具的各个零件一样还是零散分开的，这时候，就要使用装配件功能如图6所示。 图6 用得较多的是红框中两个移动功能，非常简单，就是通过在空间坐标系中将部件平移到正确的位置即可。 第五步：设置分析步 在模拟滞回曲线实验的拟静力计算中，因为需要循环加载时反复推拉形成的，所以在模拟的时候需要设定多个分析步。而单纯的静力计算则只需一步就可以完成。 分析步的设定如图7所示，本次模拟一共设定了12个分析步，即表示反复推拉的模拟次数一共12次。 在设定分析步的过程中，将分析步的属性定为“静力、通用”。除了自定义的12步以外，还有ABAQUS默认的不能删除的初始步，是调整结构受力初始状态用的。 第六步：设定约束 这一步实际上是第四步的后续，第四步中我们只是把“积木玩具”完全进行了搭接，而没有进行固定。这一步的作用就是确定“积木玩具”各个分块之间相互作用的属性与关系。 对于实验来说，加载梁是牢牢固定在气体墙之上的，因此新建一个约束，属性选择“绑定”，将加载梁与砌体牢牢固定在一起。 ABAQUS时程分析法计算地震反应得简单实例ABAQUS时程分析法计算地震反应得简单实例（在原反应谱模型上 修改） 问题描述： 悬臂柱高12m，工字型截面（图1），密度7800kg/m3，EX=2、1e11Pa，泊松比0、3，所有振型得阻尼比为2%，在3m高处有一集中质量160kg，在6m、9m、12m处分别有120kg 得集中质量。反应谱按7度多遇地震，取地震影响系数为0、08，第一组，III类场地，卓越周期Tg=0、45s。 图1 计算对象 第一部分：反应谱法 几点说明： λ本例建模过程使用CAE； λ添加反应谱必须在inp中加关键词实现，CAE不支持反应谱； λ *Spectrum不可以在keyword editor中添加，keyword editor不支持此关键词读入。 λ ABAQUS得反应谱法计算过程以及后处理要比ANSYS方便得多。 操作过程为： （1）打开ABAQUS/CAE，点击create model database。 （2）进入Part模块，点击create part，命名为column，3D、deformation、wire。continue （3）Create lines，在 分别输入0，0回车；0，3回车；0，6回车；0，9回车；0，12回车。 （4）进入property模块，create material，name：steel，general-->>density，mass density：7800 mechanical-->>elasticity-->>elastic，young‘s modulus：2、1e11，poisson’s ratio：0、3、 CAE(计算机辅助工程)是一门复杂的工程科学，涉及仿真技术、软件、产品设计和力学等众多领域。世界上几大CAE公司各自以其独到的技术占领着相应的市场。ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群，是国际上公认的最先进的大型通用非线性有限元分析软件之一。它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域。ABAQUS在技术、品质和可靠性等方面具有卓越的声誉，可以对工程中各种复杂的线性和非线性问题进行分析计算。 《ABAQUS有限元分析常见问题解答》以问答的形式，详细介绍了使用ABAQUS建模分析过程中的各种常见问题，并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法，帮助读者提高解决问题的能力。 《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。 16.1.1点对面离散与面对面离散 【常见问题16-1】 在ABAQUS/Standard分析中定义接触时，可以选择点对面离散方法(node-to-surfac e-dis-cre-tization)和面对面离散方法(surface-to-surfacediscretization)，二者有何差别? 『解答』 在点对面离散方法中，从面(slavesurface)上的每个节点与该节点在主面(mastersur face)上的投影点建立接触关系，每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点。 使用点对面离散方法时，从面节点不会穿透(penetrate)主面，但是主面节点可以穿透从面。 面对面离散方法会为整个从面(而不是单个节点)建立接触条件，在接触分析过程中同时考虑主面和从面的形状变化。可能在某些节点上出现穿透现象，但是穿透的程度不会很严重。 在如图16-l和图16-2所示的实例中，比较了两种情况。 Abaqus 使用问答 Q:abaqus的图形如何copy? A:file>print>file格式为png，可以用Acdsee打开。 Q:用Abaqus能否计算[Dep]不对称的问题? A:可以，并且在step里面的edit step对话框other里面的matrix solver有个选项。Q: 弹塑性矩阵【D】与ddsdde有何联系? A: stress=D*stran；d(stress)=ddsdde*d(stran)。 Q:在abaqus中，如果采用umat,利用自己的本构，如何让abaqus明白这种材料的弹塑性应变，也就是说，如何让程序返回弹性应变与塑性应变，好在output中输出，我曾想用最笨地方法，在uvarm中定义输出，利用getvrm获取材料点的值，但无法获取增量应力，材料常数等，研究了帮助中的例子，umatmst3.inp,umatmst3.for,他采用mises J2 流动理论，我在output history 显示他已进入塑性状态，但他的PE仍然为0！！? A: 用uvar( )勉强成功 。 Q: 本人在用umat作本构模型时， *static, 1,500,0.000001,0.1 此时要求的增量步很多，即每次增量要很小， *static 1,500 时，在弹性向塑性过度时，出现错误，增量过大，出现尖点.? A: YOU CAN TRY AS FOLLOWS: *STEP,EXTRAPOLATION=NO,INC=2000000 *STATIC 0.001,500.0,0.00001,0.1。 Q: 模型中存在两个物体的接触，计算过程中报错，怎么回事? A: 接触问题不收敛有两个方面不妨试试： 一、在*CONTACT PAIR 里调试ADJUST参数； 二、调一些模型参数，比如FRICTION等。 。 Q: 在边界条件和加载时，总是有initial这个步，然后是我们自己定义的加载步，请问这个initial步,主要作用是什么？能不能去掉？ A: 不能去掉,所有的分析都有,是默认的步。 Q: A solid extrusion base feature 这句话是什么意思？ extrusion、revolution等是什么意思？ A: 这两的是三维建模时候,在画完二维图形,如何来生成三维图形, extrusion意思是你给定一个厚度,然后二维图形第三个方向上面伸展这么多形成三维图形 revolution意思是你给定一个旋转轴,二维图形绕其旋转后形成三维轴对称图形 。 《现代机械设计方法》课程结业论文 （ 2011 级） 题目：ABAQUS实例分析 学生姓名 XXXX 学号 XXXXX 专业机械工程 学院名称机电工程与自动化学院 指导老师 XX 2013年 5 月8 日 目录 第一章Abaqus简介 (1) 一、Abaqus总体介绍 (1) 二、Abaqus基本使用方法 (2) 1.2.1 Abaqus分析步骤 (2) 1.2.2 Abaqus/CAE界面 (3) 1.2.3 Abaqus/CAE的功能模块 (3) 第二章基于Abaqus的通孔端盖分析实例 (4) 一、工作任务的明确 (6) 二、具体步骤 (6) 2.2.1 启动Abaqus/CAE (4) 2.2.2 导入零件 (5) 2.2.3 创建材料和截面属性 (6) 2.2.4 定义装配件 (7) 2.2.5 定义接触和绑定约束（tie） (10) 2.2.6 定义分析步 (14) 2.2.7 划分网格 (15) 2.2.8 施加载荷 (19) 2.2.9 定义边界条件 (20) 2.2.10 提交分析作业 (21) 2.2.11 后处理 (22) 第三章课程学习心得与作业体会 (23) 第一章： Abaqus简介 一、Abaqus总体介绍 Abaqus是功能强大的有限元分析软件，可以分析复杂的固体力学和结构力学系统，模拟非常庞大的模型，处理高度非线性问题。Abaqus不但可以做单一零件的力学和多物理场的分析，同时还可以完成系统级的分析和研究。 Abaqus使用起来十分简便，可以很容易的为复杂问题建立模型。Abaqus具备十分丰富的单元库，可以模拟任意几何形状，其丰富的材料模型库可以模拟大多数典型工程材料的性能，包括金属、橡胶、聚合物、复合材料、钢筋混泥土、可压缩的弹性泡沫以及地质材料（例如土壤、岩石）等。 Abaqus主要具有以下分析功能： 1.静态应力/位移分析 2.动态分析 3.非线性动态应力/位移分析 4.粘弹性/粘塑性响应分析 5.热传导分析 6.退火成形过程分析 7.质量扩散分析 8.准静态分析 9.耦合分析 10.海洋工程结构分析 11.瞬态温度/位移耦合分析 12.疲劳分析 13.水下冲击分析 14.设计灵敏度分析 二、Abaqus基本使用方法 1.2.1 Abaqus分析步骤 有限元分析包括以下三个步骤： 1.前处理（Abaqus/CAE）:在前期处理阶段需要定义物理问题的模型，并生 成一个Abaqus输入文件。提交给Abaqus/Standard或 Abaqus/Explicit。 2.分析计算（Abaqus/Standard或Abaqus/Explicit）:在分析计算阶段， 使用Abaqus/Standard或Abaqus/Explicit求解输入文件中所定义的 1、创建部件： Step1：执行Part/Create命令，或者单击左侧工具箱区域中的（create part）按钮，弹出如图1-1所示的Create Part对话框。在Name（部件名称）后面输入foundation，将Modeling Space（模型所在空间）设为2D Planar（二维平面），Type （类型）设为Deformable（可变形体），Base Feature（基本特征）设为Shell（壳）。单击Continue按钮退出Create Part对话框。ABAQUS/CAE自动进入绘图（Sketcher）环境。 图1-1 Step2：选择绘图工具框右上方的创建矩形工具，在窗口底部的提示区 显示“Pick a starting corner for the rectangle—or enter X,Y”，输入坐标（0,0），按下Enter键，在窗口底部的提示区显示“Pick the opposite corner for the rectangle—or enter X,Y”，输入（45.5,20），按下Enter键。单击Done，创建part 完成，如图1-2。 图1-2 Step3：单击左侧工具箱区域中的，弹出如图1-3的窗口。应用或 功能将groundwork（基础）在foundation的位置绘制出来，点击Done，返回图1-4所示窗口 图1-3 图1-4 Step4：执行Tools-Set-Create弹出如图1-5的Create Set对话框，在Name后 面输入all，点击Continue，将整个foundation模块选中如图1-6所示，点击Done，完成集合all的创建。以相同的操作，将图1-4中的小矩形区域创建Name为remove 的集合。 图1-5 图1-6 以相同的方式分别创建名称为：groundwork，retaining，backfill的part,依次如图1-7,1-8,1-9所示。并分别创建于part名称相同的集合。 图1-7 Abaqus基本操作中文教程 目录 1 Abaqus 软件基本操作 .................... 常用的快捷键 .......................... 单位的一致性 .......................... 分析流程九步走 ....................... 几何建模（Part） ..................... 属性设置（Property） ................... 建立装配体（Assembly） ................... 定义分析步（Step） ................... 相互作用（In teracti on................ ） 载荷边界（Load） ..................... 划分网格（Mesh） .................. 作业（Job） ...................... 可视化（Visualization ）................. 1 Abaqus软件基本操作 常用的快捷键 「旋转模型一Ctrl+Alt+ 鼠标左键 于平移模型一Ctrl+Alt+鼠标中键 " 缩放模型一Ctrl+Alt+ 鼠标右键 单位的一致性 CAE软件其实是数值计算软件，没有单位的概念，常用的国际单位制如下表1所示，建议采用SI （mm）进行建模。 国际单位制 SI (m) SI (mm) 「长度 m mm 力 N N 质量 kg t 时间 s s 应力 2 Pa (N/m ) 2 MPa (N/mm) 质量密度 kg/m 3 3 t/mm 加速度 m/s 2 mm/s 例如，模型的材料为钢材，采用国际单位制 SI （m ）时，弹性模量为 m,重力加速度m/s 2 ，密度为7850 kg/m 3,应力Pa;采用国际单位制SI （mm ） 时，弹性模量为 口金 重力加速度 9800 mm/s 2 ，密度为7850e-12??T/mm 5, 应力MPa 分析流程九步走 几何建模（Part 属性设置（Property ） 建立装配体（Assembly ） T 定义分析步（Step ） T 相互作用 （Interaction ）宀载荷边界（Load ） T 划分网格 （Mesh ）T 作业（Job ）T 可视化（Visualization ） ' 以上给出的是软件 ! 常规的建模和分析的流 程，用户可以根据自己 ；的建模习惯进行调整。 I 另外，草图模块可以进 !行参数化建模，建议用 」户可以参考相关资料进--- 几何建模（Part ） 关键步骤的介绍： 部件（Part ）导入 Pro/E 等CAD 软件建好的模型后，另存成 iges 、sat 、step 等格式； 然后导入Abaqus 可以直接用，实体模型的导入通常采用 sat 格式文件导 謝t fti5 忧化 fkit 可泯忧 一．创建部件 1.打开abaqus; 开始/程序/Abaqus6.10-1/Abaque CAE 2.Model/Rename/Model-1,并输入名字link4 3.单击Create part弹出Create part对话框， Name输入link-4; Modeling Space 选择2D Planar Type 选择Deformable Base Feature 选择Wire Approximate size 输入800；然后单击continue 4.单击(Create Lines:connected)通过点(0,0)、(400,0)、(400,300)、(0,300)单击(Create Lines:connected)连接(400,300)和(0,0)两点，单击提示区中的Done按钮(或者单击鼠标滚轮，也叫中键)，形成四杆桁架结构 5.单击工具栏中的(Save Model Database),保存模型为link4.cae 二.定义材料属性 6.双击模型树中的Materials(或者将Module切换到Property,单击Create Material -ε) 弹出Edit Material对话框后。 执行对话框中Mechanical/Elasticity/Elastic命令， 在对话框底部出现的Data栏中输入Young’s Module为29.5e4， 单击OK.完成材料设定。 7.单击“Create Section ”,弹出Create Section对话框， Category中选择Beam; Type中选择Truss; 单击continue按钮 弹出Edit Section对话框， 材料选择默认的Material-1,输入截面积(Cross-sectional area)为100，单击ok按钮。 1、塑性材料和接触面上都不能用C3D20R和C3D20单元，这可能是你收敛问题的主要原因。如果需要得到应力，可以使用C3D8I (在所关心的部位要让单元角度尽量接近90度)，如果只关心应变和位移，可以使用C3D8R, 几何形状复杂时，可以使用C3D10M。 2、接触对中的slave surface应该是材料较软，网格较细的面。 3、接触面之间有微小的距离，定义接触时要设定“Adjust=位置误差限度”，此误差限度要大于接触面之间的距离，否则ABAQUS会认为两个面没有接触：*Contact Pair, interaction="SOIL PILE SIDE CONTACT", small sliding, adjust=0.2. 4、定义tie时也应该设定类似的position tolerance: *Tie, name=ShaftBottom, adjust=yes, position tolerance=0.1 5、msg文件中出现zero pivot说明ABAQUS无法自动解决过约束问题，例如在桩底部的最外一圈节点上即定义了tie，又定义了contact, 出现过约束。解决方法是在选择tie或contact的slave surface时，将类型设为node region, 然后选择区域时不要包含这一圈节点（我附上的文件中没有做这样的修改）。 6、接触定义在哪个分析步取决于你模型的实际物理背景，如果从一开始两个面就是相接触的，就定义在initial或你的第一个分析步中；如果是后来才开始接触的，就定义在后面的分析步中。边界条件也是这样。 7、我在前面上传的文件里用*CONTROL设了允许的迭代次数18，意思是18次迭代不收敛时，才减小时间增量步（ABAQUS默认的值是12）。一般情况下不必设置此参数，如果在msg文件中看到opening和closure的数目不断减小（即迭代的趋势是收敛的），但12次迭代仍不足以完全达到收敛，就可以用*CONTROL来增大允许的迭代次数。 8、桩头掉在了地表下，说明接触定义得不正确。可能接触面的距离还是大于*contact pair 中的adjust=0.02, 可改为adjust=0.2 9、原则上，90度的圆弧应该划分10个单元，适当少一些可能也行。 *contact pair中的adjust=0.005，还是太小，在后处理时可以看到，接触面之间的距离大于0.005。把adjust设置大一些没关系，比如 adjust=0.1。 10、网格不好也可能产生过约束问题，不要只考虑边界条件啊！NUMERICAL SINGULARITY WHEN PROCESSING NODE 15294 打开abaqus，然后点击file——set work directory，然后选择指定文件夹，开始建模，建模完成后及时保存，在进行运算以前对已经完成的工作保存，然后点击job，修改inp文件的名称进行运算。切记切 记！！！！！！ 1、如何显示梁截面（如何显示三维梁模型） 显示梁截面：view->assembly display option->render beam profiles，自己调节系数。 2、建立几何模型草绘sketch的时候，发现画布尺寸太小了 1）这个在create part的时候就有approximate size，你可以定义合适的（比你的定性尺寸大一倍）； 2）如果你已经在sketch了，可以在edit菜单－－sketch option ——general－－grid更改 3、如何更改草图精度 可以在edit菜单－－sketch option ——dimensions－－display——decimal更改 如果想调整草图网格的疏密，可以在edit菜单－－sketch option ——general——grid spacing中可以修改。 4、想输出几何模型 part步，file，outport－－part 5、想导入几何模型？ part步，file，import－－part 6、如何定义局部坐标系 Tool－Create Datum－CSYS－－建立坐标系方式－－选择直角坐标系or柱坐标系or球坐标 7、如何在局部坐标系定义载荷 laod－－Edit load－－CSYS－Edit（在BC中同理）选用你定义的局部坐标系 8、怎么知道模型单元数目（一共有多少个单元） 在mesh步，mesh verify可以查到单元类型，数目以及单元质量一目了然，可以在下面的命令行中查看单元数。 Query---element 也可以查询的。 9、想隐藏一些part以便更清楚的看见其他part，edge等 view－Assembly Display Options——instance，打勾 10、想打印或者保存图片 File——print——file——TIFF——OK 11、如何更改CAE界面默认颜色 view->Grahphic options->viewport Background->Solid->choose the wite colour! 然后在file->save options. 12、如何施加静水压力hydrostatic load --> Pressure, 把默认的uniform 改为hydrostatic。这个仅用于standard，显式分析不支持。 13、如何检查壳单元法向 Property module/Assign/normal 14、如何输出单元体积 set步---whole model ----volume/Tickness/Corrdinate-----EVOL 15、如何显示最大、最小应力 在Visualization>Options>contour >Limits中选中Min/Max:Show Location，同样的方法可以知道具体指定值的位置。 16、如何在Visualization中显示边界条件 View——ODB display option——entity display——show boundary conditions 目录 第一章Abaqus简介 (1) 一、Abaqus总体介绍 (1) 二、Abaqus基本使用方法 (2) 1.2.1 Abaqus分析步骤 (2) 1.2.2 Abaqus/CAE界面 (3) 1.2.3 Abaqus/CAE的功能模块 (3) 第二章基于Abaqus的通孔端盖分析实例 (4) 一、工作任务的明确 (5) 二、具体步骤 (5) 2.2.1 启动Abaqus/CAE (4) 2.2.2 导入零件 (5) 2.2.3 创建材料和截面属性 (6) 2.2.4 定义装配件 (7) 2.2.5 定义接触和绑定约束（tie） (10) 2.2.6 定义分析步 (14) 2.2.7 划分网格 (15) 2.2.8 施加载荷 (19) 2.2.9 定义边界条件 (20) 2.2.10 提交分析作业 (21) 2.2.11 后处理 (22) 第三章课程学习心得与作业体会 (22) 第一章： Abaqus简介 一、Abaqus总体介绍 Abaqus是功能强大的有限元分析软件，可以分析复杂的固体力学和结构力学系统，模拟非常庞大的模型，处理高度非线性问题。Abaqus不但可以做单一零件的力学和多物理场的分析，同时还可以完成系统级的分析和研究。 Abaqus使用起来十分简便，可以很容易的为复杂问题建立模型。Abaqus具备十分丰富的单元库，可以模拟任意几何形状，其丰富的材料模型库可以模拟大多数典型工程材料的性能，包括金属、橡胶、聚合物、复合材料、钢筋混泥土、可压缩的弹性泡沫以及地质材料（例如土壤、岩石）等。 Abaqus主要具有以下分析功能： 1.静态应力/位移分析 2.动态分析 3.非线性动态应力/位移分析 4.粘弹性/粘塑性响应分析 5.热传导分析 6.退火成形过程分析 7.质量扩散分析 8.准静态分析 9.耦合分析 10.海洋工程结构分析 11.瞬态温度/位移耦合分析 12.疲劳分析 13.水下冲击分析 14.设计灵敏度分析 二、Abaqus基本使用方法 1.2.1 Abaqus分析步骤 有限元分析包括以下三个步骤： 1.前处理（Abaqus/CAE）:在前期处理阶段需要定义物理问题的模型，并生ABAQUS使用时的点点滴滴
Abaqus Explicit 接触问题
(完整word版)abaqus6.12-典型实例分析
abaqus接触动力学分析
ABAQUS使用例题
ABAQUS时程分析实例
abaqus接触分析的常见问题
ABAQUS使用解答(整理)
(完整word版)ABAQUS实例分析
(完整版)Abaqus操作说明
Abaqus基本操作中文教程
abaqus实例
ABAQUS接触分析
总结Abaqus操作技巧总结
ABAQUS实例分析论文
相关主题
文本预览