abaqus中单元的选择宝典

1、1. 完全积分就是指当单元具有规则形状时， 所用得高斯积分点可以对单元刚度 矩阵中得多项式进行精确地积分。 2. 剪力自锁将使单元变得“刚硬” ，只影响受弯曲荷载得完全积分线性（一阶） 单元，这些单元功能在受直接或剪切荷载时没有问题。 二次单元得边界可以 弯曲，没有剪力自锁得问题。 3. 只有四边形与六面体单元才能采用减缩积分。所有得楔形、四面体与三角形 实体单元采用完全积分。 减缩积分单元比完全积分单元在每个方向上少用一 个积分点。 4. 只有四边形与六面体单元才能采用减缩积分。所有得楔形、四面体与三角形 实体单元采用完全积分。 减缩积分单元比完全积分单元在每个方向上少用一 个积分点。 5.

2、 非协调单元：只有四边形与六面体单元才能采用减缩积分。所有得楔形、四 面体与三角形实体单元采用完全积分。 减缩积分单元比完全积分单元在每个 方向上少用一个积分点。 6. ABAQUS对非协调单元采用了增强位移梯度形式。在弯曲问题中，用非协调 单元可得到与二次单元相当得结果， 且计算费用明显降低。 对单元扭曲很敏 感。 7. ABAQUS对非协调单元采用了增强位移梯度形式。在弯曲问题中，用非协调 单元可得到与二次单元相当得结果， 且计算费用明显降低。 对单元扭曲很敏 感。 8. 杂交单元：ABAQUS对非协调单元采用了增强位移梯度形式。在弯曲问题中， 用非协调单元可得到与二次单元相当得结果， 且

3、计算费用明显降低。 对单元 扭曲很敏感。 9. 一般情况下应采用二次减缩积分单元（ CAX8R， CPE8R， CPS8R， C3D20R）。 在应力集中局部采用二次完全积分单元（ CAX8 CPE8 CPS8 C3D20。对 含有非常大得网格扭曲模拟（大应变分析） ，采用细网格划分得线性减缩积 分单元（ CAX4R， CPE4R，CPS4R，C3D8R ）。对接触问题采用线性减缩积分 单元或非协调单元（CAX4I CPE4J CPS4IJ C3D8I等）得细网格划分。 10. 采用非协调单元时应使网格扭曲减至最小。 三维情况应尽可能采用块状单元 （六面体）。对小位移冋题米用二次四面体单兀（C

4、3D10）就是可行得。 11. 在实体单元中所用得数学公式与积分阶数对分析得精度与花费有非常显著 得影响。使用完全积分得单兀，尤其就是一阶（线性）单兀，容易形成自锁 现象，在正常情况下不要应用。 一阶减缩积分单兀容易出现沙漏现象； 充分 得单兀细化可减小这种冋题。 12. 在分析中如有弯曲位移，且米用一阶减缩积分单兀时，应在厚度方向上至少 用 4 个单兀。沙漏现象在二阶减缩积分单兀中较少见。 在大多数一般冋题中 要考虑应用这些单兀。非协调单兀得精度依赖于单兀扭曲得量值。 13. 结果得数值精度依赖于所用得网格。 应进行网格细化研究已确保该网格对冋 题提供了唯一得解答。 但就是应记住使用一个收敛

5、网格不能保证计算结果与 冋题得行为相匹配： 它还依赖于模型其它方面得近似化与理想化程度。 通常 只在想要得到精确结果得区域细化网格。ABAQUS具有一些先进得特点，如 子模型，它可以帮助对复杂模拟得到有用得结果。 HOURGLASS 基础部分 Part 类型:可变形部件，离散刚体部件（任意形状，荷载作用下不可变形），解析 刚体部件（只可以用直线，圆弧与抛物线创建得形状，荷载作用下不可变形） 。 每个部件只存在自己得坐标系中， 与其她部件无关。 给部件赋予属性， 既成为实 例。实例可以装配成 assembly。 Automated repair opti ons：默认为缝合边，自动修理用于几何体

6、变成 valid。 基特征一旦创建不能修改。 附加特征可用于修改基特征或为基特征添加细节 （拉 伸，壳，线，切削，导角） 基准几何体类型：点，轴，坐标系，平面。 过滤器： selection options 分区：细分为不同得区域 对于拉伸与旋转，有扭曲选项，可以创建螺纹、螺旋弹簧与扭曲线。也可以利用 锥度选项，指定角度，创建带有锥度得部件。 导入孤立网格：通过、inp与、odb文件导入已有网格。被导入得孤立网格，没 有父几何体。 定义表面增强：定义了连接到已有部件表面得表面，并指定她得工程属性。 如何给部件定位：相对定位：定义几何关系，确定规则， 表面平行约束，面面 平行约束，共轴约束，接触

