abaqus中单元的选择宝典

1.完全积分就是指当单元具有规则形状时，所用得高斯积分点可以对单元刚度

矩阵中得多项式进行精确地积分。

2.剪力自锁将使单元变得“刚硬”，只影响受弯曲荷载得完全积分线性（一阶）

单元，这些单元功能在受直接或剪切荷载时没有问题。二次单元得边界可以弯曲，没有剪力自锁得问题。

3.只有四边形与六面体单元才能采用减缩积分。所有得楔形、四面体与三角形

实体单元采用完全积分。减缩积分单元比完全积分单元在每个方向上少用一个积分点。

4.只有四边形与六面体单元才能采用减缩积分。所有得楔形、四面体与三角形

实体单元采用完全积分。减缩积分单元比完全积分单元在每个方向上少用一个积分点。

5.非协调单元：只有四边形与六面体单元才能采用减缩积分。所有得楔形、四

面体与三角形实体单元采用完全积分。减缩积分单元比完全积分单元在每个方向上少用一个积分点。

6.ABAQUS对非协调单元采用了增强位移梯度形式。在弯曲问题中，用非协调

单元可得到与二次单元相当得结果，且计算费用明显降低。对单元扭曲很敏感。

7.ABAQUS对非协调单元采用了增强位移梯度形式。在弯曲问题中，用非协调

单元可得到与二次单元相当得结果，且计算费用明显降低。对单元扭曲很敏感。

8.杂交单元：ABAQUS对非协调单元采用了增强位移梯度形式。在弯曲问题中，

用非协调单元可得到与二次单元相当得结果，且计算费用明显降低。对单元扭曲很敏感。

9.一般情况下应采用二次减缩积分单元（CAX8R，CPE8R，CPS8R，C3D20R）。

在应力集中局部采用二次完全积分单元（CAX8，CPE8，CPS8，C3D20）。对含有非常大得网格扭曲模拟（大应变分析），采用细网格划分得线性减缩积分单元（CAX4R，CPE4R，CPS4R，C3D8R ）。对接触问题采用线性减缩积分单元或非协调单元（CAX4I，CPE4I，CPS4II，C3D8I等）得细网格划分。10.采用非协调单元时应使网格扭曲减至最小。三维情况应尽可能采用块状单元

（六面体）。对小位移问题采用二次四面体单元（C3D10）就是可行得。11.在实体单元中所用得数学公式与积分阶数对分析得精度与花费有非常显著

得影响。使用完全积分得单元，尤其就是一阶（线性）单元，容易形成自锁现象，在正常情况下不要应用。一阶减缩积分单元容易出现沙漏现象；充分得单元细化可减小这种问题。

12.在分析中如有弯曲位移，且采用一阶减缩积分单元时，应在厚度方向上至少

用4个单元。沙漏现象在二阶减缩积分单元中较少见。在大多数一般问题中要考虑应用这些单元。非协调单元得精度依赖于单元扭曲得量值。

13.结果得数值精度依赖于所用得网格。应进行网格细化研究已确保该网格对问

题提供了唯一得解答。但就是应记住使用一个收敛网格不能保证计算结果与问题得行为相匹配：它还依赖于模型其它方面得近似化与理想化程度。通常只在想要得到精确结果得区域细化网格。ABAQUS具有一些先进得特点，如子模型，它可以帮助对复杂模拟得到有用得结果。HOURGLASS

基础部分

Part类型:可变形部件，离散刚体部件（任意形状，荷载作用下不可变形），解析刚体部件（只可以用直线，圆弧与抛物线创建得形状，荷载作用下不可变形）。每个部件只存在自己得坐标系中，与其她部件无关。给部件赋予属性，既成为实例。实例可以装配成assembly。

Automated repair options：默认为缝合边，自动修理用于几何体变成valid。

基特征一旦创建不能修改。附加特征可用于修改基特征或为基特征添加细节（拉伸，壳，线，切削，导角）

基准几何体类型：点，轴，坐标系，平面。

过滤器：selection options

分区：细分为不同得区域

对于拉伸与旋转，有扭曲选项，可以创建螺纹、螺旋弹簧与扭曲线。也可以利用锥度选项，指定角度，创建带有锥度得部件。

导入孤立网格：通过、inp与、odb文件导入已有网格。被导入得孤立网格，没有父几何体。

定义表面增强：定义了连接到已有部件表面得表面，并指定她得工程属性。

如何给部件定位：相对定位：定义几何关系，确定规则，表面平行约束，面面平行约束，共轴约束，接触约束，重合点约束，平行坐标系约束，若定义有冲突，则将之前得相对约束转化为绝对约束。

集与表面在assembly，step，interaction与load模块中均有效。在part or property module 中创建得part集在assembly module 中有效，但不能通过set managerment修改。

