Abaqus单元类型选择

ABAQUS/analysis_单元选择标准单元类型选择标准ElementtypeA1.2ABAQUS/analysis_单元选择标准概述•ABAQUS中的单元•结构单元（壳和梁）vs.连续体单元•使用连续体单元模拟弯曲•应力集中•接触•不可压材料•网格生成•选择实体单元总结ABAQUS/analysis_单元选择标准ABAQUS中的单元A1.4ABAQUS/analysis_单元选择标准ABAQUS中的单元•ABAQUS单元库中大量的单元为不同几何体和结构建模提供了非常大的灵活性。–可以通过以下的特征为单元分类：•族•节点个数•自由度•公式•积分点A1.5ABAQUS/analysis_单元选择标准ABAQUS中的单元族•有限元族是一种广泛的分类方法。•同族的单元共享许多基本特征。•在同一族单元中又有许多变异。特殊单元，如弹簧、阻尼器和质量单元连续体(实体单元)壳单元梁单元刚体单元薄膜单元桁架单元无限单元A1.6ABAQUS/analysis_单元选择标准ABAQUS中的单元节点个数(插值)•节点的单元编号决定了单元域内节点自由度的插值方式。•ABAQUS包含一阶和二阶插值方式的单元。一次插值二次插值A1.7ABAQUS/analysis_单元选择标准ABAQUS中的单元自由度•在有限元分析过程中，单元节点的自由度是基本变量。•自由度的例子：–位移–转动–温度–电势•一些单元具有与用户定义的节点不相关的内部自由度。A1.8ABAQUS/analysis_单元选择标准ABAQUS中的单元公式•用于描述单元行为的数学公式是用于单元分类的另一种方法。•不同单元公式的例子：–平面应变–平面应力–杂交单元–非协调元–小应变壳–有限应变壳–厚壳–薄壳A1.9ABAQUS/analysis_单元选择标准积分点•在单元之内，刚度和单元质量在采样点，所谓的“积分点”，进行数值计算。•用于积分这些变量的数值算法将影响单元的行为。•ABAQUS包含“全”积分和“减缩”积分单元。ABAQUS中的单元A1.10ABAQUS/analysis_单元选择标准•全积分:–对于具有线弹性材料属性的、未扭转的单元，精确积分应变能所需的最小积分阶数。•减缩积分:–积分的阶数比全积分小一阶。ABAQUS中的单元一次插值全积分二次插值减缩积分A1.11ABAQUS/analysis_单元选择标准ABAQUS中的单元•单元命名约定：例子B21:Beam,2-D,1st-orderinterpolationCAX8R:Continuum,AXisymmetric,8-node,ReducedintegrationDC3D4:Diffusion(heattransfer),Continuum,3-D,4-nodeS8RT:Shell,8-node,Reducedintegration,TemperatureCPE8PH:Continuum,Planestrain,8-node,Porepressure,HybridDC1D2E:Diffusion(heattransfer),Continuum,1-D,2-node,ElectricalA1.12ABAQUS/analysis_单元选择标准ABAQUS中的单元•比较ABAQUS/Standard和ABAQUS/Explicit单元库–两种程序基本上具有相同的单元族：连续体、壳、梁等等。–除了应力分析，ABAQUS/Standard包括许多可以用于其它分析类型的单元：热传导、土壤固结、声学等等。•在ABAQUS/Explicit中也可以使用声学单元。–对于每个单元族，ABAQUS/Standard包含许多变种。–ABAQUS/Explicit包含几乎所有的一阶单元。•例外：二阶三角形和四面体单元、二阶梁单元–对于两种程序，许多单元选择的准则是一样的。ABAQUS/analysis_单元选择标准结构单元（壳和梁）vs.连续体单元A1.14ABAQUS/analysis_单元选择标准•连续体(实体)单元模型一般较大、并且昂贵，尤其对于三维问题。•如果合适，应该尽量使用结构单元(壳和梁)，以得到更经济的解。–与连续体单元模型相比，结构单元模型需要的单元一般少得多。•对于能得到可接受结果的结构单元，壳的厚度和梁截面的尺寸应该小于总体结构尺寸的1/10，比如：–支撑或点载荷之间的距离–尺寸变化很大的横截面之间的距离–最高振动模态的波长结构单元（壳和梁）vs.连续体单元A1.15ABAQUS/analysis_单元选择标准•壳单元–使用表面模型构成的壳单元近似模拟三维实体连续体单元。•可以有效的模拟弯曲和面内变形。–如果需要分析某个区域的细节，使用多点约束或子模型的办法可以将局部的三维实体模型加入到壳单元模型中。遭受发射冲击半球圆顶的壳模型结构单元（壳和梁）vs.连续体单元3-D实体表面模型A1.16ABAQUS/analysis_单元选择标准结构单元（壳和梁）vs.连续体单元•梁单元–使用线模型构成的梁单元近似模拟三维实体单元。•可以有效的模拟弯曲、扭转和轴力。•有许多可用的横截面形状。•还可以用工程常数的方式指定横截面属性。线模型利用梁单元建模的框架结构3-D实体ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲A1.18ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•纯弯曲的物理特征–有限元方法企图模拟的材料行为是：•在变形过程中，横截面仍然保持为平面。•沿厚度方向，轴向应变xx线性变化。•如果=0，厚度方向的应变yy为零。•没有薄膜剪切应变。–在圆弧上，线平行于梁轴。xxA1.19ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•使用二阶实体单元模拟弯曲(CPE8,C3D20R,…)–二阶全积分和减缩积分实体单元可以精确模拟弯曲：•轴向应变等于初始水平线的长度变化。•厚度方向的应变为零。•剪切应变为零。初始垂直的线长度不变(意味着yy=0)。因为假设单元边为曲线，所以变形等参线之间的夹角仍为90o(意味着xy=0)。等参线A1.20ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•使用一阶全积分实体单元模拟弯曲问题(CPS4,CPE4,C3D8)–这些单元在积分点检测到剪切应变。•不真实的；由于使用的单元公式才出现。–单元过硬的行为源于应变能用于产生剪切变形而不是产生弯曲（称为“剪切锁闭”）。因为单元的边必须仍然为直边，所以变形之后的等参线之间的夹角不再保持为90o(意味着)。0xy积分点不要在弯曲起主要作用的区域使用这些单元!A1.21ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•使用一阶减缩积分单元模拟弯曲(CPE4R,…)–这些单元去除了剪切锁闭效应。–然而，使用这些单元时会出现沙漏现象。•在中心处，只有一个积分点。•在厚度方向的一个单元不能检测到弯曲应变。•变形为零能模式（产生变形，但是没有应变；被称为“沙漏”）。长度的改变为零(意味着积分点上没有应变）。单个一阶减缩积分单元的弯曲行为。A1.22ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•沙漏在一阶减缩积分单元中可以很容易的传播，并导致不可靠的结果。•如果在厚度方向使用多个单元—最少四个单元，沙漏将不会产生问题。–每个单元捕获压缩或拉伸应变，但不能同时捕获。–轴向应变被正确的度量。–厚度方向的应变和剪切应变为零。–节省计算费用和有效的单元。在厚度方向有四个单元没有沙漏A1.23ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•检查并控制沙漏–在变形形状的绘图中，可以看到沙漏现象。•例子：带有中心点载荷简支梁的粗网格和中等网格。–ABAQUS有内建的沙漏控制方法，用以限制沙漏产生的问题。•确认用于控制沙漏的伪应变能与内能之比很小(1%)。同样的载荷和位移(1000×)。A1.24ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲–在ABAQUS/Viewer中，使用UsetheX–Y绘图功能比较能量。内能伪应变能伪应变能内能沿厚度方向有两个单元：伪应变能与内能之比为2%。沿厚度方向有四个单元：伪应变能与内能之比为0.1%。A1.25ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•使用非协调模式单元模拟弯曲(CPS4I,…)–对于以弯曲为主的问题中，这种单元可能是效率最高的实体单元。–计算费用在一阶和二阶减缩积分单元之间，兼有两种积分方法的优点。•可以正确的模拟剪切行为—在纯弯曲问题中没有剪切应变。•在厚度方向，仅用一个单元就可以模拟弯曲。•没有沙漏模式；在塑性和接触问题中，工作的很好。–如果单元严重扭曲，相对一阶减缩积分单元的优势将会减弱；然而，在严重扭曲的条件下，所有单元的精度都会下降。A1.26ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•例子：扭曲单元的悬臂梁平行扭曲梯形扭曲A1.27ABAQUS/analysis_单元选择标准使用连续体单元模拟弯曲•总结单元类型xxyyxy注释物理行为000二阶插值000OK一阶插值、全积分000剪切锁闭一阶插值、减缩积分000如果在厚度方向单元过少，将会出现沙漏现象000如果在厚度方向有足够的单元，OK非协调模式000如果不是过渡扭曲，OKABAQUS/analysis_单元选择标准应力集中A1.29ABAQUS/analysis_单元选择标准应力集中•存在应力集中的问题中，二阶单元相对一阶单元的优势非常明显；并且二阶单元非常适合于模拟（静态）裂纹。–全积分单元和减缩积分单元都工作得很好。–一般减缩积分单元的效率更高—在较低的计算费用的前提下，结果同全积分单元结果一样好，甚至优于全积分单元。A1.30ABAQUS/analysis_单元选择标准Physicalmodel一次单元模型—单元表面为直线片断二次单元模型—单元表面为二次曲线应力集中•相对一阶单元，二阶单元可以利用较少的单元捕获曲边等几何特征。A1.31ABAQUS/analysis_单元选择标准应力集中•当初始形状扭曲时，一阶和二阶四边形和六面体单元的精度都会下降。–一阶单元对扭曲的敏感性低于二阶单元，因此对出存在明显网格扭曲的问题，一阶单元是更好的选择。•相对于几乎所有的其它单元，二阶三角形和四面体单元对于初始的单元形状敏感性较低；然而，良好形状的单元将得到更好的结果。理想的好的坏的扭曲的非扭曲的A1.32ABAQUS/analysis_单元选择标准ellipticalshape应力集中•典型的应力集中问题，NAFEMS基准问题。–下面是采用不同类型单元进行分析得到的结果。PA1.33ABAQUS/analysis_单元选择标准应力集中•一阶单元(包括非协调模式单元)在应力集中问题中的效果较差。•二阶单元，如CPS6,CPS8,和CPS8R可以得到较好的结果。单元类型satDyy(Target=100.0)粗网格细网格CPS355.0676.87CPS471.9891.2CPS4I63.4584.37CPS4R43.6760.6CPS696.12101.4CPS891.2100.12CPS8R92.5697.16A1.34ABAQUS/analysis_单元选择标准应力集中•在应力集中区域，形状良好、二阶、减缩积分四边形和六面体单元将得到较好的精度。–在这些区域，扭曲单元将降低精度。ABAQUS/analysis_单元选择标准接触A1.36ABAQUS/analysis_单元选择标准接触•几乎所有的单元类型都可以很好的处理接触问题，但有以下几个例外：–二阶四边形/六面体单元–“常规的”二阶三角形/四面体（与名字以“M”结尾的、“修正的”三角形/四面体单元相对应）单元、二阶楔形单元和六节点壳和薄膜单元。•这些单元将导致收敛困难。对于二阶四面体、楔型体和六节点壳和薄膜单元，由压力载荷转化为一致节点力的方向与基于罚函数接触算法的