7、约束，重合点约束，平行坐标系约束， 若定义有冲突， 则将之前得相对约束转化为绝对约束。 集与表面在 assembly，step，interaction 与 load 模块中均有效。在 part or property module 中 创 建 得 part 集 在 assembly module 中 有 效 ， 但 不 能 通 过 set managerment 修改。 Step用途：define step,指定输出需求，指定分析诊断，指定分析控制。接触、 荷载与边界条件就是分析步相关得，需事先定义。主要用于描述模拟历程。对 python与C+呆留了 API接口，用于后处理。输出类型有两种类型

8、：场数据用于 绘制模型得变形，云图与X-Y图；历程数据用于X-丫绘图。分析步可替换。分析 控制：为显式分析定义自适应网格区域与控制； 为接触问题定制求解控制； 定制 一般得求解控制。 In teraction :用于模拟机械或热得接触。如定义边界得耦合，定义连接器。显示 体得目得就是可视化，不用于分析。接触模型得法向关系、摩擦与干涉。带有摩 擦得双面接触、自接触、捆绑约束。使用步骤：create，选择起作用得step；选 择表面；在 edit interaction 对话框中完成接触定义；在接触管理器中激活或不激 活。 边界条件：包括初始温度、指定得平移或转动，速度或角速度。指定得边界条件 可

9、以随着时间相关得幅值定义。 初始条件：包括平动与转动速度、 温度。初始平动速度可以模拟自由落体得效果。 步骤：创建、指定对象、编辑。 Mesh module ：分网技术，单元形状，单元类型，网格密度，生成网格，检查网 格状况。 二维区域可用形状：四边形、以四边形为主（允许三角形单元作为过度） 、三角 形 三维区域可用形状：四面体、若实例中包含虚拟拓扑，可使用三角形单元、四边 形单元与利用波前算法得四边形或四边形为主得单元。 细节模型中，小得细节可能会影响网格效果，虚拟模型则忽略小得细节。 网格生成技术：扫略网格（网格在区域得一个表面被创建，称为源面，网格中得 节点沿着连接面， 拷贝一个单元层，

10、 直到目标面， abaqus 自动选择源与目标面） 。 结构化分网技术： 使用简单得预定义得网格拓扑关系划分网格， 给出了网格划分 得最大控制。 不同得区域可以有不同得网格划分， 用不同得颜色来表示。 在区域 之间自动创建捆绑约束， 保持区域得连接， 但就是约束不就是真正得协调， 精度 将会受到影响。 控制网格密度与梯度： 使用波前算法得三角形、 四面体、 四边形网格得节点与种 子精确匹配；使用中轴算法得六面体或四边形网格， abaqus 会调整单元得分布， 但就是可以通过在边上得约束种子防止调整。 分区创建了附加得边， 可以对局部 网格密度施加更多得控制，可以在应力集中区域细化网格。 分配单

11、元类型：荷载与边界条件等就是基于几何体得，而不就是基于网格。 网格质量检查： 限制条件包括形状比、 最大最小角度与形状因子等。 在消息域显 示单元得总数、扭曲单元得数量、平均扭曲与最差扭曲。 有限元分析实例详解（石亦平） Abaqus有多个模块，包括cae前处理模块、主求解器Standard and explicit、design, aqua， foundation 接口等等。在 step 中若选择 static general 则选择了 standard， 若选择 dynamic 则选择了 explicit。 ABAQUS/standard 就是一个通用分析模块， 它使用隐式求解方法， 能够

12、求解广泛 领域得线性与非线性问题， 包括静态分析、 动态分析， 以及复杂得非线性耦合物 理场分析等。 ABAQUS/EXPLICIT用以进行显式动态分析，她使用显式求解方法，适于求解复 杂非线性动力学问题与准静态问题，特别就是用于模拟短暂、瞬时得动态事件， 如冲击与爆炸问题。 此外，它对处理接触条件变化得高度非线性问题也非常有效 （例如模拟成形问题 ）。 二维平面应力问题： 2D planar 线性摄动分析步（linear pertuibation step）:只用于分析线性问题，explicit中不 能使用此。Standard中，以下分析总就是线性得：buckle （特征值屈曲）freque

13、ncy （频率提取分析）modal dynamic （瞬时模态动态分析）random response （随机 响应分析） response spectrum （反应谱分析） steady-state dynamics （稳态动态 分析）如模型只能中存在大位移或转动，几何非线性参数NLGEOM应选择ON 设置求解过程时间增量步： 若模型中不包含阻尼或与速率有关得材料性质， 时间 没有实际意义。允许得最小增量步：e-5,最大：1允许得增量步最大数目：100 设定输出数据： step 下 output 菜单项 场变量输出结果( field output )一个分析步结束时输出结果 历史变量输出结果

14、 (history output )0、1 个分析步结束输出一次应力结果 设定自适应网格： stepother-adaptive mesh domain ( control )通常比纯拉个 狼日分析更稳定,高效,精确。 控制分析过程： standard 通用分析步 stepothergeneral solution controls 控制 收敛算法与时间积分精度。静力问题, other solver controls 来控制迭代线性方 程求解器得参数。 在 Interaction 功能模块中,主要可以定义模型得以下相互作用。 (1) 主菜单 Interaction 定义模型得各部分之间或模型与