Step 用途：define step，指定输出需求，指定分析诊断，指定分析控制。接触、荷载与边界条件就是分析步相关得，需事先定义。主要用于描述模拟历程。对python与c++保留了API接口，用于后处理。输出类型有两种类型：场数据用于绘制模型得变形，云图与X-Y图；历程数据用于X-Y绘图。分析步可替换。分析控制：为显式分析定义自适应网格区域与控制；为接触问题定制求解控制；定制一般得求解控制。

Interaction：用于模拟机械或热得接触。如定义边界得耦合，定义连接器。显示体得目得就是可视化，不用于分析。接触模型得法向关系、摩擦与干涉。带有摩擦得双面接触、自接触、捆绑约束。使用步骤：create ，选择起作用得step；选择表面；在edit interaction对话框中完成接触定义；在接触管理器中激活或不激活。

边界条件：包括初始温度、指定得平移或转动，速度或角速度。指定得边界条件可以随着时间相关得幅值定义。

初始条件：包括平动与转动速度、温度。初始平动速度可以模拟自由落体得效果。步骤：创建、指定对象、编辑。

Mesh module ：分网技术，单元形状，单元类型，网格密度，生成网格，检查网格状况。

二维区域可用形状：四边形、以四边形为主（允许三角形单元作为过度）、三角形

三维区域可用形状：四面体、若实例中包含虚拟拓扑，可使用三角形单元、四边形单元与利用波前算法得四边形或四边形为主得单元。

细节模型中，小得细节可能会影响网格效果，虚拟模型则忽略小得细节。

网格生成技术：扫略网格（网格在区域得一个表面被创建，称为源面，网格中得

节点沿着连接面，拷贝一个单元层，直到目标面，abaqus自动选择源与目标面）。结构化分网技术：使用简单得预定义得网格拓扑关系划分网格，给出了网格划分得最大控制。不同得区域可以有不同得网格划分，用不同得颜色来表示。在区域之间自动创建捆绑约束，保持区域得连接，但就是约束不就是真正得协调，精度将会受到影响。

控制网格密度与梯度：使用波前算法得三角形、四面体、四边形网格得节点与种子精确匹配；使用中轴算法得六面体或四边形网格，abaqus会调整单元得分布，但就是可以通过在边上得约束种子防止调整。分区创建了附加得边，可以对局部网格密度施加更多得控制，可以在应力集中区域细化网格。

分配单元类型：荷载与边界条件等就是基于几何体得，而不就是基于网格。

网格质量检查：限制条件包括形状比、最大最小角度与形状因子等。在消息域显示单元得总数、扭曲单元得数量、平均扭曲与最差扭曲。

有限元分析实例详解（石亦平）

Abaqus有多个模块，包括cae前处理模块、主求解器Standard and explicit 、design，aqua，foundation接口等等。在step中若选择static general 则选择了standard，若选择dynamic 则选择了explicit。

ABAQUS/standard 就是一个通用分析模块，它使用隐式求解方法，能够求解广泛领域得线性与非线性问题，包括静态分析、动态分析，以及复杂得非线性耦合物理场分析等。

ABAQUS/EXPLICIT ，用以进行显式动态分析，她使用显式求解方法，适于求解复杂非线性动力学问题与准静态问题，特别就是用于模拟短暂、瞬时得动态事件，如冲击与爆炸问题。此外，它对处理接触条件变化得高度非线性问题也非常有效(例如模拟成形问题)。

二维平面应力问题：2D planar

线性摄动分析步（linear pertuibation step）：只用于分析线性问题，explicit中不能使用此。Standard中，以下分析总就是线性得：buckle（特征值屈曲）frequency （频率提取分析）modal dynamic（瞬时模态动态分析）random response （随机响应分析）response spectrum （反应谱分析）steady-state dynamics （稳态动态分析）如模型只能中存在大位移或转动，几何非线性参数NLGEOM应选择ON

设置求解过程时间增量步：若模型中不包含阻尼或与速率有关得材料性质，时间没有实际意义。允许得最小增量步：e-5，最大：1允许得增量步最大数目：100 设定输出数据：step 下output 菜单项

场变量输出结果（field output）一个分析步结束时输出结果历史变量输出结果（history output）0、1个分析步结束输出一次应力结果

设定自适应网格：step—other---adaptive mesh domain （control）通常比纯拉个狼日分析更稳定，高效，精确。

控制分析过程：standard 通用分析步step—other—general solution controls控制收敛算法与时间积分精度。静力问题，other—solver controls来控制迭代线性方程求解器得参数。

在Interaction 功能模块中，主要可以定义模型得以下相互作用。

(1)主菜单Interaction 定义模型得各部分之间或模型与外部环境之间得力学或热相互作用，例如接触、弹性地基、热辐射等

(2) 主菜单Constraint 定义模型各部分之间得约束关系。

(3) 主菜单Connector 定义模型中得两点之间或模型与地面之间得连接单元( connector)，用来模拟固定连接、钱接、恒定速度连接、止动装置、内摩擦、失效条件与锁定装置等。

