专业ABAQUS有限元建模经验笔记
基于ABAQUS的有限元分析和应用
第一章绪论
1.有限元分析包括下列步骤：
2.为了将试验数据转换为输入文件，分析者必须清楚在程序中所应用的和由实验人员提供的材料数据的应力和应变的度量。
3.ABAQUS建模需注意以下内容：
4.对于许多包含过程仿真的大变形问题和破坏分析，选择合适的网格描述是非常重要的，需要认识网格畸变的影响，在选择网格时必须牢牢记住不同类型网格描述的优点。
第二章ABAQUS基础
1.一个分析模型至少要包含如下的信息：离散化的几何形体、单元截面属性、材料数据、载荷和边界条件、分析类型和输出要求。
①离散化的几何形体：模型中所有的单元和节点的集合称为网格。
②载荷和边界条件：
2.功能模块：
（1）Assembly（装配）：一个ABAQUS模型只能包含一个装配件。
（2）Interaction（相互作用）：相互作用与分析步有关，这意味着用户必须规定相互作用是在哪些分析步中起作用。
（3）Load（载荷）：载荷和边界条件与分析步有关，这意味着用户指定载荷和边界条件是在哪些分析步中起作用。
（4）Job（作业）：多个模型和运算可以同时被提交并进行监控。
3.量纲系统
ABAQUS没有固定的量纲系统，所有的输入数据必须指定一致性的量纲系统，常用的一致性量纲系统如下：
4.建模要点
（1）创建部件：设定新部件的大致尺寸的原则必须是与最终模型的最大尺寸同一量级。（2）用户应当总是以一定的时间间隔保存模型数据（例如，在每次切换功能模块时）。（3）定义装配：
在模型视区左下角的三向坐标系标出了观察模型的方位。在视区中的第2个三向坐标系标出了坐标原点和整体坐标系的方向（X，Y和Z轴）。
（4）设置分析过程：
（5）在模型上施加边界条件和荷载：
用户必须指定载荷和边界条件是在哪个或哪些分析步中起作用。
所有指定在初始步中的力学边界条件必须赋值为零，该条件是在ABAQUS/CAE中自动强加的。
在许多情况下，需要的约束方向并不一定与整体坐标方向对齐，此时用户可定义一个局部坐标系以施加边界条件。
在ABAQUS中，术语载荷通常代表从初始状态开始引起结构响应发生变化的各种因素，包括：集中力、压力、非零边界条件、体力、温度（与材料热膨胀同时定义）。
可产生运动的方向称为自由度。在ABAQUS中平移和转动自由度的标识如下图所示。
（6）模型的网格划分：
（7）材料定义：
5.比较显式与隐式过程：
6.在隐式与显式分析之间选择：对于求解光滑的非线性问题，ABAQUS／Standard更有效；对于波的传播分析，ABAQUS／Explicit更有效。
7.在隐式与显式分析中网格加密的成本：
磁盘空间与内存需求与单元数目成正比，与单元尺寸无关；应用隐式方法对于许多问题的计
算成本大致与自由度数目的平方成正比；只要网格是相对均匀的，随着网格密度的增加，显式方法比隐式方法会节省大量的计算成本。
8.小结
第三章有限单元和刚性体
有限单元是可变形的，而刚性体在空间运动不改变形状。任何物体或者物体的局部均可以定义作为刚性体；大多数单元类型都可以用于刚性体的定义。刚性体比变形体的优越性在于对刚性体运动的完全描述只需要在一个参考点上的最多6个自由度。当该变形可以忽略或者对其不感兴趣时，将模型中的一个部分作为刚性体可以极大地节省计算时间，而不影响整体结果。
3.1有限单元
3.1.1 单元的表征
每一个单元表征如下：
①单元族：单元族名字中的第一个字母或者字母串表示该单元属于哪一个单元族。
S—壳(shell)单元，B—梁(beam)单元，C—实体(continuum)单元，R—刚性(rigid)单元
②自由度：在ABAQUS中使用的关于自由度的顺序约定如下：
③节点数目与插值阶数：
④数学描述：具有其他可供选择的公式描述的单元由在单元名字末尾的附加字母来识别，如杂交单元在其名字末尾加“H”字母标识。
⑤积分：对于ABAQUS中的一些实体单元，可以选择应用完全积分或者减缩积分。ABAQUS 在单元名字末尾采用字母“R”来标识减缩积分单元（如果一个减缩积分单元同时又是杂交单元，末尾字母为RH)。
3.1.2 实体单元
在ABAQUS中，应尽可能地使用六面体单元或二阶修正的四面体单元。
平面应变单元可以用来模拟厚结构；平面应力单元可以用来模拟薄结构；无扭曲的轴对称单元（属于CAX类单元）可模拟360度的环，适合于分析具有轴对称几何形状和承受对称载荷的结构。
所有的实体单元必须赋予截面性质，它定义了与单元相关的材料和任何附加的几何数据。
3.1.3壳单元
一般性目的的壳单元和带有反对称变形的轴对称壳单元考虑了有限的膜应变和任意大转动。所有的四边形壳单元（除了S4）和三角形壳单元S3/S3R均采用减缩积分。而S4壳单元和其他三角形壳单元则采用完全积分。
3.1.4梁单元
所有梁单元必须提供梁截面性质，定义与单元有关的材料以及梁截面的轮廓(即单元截面的几何)。
所有的二维梁单元仅采用轴向的应力和应变。
3.1.5
桁架单元适合于模拟铰接框架结构。
输出轴向的应力和应变。
3.2刚性体
3.2.1确定何时使用刚性体
刚性体可以用于模拟非常坚硬的部件，这一部件或者是固定的，也可以进行任意大的刚体运动。它还可以用于模拟在变形部件之间的约束，并且提供了指定某些接触相互作用的简便放法。使模型的一部分成为刚性体有助于达到验证模型的目的。将部分模型表示为刚性体而不
是变形的有限单元体有利于提高计算效率。在ABAQUS/Explicit分析中，刚性体和部分刚性体单元并不影响整体时间增量，也不会显著影响求解的整体精度。（解析刚性表面可以十分光滑，而离散刚性体本身有很多面。）
3.2.2刚性体部件
3.2.3刚性单元
控制刚性体的规则，适合于所有组成刚性体的单元类型，包括刚性单元。
刚体单元没有单元输出变量，仅输出节点的运动。另外，可以输出在刚性体参考点处的约束反力和反力矩。
3.5小结
第四章应用实体单元
模拟的精度很大程度上依赖于在模型中采用的单元类型。
4.1单元的数学描述与积分
4.1.1完全积分（剪力自锁引起单元在弯曲时过于刚硬。）
4.1.2减缩积分
只有四边形和六面体单元才能采用减缩积分方法，而所有的楔形体、四面体和三角形实体单元可以采用完全积分。减缩积分单元比完全积分单元在每个方向少用一个积分点。
线性减缩积分单元受弯时存在沙漏数值问题，故需合理地才有细划的网格，才可以给出可接受的结果。建议当采用这类但愿模拟受弯载荷的任何结构时，沿厚度方向上至少应采用四个单元。
除了包含大应变的大位移模拟和某些类型的接触分析之外，二次减缩积分单元一般是最普遍的应力/位移模拟的最佳选择。非协调模式单元对单元的扭曲很敏感，否则应考虑应用二次减缩积分单元。
当材料行为是不可压缩（泊松比=0.5）或非常接近于不可压缩（泊松比》0.475）时，采用杂交单元。
4.2选择实体单元
建议（适用于ABAQUS/Standard和ABAQUS/Explicit):
应用ABAQUS/Standard还需考虑如下建议：
4.3建模要点
Part：大致尺寸：取部件最大尺寸的2倍。
Property：变形体的每个区域必须指定一个包含材料定义的截面属性。
Load：
Smooth Step Amplitude类型的幅值曲线经常用于准静态分析中。
Mesh：在设计网格时，首先需要考虑所采用的单元类型，其次是想从模拟中得到什么类型的结果。ABAQUS/CAE提供了三种类型的网格生成技术：
1)结构化网格划分(绿色表示区域)：应用于特定模型拓扑，将其应用于复杂模型时需将模型分割成简单的区域。
2)扫掠网格划分(黄色表示区域)：只限于具有特殊拓扑和几何体的模型。
3)自由网格划分(粉红色表示区域)：可以用于任意的模型形状。
注：橙黄色表示不能使用默认的单元类型生成网格的区域，必须被进一步分割。Interaction：定义接触面时，选择接触那一侧的箭头。无摩擦接触是ABAQUS中的默认接触属性。
后处理：研究伪应变能(artificial strain energy)室控制沙漏变形所耗散的主要能量。如果伪应变能过高(控制伪应变能与实际应变能的比率低于5%)，说明过多的应变能被用来控制沙漏变形了。
作用在单独节点上的载荷、边界条件或者接触约束易于产生沙漏现象。
4.4网格收敛性
粗糙的网格足以用来预测趋势和比较不同概念相互之间的表现如何不同。一般来说，没有必要对所分析的结构全部采用均匀的细划网格。应该在出现高应力梯度的地方采用细划网格，而在低应力梯度或应力量值不被关注的地方采用粗网格。建议：在开始时使用粗网格以识别高应力区的的位置，然后在该区域细分网格。
4.5小结
第五章应用壳单元
厚度小于典型结构尺寸的1/10，可以用壳单元模拟，否则用实体单元。
5.1单元的几何尺寸
在ABAQUS中具有两种壳单元：常规的壳单元(适用于薄壳问题)和基于连续体的壳单元(适用于模拟接触问题）。
描述壳体的横截面必须要定义壳体的厚度，此外还要选择是在分析过程中还是分析开始时计算截面的刚度。对于线弹性材料壳，选择在分析开始时计算材料刚度更为有效。（如果选择仅在模拟开始时计算横截面刚度，那么材料行为必须是线弹性的。）
5.2厚壳或薄壳
厚壳问题假设横向剪切变形对计算结果有重要的影响，而薄壳问题假设横向剪切变形非常小，可以忽略不计。判断建议见P115
5.3壳的材料方向
5.4选择壳单元
5.5建模要点后处理：
5.8小结
第六章应用梁单元
横截面的尺度小于结构典型轴向尺度的1/10，可以用梁单元模拟，否则用实体单元。
6.1梁横截面的几何形状
当梁的轮廓与梁的截面特性相关时，可以指定是在分析过程中计算界面的工程性质还是让ABAQUS预先计算它们(在分析开始时)。选择前者可以应用于线性或者非线性的材料行为；选择后者只适用于线弹性材料行为。
6.1.2横截面方向
对于三维梁单元，见P131.
6.2计算公式和积分
ABAQUS/Standard中的三次单元(假设忽略剪切变形）不能模拟剪切变形，适合模拟相对细长构件的结构。
当壁厚是典型梁横截面尺寸的1/10时，一般的薄壁假设是有效的。
ABAQUS/Standard具有剪切变形梁单元——B31OS和B32OS，它们包括了在开口薄壁横截面中翘曲的影响。
6.3选择梁单元
6.4建模要点
Property：定义梁截面：先创建截面轮廓；再创建梁截面（若在分析前指定截面的积分方式，材料性质定义就可以作为截面定义的组成部分，而不需要另外给出材料的定义；若材料数据
的定义是独立于截面属性的，那么通过编辑材料定义可以将密度包括在内。）；再定义梁截面方向。
Interaction：螺栓连接等约束，在ABAQUS中可以用多点约束、约束方程或者连接件来定义。其中，连接件允许模拟在模型装配件中任意两点之间（或者在装配件中的任意一点与地面之间的连接。每一个连接件必须提供一个连接件特性一顶一它的类型(类似于在单元与截面特性之间的关系）。
后处理：对于诸如梁这样的一维单元，应用等值线图选项可以产生的“弯矩”类图形则更为有用。
6.7小结
第七章线性动态分析
如果加载时间很短或者载荷时间在性质上是动态的，必须采用动态分析。
7.1线性动态问题简介
线性问题：小变形、线弹性材料、无接触条件。如果在模拟中存在非线性，振型叠加法将不再适用。当对非线性动态响应感兴趣时，必须对运动方程进行直接积分。在ABAQUS/Standard中的非线性动态程序采用的是隐式时间积分，适用于求解非线性结构动态问题，但计算量很大。在ABAQUS/Explicit中的非线性动态分析计算量较小，最适合于求解应力波影响非常重要并且所需模拟的事件是很短时间(典型的不超过1s)的问题，并且它能够模拟不连续的非线性问题，比采用ABAQUS/Standard更容易些。
7.2阻尼
能量耗散包括结构连接处的摩擦和局部材料的迟滞效应。
对于瞬时模态分析，在ABAQUS/Standard中可以定义一些不同类型的阻尼：
直接模态阻尼(取值范围是临界阻尼的1%~10%)、Rayleigh阻尼、复合阻尼(当结构中有多种不同的材料时)。如果使用了模态阻尼，必须指定在基于模态的动态过程中使用的特征模态。
7.3单元选择
ABAQUS的所有单元均可用于动态分析，选取单元的一般原则与静力分析相同。但在模拟冲击和爆炸荷载时，应选用一阶单元。
设计应用于动态模拟的网格时，需要考虑在响应中将被激发的振型，并且使所采用的网格能够充分地反映出这些振型。这意味着能够满足静态模拟的网格，不一定能够计算由于加载激发的高频振型的动态响应。
7.5建模要点
加载的持续时间很短意味着惯性效应，要进行动态分析。
在动态分析中，必须给定每种材料的密度，这样才能形成质量矩阵。
创建时间历史动画：(1)(动画缩放系数)：指单个增量步的变形过程；(2)(动画时间历程)表示整个时间变化过程，即所有增量步的变化过程。
7.12小结
第八章非线性
线性分析：在外加荷载与系统响应之间为线性关系。
非线性结构问题指结构的刚度随其变形而改变的问题。
8.1非线性的来源
橡胶材料可以用一种非线性、可恢复(弹性)响应的材料来模拟。
如果边界条件在分析过程中发生变化，就会产生边界非线性问题。
几何非线性发生在位移的大小影响到结构响应的情况，可能由于：大挠度或大转动；“突然
翻转”；初应力或载荷刚性化。
8.2非线性问题的求解
(1)ABAQUS/Standard将模拟划分为一定数量的载荷增量步，并在每个增量步结束时寻求近似的平衡构型。所有这些增量响应的总和就是非线性分析的近似解答。因此，为了求解非线性问题，ABAQUS/Standard组合了增量和迭代过程。
(2)ABAQUS/Explicit对于线性和非线性问题是类似的，都只需要一个小的时间增量步，只依赖于模型的最高阶自振频率，而与载荷的类型和加载时间无关，非常适合于非线性分析。
(3)分析步、增量步、迭代步的区别：
分析步：
增量步：
迭代步：
(4)几何非线性
将几何非线性的效应引入到分析中，仅需要在ABAQUS/Standard模型分析的定义中考虑几何非线性的效应，这对于ABAQUS/Explicit是默认的设置。对于一个分析步，默认的增量步数目是100。提供一个合理的初始增量步的值有利于问题的求解。ABAQUS/Standard默认的最小容许时间增量为10-5乘以分析步的总时间，且默认没有增量值的上限值。
一旦在一个分析步中包括了几何非线性，在所有的后继分析步中就都会考虑几何非线性。
8.4建模要点
非线性分析中考虑了刚性效应和压力方向的改变，而在线性分析中均未考虑。
8.7小结
第九章显示非线性动态分析
ABAQUS/Explicit适合求解的问题：高速动力学事件；复杂的接触问题；复杂的后屈曲问题；高度非线性的准静态问题；材料退化和失效。
9.2动力学显式有限元方法
所谓“显式”是指在增量步结束时的状态仅依赖于该增量步开始时的位移、速度和加速度。时间增量必须相当小，如果时间增量大于最大的时间步长，则此时间增量已经超过了稳定性限制，进而导致数值不稳定，解答不收敛。通过仔细检查能量历史和其他响应变量，用户应当确保得到了有效地响应。
稳定极限与质量(缩放质量密度能潜在的提高分析的效率)、材料(通过波速影响)、网格(稳定极限与最短的单元尺寸成比例，采用一个尽可能均匀的网格，在.sta文件中列出具有最低稳定极限的10个单元)、数值不稳定(弹簧或减振器单元在分析过程中可能成为不稳定)有关。
9.5动态振荡的阻尼
模型中加入阻尼以限制数值振荡或者为系统增加物理阻尼。ABAQUS/Explicit中加入阻尼的方法：1.体粘性(采用线性体粘性来削弱在单元最高阶频率中的振荡)；2.粘性压力(一般应用于结构问题或者准静态问题，以阻止低阶频率的动态影响)；3.材料阻尼(材料模型本身或瑞利阻尼(合理的质量比例阻尼不会明显地降低稳定极限，但刚度比例阻尼可能会明显地降低稳定极限))；4.离散减震器(能够把阻尼只施加在自认为有必要施加的节点上，且应总是与弹簧或者桁架等其他单元并行使用，因为他们不会引起稳定极限的明显下降)。
9.6能量平衡
9.7弹簧和减振器的潜在不稳定性
下面列出的单元具有使分析不稳定的潜在因素：弹簧单元、减振器单元；
识别非稳定性：
1.分析运算的输出数据(状态文件)为我们提供了分析稳定性的一些标志；
2.位移历史(不切实际的位移)；
3.能量(在ABAQUS/CAE 中建立能量的历史图；一个常值得能量平衡表明解答是稳定的；相反，一个明显的、非恒定的能量平衡则明确地提示出解答不再稳定)
消除不稳定性：1.增加质量(附加质量对结构的影响很小，只影响数值解的稳定性)；2.考虑阻尼的影响(ABAQUS/Explicit 在计算稳定时间增量时并不考虑减振器，但考虑材料阻尼)；3.控制时间增量(固定时间增量或设置时间增量的比例系数) 9.8小结
第十章 材料
通过设置含材料名称的截面属性，将模型中的不同区域与不同材料定义相联系。 10.2延性金属的塑性
在金属中的工程应力(变形前每单位面积上的力)称为名义应力，与之共轭的为名义应变(变形前每单位长度的长度变化)。当正发生颈缩时，在金属中的名义应力远低于材料的极限强度。因此，当在ABAQUS 中定义塑性数据时，必须采用真实应力和真实应变。而材料试验的数据常常是以名义应力和名义应变的形式给出的。在这种情况下，必须将塑性材料的数据从名义应力/应变的值转换为真实应力/应变的值。真实应变和名义应变的关系为：
)1ln(nom εε+=真实应力、名义应力和名义应变的关系：)1(nom nom εσσ+=
ABAQUS在提供的数据点之间进行线性插值(或者在ABAQUS/Explicit中采用规则化数据)得到材料响应，并假设在输入数据定义范围之外的响应为常数。
10.3弹塑性问题的单元选取
当模拟材料的不可压缩特性时，在ABAQUS/Standard中的完全积分二次实体单元不能用于弹塑性问题的模拟；而完全积分一次实体单元可以安全地应用于塑性问题。
减缩积分的实体单元不会发生过约束，并且可用于大多数弹塑性问题的模拟。如果应变超过了20%~40%，在使用ABAQUS/Standard中的减缩积分二次单元时需注意，因为在此量级上它们可能会承受体积自锁。这种影响可以通过加密网格来降低。
10.4建模要点
等值线图：ABAQUS/Standard是在积分店处计算单元变量。通过将积分点处的数据外推到节点，ABAQUS/CAE计算单元变量的节点值。外推算法的阶数由单元类型决定。由外推过程产生的不变量可能超过变量极限。若外推应力值与积分点处应力值的差别很大表明单元之间应力剧烈变化的事实，并且对于精确的应力计算，网格过于粗糙。如果网格细划，将会显著地减小这种外推应力的误差，但它总是会以某种程度出现。
等效塑性应变的等值线图：在材料中的等效塑性应变(PEEQ)是用来表示材料非弹性变形的标量。通过调用Contour Plot Options(等值线图选项)对话框，将等值线的最小值设置为等效塑性应变的一个小量(例如：-1E-4)，这样在ABAQUS/CAE绘制的模型图中，任何用深蓝色显示的区域仍具有弹性材料的特性。
10.5建模要点
定义压力荷载：通过定义幅值实质上定义了载荷的量值，仅需要施加一个单位压力到模型上。
10.6超弹性
典型的橡胶材料的应力-应变行为是弹性的，但是高度的非线性，这种行为称为超弹性。超弹性的变形在大应变值时(通常超过100%)仍然保持为弹性。当ABAQUS模拟超弹性材料时，作出如下假设：材料行为是弹性；材料行为是各向同性；模拟将考虑几何非线性效应。另外，ABAQUS/Standard默认地假设材料是不可压缩的；ABAQUS/Explicit假设材料是接近不可压缩的(默认的泊松比为0.475)。
可压缩性：在材料被高度限制的情况下(如用作密封垫的圆圈)，必须正确地模拟可压缩性。
10.6.2可压缩性
ABAQUS/Standard中的杂交单元必须用于模拟出现在超弹性材料中的完全不可压缩行为。
评估材料的可压缩性：材料的初始体积模量K0与其初始剪切模量μ0的比值；泊松比ν。
10.6.3应变势能
ABAQUS应用应变势能(U)（strain potential energy）来表达超弹性材料的应力-应变关系，而不是用杨氏模量和泊松比。
10.6.4应用实验数据定义超弹性行为
与塑性数据不同，对于超弹性材料的实验数据必须作为名义应力和名义应变的值提供给ABAQUS。材料的可压缩性重要时才需要给出体积压缩数据。下面的方法可以改善超弹性材料模型：从可能发生在模拟中的变形模式获得实验数据；拉伸和压缩数据均允许使用，其中压缩应力和应变作为负值键入；尽可能包含平面实验的数据；在所期望的模拟过程中材料实际承受的应变量级上，应提供更多的数据；利用在ABAQUS/CAE中的材料评估功能对实验进行模拟，并将实验数据与ABAQUS的计算结果进行比较。
10.7建模要点
10.9减少体积自锁的技术
网格细化；引入少量的可压缩性(可压缩性的存在减轻了体积自锁；通过设置材料常数D1
为一个非零值，引入压缩性，选择这个值以使初始泊松比v0接近于0.5)。
专业ABAQUS有限元建模经验笔记
基于ABAQUS的有限元分析和应用 第一章绪论 1.有限元分析包括下列步骤： 2.为了将试验数据转换为输入文件，分析者必须清楚在程序中所应用的和由实验人员提供的材料数据的应力和应变的度量。 3.ABAQUS建模需注意以下内容： 4.对于许多包含过程仿真的大变形问题和破坏分析，选择合适的网格描述是非常重要的，需要认识网格畸变的影响，在选择网格时必须牢牢记住不同类型网格描述的优点。 第二章ABAQUS基础 1.一个分析模型至少要包含如下的信息：离散化的几何形体、单元截面属性、材料数据、载荷和边界条件、分析类型和输出要求。 ①离散化的几何形体：模型中所有的单元和节点的集合称为网格。 ②载荷和边界条件： 2.功能模块： （1）Assembly（装配）：一个ABAQUS模型只能包含一个装配件。 （2）Interaction（相互作用）：相互作用与分析步有关，这意味着用户必须规定相互作用是在哪些分析步中起作用。 （3）Load（载荷）：载荷和边界条件与分析步有关，这意味着用户指定载荷和边界条件是在哪些分析步中起作用。 （4）Job（作业）：多个模型和运算可以同时被提交并进行监控。 3.量纲系统 ABAQUS没有固定的量纲系统，所有的输入数据必须指定一致性的量纲系统，常用的一致性量纲系统如下：
4.建模要点 （1）创建部件：设定新部件的大致尺寸的原则必须是与最终模型的最大尺寸同一量级。（2）用户应当总是以一定的时间间隔保存模型数据（例如，在每次切换功能模块时）。（3）定义装配： 在模型视区左下角的三向坐标系标出了观察模型的方位。在视区中的第2个三向坐标系标出了坐标原点和整体坐标系的方向（X，Y和Z轴）。 （4）设置分析过程： （5）在模型上施加边界条件和荷载： 用户必须指定载荷和边界条件是在哪个或哪些分析步中起作用。 所有指定在初始步中的力学边界条件必须赋值为零，该条件是在ABAQUS/CAE中自动强加的。 在许多情况下，需要的约束方向并不一定与整体坐标方向对齐，此时用户可定义一个局部坐标系以施加边界条件。 在ABAQUS中，术语载荷通常代表从初始状态开始引起结构响应发生变化的各种因素，包括：集中力、压力、非零边界条件、体力、温度（与材料热膨胀同时定义）。
abaqus有限元分析过程
一、有限单元法的基本原理 有限单元法（The Finite Element Method）简称有限元（FEM），它是利用电子计算机进行的一种数值分析方法。它在工程技术领域中的应用十分广泛，几乎所有的弹塑性结构静力学和动力学问题都可用它求得满意的数值结果。 有限元方法的基本思路是：化整为零，积零为整。即应用有限元法求解任意连续体时，应把连续的求解区域分割成有限个单元，并在每个单元上指定有限个结点，假设一个简单的函数（称插值函数）近似地表示其位移分布规律，再利用弹塑性理论中的变分原理或其他方法，建立单元结点的力和位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程组，从而求解结点的位移分量. 进而利用插值函数确定单元集合体上的场函数。由位移求出应变, 由应变求出应力 二、ABAQUS有限元分析过程 有限元分析过程可以分为以下几个阶段 1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型――有限元模型，从而为有限元数值计算提供必要的输入数据。有限元建模的中心任务是结构离散，即划分网格。但是还是要处理许多与之相关的工作：如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。
2.计算阶段:计算阶段的任务是完成有限元方法有关的数值计算。 由于这一步运算量非常大，所以这部分工作由有限元分析软件控制并在计算机上自动完成 3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理， 并按一定方式显示或打印出来，以便对结构性能的好坏或设计的合理性进行评估，并作为相应的改进或优化，这是惊醒结构有限元分析的目的所在。 下列的功能模块在ABAQUS/CAE操作整个过程中常常见到，这个表简明地描述了建立模型过程中要调用的每个功能模块。 “Part（部件） 用户在Part模块里生成单个部件，可以直接在ABAQUS/CAE环境下用图形工具生成部件的几何形状，也可以从其它的图形软件输入部件。 Property（特性） 截面（Section）的定义包括了部件特性或部件区域类信息，如区域的相关材料定义和横截面形状信息。在Property模块中，用户生成截面和材料定义，并把它们赋于(Assign)部件。 Assembly（装配件） 所生成的部件存在于自己的坐标系里，独立于模型中的其它部件。用户可使用Assembly模块生成部件的副本（instance），并且在整体坐标里把各部件的副本相互定位，从而生成一个装配件。 一个ABAQUS模型只包含一个装配件。
abaqus 有限元分析(齿轮轴)
Abaqus分析报告 （齿轮轴） 名称：Abaqus齿轮轴 姓名： 班级： 学号： 指导教师：
一、简介 所分析齿轮轴来自一种齿轮泵，通过用abaqus软件对齿轮轴进行有限元分析和优化。齿轮轴装配结构图如图1，分析图1中较长的齿轮轴。 图1.齿轮轴装配结构图 二、模型建立与分析 通过part、property、Assembly、step、Load、Mesh、Job等步骤建立齿轮轴模型，并对其进行分析。 1.part 针对该齿轮轴，拟定使用可变型的3D实体单元，挤压成型方式。 2.材料属性 材料为钢材，弹性模量210Gpa，泊松比0.3。
3.截面属性 截面类型定义为solid，homogeneous。 4.组装 组装时选择dependent方式。 5.建立分析步 本例用通用分析中的静态通用分析（Static，General）。 6.施加边界条件与载荷 对于齿轮轴，因为采用静力学分析，考虑到前端盖、轴套约束，而且根据理论，对受力部分和轴径突变的部分进行重点分析。 边界条件：分别在三个轴径突变处采用固定约束，如图2。 载荷：在Abaqus中约束类型为pressure，载荷类型为均布载荷,分别施加到齿轮接触面和键槽面，根据实际平衡情况，两力所产生的绕轴线的力矩方向相反，大小按比例分配。 均布载荷比计算： 矩形键槽数据： 长度：8mm、宽度：5mm、高度：3mm、键槽所在轴半径：7mm 键槽压力面积：S1 = 8x3=24mm2 平均受力半径：R1=6.5mm 齿轮数据：= 齿轮分度圆半径：R2 =14.7mm、压力角：20°、 单个齿轮受力面积：S2 ≈72mm2 通过理论计算分析，S1xR1xP1=S2xR2xP2，其中，P1为键槽均布载荷
支架的有限元分析ABAQUS
支架的线性静力学分析实例：建模和分析计算 在此实例中读者将学习ABAQUS/CAE的以下功能。 1) Sketch功能模块:导人CAD二维图形，绘制线段、圆弧和倒角，添加尺寸，修改平面图，输出平面图。 2) Part功能模块:通过拉伸来创建几何部件，通过切割和倒角未定义几何形状。 3) Property功能模块:定义材料和截面属性。 4) Mesh功能模块:布置种子，分割实体和面，选择单元形状、单元类型、网格划分 技术和算法，生成网格，检验网格质量，通过分割来定义承受载荷的面。 5) Assembly功能模块:创建非独立实体。 6) Step功能模块:创建分析步，设置时间增量步和场变量输出结果。 7) Interaction功能模块:定义分布榈合约束(distributing coupling constraint)。 8) Load功能模块:定义幅值，在不同的分析步中分别施加面载荷和随时间变化的集中力，定义边界条件。 9) Job功能模块:创建分析作业，设置分析作业的参数，提交和运行分析作业，监控运行状态。 10) Visualization功能模块:后处理的各种常用功能。 结构静力学分析（static analysis）是有限元法的基本应用领域，适用于求解惯性及阻尼对结构响应不显著的问题。主要用来分析由于稳态外载荷引起的位移，应力和应变等。本章的静力学分析实例按照ABAQUS工程分析的流程对支架进行线性静力学分析，通过实例基本掌握了分析的流程，同时了解接触的定义。 1.问题描述 所示的支架，一端牢固地焊接在一个大型结构上，支架的圆孔中穿过一个相对较软的杆件，圆孔和杆件用螺纹连接。材料的弹性模量E=2100000MPa,泊松比为0.3。
ABAQUS有限元接触分析的基本概念
ABAQUS有限元接触分析的基本概念2009-11-24 00:06:28 作者：jiangnanxue 来源：智造网—助力中国制造业创新—https://www.sodocs.net/doc/0a4889831.html, CAE(计算机辅助工程)是一门复杂的工程科学，涉及仿真技术、软件、产品设计和力学等众多领域。世界上几大CAE公司各自以其独到的技术占领着相应的市场。ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群，是国际上公认的最先进的大型通用非线性有限元分析软件之一。它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域。ABAQUS在技术、品质和可靠性等方面具有卓越的声誉，可以对工程中各种复杂的线性和非线性问题进行分析计算。 《ABAQUS有限元分析常见问题解答》以问答的形式，详细介绍了使用ABAQUS建模分析过程中的各种常见问题，并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法，帮助读者提高解决问题的能力。 《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。 16.1.1 点对面离散与面对面离散 【常见问题16-1】 在ABAQUS/Standard分析中定义接触时，可以选择点对面离散方法(node-to-surface-dis - cre-tization)和面对面离散方法(surface-to-surface discretization)，二者有何差别? 『解答』 在点对面离散方法中，从面(slave surface)上的每个节点与该节点在主面(master surface)上的投影点建立接触关系，每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点。 使用点对面离散方法时，从面节点不会穿透(penetrate)主面，但是主面节点可以穿透从面。 面对面离散方法会为整个从面(而不是单个节点)建立接触条件，在接触分析过程中同时考虑主面和从面的形状变化。可能在某些节点上出现穿透现象，但是穿透的程度不会很严重。 在如图16-l和图16-2所示的实例中，比较了两种情况。
Abaqus有限元分析中的沙漏效应
Abaqus有限元分析中的沙漏效应[转] 2011-09-21 17:34:27| 分类：有限元 | 标签： |字号大中小订阅 1. 沙漏的定义 沙漏hourglassing一般出现在采用缩减积分单元的情况下： 比如一阶四边形缩减积分单元，该单元有四个节点“o”，但只有一个积 分点“*”。而且该积分点位于单元中心位置，此时如果单元受弯或者受剪，则必然会发生变形，如下图a所示。 关于沙漏问题，建议看看abaqus的帮助文档，感觉讲的非常好，由浅入深，把深奥的东西讲的很容易理解。 沙漏的产生是一种数值问题，单元自身存在的一种数值问题，举个例子，对于单积分点线性单元，单元受力变形没有产生应变能--也叫0能量模式，在 这种情况下，单元没有刚度，所以不能抵抗变形，不合理，所以必须避免这种情况的出现，需要加以控制，既然没有刚度，就要施加虚拟的刚度以限制沙漏 模式的扩展---人为加的沙漏刚度就是这么来的。 关于沙漏现象的判别，也就是出现0能模式的方法最简单的是察看单元变 形情况，就像刚才所说的单点积分单元，如果单元变成交替出现的梯形形状， 如果多个这样的单元叠加起来，是不是象我们windows中的沙漏图标呢？ ABAQUS中沙漏的控制： *SECTION CONTROLS：指定截面控制 警告：对于沙漏控制，使用大于默认值会产生额外的刚度响应，甚至当值 太大时有时导致不稳定。默认沙漏控制参数下出现沙漏问题表明网格太粗糙， 因此，更好的解决办法是细化网格而不是施加更大的沙漏控制。 该选项用来为减缩积分单元选择非默认的沙漏控制方法，和standard中的修正的四面体或三角形单元或缩放沙漏控制的默认系数；在explicit中，也 为8节点块体单元选择非默认的运动方程：为实体和壳选择二阶方程、为实体 单元激活扭曲控制、缩放线性和二次体积粘度、设置当单元破损时是否删除他们、或为上述完全破损的单元指定一标量退化参数。等 必需参数： NAME：名字 可选参数： DISTORTION CONTROL：只用于explicit分析。=YES激活约束防止负体积 单元出现或其他可压缩材料的过度变形，这对超弹材料是默认的。DISTORTION
abaqus有限元分析简支梁
1.梁C 的主要参数： 其中：梁长3000mm ，高为406mm ，上下部保护层厚度为38mm ，纵筋端部保护层厚度为25mm 抗压强度：35.1MPa 抗拉强度：2.721MPa 受拉钢筋为2Y16，受压钢筋为2Y9.5，屈服强度均为440MPa 箍筋：Y7@102，屈服强度为596MPa 2.混凝土及钢筋的本构关系 1、运用陈光明老师的论文（Chen et al. 2011）来确定混凝土的本构关系： 受压强度： 其中C a E ==28020，c f ρσ'=，0.002ρε= 2、受压强度与开裂位移的相互关系：
其中123.0, 6.93c c == 3、损伤因子： 其中c h = e=10（选取网格为10mm ） 4、钢筋取理想弹塑性 5、名义应力应变和真实应力及对数应变的转换： ln (1)ln(1)true nom nom Pl true nom E σσεσε ε=+=+- 6、混凝土最终输入的本构关系如下： compressive behavior tensile behavior tension damage yield stress inelastic strain yield stress displacement parameter displacement 21.50274036 2.721 25.56359281 2.72247E-05 2.683556882 0.0003129 0.18766492 0.0003129 28.88477336 8.85105E-05 2.646628319 0.0006258 0.31902609 0.0006258 31.43501884 0.000177278 2.610210508 0.0009387 0.41606933 0.0009387 33.24951537 0.000292271 2.574299562 0.0012516 0.49065237 0.0012516 34.40787673 0.000430648 2.538891515 0.0015645 0.54973463 0.0015645 35.01203181 0.000588772 2.503982327 0.0018774 0.5976698 0.0018774 35.16872106 0.000762833 2.46956789 0.0021903 0.63732097 0.0021903 34.97805548 0.000949259 2.435644029 0.0025032 0.67064827 0.0025032 34.52749204 0.001144928 2.402206512 0.0028161 0.69903885 0.0028161 33.88973649 0.001347245 2.369251048 0.003129 0.72350194 0.003129 33.17350898 0.001541185 2.336773294 0.0034419 0.74478941 0.0034419 32.38173508 0.001737792 2.30476886 0.0037548 0.76347284 0.0037548 31.54367693 30.68161799 0.001936023 0.002135082 2.27323331 2.242162167 0.0040677 0.0043806 0.77999451 0.79470205 0.0040677 0.0043806
abaqus有限元建模小例子
问题一: 工字梁弯曲 1.1 问题描述: 在<<材料力学实验>>中，弯曲实验測定了工字梁弯曲应变大小及其分布，以验证弯曲正应力公式。在这里，採用ABAQUS/CAE建立试验件的有限元模型，ABAQUS/Standard模块进行分析求解，得到应力、应变分布，对比其与理论公式计算值及实验測量值的差別。 弯曲实验的相关数据： 材料：铝合金E=70GPa 泊松比0.3 实验装置结构简图如图所示： 结构尺寸测量值：H=50（+/-0.5mm） h=46(+/-0.5mm) B=40（+/-0.5mm） b=2(+/-0.02mm) a=300(+/-1mm) F1=30N Fmax=300N N ? F100 = 1.2 ABAQUS有限元建模及分析 一对象: 工字型截面铝合金梁 梁的结构简图如图1所示，結构尺寸、载荷、約束根据1.1设定，L取1600mm，两端各伸出100mm。 二用ABAQUS/CAE建立实验件的有限元模型,效果图如下： 边界条件简化： 左侧固定铰支座简化为下表面左参考点处的约束U1=U2=U3=0
右侧活动铰支座简化为下表面右参考点处的约束U1=U2=UR3=0 几何模型
有限元模型 三ABAQUS有限元分析結果 ①应力云图（Z方向正应力分量）：施加载荷前 F=300N
②应变（Z方向分量）： 中间竖直平面的厚度方向应变分布图： F=100N F=200N
F=300N 由上图可以看出应变沿着厚度方向呈线性比例趋势变化，与实验测得的应变值变化趋势相同。中性轴处应变均接近零值，应变与距离中性轴位移基本为正比关系。 1.3分析结果: 中间竖直截面上下边缘轴向应力数值对比：*10^-6 MPa 距中性轴距ABAQUS模拟实验测量值平均理论值 1/2H -96.182*70000 -97*70000 -6.9165=-70000*98.807 -1/2H 95.789*70000 92*70000 6.9165
abaqus有限元分析报告开裂梁要点
Abaqus梁的开裂模拟计算报告 1.问题描述 利用ABAQUS有限元软件分析如图1.1所示的钢筋混凝土梁的裂缝开展。参考文献Brena et al.（2003）得到梁的基本数据： 图1.1 Brena et al.（2003）中梁C尺寸 几何尺寸：跨度3000mm，截面宽203mm，高406mm的钢筋混凝土梁 由文献Chen et al. 2011得材料特性： 1.混凝土：抗压强度f c’=35.1MPa，抗拉强度f t= 2.721MPa，泊松比ν=0.2，弹性模量 E c=28020MPa； 2.钢筋：弹性模量为E c=200GPa，屈服强度f ys=f yc=440MPa，f yv=596MPa 3.混凝土垫块：弹性模量为E c=28020MPa，泊松比ν=0.2 2.建模过程 1)Part 打开ABAQUS使用功能模块，弹出窗口Create Part，参数为：Name：beam；Modeling Space：2D；Type：Deformable；Base Feature─Shell；Approximate size：2000。点击Continue 进入Sketch二维绘图区。由于该梁关于Y轴对称，建模的时候取沿X轴的一半作为模拟对象。 使用功能模块，分别键入独立点(0,0),(1600,0),(1600,406),(406,0),(0,0)并按下下方提 示区的Done，完成草图。 图2.1 beam 部件二维几何模型
相同的方法建立混凝土垫块： 图2.2 plate 部件二维几何模型 所选用的点有(0,0),(40,0),(40,10),(0,10) 受压区钢筋： 在选择钢筋的base feature的时候选择wire，即线模型。 图2.3 compression bar 部件二维几何模型 选取的点(0,0),(1575,0) 受拉区钢筋： 图2.4 tension bar 部件二维几何模型 选取的点(0,0),(1575,0) 箍筋： 图2.5 stirrup 部件二维几何模型 选取的点为(0,0),(0,330) 另外，此文里面为了作对比，部分的模型输入尺寸的时候为m，下面无特别说明尺寸都为mm。
石亦平ABAQUS有限元分析实例详解之读后小结-完整版
目录 第一章ABAQUS简介 (1) 第二章ABAQUS基本使用方法 (1) 第三章线性静力分析实例 (6) 第四章 ABAQUS的主要文件类型 (8) 第五章接触分析实例 (9) 第六章弹塑性分析实例 (13) 第七章热应力分析实例 (15) 第八章多体分析实例 (16) 第九章动态分析实例 (17) 第十章复杂工程分析综合实例 (20)
第一章ABAQUS简介 [1] (pp7) 在[开始] →[程序] →[ABAQUS 6.5-1]→[ABAQUS COMMAND]，DOS提示符下输入命令 Abaqus fetch job =
Abaqus螺栓有限元分析
1.分析过程 1.1.理论分析 1.2.简化过程 如果将Pro/E中的3D造型直接导入Abaqus中进行计算，则会出现裂纹缝隙无法修补，给后期的有限元分析过程造成不必要的麻烦，因此，在Abaqs中进行计算之前，对原来的零件模型进行一些简化和修整。 A.法兰部分不是分析研究的重点，因此将其简化掉； B.经计算，M24×3的螺纹的升角很小，在度，因此可以假设螺旋升角为0； C.忽略螺栓和螺母的圆角等细节； 1.3.Abaqus中建模 查阅机械设计手册，得到牙型如下图所示，在Abaqus中按照下图所示创建出3D模型，如错误!未找到引用源。所示。同样的方式，我们建立螺母的3D模型nut，如错误!未找到引用源。所示。
图 1-1 图 1-2 建立材料属性并将其赋予模型。在Abaqus的Property模块中，选择Material->Manager->Create，创建一个名为Bolt&Nut的新材料，首先设置其弹性系数。在Mechanical->Elastic中设置其杨氏模量为193000Mpa，设置其泊松比为，如错误!未找到引用源。所示。 建立截面。点击Section->Manager->Creat，建立Solid，Homogeneous的
各向同性的截面，选择材料为Bolt&Nut，如错误!未找到引用源。所示。 将截面属性赋予模型。选择Assign->Section，选择Bolt模型，然后将刚刚建立的截面属性赋予它。如错误!未找到引用源。所示。同样，给螺母nut赋予截面属性。 图 1-3 图 1-4
图 1-5 然后，我们对建立的3D模型进行装配，在Abaqus中的Assembly模块中，我们同时调入两个模型，然后使用Constraint->Coaxial命令和Translate和Instance命令对模型进行移动，最终的装配结果如错误!未找到引用源。所示。 图 1-6 第四步，对模型进行网格划分。进入Abaqus中的Mesh模块，然后选择Bolt 零件，使用按边布种的方式对其进行布种，布种结果如错误!未找到引用源。所示。在菜单Mesh->Control中进行如错误!未找到引用源。所示的设置使用自由网格划分，其余设置使用默认。在菜单Mesh->Element type中选用如错误!未找到引用源。所示的设置。按下Mesh图标，对工件进行网格划分，最终的结果如错误!未找到引用源。所示。同样的方式对螺母模型nut进行网格划分，最终结
石亦平ABAQUS有限元分析实例详解之读后小结 (Part 4)
石亦平《ABAQUS有限元分析实例详解》之读后小结 第九章动态分析实例 [95] (pp280) ABAQUS包括两大类方法： 振型叠加法（modal superposition procedure）：用于求解线性动态问题； 直接解法（direct-solution dynamic analysis procedure）：主要用于求解非线性动态问题。 提示：ABAQUS的所有单元均可用于动态分析，选取单元的一般原则与静力分析相同。但在模拟冲击和爆炸载荷时，应选用一阶单元，因为它们具有集中质量公式，模拟应力波的效果优于 二次单元所采用的一致质量公式。 [96] (pp281) 振型叠加法的基础是结构的各阶特征模态（eigenmode），因此在建模时要首先定义一个 频率提取分析步（frequency extraction），从而得到结构的振型（mode shape）和固有频率（natural frequency），然后才能定义振型叠加法的各种分析步。振型叠加法包括4种分析类型： （1）瞬时模态动态分析（transient modal dynamic analysis）计算线性问题在时域（time domain）上的动态响应。用此分析要满足如下5个基本条件： (a) 系统是线性的（线性材料特性，无接触行为，不考虑几何非线性）。 (b) 响应只受相对较少的频率支配。当在响应中频率的成分增加时（例如打击和碰撞问题），振 型叠加法的效率将会降低。 (c) 载荷的主要频率应该在所提取的频率范围之内，以确保对载荷的描述足够精确。 (d) 特征模态应该能精确地描述任何突然加载所产生的初始加速度。 (e) 系统的阻尼不能过大。 （2）基于模态的稳态动态分析（mode-based steady-state dynamic analysis）在用户指定频率内的谐波激励下，计算引起结构响应的振幅和相位，得到的结果是在频域（frequency domain）上的。其典型分析对象包括发动机的零部件和建筑物中的旋转机械等。 （3）反应谱分析（response spectrum analysis）当结构的固定点处发生动态运动时，计算其峰值响应（位移、应力等），得到的结果是在频域上的。其典型应用是计算在发生地震时建筑物 的峰值响应。 （4）随机响应分析（random response analysis）当结构随机连续的激励时，计算其动态响应，
Abaqus螺栓有限元分析
Abaqus螺栓有限元分析
1.分析过程 1.1.理论分析 1.2.简化过程 如果将Pro/E中的3D造型直接导入Abaqus中进行计算，则会出现裂纹缝隙无法修补，给后期的有限元分析过程造成不必要的麻烦，因此，在Abaqs中进行计算之前，对原来的零件模型进行一些简化和修整。 A.法兰部分不是分析研究的重点，因此将其简化掉； B.经计算，M24×3的螺纹的升角很小，在度，因此可以假设螺旋升角为0； C.忽略螺栓和螺母的圆角等细节； 1.3.Abaqus中建模 查阅机械设计手册，得到牙型如下图所示，在Abaqus中按照下图所示创建出3D模型，如图错误！文档中没有指定样式的文字。-1所示。同样的方式，我们建立螺母的3D模型nut，如图错误！文档中没有指定样式的文字。-2所示。
图错误！文档中没有指定样式的文字。-1 图错误！文档中没有指定样式的文字。-2 建立材料属性并将其赋予模型。在Abaqus的Property模块中，选择Material->Manager->Create，创建一个名为Bolt&Nut的新材料，首先设置其弹性系数。在Mechanical->Elastic中设置其杨氏模量为193000Mpa，设置其泊松比为0.3，如图错误！文档中没有指定样式的文字。-4所示。 建立截面。点击Section->Manager->Creat，建立Solid，Homogeneous的各向同性的截面，选择材料为Bolt&Nut，如图错误！文档中没有指定样式的文字。-5所示。
将截面属性赋予模型。选择Assign->Section，选择Bolt模型，然后将刚刚建立的截面属性赋予它。如图错误！文档中没有指定样式的文字。-3所示。同样，给螺母nut赋予截面属性。 图错误！文档中没有指定样式的文字。-3 图错误！文档中没有指定样式的文字。-4
相关文档
- 《ABAQUS动力学有限元分析指南》张文元3
- ABAQUS有限元分析从入门到精通(第2版)
- Abaqus螺栓有限元分析
- Abaqus螺栓有限元分析
- 专业ABAQUS有限元建模经验笔记
- 石亦平ABAQUS有限元分析实例详解之读后小结-完整版
- ABAQUS有限元接触分析的基本概念
- abaqus有限元分析报告开裂梁要点
- ABAQUS有限元分析实例详解 3
- ABAQUS有限元分析用的材料属性表material
- 支架的有限元分析ABAQUS
- abaqus有限元分析简支梁
- 石亦平ABAQUS有限元分析实例详解之读后小结-完整版
- 石亦平ABAQUS有限元分析实例详解之读后小结 (Part 4)
- (完整版)ABAQUS有限元分析实例详解3
- abaqus有限元分析
- abaqus有限元建模小例子
- (完整版)ABAQUS有限元分析实例详解
- abaqus有限元分析过程
- abaqus 有限元分析(齿轮轴)