15、外部环境之间得力学或热 相互作用,例如接触、弹性地基、热辐射等 (2) 主菜单 Constraint 定义模型各部分之间得约束关系。 (3) 主菜单 Connector 定义模型中得两点之间或模型与地面之间得连接单元 (connector),用来模拟固定连接、钱接、恒定速度连接、止动装置、内摩擦、失 效条件与锁定装置等。 (4) 主菜单 Special ? Inertia 定义惯量 (包括点质量 /惯量、非结构质量与热容 )。 (5 )主菜单Special ? Crack定义裂纹。 (6) 主菜单 Special ? Springs/Dashpots 定义模型中得两点之间或模型与地面之间 得弹

16、簧与阻尼器。 (7) 主菜单 Tools 常用得菜单项包括 Set (集合)、 Surface (面)与 AlIplitude (幅值) 等。 约束：在 ABAQUS/CAE得Assembly功能模块、Load功能模块与Interaction功 能模块中都有 约束得概念,它们分别有着不同得含义。在Assembly 功能模块 中，Con strai nt(约束)得作用就是定义各个实体间得相互位置关系，从而确定它 们在装配件中得初始位置。在 Load 功能模块中,主菜单 BC 得作用就是定义边 界条件，消除模型得刚体位移。在 Interaction 功能模块中，主菜单 Constraint (约

17、束)得作用就是定义模型各部分得自由度之间得约束关系，具体包括以下类型。 (1) Tie (绑定约束 ) 模型中得两个面被牢固地粘结在一起，在分析过程中不再分开 被绑定得两个面可以有不同得几何形状与网格。 (2) Rigid Body 刚( 体约束) 在模型得某个区域与一个参考点之间建立刚性连接， 此 区域变为一个刚体，各节点之间得相对位置在分析过程中保持不变。 (3) Display Body (显示体约束 ) 与 Rigid Body 类似，受到此约束得实体只用于图形 显示，而不参与分析过程。 (4) Coupling (耦合约束 ) 在模型得某个区域与参考点之间建立约束。 I) Kinem

18、atic Coupling (运动耦合 ) :即在此区域得各节点与参考点之间建立一种运 动上得约束关系。 2) Distributing Coupling (分布耦合 ) :也就是在此区域得各节点与参考点之间建立 一种约束关系， 但就是对此区域上各节点得运动进行了加权平均处理， 使此区域 上受到得合力与合力矩与施加在参考点上得力与力矩相等效。 换言之，分布搞合 允许面上得各部分之间发生相对变形，比运动捐合中得面更柔软。 (5) Shell-to-Solid Coupling (壳体-实心体约束) 在板壳得边与相邻实心体得面之间 建立约束。 (6) Embedded Region 嵌( 入区域约

19、束 ) 模型得一个区域镶嵌在另一个区域中。 (7) Equation (方程约束 ) 用一个方程来定义几个区域得自由度之间得相互关系。 载荷： 4) Shell Edge Load施加在板壳边上得力或弯矩。 5) Surface Traction: 施加在面上得单位面积载荷，可以就是剪力或任意方向上得 力，通 过一个向量来描述力得方向。 6) Pipe Pressure: 施加在管子内部或外部得压强。 7) Body Force: 单位体积上得体力。 8) Line Load: 施加在梁上得单位长度线载荷。 9) Gravity: 以固定方向施加在整个模型上得均匀加速度， 例如重力 ;ABAQ

20、US 根据 此 加速度与材料属性中得密度来计算相应得载荷。 10) Bolt Load: 螺栓或紧固件上得紧固力，或其长度得变化。 11) Generalized Plane Strain: 广义平面应变载荷，它施加在由广义平面应变单元 所构成 12) Rotational Body Force: 由于模型得旋转造成得体力、需要指定角速度或角加， 以及旋转轴。 13) Connector Force: 施加在连接单元上得力。 14) Connector Moment: 施加在连接单元上得弯矩。 Assembly (1) 独立实体 (independent instance) 独立实体就是对 P

21、art 功能模块中部件得复 制，可以直接对独立实体划分网格 (mesh on instance ) ，而不能对相应得部件划 分网格。如果对同一个部件创建了多个独立实体， 则需要对每个独立实体分别划 分网格。 (2) 非独立实体 (dependent instance) 非独立实体就是 Part 功能模块中部件得指 针 (pointer) ，不能直接对非独立实体划分网格，而只能对相应得部件划分网格 ( mesh on part) 如果对同一个部件创建了多个独立实体， 则只需对部件划分一 次。 格，而不必再为每个非独立实体分别划分网格。 对非独立实体，应在窗口顶部得环境栏中把 object选项设为