(4) 主菜单Special ? Inertia 定义惯量(包括点质量/惯量、非结构质量与热容)。

(5 )主菜单Special ? Crack 定义裂纹。

(6) 主菜单Special ? Springs/Dashpots 定义模型中得两点之间或模型与地面之间得弹簧与阻尼器。

(7) 主菜单Tools 常用得菜单项包括Set (集合)、Surface (面)与AlI\plitude (幅值)等。

约束：在ABAQUS/CAE 得Assembly 功能模块、Load 功能模块与Interaction 功能模块中都有"约束"得概念，它们分别有着不同得含义。在Assembly 功能模块中，Constraint(约束)得作用就是定义各个实体间得相互位置关系，从而确定它们在装配件中得初始位置。在Load 功能模块中，主菜单BC 得作用就是定义边界条件，消除模型得刚体位移。在Interaction功能模块中，主菜单Constraint (约束)得作用就是定义模型各部分得自由度之间得约束关系，具体包括以下类型。(1)

Tie (绑定约束) 模型中得两个面被牢固地粘结在一起，在分析过程中不再分开。被绑定得两个面可以有不同得几何形状与网格。

(2) Rigid Body (刚体约束) 在模型得某个区域与一个参考点之间建立刚性连接，此区域变为一个刚体，各节点之间得相对位置在分析过程中保持不变。

(3) Display Body (显示体约束) 与Rigid Body 类似，受到此约束得实体只用于图形显示，而不参与分析过程。

(4) Coupling (耦合约束) 在模型得某个区域与参考点之间建立约束。

I) Kinematic Coupling (运动耦合) :即在此区域得各节点与参考点之间建立一种运动上得约束关系。

2) Distributing Coupling (分布耦合) :也就是在此区域得各节点与参考点之间建立一种约束关系，但就是对此区域上各节点得运动进行了加权平均处理，使此区域上受到得合力与合力矩与施加在参考点上得力与力矩相等效。换言之，分布搞合允许面上得各部分之间发生相对变形，比运动捐合中得面更柔软。

(5) Shell-to-Solid Coupling (壳体-实心体约束) 在板壳得边与相邻实心体得面之间建立约束。

(6) Embedded Region (嵌入区域约束) 模型得一个区域镶嵌在另一个区域中。

(7) Equation (方程约束) 用一个方程来定义几个区域得自由度之间得相互关系。载荷：

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) Connector Force: 施加在连接单元上得力。

14) Connector Moment: 施加在连接单元上得弯矩。

Assembly

(1)独立实体(independent instance) 独立实体就是对Part 功能模块中部件得复制，可以直接对独立实体划分网格(mesh on instance ) ，而不能对相应得部件划分网格。如果对同一个部件创建了多个独立实体，则需要对每个独立实体分别划分网格。

(2) 非独立实体(dependent instance) 非独立实体就是Part 功能模块中部件得指针(pointer) ，不能直接对非独立实体划分网格，而只能对相应得部件划分网格（mesh on part）如果对同一个部件创建了多个独立实体，则只需对部件划分一次。

格，而不必再为每个非独立实体分别划分网格。

对非独立实体，应在窗口顶部得环境栏中把object选项设为part，即对部件划分网格；反之，对独立实体划分网格，应设为assembly，对整个装配件划分网格。设置边上得种子，可以点击窗口右下角得constraints选择约束条件

无约束：节点数目可以超出或者少于种子；部分约束：只能超出，不能少于；完全约束

单元形状选择：二维（quad：完全使用四边形，quad-dominated：过渡区允许出现三角形单元，tri：完全使用三角形）

三维（hex：完全使用六面体，hex-dominated：过渡区允许出现楔形，tet：完全使用四面体，wedge：完全使用楔形）

网格颜色：structured 绿色sweep 黄色free 粉红色自由网格划分采用tri 与tet得二次单元来保证精度，structured and sweep 一般采用quad and hex ，如果定义seeds完全约束，可能划分不成功，可去除种子。

如果某个区域显示为橙色表明无法使用目前赋予它得网格划分技术来生成网格。

可把实体分割( partition) 为几个简单得区域，再划分网格。

Medial axis 算法：首先把要划分网格得区域分为一些简单得区域，然后使用结构化网格划

分技术来为这些简单得区域划分同格。

1 )使用Medial Axis 算法更容易得到单元形状规则得网格，但网格与种子得位置吻合较差。

2) 在二维模型中使用Medial Axis 算法时，选择Minimize the mesh transition (最小化网格得过渡)可以提高网格得质量，但用这种方法生成得网格更容易偏离种子。

3）如果在某些边设置了受完全约束得seeds，则该算法会自动会其她边设置最佳得种子分布。

4）不支持由cad导入得粗糙模型与虚拟拓扑（virtual topology）

Advancing Front 算法：首先在边界上生成四边形网格，然后再向区域内部扩展。

