Abaqus 建模流程
Abaqus标准版共有“部件（part）”、“材料特性（propoterty）”、“装配（assemble）”、“计算步骤（step）”、“交互（interaction）”、“加载（load）”、“单元划分（mesh）”、“计算（job）”、“后处理（visualization）”、“草图（sketch）”十大模块组成。 建模方法：
1 首先建立“部件”
（1）根据实际模型的尺寸决定绘图区的大小，一般为模型的1.5倍，间距大小可以在edit菜单sketcher options选项里调整。
（2）在绘图区分别建立部件中的各个特征体，建立特征体的方法主要有挤压、旋转、平扫三种。同一个模型中两个不同的部件可以有同名的特征体组成，也就是说不同部件中可以有同名的特征体，同名特征体可以相同也可以不同。部件的特征体包括用各种方法建立的基本特征体、数据点（datum point）、数据轴（datum axis）、数据平面（datum plane）等等。选择多个元素时，可以同时按住shift键，或者按住鼠标左键进行窗选；如果取消对某个元素的选择可以同时按住ctrl键。同时按住ctrl、shift和鼠标左键（中键、右键）然后平移鼠标可以进行旋转（平移、缩放）。如果想修改或撤销已经完成的操作，可以在窗口左侧的模
型树中找到此项操作，在上面点击右键，选择Edit或delete。
（3）编辑部件可以用部件管理器进行部件复制，重命名，删除等，部件中的特征体可以是直接建立的特征体，还可以间接手段建立，如首先建立一个数据点特征体，通过数据点建立数据轴特征体，然后建立数据平面特征体，再由此基础上建立某一特征体，最先建立的数据点特征体就是父特征体，依次往下分别为子特征体，删除或隐藏父特征体其下级所有子特征体都将被删除或隐藏。 （4）部件类型：
? 可变形体：任意形状的，可以包含不同维数的特征（实体、表面、线）；在荷载作用下可以变形。 ? 不连续介质刚体：任意形状的；在荷载作用下不可变形。
? 解析刚体：只可以用直线、圆弧和抛物线创建的形状；在荷载作用下不可变形 ? 欧拉部件：实体区域；定义在欧拉分析中材料可以流动的区域
刚体是不能够施加质量、惯性轴等特性的，建立刚体后必须给刚体指定一个参考点（reference point），在加载模块里对参考点施加约束和定义其运动，对参考点施加的荷载或运动就相当于施加给了整个刚体。除了刚形体有旋转的情况或者要求绕刚体中的某一轴的反力矩情况外，参考点的位置并不重要，上述两种情况，参考点应该位于绕其转动的轴上。在创建部件时需要指定部件的类型，一旦建立后就不能更改其类型。对于形状简单的刚性部件，使用解析刚体可以精确模拟部件的几何形状，而且可以减小计算代价，但如果刚性部件的几何形状较复杂，无法用解析刚体来建模，就需要使用离散刚体。
解析刚体不需要画网格，离散刚体需要画网格（边界由网格节点控制），且要在发生接触的部位划分足够细的网格，以保证不出现大的尖角。
创建刚体的三种方法：1.离散刚体和解析刚体；2.Interaction模块中的刚体约束和显示体约束，可以将变形体变为刚体；3.定义一块钢板，其属性定义弹模无限大、泊松比无限小，可以模拟刚体。
（5）分区将部件再细分为不同的区域，区域可以用于创建几何集，还可以用于划分网格，一般在Assembly和Mesh模块创建可划分网格的分区效果更好。
（6）在修改部件几何形状时，尽量修改顶点位置或编辑尺寸，而不要创建或删除线段，这样可以减少对已定义的部件特征、集合和面的影响。在修改几何模型后，必须对原模型的截面属性、面、集合、载荷、边
界条件和约束进行全面检查，以便确定原模型是否受到影响。
（7）在创建轴对称部件时，ABAQUS/CAE要求旋转轴必须是竖直方向的辅助线，而且轴对称部件的整个平面图都要位于旋转轴的右侧。
（8）ABAQUS/CAE 推荐的建模方法是把整个数值模型（如材料、边界条件、荷载等）都直接定义在几何模型上，而不是像其他前处理器那样定义在单元和节点上，这样在修改网格时不必重新定义材料和边界条件等模型参数。在处理复杂问题时，可以先简单地划分粗网格，得到初步的模拟结果，然后再在适当的区域细化网格。
一般先划分网格，这样做的好处是，往往在划分网格的过程中，会发现部件的几何模型需要进一步修改，例如存在过小的圆角或线段，导致不必要的细化网格；而经过这些修改后，已经定义好的边界条件、载荷和接触等可能变为无效的，需要再重新定义。
（9）利用Sketch模块创建独立的草图。该种方法创建的草图不与任何的部件相关联，可以保留，作后继使用。Sketch 约束定义了几何实体之间的逻辑关系，如平行、垂直、切线、一致、同心等等。
2 建立材料特性
（1）输入材料特性参数如弹性模量、泊松比等
大多数实验数据常常是用名义应力和名义应变的值给出的。这时，必须把塑性材料的数据从名义应力/应变的值转换为真实应力/应变的值。当应变很小时，真实值和名义值之间差别很小，而当应变很大时，二者之间就会有明显的差别；因此，如果模拟的应变比较大，就一定要向ABAQUS提供合适的应力－应变数据，这是极为重要的。对一般多维应力状态，用屈服准则确定应力属于弹性还是弹塑性范围。
定义截面属性时，平面应力单元、平面应变单元和轴对称单元都应该定义为实体截面属性(*SOLID SECTION)，而不是壳截面属性(*SHELL SECTION) 。
在进行弹塑性分析时，同样可以使用分区的方法，将部件中重要的、塑性变形较大的区域定义为弹塑性材料，将不重要的、几乎不发生塑性变形的区域定义为弹性材料，以便使分析更容易收敛，缩短计算时间。
尽量不要对塑性材料施加点载荷，而是根据实际情况来使用面载荷或线载荷。如果必须在某个节点上施加点载荷，可以使用耦合约束来为载荷作用点附近的几个节点建立刚性连接，这样这些节点就会共同承担点载荷。
材料方向：对于壳、梁和桁架单元，局部的材料方向总是随着变形而转动。对于实体单元，仅当单元中提供了非默认的局部材料方向时，它的局部材料方向才随着变形而转动，否则，默认的局部材料方向在整个分析中将始终保持不变。
（2）建立截面（section）特性，如均质的、各项同性、平面应力平面应变等等，截面特性管理器依赖于材料参数管理器
（3）分配截面特性给特征体，把截面特性分配给部件的某一区域就表示该区域已经和该截面特性相关联
3 模型装配
在装配（assemble）模块里首先建立部件实例（part instance），一个部件实例可以看作部件的代表，但并不是原部件的拷贝。实例一直和原部件保持关联，当原部件几何形状发生变化时，实例也发生相应变化。一个装配模型可以包含一个部件的多个实例，在创建第一个实例时所生成的装配模型总体坐标系是该装配模型的一个实例。
同一个部件中所有特征体在装配模块中对该部件建立实例时会形成一个整体，选择该实例时，该实例在装配之前原部件中所有特征体都被选择了。后续所有模块的操作对象就是所生成的部件实例，也即装配模型中的特征体，而不是原来的部件。
对于各部件的实例，可以在view菜单assembly display options选项里选择instance标签对现有的各实例决定其是否显示在当前视窗中，这一功能对选择视窗中的对象很有帮助。
（1）部件实例有独立的和非独立的两种，缺省状态是非独立实例。
（2）在交互模块、加载模块和单元划分模块里操作的对象都是装配模型中各个部件实例。
（3）创建了一个部件实例后，ABAQUS需要生成一个装配体的总体坐标系定位该实例，该装配体的总体坐标系与部件的总体坐标系是两个不同的坐标系。创建部件基特征体时的绘图（sketch）坐标原点与装配体的总体坐标系原点重合，并且xy坐标平面和装配体总体坐标系xy平面平行。创建了第一个实例后，ABAQUS定位该实例的方法就是将该实例基特征体的坐标原点（绘制平面草图的坐标原点）与装配体总体坐标系原点重合。
（4）定位各个部件实例
常见的定位标准包括：平行面、面对面、平行边、边对边、共轴、点重合、坐标系平行、接触。各定位标准之间互不影响，可以用新的定位标准替换原定位标准。箭头指向相同的方向。 每一个定位标准都作为装配模型的特征体而保存，可以在特征体管理器里进行编辑。 （5）集和面
如果当前的功能模块是Assembly、Interaction、Load或Mesh（处在为装配件划分网格的状态下），则使用主菜单Tools定义的面或集合是属于整个装配间的；而如果当前的功能模块式Part或Mesh（处在为部件划分网格的状态下），则使用主菜单Tools定义的面或集合只是属于此部件，不能在Assembly、Interaction 或Load 功能模块中使用。因此，创建集合或面时，要注意首先选择正确的功能模块（恰当的做法是在需要的模块中建立集和面）。
在定义约束、边界条件、载荷、接触或场变量等模型参数时，都应事先定义相应的集合和面，并给出容易识别的名称，这样在建立复杂模型时，会大大降低出错的可能性。
4 设置分析步（step）
（1）对模型施加荷载和边界条件之前或者定义模型的接触问题之前，必须定义不同的分析步骤。然后可以指定在哪一步施加荷载，在哪一步施加边界条件，哪一步确定相互关联。 (2)CAE缺省地创建初始步（initial）
分析步创建完成后会自动生成输出结果管理器 (3)输出结果要求
ABAQUS求解器通常计算每一个增量步的许多变量值，而往往我们只对其中某一小部分计算数据感兴趣，软件提供了指定要输出到计算结果数据库中的某些变量结果的功能。输出要求包括以下一些信息： （a）所需要的变量或者变量分量；
（b）模型中某一特定区域和积分点的计算结果； （c）写到计算结果数据库中各变量值的写入频率；
建立了第一分析步后，CAE缺省地选择和相应的分析过程中输出变量集。缺省的情况下，CAE输出模型中每个节点或积分点的计算值。
在一般分析步中，载荷必须以总量而不是以增量的形式给定。例如，如果在分析步1中有一个10kN的集中载荷，而在分析步2中此载荷变为40kN，那么在这两个分析步中，对载荷的定义应该分别是10kN和40kN，而不是10kN和30kN。
? 场变量输出（field）和历程输出（history）
（a）场变量输出：
在通常情况下，用于绘制模型的变形、云图和X–Y图，由于ABAQUS生成的实时输出结果数据库文件都很大，因此可以通过修改输出要求来限制结果数据库的大小。 （b）历程输出：
ABAQYUS对模型中指定点产生历程输出数据。使用后处理模块在XY坐标系中查看历史输出结果。结果的输出频率依赖于如何使用计算生成的各种数据，输出频率可以很高。可以建立历史输出要求，通过该要求限制历史输出频率。在建立历史输出要求时可以指定某一个独立的变量写入输出结果数据库。 ? 通用分析步（general step）和线性摄动分析步（linear perturbation step）
分析步包括通用步和线性摄动步两大类，当在已有的分析步中插入新的通用分析步或者线性摄动分析步时，其上一个分析步相应的输出结果要求会自动传递给该分析步。如果删除一个分析步，相应的结果输出要求以及其后由该步传递的各分析步的输出结果要求都将被删除。如果某一个分析步没有相应的结果输出要求，在计算模块（job）里生成输入文件时将会给出警告。 （1）通用分析步
定义的是一个接一个顺序的分析流程，可以用于线性和非线性分析，主要有以下类型： -static，general 使用ABAQUS/standard进行静力分析 -dynamics，implicit 使用ABAQUS/standard进行隐式动力分析 -dynamics，explicit 使用ABAQUS/explicit进行显式动态分析 （2）线性摄动分析
分析“基础状态”基础上的线性响应，而基础状态是前溯最近的general step(通用分析步)，下一个分析步和Linear perturbation steps是没有关系的。只能用于分析线性问题，在ABAQUS/explicit不能用线性摄动分析，以下类型总是采用线性摄动分析步： -buckle 线性特征值屈曲 -frequency 频率提取分析 -modal dynamics 瞬时模态动力分析 -random response 随机响应分析 -response spectrum 反应谱分析
-steady-state dynamic 谐波激励稳态动力分析
线性分析是基状态（初始构型或当前构型）的线性摄动，基状态之前的响应可以是非线性的。但是，模型必须是静态平衡的(在进行线性摄动分析之前，只有先利用*STATIC分析步达到静力平衡，才可以应用?DYNAMIC选项）。在摄动分析步之后，可以继续进行非线性分析步。在Abaqus/Explicit中，只有通用分析步。
? 时间增量步的设置 （1） 增量步的类型：
ABAQUS/Standard使用Newton-Raphson算法来求解非线性问题，把所有载荷按一定的要求分成若干
载荷步step，每一步step根据ABAQUS自动载荷增量，分成若干增量increments，每一增量施加一定的载荷，然后每一增量通过若干迭代步iteration 进行迭代，当系统达到平衡时，迭代结束，完成一个增量。当所有的增量都完成后，计算结束，所有增量响应的总和就是非线性分析的近似解；反之，计算可能出现发散。这时，可以通过采用多钟方法（如调整放大质量系数，单元网格优化等）调整增量大小，使计算继续进行。
ABAQUS/Explicit在求解非线性问题时不需要进行迭代，而是显示地从上一个增量步的静力学状态来推出动力学平衡方程的解。ABAQUS/Explicit 的求解过程需要大量的增量步，但由于不进行迭代，也不需要求解全体方程组，其每个增量步的计算成本很小，可以很高效地求解复杂的非线性问题。
Automatic即增量步的大小由ABAQUS自动控制，根据分析结果的收敛情况自动增大或减小增量步。在默认情况下，如果经过16次迭代的解仍不能收敛或者结果显示出发散，ABAQUS/Standard就放弃当前增量步，并将增量步的值设置为原来值的25%，重新开始计算。利用比较小的载荷增量来尝试找到收敛的解答。若此增量仍不能使其收敛，ABAQUS/Standard将再次减小增量步的值。在中止分析之前，ABAQUS/Standard 默认地允许至多5次减小增量步的值。如果连续两个增量步都只需少于5次的迭代就可以得到收敛解，ABAQUS/Standard 会自动地将增量步的值提高50%。
（2）允许的增量步最大数目：100，即如果经过100个增量步后结果还不收敛，则分析中止。
（3）初始增量步大小：0.1。用户只需在每个分析步模拟中给出第1个增量步的值，然后，ABAQUS/Standard 自动地调整后续增量步的值。对于简单的问题，可以直接令初始增量步等于分析步时间（例如令初始增量步等1）。对于复杂的非线性问题（例如模型中有复杂的接触或大的塑性变形），ABAQUS/Standard不得不反复减小增量步，从而导致占用了CPU时间以及甚至不能收敛，可以尝试减小初始增量步。 （4）允许的最小增量步：10-5允许的最大增量步：1
（5）在静态分析中，如果模型中不包含阻尼或与速率相关的材料性，“时间”就没有实际的物理意义。方便起见，一般都把分析步时间设为默认的1。
（6）对于复杂的三维问题，如果出现收敛困难，可以使用额外的分析步和边界条件，将荷载逐步施加到模型上。即在接触分析中，如果在第一个分析步中就把全部载荷施加到模型上，有可能分析无法收敛，建议先定义一个只有很小载荷（或位移）的分析步，让接触关系平稳地建立起来，然后在下一个分析步中再施加真实的载荷。这样虽然分析步的数目增多了，但减小了收敛的困难，计算时间可能反而会缩短。 ? 设定自适应网格
分析锻压、拉拔和轧制等大变形问题时，模型的几何形状发生显著变化，网格会产生严重的扭曲变形，导致分析精度下降，稳定步长缩短，甚至无法达到收敛。ABAQUS的自适应网格功能允许单元网格独立于材料移动，从而在大变形分析过程中也能始终保证高质置的网格。
自适应网格主要用于ABAQUS/Explicit, 以及ABAQUS/Standard中的表面磨损过程模拟。在一般的 ABAQUS/Standard分析中尽管也可以设定自适应网格，但不会起到明显的作用。
点击Step 功能模块的主菜单Other——Adaptive Mesh Domain可以设定自适应网格的有效区域，点击主菜单Other——Adaptive Mesh Controls可以设置自适应网格的参数。ABAQUS的自适应网格不改变网格的拓扑结构（单元和连接关系），它结合了纯拉格朗日分析（网格跟随材料终动）和欧拉分析（网格位置固定，材料在网格中流动），被称为“任意拉格朗日- 欧拉（ALE) 分析”。它通常比纯拉格朗日分析更有效、
更精确和更稳定。
对于ABAQUS/Standard 的通用分析步，可以点击Step功能模块的主菜单 Other—General Solution Controls来控制收敛算法和时间积分精度。对于静力问题的通用分析步和线性摄动分析步，以及稳态传热问题，可以点击主菜单Other->Solver Controls 来控制迭代线性方程求解器的参数。 ? 设定几何非线性（Nlgeom）
进行弹塑性分析时，如果模型的位移较大，则设定几何非线性为on。当然弹塑性分析中并不一定要考虑几何非线性，几何非线性的含义是位移的大小对结构的响应发生影响，例如大位移、大转动、初始应力、几何刚性化和突然翻转等。 ? 分析控制
?为Abaqus/Explicit分析定义自适应网格区域和自适应网格控制。 ?为接触问题定制求解控制。
?定制一般求解控制，用于控制Abaqus中的收敛控制参数和时间积分精度算法。
5 选择监视自由度
在分析过程中，可以有效的显示求解历程，为求解过程提供简单的指示。选定某个自由度，指示当前解的位置。比如在结构突变分析过程中，监控选定薄板拱形结构的中点。
6 建立交互作用（接触、约束）
ABAQUS/CAE中的接触分析主要包括以下建模步骤：
1) 在Interaction功能模块、Assembly功能模块或 Load功能模块中定义各个接触面。 2) 在Interaction功能模块中定义接触属性（包括法向接触属性和切向的摩擦属性）。
3) 在Interaction功能模块中定义接触（包括主面、从面、滑动公式、从面位置调整、接触属性、接触面距离和接触控制等）。
4) 在Load功能模块中定义边界条件，保证消除模型的刚体位移。 在Interaction功能模块中，主要可以定义模型的以下相互作用：
(1) 主菜单Interaction定义模型的各部分之间或模型与外部环境之间的力学或热相互作用，例如接触、弹性地基、热辐射等。
(2) 主菜单Constraint 定义模型各部分之间的约束关系。
(3) 主菜单 Connector 定义模型中的两点之间或模型与地面之间的连接单元，用来模拟固定连接、铰接、恒定速度连接、止动装置、内摩擦、失效条件和锁定装置等。
(4) 主菜单Special—Inertia 定义惯量（包括点质量/ 惯量、非结构质量和热容）。 (5) 主菜单Special—Crack 定义裂纹。
(6) 主菜单Special—Springs/Dashpots 定义模型中的两点之间或模型与地面之间的弹簧和阻尼器。 (7) 主菜单Tools常用的菜单项包括Set (集合）、Surface (面)和Amplitude (幅值)等。 ?
接触
接触分析中的关键问题是定义接触属性、接触面和接触关系。
即使两个实体之间或一个装配件的两个区域之间在空间位置上是互相接触的，ABAQUS/CAE也不会自动认为它们之间存在着接触关系，需要使用 Interaction 模块中的主菜单Interaction来定义这种接触关系。
相互作用与分析步有关，必须规定相互作用是在哪些分析步中起作用。
在必要的时候，利用接触管理器激活/不激活接触，以分析其区别。 在三维模型中可以使用自动约束探测快捷方便地定义接触和绑定约束。 接触对中的 slave surface 应该是材料较软，网格较细的面。
接触面之间有微小的距离，定义接触时要设定“Adjust=位置误差限度”，此误差限度要大于两触面之 间的距离，否则ABAQUS 会认为个面没有接触。由于模型中存在数值误差，所以一般要设置这个位置误差限度（0.02）。
小滑移问题的接触压强总是根据未变形时的接触面积来计算的，有限滑移问题的接触压强则是根据变化的接触面积来计算。
如果模型中有塑性材料，或分析过程中会发生很大的位移或局部变形，或施加载荷后会使接触状态发生很大的变化，则应设置较小的初始时间增量步。
在对分析步的定义中可以使用下面关键词*CONTACT PRINT将接触信息输出到DAT文件（ABAQUS/CAE不支持）。CPRESS和CFN的区别是：CPRESS是从面各个节点上各自的接触压强，而CFN代表接触面所有节点接触力的合力，它包含四个变量：CFNM、CFN1、CFN2和CFN3。接触面所有节点在垂直于接触面方向上接触力的合力称为法向接触力。如果接触面是曲面，就无法由CFN直接得到法向接触力，这时可以通过各个从面节点的CPRESS来计算法向接触力
法向接触力=从面上所有节点的CPRESS之和X 从面的面积/从面上的节点数 摩擦力=法向接触力X摩擦系数
利用MSG文件可以查看分析迭代的详细过程，从面节点有开放和闭合两种接触状态。如果在一次迭代中节点的接触状态发生了变化，称为“严重不连续迭代（SDI）”。如果分析能够收敛，每次严重不连续迭代中CLOSURES和OPENINGS的数目会逐渐减少，最终所有从面节点的接触状态都不再发生变化，就进入平衡迭代，直至收敛。如果CLOSURES和OPENINGS的数目逐渐减少，但最终不断重复出现“0 CLOSURES，1OPENINGS”和“1 CLOSURES，0 OPENINGS”（此处的数字也可以大于1），即所谓“振颤”。如果CLOSURES和OPENINGS的数目逐渐减少，但减小的速度很慢，达到第12次严重不连续迭代后，ABAQUS就自动减小增量步长，重新开始迭代。如果增大这个最大次数，允许ABAQUS多进行几次迭代，就有可能达到收敛。操作方法：进入Step模块，主菜单Other→General Solution Controls→Edit，选择相应的分析步，点击Continue，选中Specify，点击Time Incrementation标签页，点击第一个More，把Is由默认的12改为适当的值，然后点击OK。
如果希望在MSG文件中看到更详细的接触分析信息，可以在Step模块中选择菜单Output→Diagnostic Print然后选中Contact。其相应的关键词是*PRINT, CONTACT=YES。 定义主面和从面的一般规则为：
（1）选取刚度大的面作为主面。这里的“刚度”指材料特性和结构刚度。解析面或由刚性单元构成的面必须作为主面，从面则必须是柔体上的面（可以是施加了刚性约束的柔体）。
（2）若两接触面刚度相似，则选取粗糙网格的面作为主面。
（3）如果能使两接触面的网格节点位置一一对应，则能使结果更精确。
（4）主面必须是连续的，由节点构成的面不能作为主面。如果是有限滑移，主面在发生接触的部位必须是光滑的，即不能有尖角。
（5）若主面在发生接触的部位存在尖锐的凹角或凸角，应该在此尖角处把主面分为两部分来分别定义，即定义为两个面。对于有单元构成的主面，ABAQUS会自动进行平滑处理。
（6）若是有限滑移，则在整个分析过程中，都尽量不要让从面节点落到主面之外（尤其不要落在主面的背面），否则容易出现收敛问题。
（7）一对接触面的法线方向应该相反，都指向实体的外部。一般来说，对于三维柔性实体，ABAQUS会自动选择正确的法线方向，而在使用梁单元、壳单元、膜单元、绗架单元或刚体单元来定义接触面时，
用户往往需要自己制订法线方向，就容易出现错误。
解决接触分析中的收敛问题：
（1）检查接触关系、边界条件和约束。
（2）消除刚体位移（Numerical Singularity（数值奇异），有些情况下，还会显示Negative Eigenvalue（负特征值）警告信息）。
（3）一般来说，如果从面上有90o的圆角，建议在此圆角处至少划分10个单元。
（4）如果接触属性为“硬接触”，应尽可能使用六面体一阶单元（C3D8）。如果无法划分六面体单元网格，可以使用修正的四面体二次单元(C3D10M)。
（5）避免过约束，即节点的某个自由度上同时定义了两个以上的约束条件。可能造成过约束的因素有：(a)接触：从面节点会受到沿主面法线方向的约束；(b)边界条件；(c)连接单元；(d)子模型边界；(e)各种约束。
（6）摩擦系数越大，接触分析就越不容易达到收敛。 ?
约束 .
在 ABAQUS/CAE的Assembly 功能模块、Load 功能模块和Interaction 功能模块中都有“约束”的概念，它们分别有不同的含义。在 Assembly 功能模块中，主菜单 Constraint(约束)的作用是定义各个实体间的相互位置关系，从而确定它们在装配件中的初始位置。在Load功能模块中，主菜单 BC 的作用是定义边界条件，消除模型的刚体位移。在Interaction功能模块中，主菜单 Constraint (约束)的作用是定义模型各部分的自由度之间的约束关系，具体包括以下类型。
(1) Tie (绑定约束）模型中的两个面被牢固地粘结在一起，在分析过程中不再分开。被绑定的两个面可以有不同的几何形状和网格。
(2) Rigid Body(刚体约束）在模型的某个区域和一个参考点之间建立刚性连接，此区域变为一个刚体，各节点之间的相对位置在分析过程中保持不变（为整个实体施加刚体约束时，无论实体的类型是Solid、Shell或Wire，都应将刚体约束施加在实体的单元上，即选择Body(elements)）。
(3) Display Body (显示体约束）与 Rigid Body类似，受到此约束的实体只用于图形显示，而不参与分析过程。
*刚体约束和显示体约束本质上是一样的，其共同优点是只要去掉约束，部件就恢复成柔体，可以进行多柔体分析。
(4) Coupling (耦合约束）在模型的某个区域和参考点之间建立约束。
1) Kinematic Coupling (运动耦合)：即在此区域的各节点与参考点之间建立一种运动上的约束关系。 2) Distributing Coupling (分布耦合）：也是在此区域的各节点与参考点之间建立一种约束关系，但是对此区域上各节点的运动进行了加权平均处理，使此区域上受到的合力和合力矩与施加在参考点上的力和力矩相等效。换言之，分布耦合允许面上的各部分之间发生相对变形，比运动耦合中的面更柔软。 . (5) Shell-to-Solid Coupling (壳体-实心体约束）在板壳的边和相邻实心体的面之间建立约束。 (6) Embedded Region (嵌入区域约束）模型的一个区域镶嵌在另一个区域中。 (7) Equation (方程约束）用一个方程来定义几个区域的自由度之间的相互关系。 ?
连接
ABAQUS模拟多体系统的基本思路是：使用2节点的连接单元在模型各部分之间建立连接，并通过定义连接属性来描述各部分之间的相对运动约束关系。
连接点可以是模型中的参考点、网格实体的节点、几何实体的顶点或地面。应尽量选择参考点作为连接单元的连接点，而不要直接使用Solid实体的节点，因为具有旋转属性的连接单元会激活Solid实体节点上的旋转自由度，如果这些旋转自由度没有得到充分的约束，就会造成收敛问题。
连接单元的作用不仅仅是在两个连接点之间施加运动约束，它还有另外一个重要的作用：度量两个连接点的相对运动、力和力矩。
在多体分析中，如果连接属性或边界条件选择不正确，很容易出现过约束。如果ABAQUS无法自动解决过约束问题，则可能出现以下结果：（1）分析过程无法达到收敛；（2）虽然能够达到收敛，但出现远远超过正常数量级的刚体位移；（3）虽然能够达到收敛，位移结果也正确，但某个连接单元反作用力或力矩远远大于应有的值。出现过约束时，在MSG文件中会显示Overconstraint Check和Zero Pivot等警告信息。提示：ABAQUS/Explicit不会显示Zero Pivot等警告信息，因此在进行显示分析前，应首先使用ABAQUS/Standard进行分析，确保模型中没有过约束。
一个正确的多体分析模型应满足如下关系：
实体总数x 6 = 位移边界条件所约束的自由度总数+ 连接单元中受约束的相对运动分量总数 基准坐标系的原点不一定要在连接单元的连接点所在的位置上，只要坐标轴的方向正确即可。 ?
过约束
msg文件中出现 zero pivot 说明 ABAQUS 无法自动解决过约束问题，例如在桩底部的最外一圈节点上即定义了 tie，又定义了 contact, 出现过约束。解决方法是在选择 tie 或 contact 的 slave surface 时，将类型设为 node region, 然后选择区域时不要包含这一圈节点。
7 建立交互作用特性
交互作用是用来建立模型中接触表面或相距很近的表面之间力学关系的对象。可以建立一系列交互作用特性，它和交互作用相互独立，每个交互作用都可以被分配到交互作用特性。交互作用特性共有三种：接触特性（contact）、膜条件特性（file condition）、激励和传导特性（actuator/sensor）
接触交互作用特性可以是切向接触和法向接触，接触面间可以是有摩擦、无摩擦和阻尼接触，还可以相互间分离。接触交互作用特性中通常包含阻尼、热传导、热辐射、摩擦生热等信息。接触交互作用特性可以被通用接触、面对面接触或自我接触等交互作用引用。
膜条件交互作用特性定义膜层传热系数为温度的函数。膜条件特性只能被膜条件交互作用引用。
8 施加边界条件和荷载
实体单元（solid element）只有平动自由度，没有转动自由度，所以施加边界条件时只需约束起平动自由度即可。对于分析刚体来说，约束只能施加给分析刚体的参考点。缺省的情况下，ABAQUS将边界条件传递给其后的每一个分析步。对每一个分析步中的边界条件可以进行编辑和修改。
指定载荷和边界条件可以随着时间相关的幅值定义变化，而且幅值定义既可以参考分析步时间也可以参考总时间。默认情况下，一般静态分析步中使用斜坡幅值曲线（0→1）。在一般分析步中，载荷必须以总量而不是以增量的形式给定。 ?
荷载
1) Concentrated Force ：施加在节点或几何实体顶点上的集中力，表示为力在三个方向上的分量。 2)Moment：施加在节点或几何实体顶点上的弯矩，表示为力矩在三个方向上的分量。 3)Pressure：单位面积荷载（荷载的方向总是与面或边垂直，正值为压力，负值为拉力）。 4) Shell Edge Load：施加在板壳边上的力或弯矩。
5) Surface Traction：施加在面上的单位面积荷载，可以是剪力或任意方向上的力，通过一个向量来描述力的方向。
6) Pipe Pressure：施加在管子内部或外部的压强。 7) Body Force：单位体积上的体力。
8) Line Load：施加在梁上的单位长度线荷载。
9) Gravity：以固定方向施加在整个模型上的均匀加速度，例如重力；ABAQUS根据此加速度和材料属性中的密度来计算相应的荷载。
10) Bolt Load：螺栓或紧固件上的紧固力，或其长度的变化。
11）Generalized Plane Strain：广义平面应变荷载，它施加在由广义平面应变单元所构成区域的参考点上。 12) Rotational body force：由于模型的旋转造成的体力，需要指定角速度或角加速度，以及旋转轴。 13) ConnectorForce：施加在连接单元上的力。 14) Connector Moment：施加在连接单元上的弯矩。 15)温度和电场变量 ?
边界条件
使用主菜中 BC可以定义以下类型的边界条件：（约束）对称/反对称/固支、（施加或约束）位移/转角、速度/角速度、加速度/角加速度、连接单元位移/速度/加速度、温度、声音压力、孔隙压力、电势、质量集中。
XSYMM：对称边界条件，对称面为与坐标轴1垂直的平面，即U1＝UR2＝UR3=0; YSYMM：对称边界条件，对称面为与坐标轴2垂直的平面，即U2＝UR1＝UR3=0; ZSYMM：对称边界条件，对称面为与坐标轴3垂直的平面，即U3＝UR1＝UR2=0; XASYMM：反对称边界条件，对称面为与坐标轴1垂直的平面，即U2＝U3＝UR1=0; YASYMM：反对称边界条件，对称面为与坐标轴2垂直的平面，即U1＝U3＝UR2=0; ZASYMM：反对称边界条件，对称面为与坐标轴3垂直的平面，即U1＝U2＝UR3=0; PINNED：约束所有平移自由度，即U1=U2=U3=0;
ENCASTRE：约束所有自由度（固支边界条件），即U1=U2=U3=UR1=UR2=UR3=0.
在边界条件中给出的位移值是相对于模型初始状态的绝对位移值，而不是当前分析步中的增量值。 ?
场变量和荷载状况
使用主菜单Field可以定义场变量（包括初始速度场和温度场变量）。有些场变量与分析步有关，也有些仅仅作用于分析的开始阶段。使用主菜单Load Case 可以定义荷载状况，荷载状况由一系列的荷载和边界条件组成，用于静力摄动分析和稳态动力分析。
9 网格划分
什么是网格？物理部件模型的几何近似，包含许多几何上简单的节点和单元的离散几何体。 （1）进入单元划分模块后，ABAQUS的颜色代表该模型中不同区域适合用哪种方法就行单元划分。绿色表示可以采用结构法划分，黄色表示可以用扫掠法划分，橙色表示该区域不能用缺省的单元（实体单元缺省的单元为六面体单元hexahedral）形状进行单元划分，必须对该区域进行分解后才能用缺省的单元形状进
行单元划分。当然，可以采用四面体单元（tetrahedral）利用自由网格技术对任何形状的模型区域进行单元划分。 ? ? ? ?
自由划分：可以应用到任意平面和曲面
结构划分：结构化的网格划分通常给出了对网格的最大的控制 扫掠划分：在源面中，相邻面之间的二面角不能和180°相差太远。
虚拟拓扑：在某些情况下，装配件的部件实例可能包含一些小的细节，比如表面和边。虚拟拓扑可以忽略这些不重要的细节。
如果部件实例中包含虚拟拓扑，那么它只能使用以下单元通过自由网格技术划分网格： ?自由网格：三角形和四面体单元、用波前法划分的四边形为四边形为主单元网格 ?扫略网格：六面体或者楔形单元
?映射网格：四边形, 三角形, 或者六面体单元 （2） 分解模型（partition）
可以对模型中的边（edge）、面（face）和体（cell）进行分解。用来将边、面、体分解成更小部分的点、边、面都成为模型中的特征体，这些特征体和其他特征体一样可以在特征体管理器中查看。（如：将一个体分解成两部分需要用一个面将体切割成两部分，这个面就成了模型中一个新的特征体。）。
有五种分割特征体可以将一个特征体分解：定义切割面（define cutting plane）、使用基准平面（use datum plane）、延伸平面（extended face）、挤压或旋转边（extrude/sweep edges）、N-sided patch。
一次分解操作仅仅只是将被分解的对象分解成两部分，并不能改变被分解对象所在特征体（部件实例）的整体性。不能对刚体进行单元划分。 （3）设置网格种子 选择Seed part或SeedEdge ? ?
设置全局种子 设置边上的种子 边上的种子无约束圆圈 边上的种子受部分约束三角形 边上的种子受完全约束方形
对于可以用扫掠方法进行网格划分的区域，边的网格种子从选定的边到匹配边自动传播，由分区创建的新边自动继承总体网格种子。 （4）单元划分控制
在 Assign mesh control中指定单元划分方法（结构划分法（structured）、自由划分法（free）、扫掠法（sweep）等等）。操作的对象是被分解后的边、面、和体，可以对同一实例（装配模型的特征体）分解后产生的不同边、面、体分别采用不同的单元划分方法，指定不同的单元类型。
在Assign element type中指定单元类型（六面体单元、四面体单元等等），选择单元库（standard、explicit）、确定线性单元（linear）或者二次单元（quadratic）、确定这两种单元的特性：杂交元（hybrid formulation）、缩减积分（reduced integration）、非协调单元模式（incompatible modes）。 1)
单元形状 二维问题
a) Quad：网格中完全使用四边形单元。
b) Quad-dominated：网格中主要使用四边形单元，但在过渡区域允许出现三角形单元，更容易实现从粗网
格向细网格过渡。
c) Tri: 网格中完全使用三角形单元。 三维问题
a) Hex：网格中完全使用六面体单元。
b) Hex-dominated：网格中主要使用六面体单元，但在过渡区域允许出现锲形（三棱柱）单元。
c) Tet：网格中完全使用四面体单元。 d) Wedge：网格中完全使用楔形单元。
Quad单元（二维）和Hex单元（三维）可以用较小的计算代价得到较高的精度。自由网格采用Tri单元（二维）和Tet单元（三维），一般应选择带内部节点的二次单元来保证精度。结构化网格和扫掠网格一般采用quad单元（二维单元）和hex单元（三维），分析精度较高，因此在划分单元时选择后两种。 2)
三维实体单元类型 ? ? ?
线性（linear)单元又称一阶单元，仅在单元的角点处布置节点，在各方向都采用线性插值； 二次（quadratic)单元又称二阶单元，在每条边上有中间节点，采用二次插值；
修正的（modified)二次单元只有Tri或Tet单元才有这种类型，即在每条边上有中间节点，并采用修正的二次插值
b) 连续体单元
ABAQUS/Standard 的连续体单元库包括二维和三维的线性单元和二次单元，分别可以采用完全积分或减缩积分，另外还有修正的二次Tri 和Tet单元，以及非协调模式单元和杂交单元。
ABAQUS/Explicit 的连续体单元库包括二维和三维的线性减缩积分单元，以及修正的二次Tri和Tet单元。ABAQUS/Explicit 中没有二次完全积分的连续体单元。 c) 线性完全积分（ linear full-integration ) 单元
保持默认的Linear参数，取消对Reduced integration(减缩积分)的选择，例如CPS4单元（4节点四边形双线形平面应力完全积分单元）和C3D8单元（8节点六面体线性完全积分单元）。所谓“完全积分”是指当单元具有规则形状时，所用的离斯积分点的数目足以对单元刚度矩阵中的多项式进行精确积分。承受弯曲荷载时，线性完全积分单元会出现剪切自锁(shear locking)问题，造成单元过于刚硬，即使划分很细的网格，计算精度仍然很差。
d) 二次完全积分（quadratic-full-integration ) 单元
在Element Type对话框选择 Quadratic 参数，取消对Reduced integration的选择，例如CPS8 (8节点四边形二次平面应力完全积分单元）和C3D20 (20节点六面体二次完全积分单元）。
二次完全积分单元的优点如下。
? 对应力的计算结果很精确，适于模拟应力集中问题； ? —般情况下没有剪切自锁。
但使用这种单元时要注意以下问题 ? 不能用于接触分析；
? 对于弹塑性分析，如果材料是不可压缩性的（例如金属材料），则容易产生体积自锁 (volumetric locking)； ? 当单元发生扭曲或弯曲应力有梯度时，有可能出现某种程度的自锁。 e) 线性减縮积分（linear reduced-integration) 单元
对于Quad单元和Hex单元ABAQUS/CAE 默认的单元类型是线性减缩积分单元，例如CPS4R (4节点四边形双线形平面应力减缩积分单元)和C3D8R（8节点六面体线性减缩积分单元）。减缩积分单元比通常的完全积分单元在每个方向少用一个积分点。线性减缩积分单元只在单元的中心有一个积分点，由于存在所谓“沙漏”（hourglass) 数值问题而过于柔软。ABAQUS在线性减缩积分单元中引人了“沙漏刚度”以限制沙漏模式的扩展。模型中的单元越多，这种刚度对沙漏模式的限制越有效。可以选择沙漏控制参数为Enhanced , RelaxStiffness ,Stiffness ,Viscous或combined。采用线性减缩积分单元模拟承受弯曲荷载的结构时，沿厚度方向上至少应划分四个单元。
线性缩减积分单元有以下优点： ? 对位移的求解结果较精确。
? 网格存在扭曲变形时（例如 Quad单元的角度远远大于或小于 90°），分析精度不会受到大的影响。 ? 在弯曲荷载下不容易发生剪切自锁。 a) 节点数目和插值阶数
其缺点如下：
? 需要划分较细的网格来克服沙漏问题。
? 如果希望以应力集中部位的节点应力作为分析指标，则尽量不要使用线性减缩积分单元，而应使用二次单元，因为线性缩减积分单元只在单元的中心有一个积分点，相当于常应力单元，它在积分点上的应力结果是相对精确的，而经过外插值和平均后得到的节点应力则不精确。如果在应力集中部位进行了网格细化，使用二次减缩积分单元与二次完全积分单元得到的应力结果相差不大，而二次减缩积分单元的计算时间相对较短。
Note：在查看模型的应力结果时有两种选择：
? 查看节点上的应力：这是通常用的方法，其优点是简单方便，但事实上，后处理中得到的节点应力是对单元积分点上的应力进行外插值和平均后得到的并不精确；
? 查看单元积分点上的应力：这是ABAQUS所推荐的方法。线性缩减积分单元只有一个积分点，可以很方便地查看积分点上的分析结果，但其他类型的单元有多个积分点，就需要详细了解节点的编号顺序，并根据模型的实际情况来决定查看哪个积分点，这一过程很烦琐。另外要注意，单元积分点上的应力值往往不是应力集中区城的最大应力。
(用户可以自己在上述两种方法中做出选择，需要注意的是，如果希望查看节点上的应力，就尽量不要使用线性缩减积分单元；如果使用了线性缩减积分单元，就应该查看单元积分点上的分析结果，并且要在应力变化剧烈的部位划分足够细的网格。如果外推应力值与积分点应力值差别很大，说明单元间应力变化剧烈，单元网格过于粗糙，计算的应力不够精确。这种外推应力误差会因为单元网格细化而减小，但总是存在。所以在使用单元变量的节点值时要谨慎。默认的应力不变量的计算方法是computer scalars before averaging（先不变量再平均），得到的节点应力偏大，作为工程分析的结果会更安全；也可以在result options中设置为computer scalars after averaging（先平均再不变量），得到的节点应力偏小。事实上只有单元积分点上的应力结果是相对精确的，一般在所关心的高应力部位细化网格，使用默认设置即可。更广义地说，有限元作为一种数值计算方法，其本身就是不精确的。 f) 二次缩减积分（quadratic-reduced-integration）单元
对于Quad单元或Hex单元，可以在 Element Type 对话框中将单元类型设置为二次减缩积分单元，如 CPS8R (8节点四边形二次平面应力减缩积分单元）和C3D20R(20节点六面体二次缩减积分中元），这种单元不但保待了前面介绍的线性缩减积分单元的优点，而且还具有以下特性： ? 即使不划分很细的网格也不会出现严重的沙漏问题 ? 即使在复杂应力状态下，对自锁问题也不敏感 但使用这种单元时需要注意以下问题： ? 不能在接触分析中使用 ? 不适于大应变问题
? 存在与线性减缩积分单元相类似的问题，由于积分点少，得到的节点应力的精度往往低于二次完全积分单元。
g) 非协调模式（ incompatible modes) 单元
对于Quad单元或Hex单元，可以在Element Type 对话框中将单元类型设为非协调模式单元。例如CPS4（4节点四边形双线形平曲应力非协调单元）和C3D8I（8节点六面体线性非协调模式单元）。仅在 ABAQUS/Standard中有非协调模式单元。其目的是克服在线性完全积分单元中的剪切自锁问题。
非协调摸式单元的优点如下：
? 克服了剪切自锁问题，在单元扭曲比较小的情况下，得到的位移和应力结果很精确；
? 在弯曲问题中，在厚度方向上只需很少的单元，就可以得到与二次单元相当的结果，而计算成本却明显降低；
? 使用了增强变形梯度的非协调模式，单元交界处不会重叠和开洞，因此很容易扩展到非线性、有限应变的位移。
但使用这种单元时需注意，如果所关心部位的单元扭曲比较大，尤其是出现交错扭曲时，分析精度会降低。
h) Tri单元和Tet单元
使用Tri单元和Tet单元时应注意以下问题：
? 线性Tri单元和Tet单元的精度很差，所以不要在模型中所关心的部位及其附近区域使用。
? 二次 Tri 单元和 Trt 单元精度较高，而且能模拟任意的几何形状，但计算代价比Quad单元或Hex单元大，因此如果模型中能够使用 Quad单元或Hex单元，尽量不要使用Tri 单元或Trt单元。 ? 二次Trt单元（C3D10) 适于ABAQUS/standard中的小位移无接触问题；修正的二次Tet单元（C3D10M) 适于ABAQUS/Explicit，以及 ABAQUS/standard中的大变形和接触问题。 ? 使用自由网格不易通过布置种子来控制实体内部的单元大小。 i) 杂交(hybrid) 单元
在ABAQUS/standard中，每一种实体单元（包括所有缩减积分和非协调模式单元）都有其相应的杂交单元，用于不可压縮材料（泊松比为0.5) 或近似不可压缩材料（泊松比大于0.475)，橡胶就是一种典型的不可压缩材料。除了平面应力问题之外，不能用普通单元来模拟不可压缩材料的响应，因为此时单元中的压应力是不确定的。杂交单元在它的名字中含有字母H。ABAQUS/Explicit 中没有杂交单元。 j) 混合使用不同类型的单元
当三维实体几何形状较复杂时，无法在整个实体上使用结构化网格或扫掠网格划分技术，得到Hex单元网格，这时一种常用的做法是对于实体不重要的部分使用自由网格划分技术，生成Tet单元网格，而对于所关心的部分采用结构化网格或扫掠网格生成 Hex 单元网格，在生成这样的网格时，ABAQUS会给出提示将生成非协调的网格，在不同单元类型的交界处将自动创建绑定（tie）约束。
需要注意的是，在不同单元类型网格的交界处.，即使单元角部节点是重合的，仍然有可能出现不连续的应力场，而且在交界处的应力可能大幅度地增大。如果在实体中混合使用线性和二次单元，也会出现类似的问题。因此在混合使用不同类型的单元时，应确保其交界处远离所关心的区域，并仔细检查分析结果 是否正确。
对于无法完全采用Hex单元网格的实体，还可以使用以下方法：
? 对整个实体划分Tet 单元网格，使用二次单元C3D10或修正的二次单元 C3D10M，同样可以达到所需的精度，只是计算时间较长。
? 改变实体中不重要部位的几何形状，然后对整个实体采用Hex单元网格。 k)
选择三维实体单元类型的基本原则
计算代价，提高计算精度。当几何形状复杂时，也可以在不重要的区域使用少量楔形（Wedge) 单元。 ? 如果使用了自由网格划分技术，Tet单元的类型应选择二次单元。在ABAQUS/Explicit中应选择修正的 Tet单元C3D10M，在,ABAQUS/Standard中可以选择 C3D10，但如果有大的塑性变形，或模型中存在接触，而且使用的是默认的“硬”接触关系，则也应选择修正的Tet单元 C3D10M。
? ABAQUS的所有单元均可用于动态分析，选取单元的一般原则与静力分析相同。但在使用 ABAQUS/Explicit 模拟冲击或爆炸荷载时，应选用线性单元，因为它们具有集中质量公式，模拟应力波的效果优于二次单元所采用的一致质量公式。
如果使用的求解器是 ABAQUS/Standard, 在选择单元类型时还应注意以下方面：
? 对于应力集中问题，尽量不要使用线性减缩积分单元，可使用二次单元来提高精度。如果在应力集中部位进行了网格细化，使用二次减缩积分单元与二次完全积分单元得到的应力结果相差不大，而二次减缩积分单元的计算时间相对较短。
? 对于弹塑性分析，如果材料是不可压缩性的（例如金属材料），则不能使用二次完全积分单元（C3D20），否则会出现体积自锁问题，也不要使用二次Tri 单元或Tet单元。推荐使用的是修正的二次Tri单元或Tet单元(C3D10M)、非协调单元（C3D8I），以及线性减缩积分单元（C3D8R）。如果使用二次减缩积分? 对于三维区域，尽可能采用结构化网格划分技术或扫掠网格划分技术，从而得到Hex 单元网格，减小
单元，当应变超过20%-40%时要划分足够密的网格。
? 塑性材料和接触面上都不能用 C3D20R 和 C3D20 单元，这可能是你收敛问题的主要原因。如果需要得到应力，可以使用 C3D8I (在所关心的部位要让单元角度尽量近 90 度)，如果只关心应变和位移，可以使用 C3D8R, 几何形状复杂时，可以使用 C3D10M。
? 如果模型中存在接触或大的扭曲变形，则应使用线性Quad或Hex单元，以及修正的二次Tri单元或Tet 单元，而不能使用其他的二次单元。
? 对于以弯曲为主的问题，如果能够保证在所关心部位的单元扭曲较小，使用非协调单元（例如C3D8I单元）可以得到非常精确的结果。
? 除了平面应力问题之外，如果材料是完全不可压缩的（例如橡胶材料），则应使用杂交单元；在某些情况下，对于近似不可压缩材料也应使用杂交单元。 3)
选择壳单元的类型
如果一个薄壁构件的厚度远小于其典型整体结构尺寸（一般为小于 1/10)，并且可以忽略厚度方向的应力，就可以用壳单元来模拟此结构。壳体问题可以分为两类：薄壳问题(忽略横向剪切变形）和厚壳问题（考虑横向剪切变形）。对于单一各向同性材料，一般当厚度和跨度的比值小于1/15 时，可以认为是薄壳；大于 1/15 时，则可以认为是厚壳。对于复合材料，这个比值需要更小一些。 ABAQUS 的壳单元可以有多种分类方法，按照薄壳和厚壳可划分为： ? 通用目的（general-purpose) 壳单元：此类单元对薄壳和厚壳问题均有效。
? 特殊用途(special-purpose) 壳单元：包括纯薄壳 (thin-only) 单元和纯厚壳 ( thick-only) 单元。 根据单元的定义方式，还可以将 ABAQUS 壳单元划分为：
? 常规（conventional) 壳单元：通过定义单元的平面尺寸、表面法向和初始曲率来对参考面进行离散，只能在截面厲性中定义壳的厚度，而不能通过节点来定义壳的厚度。
? 连续体（continuum) 壳单元：类似于三维实体单元，对整个三维结构进行离散。 选择壳单元的类型时可以遵循以下原则：
? 对于薄壳问题，常规壳单元的性能优于连续体壳单元；而对于接触问题，连续体壳单元的计算结果更加精确，因为它能在双面接触中考虑厚度的变化。
? 如果需要考虑薄膜模式或弯曲模式的沙漏问题，或模型中有面内弯曲，在 ABAQUS/Standard 中使用S4单元（4节点四边形有限薄膜应变线性完全积分壳单元）可以获得很高的精度。 ? S4R单元（4节点四边形有限薄膜应变线性减缩积分壳单元）性能稳定，适用范围很广。
? S3/S3R 单元（3节点三角形有限薄膜应变线性壳单元）可以作为通用壳单元使用。由于单元中的常应变近似，需要划分较细的网格来模拟弯曲变形或高应变梯度。
? 对于复合材料，为模拟剪切变形的影响，应使用适于厚壳的单元（例如 S4、S4R、S3、S3R、S8R），并要注意检查截面是否保持平面。
? 四边形或三角形的二次売单元对剪切自锁或薄膜自锁都不敏感，适用于一般的小应变薄壳。
? 在接触模拟中，如果必须使用二次单元，不要选择 STRI65单元（三角形二次壳单元），而应使用 S9R5 单元（9节点四边形壳单元）。
? 如果模型规模很大且只表现几何线性，使用 S4R5单元（线性薄壳单元）比通用壳单元更节约计算成本。 ? 在ABAQUS/Explicit中，如果包含任意大转动和小薄膜应变，应选用小薄膜应变单元。 4)
选择梁单元的类型
如果一个构件横截面的尺寸远小于其轴向尺度（一般的判据为小于1/10)，并且沿长度方向的应力是最重要的因素，就可以用梁单元来模拟此结构。ABAQUS中的所有梁单元都是梁柱类单元，即可以产生轴向变形、弯曲变形和扭转变形。Timoshenko 梁单元还考虑了横向剪切变形的影响。B21 和 B31 单元（线性梁单元）以及B22和 B32单元（二次梁单元）是考虑剪切变形的Timoshenko梁单元，它们既适用于模拟剪切变形起重要作用的深梁，又适用于模拟剪切变形不太重要的细长梁。这些单元的横截面特性与厚壳单元的横截面特性相同。
ABAQUS/Standard 中的三次单元B23和B33被称为Euler-Bernoulli 梁单元，它们不能模拟剪切变形，但适合于模拟细长的构件（横截面的尺寸小于轴向尺度的1/10 )。由于三次单元可以模拟沿长度方向的三阶变量，所以只需划分很少的单元就可以得到很精确的结果。
选择梁单元的类型可以遵循以下原则。
? 在任何包含接触的问题中，应使用B21或B31单元（线性剪切变形梁单元）。 ? 如果横向剪切变形很重要，则应采用B22和B32单元（二次Timoshenko梁单元）。
? 在ABAQUS/Standard 的几何非线性模拟中，如果结构非常刚硬或非常柔软，应使用杂交单元，例如B21H和B32H单元。
? 如果在 ABAQUS/Standard中模拟具有开口薄壁横截面的结构，应使用基于横截面翘曲理论的梁单元，例如B310S、B320S单元。 5) ?
划分网格的算法： Medial Axis （中轴）算法
首先把要划分网格的区域分为一些简单的区域，然后使用结构化网格划分技术来为这些简单的区域划分网格。Medial Axis 算法具有以下特性：
a) 使用 Medial Axis 法更容易得到单元形状规则的网格，但网格与种子的位置吻合得较差。
b) 在二维模型中使用Medial Axis算法时，选择 Minimize the mesh transition（最小化网格的过渡）可以
提髙网格的质量，但是生成的网格更容易偏离种子。
c) 如果在模型的一部分边上定义了受完全约束的种子，Medial Axis算法会自动为其他的边选择最佳的种
子分布。
d) MedialAxis算法不支持由CAD模型导人的不精确模型 (imprecise part) 和虚拟拓扑( virtual topology ) ?
Advancing Front（波前）算法
Advancing Front算法首先在边界上生成四边形网格，然后再向区域内部扩展，它具有以下特性： a) 使用Advancing Front算法得到的网格可以与种子的位置很好地吻合，但在较窄的区域内，精确匹配每
粒种子可能会使网格歪斜。
b) 使用 Advancing From算法更容易得到单元大小均匀的网格。有些情况下，单元尺寸均匀很重要，例如
在ABAQUS/Explicit中，网格中的小单元会限制增量步长。 c) 使用Advancing From算法容易实现从粗网格到细网格的过渡。 d) Advancing From算法支持不精确模型和二维模型的虚拟拓扑。
在实际应用中，具体选择哪种算法更好，往往需要自己去尝试。一般情况下，种子布置得较密，使用Advancing From算法得到的单元大小更均匀，而且能与种子的位置精确地匹配；当种子布置得较稀疏时，使用MedialAxis算法得到的单元形状更加规则，但没有准确地匹配种子的位置，Advancing From算法准确地匹配了种子的位置，但因此导致单元形状歪斜。 （5）划分单元
划分网格失败可能有多种原因，例如：
1) 几何模型有问题，例如模型中有自由边或很小的边、面、尖角、缝隙等 2) 种子布置得太稀疏。 （6）检查网格质量
在Mesh功能模块中点击左侧工具区中的 (Verify Mesh)，可以选择部件、实体、几何区域或单元，检查其网格的质量，获得节点和单元信息。选择Analysis Checks (分析检查）可以检查分析过程中会导致错误或警告信息的单元。
在开始模拟时，可以先简单地划分几何粗网格，通过理解粗网格模拟的结果（高应力区），然后在适当的区域细分网格。预测准确的应力比计算准确的位移需要更加细化的网格。
10 提交工作job
?在Abaqus/CAE中可以使用Keywords Editor直接编辑输入文件，或者用其它的文本编辑器对输入文件进行编辑。Abaqus/Standard或Abaqus/Explicit中的一些功能在Abaqus/CAE中不支持，用户可以利用Keywords Editor添加这些功能。
?利用Keywords Editor对模型所作的修改将作为模型数据库的一部分。所作的修改将被保留。 ?如果利用文本编辑器编辑输入文件，所作的修改将不影响模型数据库。 监控求解过程：
引起警告和错误的区域将自动被创建为单元集和节点集。可视化模块可以利用这些集，以帮助调试模型。被监控自由度的X–Y曲线在新的视图中被自动创建，并随着程序的运行自动更新。
11 画布对象
画布可以看作是一个无限的屏幕或黑板，在上面可以布置各种对象。 画布对象包括三大类：视窗、画布文字注释、画布箭头注释。
（1）视窗是画布上显示模型和分析结果的对象。可以画布上随意建立和删除视窗，控制其尺寸、位置和外观，但是画布上至少有一个视窗对象，不能全部删除所有的视窗对象。
（2）文字注释和箭头注释只能根据画布定位，与视窗无关，可以在视窗之内也可以在视窗之外，移动视窗对文字和箭头注释的位置没有任何影响，但是可以调整他们的位置使得他们处于视窗之中。
12 草图模块（sketch）
草图是二维的剖面图，可以用于生成三维部件。在草图模块中可以定义平面部件、梁、或者分割体用于挤压、平扫、旋转等方法形成三维部件。在草图模块中也可以定义与特征体无关的独立的二维平面断面图。
13 后处理（Visualization ）文件输出
(1) *File Output：定义输出到结果文件
File Output选项可以输出节点、单元、整体数据到选定的文件。*EL FILE、*ENGERGY FILE和*NODE FILE 选项必须和*FILE OUTPUT选相联使用。
ABAQUS输入文件（input file）
ABAQUS输入文件包含模型数据和历史数据。模型数据定义有限元模型：单元、节点、单元特性、材料定义等等。模型数据用来组织生成部件，部件经过装配后生成各种模型。
历史数据定义对模型的操作，即求解模型响应所需要的时间顺序或加载情况等。在ABAQUS里将这个历史过程分解为不同的分析步。每一个分析步都是某一特定类型的响应，如静载、动力响应，土体瞬时固结等等。分析步的定义必须包括过程类型（静态应力分析、热传导分析等）、时间积分和非线性求解控制参数、荷载和输出控制。 非线性求解步和线性慑动分析步
ABAQUS中的非线性求解步和线性慑动分析步有着明显的差别。非线性分析步定义一系列事件，上一个非线性步必须为下一个非线性步提供初始条件。线性慑动分析步提供了系统基本状态（BASE STATE）的线性响应，基本状态也就是优先于线性慑动分析步的最后一个非线性分析步。
每一个非线性分析步都必须把前一个非线性分析步的状态作为自己的初始条件。例如，动力分析可以不加载，动力响应主要来自静力分析步中所储存应变能的释放。 （2）计算结果输出到data file或者results file
所给定的场变量或历程变量可以通过下面Keyword写入.dat文件，但是不能在CAE中实现。
*CONTACT PRINT *EL PRINT *ENERGY PRINT *INTERACTION PRINT *MODAL PRINT *NODE PRINT *SECTION PRINT
14 查看诊断信息
具体操作方法： 在Visualization功能模块的主菜单中选择Tools→Job Diagnostics，在弹出的Job Diagnostics对话框中，点击左侧区域中的加号，可以看到整个分析过程。如果警告信息（如严重不连续迭代SDI）没有出现在各个增量步的最后一次迭代，那么分析结果就是正确的。
15 子模型
子模型是在全局模型的基础上，对局部进行网格细化，作进一步分析，子模型是从全局模型上切分下来的一部分；子结构是将模型的局部作为一个整体来处理，缩聚其内部自由度，只保留与外部有连接关系的自由度，从而减小刚度矩阵和质量矩阵的规模和计算量。子结构往往用于具有相同特征和性质的重复性局部结构。
子模型的驱动变量（driven variable）一般是位移。全局模型在子模型边界上的位移结果， 被作为边界条件来引入子模型。如果全局模型和子模型在子模型边界上的节点分布不同，ABAQUS 会对全局模型在此处的位移结果进行插值处理。 子模型分析的基本步骤：
（1）完成对全局模型的分析，并保存子模型边界附近的分析结果。
提示：全局模型在子模型边界上的位移结果是否准确，会在很大程度上影响子模型的分析结果精度。因此要保证全局模型在子模型边界上有足够细化的网格，另外还要尽量选择位移变化不剧烈的位置作为子模型边界。
（2）创建子模型，定义子模型边界。 （3）设置各个分析步中的驱动变量。
（4）设置子模型的边界条件、载荷、接触和约束。 （5）提交对子模型的分析，检查分析结果。
对于同一个分析步，全局模型和子模型的增量步长可以不同，ABAQUS会自动对其进行插值处理（对于大变形分析也没有问题）。此外，选择子模型边界时，要注意避免发生过约束。
ABAQUS/Standard 和 ABAQUS/Explicit 之间的主要区别 参 量 单元库 分析过程 材料模型 ABAQUS/Standard 提供了丰富的单元库 单元是在ABAQUS/Standard 中单元的子集 一般过程和线性摄动过程 提供了广泛的材料模型 个显著的E别是提供了允许材料失效的模型 对于求解接触问题 ，具有强健的 接触公式 能力 应用基于刚度的求解技术 ，具有无条求解技术 件稳定性 由于在增量步中大量的迭代，可占用磁盘空间和内存 大量的磁盘空间和内存 ABAQUS/Standard 要小很多 使用显式方法，机时消耗与单元数量成正比，并且大致与最小单元的尺寸成反比；磁盘空间和内存需求与单元数目成正比，与单元尺寸无关。隐式方法，经验表明，对于许多问题的计算成本大致与自由度数目的平方成正比。
磁盘空间和内存的占用量相对于应用显式积分求解技术，具有条件稳定性 接触模拟 具有强健的接触功能，甚至能够解决最复杂的一般过程 类似于在 ABAQUS/Standard 中的材料模型，一ABAQUS/Explicit 提供了适用于显式分析的丰富的单元库 ，这些