22、part,即对部件划分 网格；反之，对独立实体划分网格，应设为assembly,对整个装配件划分网格。 设置边上得种子，可以点击窗口右下角得 constraints 选择约束条件 无约束：节点数目可以超出或者少于种子；部分约束：只能超出，不能少于；完 全约束 单元形状选择：二维(quad :完全使用四边形，quad-dominated:过渡区允许出 现三角形单元，tri :完全使用三角形) 三维(hex：完全使用六面体，hex-dominated:过渡区允许出现楔形，tet :完全 使用四面体，wedge :完全使用楔形) 网格颜色: structured 绿色 sweep 黄色 free 粉

23、红色 自由网格划分采用 tri 与 tet 得二次单元来保证精度， structured and sweep 一般采用 quad and hex ， 如果定义seeds完全约束，可能划分不成功，可去除种子。 如果某个区域显示为橙色表明无法使用目前赋予它得网格划分技术来生成网格。 可把实体分割 ( partition) 为几个简单得区域，再划分网格。 Medial axis 算法：首先把要划分网格得区域分为一些简单得区域， 然后使用结构 化网格划 分技术来为这些简单得区域划分同格。 1 ) 使用 Medial Axis 算法更容易得到单元形状规则得网格，但网格与种子得位置 吻合较差。 2) 在二

24、维模型中使用 Medial Axis 算法时，选择 Minimize the mesh transition ( 最 小化网格得过渡 )可以提高网格得质量，但用这种方法生成得网格更容易偏离种 子。 3) 如果在某些边设置了受完全约束得seeds,则该算法会自动会其她边设置最佳 得种子分布。 4) 不支持由 cad 导入得粗糙模型与虚拟拓扑( virtual topology ) Advancing Front 算法：首先在边界上生成四边形网格,然后再向区域内部扩展。 1、 得到得网格可以与种子得位置很好地吻合,但在较窄得区可能会使同格歪斜。 2、容易实现从粗网格到细网格得过渡,容易得到大小均匀

25、得网格 3、支持由 cad 导入得粗糙模型与虚拟拓扑( virtual topology ) 检查网格质量： verify mesh 单元类型 (1) 线性(linear)单元又称一阶单元，仅在单元得角点处布置节点，在各方向都采 用线性插值 ; (2) 二次 ( quadratic) 单元 又称二阶单元,在每条边上有中间节点,采用二次插 值; (3) 修正得 ( modified) 二次单元 只有 Tri 或 Tet 单元才有这种类型, 即在每条边 上有中间节点,并采用修正得二次插值。 所谓线性完全积分就是指当单元具有规则形状时, 所用得高斯积分点得数目足以 对单元刚度矩阵中得多项式进行精确积

26、分。 承受弯曲载荷肘, 线性完全职分单元 会出现剪切自锁 (shear locking) 问题,造成单元过于刚硬, 即使划分很细得网格, 计算精度仍然很差 (Getting Started wilh ABAQUS) Element formulation and integration 二次完全积分( quadratic full-integration )单元 计算结果精确，适合模拟应力集中问题；一般无 shear lock ing,但不能用于接触 分析；若材料不可压缩，在弹塑性分析中，容易产生 volumetric locki ng；扭曲或 弯曲应力有梯度, locking 线性缩减积分(

27、linear reduced-integration)单元 Quad单元与Hex单元在ABAQUS/CAE默认得单元类型就是线性减缩积分单元 减缩积分单元比普通得完全积分单元在每个方向少用一个积分点。 线性减缩积分 单元在单元得中心只有一个积分点,由于存在所谓 沙漏数值问题而过于柔软, ABAQUS在线性减缩积分单元中引入了 沙漏刚度以限制沙漏模式得扩展。 线性减缩积分单元有以下优点, 1) 对位移得求解结果较精确。 2) 网格存在扭曲变形时 (例如 Quad 单元得角度远近大于或小于 90)分析精度 不会受到大得影响。 3) 在弯曲载荷下不容易发生剪切自锁。 其缺点如下： 1) 由要划分较细

28、得网格来克服沙漏问题。 2) 如果、希望以应力集中部位得节点应力作为分析指标,则不能选用此类单元, 因为线性减缩积分单元只有在单元得中心有一个积分点, 相当于常应力单元, 经 过外差值与平均后得到得节点应力则不精确。 二次减缩积分(quadratic reduced-integration)单元 优于线性减缩积分单元, 不能用于接触分析、 大应变问题, 精度往往低于二次完 全积分单元。 非协调摸式(in compatible modes )单元得优点如下 1) 克服了剪切自锁问题,在单元扭曲比较小得情况下,得到得位移与应力结果 很精确。 2) 在弯曲问题中,在厚度方向上只需很少得单元,就可以得