1、得到得网格可以与种子得位置很好地吻合，但在较窄得区可能会使同格歪斜。

2、容易实现从粗网格到细网格得过渡，容易得到大小均匀得网格

3、支持由cad导入得粗糙模型与虚拟拓扑（virtual topology）

检查网格质量：verify mesh

单元类型

(1)线性( linear)单元又称一阶单元，仅在单元得角点处布置节点，在各方向都采用线性插值;

(2) 二次( quadratic) 单元又称二阶单元，在每条边上有中间节点，采用二次插值;

(3) 修正得( modified) 二次单元只有Tri 或Tet 单元才有这种类型，即在每条边上有中间节点，并采用修正得二次插值。

所谓线性完全积分就是指当单元具有规则形状时，所用得高斯积分点得数目足以对单元刚度矩阵中得多项式进行精确积分。承受弯曲载荷肘，线性完全职分单元会出现剪切自锁(shear locking) 问题，造成单元过于刚硬，即使划分很细得网格，计算精度仍然很差(Getting Started wilh ABAQUS) "Element formulation and integration"

二次完全积分（quadratic full-integration）单元

计算结果精确，适合模拟应力集中问题；一般无shear locking，但不能用于接触分析；若材料不可压缩，在弹塑性分析中，容易产生volumetric locking；扭曲或弯曲应力有梯度，locking

线性缩减积分（linear reduced-integration）单元

Quad 单元与Hex 单元在ABAQUS/CAE默认得单元类型就是线性减缩积分单元减缩积分单元比普通得完全积分单元在每个方向少用一个积分点。线性减缩积分单元在单元得中心只有一个积分点，由于存在所谓"沙漏"数值问题而过于柔软，ABAQUS 在线性减缩积分单元中引入了"沙漏刚度"以限制沙漏模式得扩展。

线性减缩积分单元有以下优点，

1)对位移得求解结果较精确。

2) 网格存在扭曲变形时(例如Quad 单元得角度远近大于或小于90°)分析精度不会受到大得影响。

3)在弯曲载荷下不容易发生剪切自锁。

其缺点如下：

1)由要划分较细得网格来克服沙漏问题。

2) 如果、希望以应力集中部位得节点应力作为分析指标，则不能选用此类单元，因为线性减缩积分单元只有在单元得中心有一个积分点，相当于常应力单元，经过外差值与平均后得到得节点应力则不精确。

二次减缩积分（quadratic reduced-integration）单元

优于线性减缩积分单元，不能用于接触分析、大应变问题，精度往往低于二次完全积分单元。

非协调摸式（incompatible modes）单元得优点如下

1) 克服了剪切自锁问题，在单元扭曲比较小得情况下，得到得位移与应力结果很精确。

2) 在弯曲问题中，在厚度方向上只需很少得单元，就可以得到与二次单元相当得结而计算成本明显降低。

3)使用了增强变形梯度得非协调模式，单元交界处不会重叠或开洞，因此很容易扩展到非线性、有限应变得位移。

注意，如果所关心部位得单元扭曲比较大，尤其就是出现交错扭曲时，分析精度会降低。

综上所述，选择三维实体单元类型时应遵循以下原则。

1)对于三维区域，尽可能采用结构化网格划分技术或扫掠网格划分技术，从而得到Hex 单元网格，减小计算代价，提高计算精度。当几何形状复杂时，也可以在不重要得区域使用少量模形(Wedge) 单元。

2) 如果使用了自由网格划分技术，Tet 单元得类型应选择二次单元。在ABAQUS/Explicit

中应选择修正得Tet 单元。C3D10M ，在ABAQUS/Standard 中可以选择C3D10，但如果有

大得塑性变形，或模型中存在接触，而且使用得就是默认得"硬"接触关系( " hard" contact relationship) ，则也应选择修正得Tet 单元C3Dl0M 。

3) ABAQUS 得所有单元均可用于动态分析，选取单元得一般原则与精力分析相同。但在使用ABAQUS/Explicit模拟冲击或爆炸载荷时，应选用线性单元，因为它们具有集中质量公式，模拟应力波得效果优于二次单元所采用得一致质量公式。

如果使用得求解器就是ABAQUS/Standard，在选择单元类型时还应注意以下方面。

1) 对于应力集中问题，尽量不要使用线性减缩积分单元，可使用二次单元来提高精度。如果在应力集中部位进行了网格细化，使用二次减缩积分单元与二次完全积分单元得到得应力结果相差不大，而二次减缩积分单元得计算时间相对较短。

2) 对于弹塑性分析，如果材料就是不可压缩性得(例如金属材料)，则不能使用二次完全积分单元，否则会出现体积自锁问题，也不要使用二次Tri 单元或Tet 单元。推荐使用得就是修正得二次Tri 单元或Tet 单元、非协调单元，以及线性减缩积分单元。如果使用二次减缩积分单元，当应变超过20% -40% 时要划分足够密得网格。