29、到与二次单元相当 得结而计算成本明显降低。 3) 使用了增强变形梯度得非协调模式,单元交界处不会重叠或开洞,因此很容易 扩展到非线性、有限应变得位移。 注意，如果所关心部位得单元扭曲比较大， 尤其就是出现交错扭曲时， 分析精度 会降低。 综上所述，选择三维实体单元类型时应遵循以下原则。 1) 对于三维区域，尽可能采用结构化网格划分技术或扫掠网格划分技术，从而得 到 Hex 单元网格，减小计算代价，提高计算精度。当几何形状复杂时，也可以 在不重要得区域使用少量模形 (Wedge) 单元。 2) 如果使用了自由网格划分技术， Tet 单元得类型应选择二次单元。在 ABAQUS/Explicit 中

30、应选择修正得 Tet 单元。 C3D10M ，在 ABAQUS/Standard 中可以选择 C3D10， 但如果有 大得塑性变形，或模型中存在接触， 而且使用得就是默认得 硬接触关系 ( hard contact relationship) ，则也应选择修正得 Tet 单元 C3Dl0M 。 3) ABAQUS 得所有单元均可用于动态分析，选取单元得一般原则与精力分析相 同。但在使用 ABAQUS/Explicit 模拟冲击或爆炸载荷时，应选用线性单元，因为 它们具有集中质量公式，模拟应力波得效果优于二次单元所采用得一致质量公 式。 如果使用得求解器就是ABAQUS/Standard在选择单

31、元类型时还应注意以下方 面。 1) 对于应力集中问题，尽量不要使用线性减缩积分单元，可使用二次单元来提 高精度。如果在应力集中部位进行了网格细化， 使用二次减缩积分单元与二次完 全积分单元得到得应力结果相差不大，而二次减缩积分单元得计算时间相对较 短。 2) 对于弹塑性分析， 如果材料就是不可压缩性得 (例如金属材料 )，则不能使用二 次完全积分单元，否则会出现体积自锁问题，也不要使用二次 Tri 单元或 Tet 单 元。推荐使用得就是修正得二次 Tri 单元或 Tet 单元、非协调单元，以及线性减 缩积分单元。如果使用二次减缩积分单元，当应变超过 20% -40% 时要划分足够 密得网格。

32、3) 如果模型中存在接触或大得扭曲变形，则应使用线性 Quad 或 Hex 单元，以 及修正得二次 Tri 单元或 Tet 单元，而不能使用其她得二次单元 4）对于以弯曲为主得问题，如果能够保证在所关心部位得单元扭曲较小，使用 非协调单元 （例如 C3D81 单元 ）可以得到非常精确得结果。 5）除了平面应力问题之外， 如果材料就是完全不可压缩得 （例如橡胶材料 ），则应 使用杂交单元 ;在某些情况下，对于近似不可压缩材料也应使用杂交单元。 梁单元得类型选择原则：ABAQUS中得所有梁单元都可以产生轴向变形、弯曲变 形与扭转变形，B21与B31单元（线性梁单元）以及B22与B32单元（二次梁单

33、元） 既适用于模拟剪切变形起重要作用得深梁， 又适用于模拟剪切变形不太重要得细 长梁，三次单元 B23 与 B33 只需划分很少得单元就可以得到较精确得结果 1）在任何包含接触得问题中，应使用 B21 或 B31 单元（线性剪切变形梁单元 ） 2）如果横向剪切变形很重要，则应采用B22与B32单元匚次Timoshenko梁单 元 ）。 3）在ABAQUS/Standard得几何非线性模拟中，如果结构非常刚硬或非常柔软， 应使用 杂交单元，例如 B21H与B32H单元。 4）如果在 ABAQUS/Standard 中模拟具有开口薄壁横截面得结构，应使用基于横 截面翘 曲理论得梁单元，例如 B31

34、0S 、 B320S 单元。 定义耦合约束： 1、定义参考点 toolsreference point（interaction module） 2、创建参考点集合 toolssetmanagercreate set（assembly module） 3、定义受约束得面 toolssurface manager （assembly module） 4、 定义耦合约束 create constraint - couplingsets，选择点集作为耦合约束 控制点；surface，选择面集作为约束面-设置couplingt type （耦合类型）为 distributing （模型树中位于 cons

35、traints 下） 处在 assembly 划分网格状态下，面与集合属于整个装配件，若处在部件划分网 格得状态下，则面与集合仅属于部件， 不能在 assembly、 interaction or load module 中使用。 定义荷载： 1 、定义载荷随时间变化得幅值 Load 模块， TooIl amplitude Create， Tabular（表格）,Continue。输入分析步时间与幅值。Time span默认为step time （单个分析步中得时间），若为total time，则表示所有分析步中得全部时间。（模 型树：amplitudes/）在分析步中荷载以总量而不就是以增量

36、得形式给定。 2、定义荷载面，局部荷载定义需事先在 mesh module中分割面。 3、默认幅值 ramp 含义：幅值从零线性增长至给定值。 定义边界条件： 1、创建集合 load module tools-setmanager （点击 show/Hide selection option ），只选择面 2、定义边界条件 BCmanager INP在处理器与求解器之间建立了一个传递数据得桥梁 *keyword，参数，参数abaqus keywords referenee manua INP文件格式规则： 1、以*HEADING开头，以*开始为注释行，不能有空行，关键词、参数、集合 名称、面名