3) 如果模型中存在接触或大得扭曲变形，则应使用线性Quad 或Hex 单元，以及修正得二次Tri 单元或Tet 单元，而不能使用其她得二次单元。

4) 对于以弯曲为主得问题，如果能够保证在所关心部位得单元扭曲较小，使用非协调单元(例如C3D81 单元)可以得到非常精确得结果。

5) 除了平面应力问题之外，如果材料就是完全不可压缩得(例如橡胶材料)，则应使用杂交单元;在某些情况下，对于近似不可压缩材料也应使用杂交单元。

梁单元得类型选择原则：ABAQUS 中得所有梁单元都可以产生轴向变形、弯曲变形与扭转变形，B21 与B31单元(线性梁单元)以及B22 与B32 单元(二次梁单元)既适用于模拟剪切变形起重要作用得深梁，又适用于模拟剪切变形不太重要得细长梁，三次单元B23 与B33 只需划分很少得单元就可以得到较精确得结果

1)在任何包含接触得问题中，应使用B21 或B31 单元(线性剪切变形梁单元)

2) 如果横向剪切变形很重要，则应采用B22 与B32 单元(二次Timoshenko 梁单元)。

3) 在ABAQUS/Standard 得几何非线性模拟中，如果结构非常刚硬或非常柔软，应使用

杂交单元，例如B21H 与B32H 单元。

4) 如果在ABAQUS/Standard 中模拟具有开口薄壁横截面得结构，应使用基于横截面翘

曲理论得梁单元，例如B310S 、B320S 单元。

定义耦合约束：

1、定义参考点tools—reference point（interaction module）

2、创建参考点集合tools—set—manager—create set（assembly module）

3、定义受约束得面tools—surface—manager （assembly module）

4、定义耦合约束create constraint –coupling—sets，选择点集作为耦合约束控制点；surface，选择面集作为约束面---设置couplingt type（耦合类型）为distributing （模型树中位于constraints下）

处在assembly划分网格状态下，面与集合属于整个装配件，若处在部件划分网格得状态下，则面与集合仅属于部件，不能在assembly、interaction or load module 中使用。

定义荷载：

1、定义载荷随时间变化得幅值Load模块，Tools→amplitude→Create ，

Tabular（表格) ，Continue。输入分析步时间与幅值。Time span 默认为step time （单个分析步中得时间），若为total time ，则表示所有分析步中得全部时间。（模型树：amplitudes/）在分析步中荷载以总量而不就是以增量得形式给定。

2、定义荷载面，局部荷载定义需事先在mesh module中分割面。

3、默认幅值ramp 含义：幅值从零线性增长至给定值。

定义边界条件：

1、创建集合load module tools---set—manager（点击show/Hide selection option），只选择面

2、定义边界条件BC—manager

INP在处理器与求解器之间建立了一个传递数据得桥梁

*keyword，参数，参数《abaqus keywords reference manual》

INP文件格式规则：

1、以*HEADING开头，以**开始为注释行，不能有空行，关键词、参数、集合名称、面名称不区分大小写

2、没行不超过256个字符，*ELEMENT 包含节点不超过15个，*ELSET与*NSET 数据行包含数据不超过16个，超出部分被忽略。

3、一行未结束需换行时，加逗号

4、关键词与各参数之间、数据之间均要加逗号，表明下一行就是这行延续INP 带孔方板实例

1、*PREPRINT,,echo=no，model=no，history=no，contact=no；设置dat文件中记录得内容

2、*PART,NAME=名字；非独立实体（网格划分在部件上），part数据块包含节点、单元、集合与截面属性等数据，若为独立实体，则该数据出现在*instance中。

3、*NODE，编号，坐标1，坐标2，，，不同部件与实体可以有相同得编号，在引用时需加上实体名，如PART-A-1、5，表示PART-A-1中得节点5。

4、*ELEMENT,TYPE=单元类型

单元编号，节点1编号，节点2编号。。。。。

5、节点集合*NSET单元集合*ELSET分两类

（1）、定义在part或instance数据块中，出现在*part与*end part之间，用于

定义截面属性

（2）、定义在assembly数据块中，出现在*end instance之后、*end assembly之前，用于定荷载、边界条件、面、接触或约束等。节点集合与单元集合得名称不得超过80个字符

*NEST,NSET=名称，GENERATE

起点编号，终点编号，编号增量

6、*solid section，ELSET=单元集合名称，material=材料名称

<截面参数>截面参数可以就是二维模型得厚度或一维模型得截面面积材料名称不超过80个字符，必须以字母开头

7、*assembly，name=名称9、assembly中得集合与part中基本相同，只

……………需加上INSTANCE=名称*end assembly

8、*instance,name=名称，part=名称10、*surface，type=面类型，name=面名称

……………、构成此面得集合1，名称1