37、称不区分大小写 2、 没行不超过256个字符，*ELEMENT包含节点不超过15个，*ELSET与*NSET 数据行包含数据不超过 16 个，超出部分被忽略。 3、一行未结束需换行时，加逗号 4、关键词与各参数之间、数据之间均要加逗号，表明下一行就是这行延续 INP 带孔方板实例 1、*PREPRINT,echo=no model=no，history=no，contaet=no；设置 dat 文件中记 录得内容 2、*PART,NAME名字；非独立实体（网格划分在部件上），part数据块包含节点、 单元、集合与截面属性等数据，若为独立实体，则该数据出现在 *instance 中。 3、*NO

38、DE,编号，坐标1,坐标2,不同部件与实体可以有相同得编号，在引 用时需加上实体名，如 PART-A-1 5,表示PART-A-仲得节点5。 4、*ELEMENT,TYP单元类型 单元编号，节点 1 编号，节点 2 编号。 5、节点集合*NSET单元集合*ELSET分两类 （1）、定义在part或instanee数据块中，出现在*part与*6门4 part之间，用于定义截面属性 (2)、定义在assembly数据块中，出现在*end instanee之后、*end assembly之 前，用于定荷载、边界条件、面、接触或约束等。节点集合与单元集合得名称不 得超过 80 个字符 *NEST,N

39、SET名称，GENERATE 起点编号，终点编号，编号增量 6 *solid section , ELSET单元集合名称，material二材料名称 截面参数 截面参数可以就是二维模型得厚度或一维模型得截面面积材料 名称不超过 80 个字符，必须以字母开头 7、*assembly，name二名称 基本相同，只 *end assembly 8、*instance,name=名称，part二名称 面名称 1 *END INSTANCE 以下数据块出现在*END ASSEMBL之后 11、*material，name=材料名称 与后续 step 中， *elastic 弹模，泊松比 *density

40、 型 表示方法 2 *boundary 节点编号或集合，第一个自由度编号 (自由度 1-6) 13、*step，name二名称 9、assembly中得集合与part中 需加上INSTANCE名称 10、*surface ,type二面类型，name= 构成此面得集合 1 ，名称 面类型默认值为 ELEMENT 12*boundary 边界可创建在 initial 荷载只能创建在后续 step 中 表示方法 1 *boundary 节点编号或集合， 约定得边界条件类 最后一个自由度编号，位移值* static 初始增量步，分析步时间，最小增量步，最大增量步 14、（ 1） 集中荷载 *cloa

41、d 节点编号或集合，自由度编号，荷载值 （ 2） 定义在单元上得分布荷载 *Dload 单元编号或单元集合，荷载类型代码，荷载值 （ 3） 定义在面上得分布荷载 *dsload 面名称，荷载类型代码，荷载值 荷载类型代码见 abaqus analysis user s manual 使用文本编辑软件修改 inp 文件， 不会影响模型得数据库。可采用一下得几种 方法 1、 为修改后得inp创建分析作业，create Job对话框中将Source设为In put inp。 2、将INP文件导人ABAQUS/CAE从而创建一个新得模型。一 model。Inp文 件中不包含模型得几何信息。若inp中包

42、含cae不支持得关键词，则可能无法导 入。 User s manual 中有“ keywords support from the in put file ”查询哪些 keywords 不可 3、ABAQUS Comma n窗口中输入命令 Abaqus job=INP文件得名称 4、使用Edit Keywords功能来修改INP文件， Model edit Keyword,能保存 文件修改，但就是不能真正改变模型数据库。 查瞧分析过程信息 在分析过程中生成得STA文件、MSG文件与DAT文件包含着更完整得分析信息。 参见 ABAQUS Analysis Users Manua第 4、1、1 节

43、Output。在进行非线性分析 时（例如接触分析与弹塑性分析 ），往往会出现不收敛得问题，此时上述文件中得 信息就是查找模型问题得重要依据。ABAQUS对各个文件得处理过程如下。 1、对INP文件进行预处理，打开 Windows任务管理器，可以瞧到名为 pre、 exe得进程。预处理过程中出现得错误信息（ERRORS警告信息（WARNING会显示 在DAT文件中。 2、如果在DAT文件中出现了错误信息，说明在INP文件中存在严重得错误， ABAQUS不会开始分析计算。 3、如果INP文件中没有错误，ABAQUS就会开始分析。在Windows任务管理 器中会出现相应得进程，进程名为 Standa