*END INSTANCE 面类型默认值为ELEMENT

以下数据块出现在*END ASSEMBLY之后

11、*material，name=材料名称12*boundary 边界可创建在initial 与后续step中，

*elastic 荷载只能创建在后续step中

弹模，泊松比表示方法1 *boundary

*density 节点编号或集合，约定得边界条件类型

表示方法2 *boundary

节点编号或集合，第一个自由度编号，最后一个自由度编号，位移值**（自由度1-6）

13、*step，name=名称

*static

初始增量步，分析步时间，最小增量步，最大增量步

14、（1）集中荷载

*cload

节点编号或集合，自由度编号，荷载值

（2）定义在单元上得分布荷载

*Dload

单元编号或单元集合，荷载类型代码，荷载值

（3）定义在面上得分布荷载

*dsload

面名称，荷载类型代码，荷载值

荷载类型代码见《abaqus analysis user’s manual》

使用文本编辑软件修改inp文件，不会影响模型得数据库。可采用一下得几种方法

1、为修改后得inp创建分析作业，create Job对话框中将Source设为Input inp。

2、将INP文件导人ABAQUS/CAE、从而创建一个新得模型。—model。Inp文件中不包含模型得几何信息。若inp中包含cae不支持得关键词，则可能无法导入。

User’s manual中有“keywords support from the input file”查询哪些keywords 不可

3、ABAQUS Command窗口中输入命令Abaqus job=

4、使用Edit Keywords功能来修改INP文件，Model→edit Keyword，能保存文件修改，但就是不能真正改变模型数据库。

查瞧分析过程信息

在分析过程中生成得STA文件、MSG文件与DAT文件包含着更完整得分析信息。参见ABAQUS Analysis User's Manual~第4、1、1节"Output"。在进行非线性分析时(例如接触分析与弹塑性分析)，往往会出现不收敛得问题，此时上述文件中得信息就是查找模型问题得重要依据。ABAQUS对各个文件得处理过程如下。

1、对INP文件进行预处理，打开Windows任务管理器，可以瞧到名为pre、exe得进程。预处理过程中出现得错误信息(ERROR)与警告信息(WARNING)会显示

在DAT文件中。

2、如果在DAT文件中出现了错误信息，说明在INP文件中存在严重得错误，ABAQUS不会开始分析计算。

3、如果INP文件中没有错误，ABAQUS就会开始分析。在Windows任务管理器中会出现相应得进程，进程名为Standard、exe；Explicit、exe。

4、如果ABAQUS/Standard在分析过程中发现问题，会在MSG文件中显示相应得错误信息或警告信息。另外各个时间增量步得迭代过程也将显示在MSG文件中。

5、ABAQUS/Explicit会在STA文件中详细地列出分析过程信息、ABAQUS/Standard 只就是在STA文件中简要列出已完成得分析步与迭代收敛悄况

6、ABAQUS/Standard在MSG文件中详细列出与迭代收敛有关得参数设置与分析过程

7、ABAQUS/Standard会在DAT文件得后半部分显示用户所要求输出得分析结果，以及模型得规模、求解所占用得内存与磁盘空间、分析所周时间等内容。如果在INP文件得Step数据块中使用*NODE PRINT或*EL PRINT等关键词，就可以将节点或单元得分析结果输出到DAT文件中，*NODE PRINT，NSET=Set-Point u，8、abaqus运行环境得设置，安装目录下site中环境文件abaqus-v6、env，可修改参数，详细见abaqus installation and license guide----4、1与user’s manual ------3、4、1

接触分析

1、非线性问题三种类型：材料非线性material nonlinearity（应力应变关系），几何非线性geometric nonlinearity（位移得大小对结构得响应发生影响，如大位移、大转动、初始应力、几何刚性化与突然翻转(snap through)），边界条件非线性boundary nonlinearity（边界条件在分析过程中发生变化，如接触问题）

2、ABAQUS/Standard使用Newton-Raphson算法来求解非线性问题，它把分析过程划分一系列得载荷增量步，在每个增量步内进行若干次迭代(iteration)，得到可接受得解后，再求解下一个增量步，所有增量响应得总与就就是非线性分析得近似解

ABAQUS/Explicit不需要进行迭代，而就是显示地从上个增量步得静力学状态来推

出动力学平衡方程得解。ABAQUS/Explicit得求解过程需要大量得增量步，但由于不进行迭代，也不需要求解全体方程组，其每个增量步得计算成本很小，可以很高效地求解复杂得非线性问题。

3、若部件刚度大，且变形、应力不就是重点，则可将其当作刚体，减小模型规模。

（1）、创建刚体：part---analytical rigid，，，

（2）、tools--- reference point，边界条件与荷载均施加在reference point上，RP 黄色

（3）ASSEMBLY---instance part mesh----设置单元类型，，刚体部件不需划分网格与设置单元类型，也不需要材料与截面属性

（4）step---initial：定义边界（仅定义，加载在load模块中进行），接触分析中，建议先施加较小荷载，建立接触关系，再施加真实荷载，容易收敛。每个step 后均可设置场变量与历史变量变量输出

（5）定义接触面interaction---tools---surface—manager，，接触对得法线方向相反，指向实体外部。

（6）定义接触属性creat interaction property

（7）定义接触creat interaction，选择主面master surface，从面slave surface。刚体单元得面必须就是主面。Sliding formulation：finite sliding （有限滑移）small sliding

（8）边界与荷载：load---tools---set创建参考点集合，在参考点上定义边界BC，加载,此时需选择step

4、主要问题

1、三类接触面：由单元构成得柔体接触面或刚体接触面，由节点构成得接触面，解析刚体截面面，，，一个接触对（contact pair）最多只能有一个由节点构成得接触面。若只有一个接触面，称为自接触（self contact）

2、两种算法：通用接触算法，接触对算法（需指定接触面）

3、接触方向总就是主面得法线方向，从面上得节点不会穿越主面，但主面上得节点可以穿越从面。主面选择原则：刚度大，网格较粗，主面不能就是由节点构成得面，并且必须就是连续得，如果就是有限滑移(finite sliding)，主面

在发生接触得部位必须就是光滑得(即不能有尖角)。如果接触面在发生接触得部位有很大得凹角或尖角，应该将其分别定义为两个面。如果就是有限滑移(finite sliding)，则在整个分析过程中，都尽量不要让从面节点落到主面之外(尤其就是不要落到主面得背面)，否则容易出现收敛问题。一对接触面得法线方向应该相反，如果法线方向错误，ABAQUS往往会将其理解为具有很大过盈量得过盈接触，因而无法达到收敛。*CONTACT PAIR，INTERACTION=接触属性

从面名称，主面名称

4、有限滑移（finite sliding）两个接触面之间可以有任意得相对滑动。这就是定义接时得默认特性。其关键词为*CONTACT PAIR，INTERCTION=<接触属性〉

〈从面名称>，<主面名称〉

Standard需不断判定从面节点与主面哪部分发生接触，要求主面就是光滑得，否则不容易收敛

5、小滑移（small sliding）滑动量大小只就是单元尺寸得一小部分

*CONTACT PAIR，INTERCTION=<接触属性〉,SMALL SLLIDING

〈从面名称>，<主面名称〉

在分析开始就确定了从面节点与主面得哪些部分发生了接触，接触关系不会改变，小滑静也可以用于儿何非线性问题(即使用、STEP、NLGEOM 定义得分析步)，并考虑面得大转动与大变形，更新解除力得传递路径。如果在模型中没有几何非线性，则忽略面得转动与变形，载荷得路程保持不变。

小滑移问题得接触压强根据未变形时得接触面积来计算，有限滑移则就是根据变化得接触面积来计算。

6、不做设置，abaqus自动根据模型中主面与从面得距离判断接触状态

*CONTACT PAIR，INTERCTION=<接触属性名称〉，adjust=位移误差限度

〈从面名称>，<主面名称〉

位移误差限度：若从面节点与主面得距离小于此限度，则调整节点初始坐标，使其与主面得距离为0

7、Contact property 两部分：接触面之间得法向作用与切向作用。对

于法向作用，ABAQUS中接触压力与间隙得默认关系就是硬接触(hard contact)，其接触面之间能够传递得接触压力得大小不受限制；当接触压力变为0或负值时，两个接触面分离，并且去掉相应节点上得接触约束。另外，ABAQUS还提供了多种软接触(soft contact)，包括指数模型、表格模型、线性模型等。对于切向作用，ABAQUS中常用得摩擦模型为库伦摩擦，默认得摩擦系数为0。切向力达到切应力之前，摩擦面之间不会相对滑动。τ=μ×P，P为法向接触压强(CPRESS)摩擦类型参见：abaqus analysis user’s manual-----22、1、4—friction behavior 接触属性定义*surface interaction，name=接触属性名称

*friction

<摩擦系数> Cae—interaction—property—create，machanical—tangential behavior设定摩擦，machanical—normal behavior设定法向作用类型

8、接触信息得输出：*contact print，slave=。。，master=。。，nset=。。

<结果变量名称>。

CFN：接融压力得合力CFS：摩擦应力得合力CAREA：接触面积CMN：接触压力相对于原点得合力矩CMS：摩擦应力相对于原点得合力矩CFT：接触压力与摩擦应力得合力CMT：接触压力与摩擦应力相对于原点得合力矩。

若未给定结果变量名称，则输出以下从面节点得变量结果于dat中

Status（接触状态），cpress（接触压强），cshear1（在局部方向1上得摩擦剪应力）

copen（从面节点与主面得距离），cslip1（在局部方向1上得相对切向滑移：各增量步中滑移得总与）