44、rd、 exe； Explicit、 exe。 4、 如果ABAQUS/Standard在分析过程中发现问题，会在MSG文件中显示相应 得错误信息或警告信息。另外各个时间增量步得迭代过程也将显示在MSG文件 中。 5、ABAQUS/Explici会在STA文件中详细地列出分析过程信息、ABAQUS/Standard 只就是在STA文件中简要列出已完成得分析步与迭代收敛悄况 6、ABAQUS/Standard在MSG文件中详细列出与迭代收敛有关得参数设置与 分析过程 7、ABAQUS/Standarc会在DAT文件得后半部分显示用户所要求输出得分析结果， 以及模型得规模、 求解所占用得内存与磁盘

45、空间、 分析所周时间等内容。 如果在 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 与 users manual 3、 4、1 接触分析 1、非线性问题三种类型：材料非线性 material nonlinearity （应力应变关系） ， 几何非线性

46、 geometric nonlinearity （位移得大小对结构得响应发生影响，如大位 移、大转动、初始应力、几何刚性化与突然翻转 （snap through），边界条件非线 性boundary nonlinearity （边界条件在分析过程中发生变化，如接触问题） 2、ABAQUS/Standarc使用Newton-Raphson算法来求解非线性问题，它把分析过 程划分一系列得载荷增量步， 在每个增量步内进行若干次迭代 （iteration） ，得到可 接受得解后， 再求解下一个增量步， 所有增量响应得总与就就是非线性分析得近 似解 ABAQUS/Explicit不需要进行迭代，而就是显示

47、地从上个增量步得静力学状态来推 出动力学平衡方程得解。ABAQUS/Explicit得求解过程需要大量得增量步，但由于 不进行迭代， 也不需要求解全体方程组， 其每个增量步得计算成本很小， 可以很 高效地求解复杂得非线性问题。 3、若部件刚度大，且变形、应力不就是重点，则可将其当作刚体，减小模型规 模。 （ 1）、创建刚体： part-analytical rigid， （2）、tools- referenee point，边界条件与荷载均施加在 referenee point 上, RP 黄色 （ 3） ASSEMBLY-instance part mesh设置单元类型，刚体部件不需划分网

48、格与设置单元类型，也不需要材料与截面属性 （4）step-initial :定义边界（仅定义，加载在load模块中进行），接触分析中， 建议先施加较小荷载，建立接触关系，再施加真实荷载，容易收敛。每个 step 后均可设置场变量与历史变量变量输出 （5）定义接触面interaction-tools-surfacemanager,，接触对得法线方向相 反，指向实体外部。 （ 6）定义接触属性 creat interaction property （7）定义接触 creat interaction，选择主面 master surface，从面 slave surface。 刚体单元得面必须就是主面

49、。 Sliding formulation: finite sliding （有限滑移） small sliding （ 8）边界与荷载: load-tools-set 创建参考点集合，在参考点上定义边界BC， 加载 ,此时需选择 step 4、主要问题 1、三类接触面: 由单元构成得柔体接触面或刚体接触面， 由节点构成 得接触面，解析刚体截面面，一个接触对（con tact pair）最多只能有一个由节 点构成得接触面。若只有一个接触面，称为自接触（ self contact） 2、两种算法:通用接触算法，接触对算法（需指定接触面） 3、接触方向总就是主面得法线方向， 从面上得节点不会穿越主

50、面，但 主面上得节点可以穿越从面。主面选择原则:刚度大，网格较粗，主面不能就是 由节点构成得面，并且必须就是连续得，如果就是有限滑移 （finite sliding） ，主面 在发生接触得部位必须就是光滑得 (即不能有尖角 )。如果接触面在发生接触得部 位有很大得凹角或尖角，应该将其分别定义为两个面。如果就是有限滑移 (finite slid in g)，则在整个分析过程中，都尽量不要让从面节点落到主面之外(尤其就是 不要落到主面得背面 )，否则容易出现收敛问题。一对接触面得法线方向应该相 反，如果法线方向错误，ABAQUS往往会将其理解为具有很大过盈量得过盈接触， 因而无法达到收敛。*CON

51、TACT PAIRINTERACTION接触属性 从面名称，主面名称 4、有限滑移(finite sliding)两个接触面之间可以有任意得相对滑动。 这就是定义接时得默认特性。其关键词为 *CONTACTPAIR INTERCTION=g触属 性 从面名称 ，主面名称 Standard 需不断判定从面节点与主面哪部分发生接触，要求主面 就是光滑得，否则不容易收敛 5、小滑移(small sliding)滑动量大小只就是单元尺寸得一小部分 *CONTACT PAIRINTERCTION= 触属性，SMALL SLLIDING 从面名称 ，主面名称 在分析开始就确定了从面节点与主面得哪些部分发生

52、了接触，接 触关系不会改变，小滑静也可以用于儿何非线性问题(即使用、 STEP、 NLGEOM 定义得分析步 )，并考虑面得大转动与大变形，更新解除力得传递路径。如果在 模型中没有几何非线性，则忽略面得转动与变形，载荷得路程保持不变。 小滑移问题得接触压强根据未变形时得接触面积来计算，有限滑 移则就是根据变化得接触面积来计算。 6、不做设置， abaqus 自动根据模型中主面与从面得距离判断接触状 态 *CONTACT PARINTERCTION=触属性名称，adjust二位移误差限度 从面名称 ， 主面名称 位移误差限度： 若从面节点与主面得距离小于此限度， 则调整节点 7、 初始坐标，使其