CPRESS就是从面上各个节点上各自得接触压强，CFN代表接触面所有节点接触力得合力。接触面上所有节点在垂直于接触面方向上接触力得合力称为法向接触力，若接触面就是曲面，无法由CFN直接得到法向接触力，这时可以通过各个从面节点得CPRESS来计算向接触力。

法向接触力=从面节点上得CPRESS之与×从面面积/从面上得节点数

摩擦力=法向接触力×摩擦系数

9、迭代过程问题：如果当前得时间增量步无法在规定得迭代次数内达到收敛ABAQUS会自动减小时间步，重新开始迭代即Cutback。如果仍不能收敛，则会继续减小时间增量步，如果达到了现定得Cutback最大次数(默认值为5次)或时间增量步长减小到所规定得最小限度(默认值为10-5)ABAQUS就会中止分析。参见Analysis User's Manual—8、3--Analysis convergence controls。visualization module—tools—job diagnostics，可查瞧收敛过程得诊断信息。

从面节点有开放（open）与闭合（closure）两种接触状态。如果在一次迭代中节点得接触状态发生了变化，则称之为"严重不连续迭代"(severe discontinuity iteration)。在MSG文件中显示了接触状态发生变化得节点数日(例如"0 CLOSURES,10 OPENINGS")

如果分析能够收敛，每次严重不连续迭代中CLOSURES与OPENINGS得数目会逐渐减少。当所有从面节点得接触状态都不再发生变化，就进入平衡迭代，最终达到收敛。当closures与openings数目时而减小时而增大，则可以尝试减小时间增量步；当closures与openings数目在0与1之间不断变换，意味着一个从面节点得接触状态不断在打开与闭合，所谓得振颤（chattering），这时无法通过减小时间步来达到收敛。

当closures与openings数目减小速度很慢时，达到第12次严重不连续迭代之后，abaqus会自动缩小时间步长，重新开始迭代。此时，可以改变迭代参数得设置*controls，parameters=time incrementation

，，，，，，25，，，，，

Cae：step module---other---general solution controls---edit，specify--time incrementation---more

若想瞧更详细信息，step module—output---diagnostic print

*PRINT,CONTACT=YES

10、收敛问题解决办法：Analysis User's Manual—common difficulties associated with contact modeling in abaqus

A、检查接触关系、边界条件与约束打开ODB文件或者进入visualization模块，view---ODB display options，entity display—show boundary conditions与show coupling constraints。

B、消除刚体位移：查瞧ODB文件，visualization module---tools---job diagnostic，选中highlight selection in viewport，可以显示出现numerical singularity得节点。也可以利用接触或摩擦来约束刚体位移，可以在接触对上设置微小过盈量，保证在分析一开始就已经建立起接触关系，也可以施加临时边界条件，还可以在实体上得任意一点与地面之间定义一个软弹簧，约束刚体位移。Interaction module，special---spring/dashpots—creat，弹簧类型为connect points to ground，选节点，将degree of freedom设为出现了刚体位移得自由度，将spring stiffness设为一个较小值。设置完之后继续求解，若刚开始出现warings，后面不再出现，则ok

C、使用绑定约束：如果某一接触对得接触状态对整个模型得影响不大，或者这一对接触面在整个分析过程中都就是始终紧密接触得，可以考虑将它们之间得接触关系改为绑定约束，这样有助于消除刚体位移，并且大大减少计算接触状态所需要得迭代。CAE:interaction module，interaction—constraints，将type设为tie。若使用了绑定约束或者定义了过盈约束，则必须让位置误差限度略大于主面与从面在模型中得距离。

D、平稳建立接触关系，先定义一个很小荷载得分析步，在下个分析步中加真实荷载

E、细化网格：粗糙得网格会使ABAQUS难以确定接触状态，例如，如果在接触面得宽度方向上只有一个单元，则常常会出现收敛问题。一般来说，如果从面上有90°得圆角，建议在此圆角处至少划分10个单元。

F、如果接触属性为默认得"硬"接触，则不能使用六面体二次单元(C3D20与C3D20R)，以及四面体二次单元(C3D10)，而应尽可能使用六面体一阶单元。如果无法划分六面体单元网格，可以使用修正得四面体二次单元(C3D10M) ABAQUS Analysis User's Manual--Defining contact pairs in ABAQUS/Standard。

G、避免过约束（overconstraint）：接触:从面节点会受到沿主面法线方向得约束边界条件。连接单元(connector)子模型边界(*SUBMODEL)各种约束，例如耦合约束(*COUPLING)、刚体约束(*RIGID BODY)、绑定约束(*TIE)、旋转周期对称约束(*TIE,CYCLIC SYMMETRY)、多点约束(*MPC)、线性方程约束(*EQUATION)等。例如，如果在节点上同时定义了绑定约束与边界条件，或者既约束了沿切向得位移，又定义了使用Lagrange摩擦或粗糙摩擦得接触关系，都会造成过约束。对于一些