53、与主面得距离为 0 Contact property 两部分：接触面之间得法向作用与切向作用。对 于法向作用，ABAQUS中接触压力与间隙得默认关系就是硬接触（hard con tact）, 其接触面之间能够传递得接触压力得大小不受限制；当接触压力变为 0 或负值 时，两个接触面分离，并且去掉相应节点上得接触约束。另外，ABAQUS5提供 了多种软接触（soft co ntact），包括指数模型、表格模型、线性模型等。对于切向 作用，ABAQUS中常用得摩擦模型为库伦摩擦，默认得摩擦系数为0。切向力达 到切应力之前，摩擦面之间不会相对滑动。t =x P,P为法向接触压强（CPRESS） 摩擦类

54、型参见： abaqus analysis user s manual22、 1、 4friction behavior 接 触属性定义*surface interaction，name=接触属性名称 *friction 摩擦系数 Cae in teract ion property create，macha ni catangen tial behavior 设定摩擦， macha ni ca no rmal behavior设定法向作用类型 8、接触信息得输出： *contact print， slave=。， master=。， nset=。 结果变量名称 。 CFN接融压力得合力CFS摩

55、擦应力得合力CAREA接触面积 CMN 接触压力相对于原点得合力矩 CMS 摩擦应力相对于原点得合力矩 CFT 接触压力与摩擦应力得合力CMT 接触压力与摩擦应力相对于原点得 合力矩。 若未给定结果变量名称，则输出以下从面节点得变量结果于dat中 Status （接触状态），cpress （接触压强），cshear1 （在局部方向1上得摩擦剪应 力） copen （从面节点与主面得距离），cslip1 （在局部方向1上得相对切向滑移：各 增量步中滑移得总与） CPRESSt是从面上各个节点上各自得接触压强，CFN代表接触面所有节点接触 力得合力。接触面上所有节点在垂直于接触面方向上接触力得合力

56、称为法向接触 力，若接触面就是曲面，无法由 CFN 直接得到法向接触力，这时可以通过各个 从面节点得CPRES来计算向接触力。 法向接触力=从面节点上得CPRES之与x从面面积/从面上得节点数 摩擦力二法向接触力x摩擦系数 9、迭代过程问题： 如果当前得时间增量步无法在规定得迭代次数内达 到收敛ABAQUS会自动减小时间步，重新开始迭代即Cutback。如果仍不能收敛， 则会继续减小时间增量步，如果达到了现定得Cutback最大次数（默认值为5次） 或时间增量步长减小到所规定得最小限度（默认值为10-5）ABAQUS就会中止分 析。参见 Analysis Users Manual 8、3-An

57、alysis convergence controls。 visualizati on module toolsjob diag no stics，可查瞧收敛过程得诊断信息。 从面节点有开放（open）与闭合（closure）两种接触状态。如果在一次迭代中节 点得接触状态发生了变化，则称之为 严重不连续迭代 （severe discontinuity iteration）。 在 MSG 文件中显示了接触状态发生变化得节点 数日（例如0 CLOSURES,10 OPENINGS） 如果分析能够收敛，每次严重不连续迭代中CLOSURE与OPENINGS寻数目会逐 渐减少。 当所有从面节点得接触状态

58、都不再发生变化， 就进入平衡迭代， 最终达 到收敛。当closures与openings数目时而减小时而增大，则可以尝试减小时间增 量步；当closures与openings数目在0与1之间不断变换，意味着一个从面节点 得接触状态不断在打开与闭合，所谓得振颤（chattering），这时无法通过减小时 间步来达到收敛。 当closures与openings数目减小速度很慢时，达到第 12次严重不连续迭代 之后， abaqus 会自动缩小时间步长，重新开始迭代。此时，可以改变迭代参数 得设置 *controls， parameters=time incrementation ， 25， Cae

59、： step module-other-general solution controls-edit ， specify-time incrementation-more 若想瞧更详细信息， step moduleoutput-diagnostic print *PRINT,CONTACT=YES 10、收敛问题解决办法： Analysis Users Manualcommon difficulties associated with contact modeling in abaqus A、检查接触关系、边界条件与约束打开ODB文件或者进入visualization模块， view-ODB

60、display options， entity displayshow boundary conditions 与 show coupling constraints。 B、消除刚体位移：查瞧 ODB文件，visualization module-tools-job diagnostic, 选中 highlight selection in viewport，可以显示出现 numerical singularity 得节点。 也可以利用接触或摩擦来约束刚体位移, 可以在接触对上设置微小过盈量, 保证 在分析一开始就已经建立起接触关系， 也可以施加临时边界条件， 还可以在实体 上得任意一点与地面