(完整版)ABAQUS实体单元类型总结
Abaqus单元类型选择 ppt课件

A1.7
Abaqus单元类型选择
ABAQUS/analysis_单元选择 标准
公式 • 用于描述单元行为的数学公式是用于单元分类的另一种方法。 • 不同单元公式的例子:
– 平面应变 – 平面应力 – 杂交单元 – 非协调元 – 小应变壳 – 有限应变壳 – 厚壳 – 薄壳
构尺寸的1/10,比如: – 支撑或点载荷之间的距离 – 尺寸变化很大的横截面之间的距离 – 最高振动模态的波长
A1.14
Abaqus单元类型选择
• 壳单元 – 使用表面模型构成的壳单元近似 模拟三维实体连续体单元。
• 可以有效的模拟弯曲和面内 变形。
– 如果需要分析某个区域的细节, 使用多点约束或子模型的办法可 以将局部的三维实体模型加入到 壳单元模型中。
– 对于具有线弹性材料属性的、 未扭转的单元,精确积分应变 能所需的最小积分阶数。
• 减缩积分:
– 积分的阶数比全积分小一阶。
一次插值
二次插值
ABAQUS/analysis_单元选择 标准
全积分
减缩积分
A1.10
Abaqus单元类型选择
• 单元命名约定:例子
B21: Beam, 2-D, 1st-order interpolation
ABAQUS/analysis_单元选择 标准
3-D 实体
线模型
利用梁单元建模的框架结构
A1.16
ABAQUS/analysis_单元选择标准
Abaqus单元类型选择
Abaqus单元类型选择
• 纯弯曲的物理特征 – 有限元方法企图模拟的材料行为是: • 在变形过程中,横截面仍然保持为 平面。 • 沿厚度方向,轴向应变xx 线性变化。
abaqus中实体与实体壳单元的区别

abaqus中实体与实体壳单元的区别
在Abaqus中,实体单元和实体壳单元是两种不同的有限元单元类型。
实体单元是一个三维单元,用于建模实体结构,例如实心体或复杂的几何形状。
它具有六个自由度(三个平移和三个旋转自由度),可以用来模拟实体物体的力学行为,例如固体的应力和应变分析。
实体壳单元是一种专用于建模薄壳结构的单元。
实质上,它将薄壳结构简化为一个二维平面,可以有效地模拟薄壳结构的行为。
它通常用于模拟平面板、屋盖、车身等具有表面功能的结构。
实体壳单元只有三个平移自由度(在任意平移方向上),并且不具有旋转自由度。
这意味着实体壳单元不能够准确地模拟薄壳结构中的扭转和转动效应。
综上所述,实体单元和实体壳单元在应用和模拟能力上有所不同。
实体单元更适用于建模实心体和复杂几何结构,而实体壳单元则适用于模拟薄壳结构的应力和变形行为。
如何选择ABAQUS单元类型

1、按照节点位移插值的阶数,可以将ABAQUS单元分为线性单元、二次单元和修正的二次单元2、线性完全积分单元在承受弯曲载荷时会出现剪切自锁,造成单元过于刚硬,即使划分很细的网格,计算精度仍然很差3、二次完全积分单元适于模拟应力集中问题,一般情况下不会出现剪切自锁,但不能在接触分析和弹塑性分析中使用4、线性减缩积分单元对位移的求解结果较精确,在弯曲载荷下不容易发生剪切自锁,网格的扭曲变形(例如Quad单元的角度远远大于或小于90°)对其分析精度影响不大,但这种单元需要划分较细的网格来克服沙漏问题,且不适于求解应力集中部位的节点应力5、二次减缩积分单元不但支持了线性减缩积分单元的优点,而且不划分很细的网格也不会出现严重的沙漏问题,即使在复杂应力状态下,对自锁问题也不敏感,但它不适于接触分析和大应变问题6、非协调模式单元克服了剪切自锁问题,在单元扭曲比较小的情况下得到的位移和应力结果很精确,但如果所关心部位的单元扭曲比较大,其分析精度会降低7、线性Tri单元和Tet单元的精度很差,二次Tet单元(C3D10)适于ABAQUS/Standand中的小位移无接触问题,修正的二次Tet单元(C3D10M)适于ABAQUS/Explicit,以及ABAQUS/Standand中的大变形和接触问题8、ABAQUS的壳单元可以有多种分类方法,按照薄壳和厚壳来划分,可以分为通用目的(general-purpose)壳单元和特殊用途(special-purpose)壳单元;按照单元的定义方式,可以分为常规(conventional)壳单元和连续体(continuum)壳单元9、ABAQUS中的所有梁单元都可以产生轴向变形、弯曲变形和扭转变形,B21和B31单元(线性梁单元)以及B22和B32单元(二次梁单元)即适用于模拟剪切变形引起重要作用的深梁,又适用于模拟剪切变形不太重要的细长梁,三次单元B23和B33只需划分很少的单元就可以得到较精确的结果1、对于应力集中问题,尽量不要使用线性减缩积分单元,可使用二次单元来提高精度。
abaqus 实体单元剪应力

abaqus 实体单元剪应力
Abaqus是一款工程仿真软件,广泛应用于材料、结构和流体动力学等领域。
在Abaqus中,实体单元剪应力是指在不同方向上的剪切应力分量。
实体单元类型:在Abaqus中,实体单元通常包括三维实体单元、二维实体单元和一维实体单元。
这些实体单元可以用来模拟各种材料和结构的力学行为。
剪应力计算:剪应力是三维空间中的应力分量之一,它与正应力(法向应力)和切应力(切向应力)一起描述了物体的应力状态。
在Abaqus中,可以通过在应力分析模块中查看剪应力分量。
剪应力应用:剪应力在工程中的应用非常广泛,例如在结构稳定性分析、材料屈服强度分析和疲劳寿命预测等方面都有重要应用。
通过在Abaqus中进行剪应力分析,可以更好地了解物体在不同受力条件下的行为。
剪应力影响因素:剪应力的影响因素很多,包括材料性质、几何形状和边界条件等。
例如,对于金属材料,剪应力通常会随着应变的增加而逐渐增大;而对于陶瓷材料,剪应力则可能在应变较小时就达到饱和。
剪应力控制方法:在实际工程中,为了控制结构的力学性能和稳定性,通常需要对剪应力进行合理控制。
例如,可以通过优化结构设计、选择合适的材料和加工方法等措施来降低和分散剪应力集中部位的力量。
总之,在Abaqus中,实体单元剪应力是指在不同方向上的剪切应力分量,它对于理解和预测材料的力学行为具有重要意义。
ABAQUS简支梁分析报告(梁单元和实体单元)

基于ABAQUS简支梁受力和弯矩的相关分析(梁单元和实体单元)对于简支梁,基于 ABAQUS2016,首先用梁单元分析了梁受力作用下的应力,变形,剪力和力矩;对同一模型,并用实体单元进行了相应的分析。
另外,还分析了梁结构受力和弯矩作用下的剪力及力矩分析。
对于CAE仿真分析具体细节操作并没有给出详细的操作,不过在后面上传了对应的cae,odb,inp文件。
不过要注意的是本文采用的是ABAQUS2016进行计算,低版本可能打不开,可以自己提交inp文件自己计算即可。
可以到小木虫搜索:“基于ABAQUS简支梁受力和弯矩的相关分析”进行相应文件下载。
对于一简支梁,其结构简图如下所示,梁的一段受固支,一段受简支,在梁的两端受集中载荷,梁的大直径D=180mm,小直径d=150mm,a=200mm,b=300mm,l=1600mm,F=300000N。
现通过梁单元和实体单元分析简支梁的受力情况,变形情况,以及分析其剪力和弯矩等。
材料采用45#钢,弹性模量E=2.1e6MPa,泊松比v=0.28。
图1 简支梁结构简图1.梁单元分析ABAQUS2016中对应的文件为beam-shaft.cae ,beam-shaft.odb,beam-shaft.inp。
在建立梁part的时候,采用三维线性实体,按照图1所示尺寸建立,然后在台阶及支撑梁处进行分割,结果如图2所示。
图2 建立part并分割接下来为梁结构分配材料,创建材料,定义弹性模量和泊松比,创建梁截面形状,如图3,非别定义两个圆,圆的直接分别为180和150mm。
然后创建两个截面,截面选择梁截面,再选择图2中的所有梁,定义梁的方向矢量为(0,0,-1)(点击图3中的n2,n1,t那个图标即可创建梁的方向矢量),最后把创建好的梁赋给梁结构。
图3 创建梁截面形状接下来装配实体,再创建分析步,在创建分析步的时候,点击主菜单栏的Output,编辑Edit Field Output Request,在SF前面打钩,这样就可以在结果后处理中输出截面剪力和力矩,如图4所示。
ABAQUS单元类型

ABAQUS单元类型Advanced Finite Element Analysis–And ApplicationsDaming Zhang, Ph.D.Associate Professor of Transportation SystemsDepartment of Industrial TechnologyCollege of Agricultural Sciences and TechnologyCalifornia State University, FresnoMay 27, 2009Dr. Daming Zhang -Cal State Univ Fresno1Advanced Finite Element Analysis -And ApplicationsLecture 4:ABAQUS Element LibraryDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications2ABAQUS Elements OverviewA wide range of elements available for solving different problemsFive characteristics of an element:–Family–Degree of freedom–Number of nodes–Formulation–IntegrationUnique name: T2D2, S4R, C3D8Iused on the *ELEMENT option, TYPE parameterDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications3FamilyUsed to distinguish the geometryIndicated by first letter or letters: S4R, C3D8I, CINPE4 Dr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications4Degree of FreedomDegree of Freedom (DOF) are the fundamental variables 1Translation in direction 12Translation in direction 23Translation in direction 34Rotation about the 1-axis5Rotation about the 2-axis6Rotation about the 3-axis7Warping in open-section beam elements8Acoustic pressure or pore pressure9Electric potential11Temperature12+Temperature at other points through the thickness of beamsand shellsDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications5Number of Nodes ?Determines the interpolation orderfirst order, second order, …Clearly identified in the name: C3D8, S8RBeam family indicating order of interpolation: B31, B32 Dr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications6FormulationRefers to the mathematical theoryLagrangian:material descriptionEulerian:Spatial descriptionShell family has 3 classes:General purposethin-onlythick-onlyAlternative formulations (end of element name)Hybrid formulation: C3D8H, B31HIncompatible formulation: C3D8IDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications7IntegrationGaussian Quadrature to integrate quantities over the volume of each elementFull or Reduced integrationUse “R” at the end of element name to distinguish the reduced-integration elements: CAX4, CAX4RWill significantly affect the accuracyDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications8Continuum Elements (1)Used to model the widest variety of components ?Element names begin with “C”Next two indicate the dimensionality: 3D, PE, PS, AX ?The last shows the degree of freedom3D continuum elements:hexa, penta, tetra2D continuum elements:plane strainplane stressaxisymmetricshape: quadrilateralor triangularDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications 9Continuum Elements (2)2D continuum elements must be defined in 1-2 plane Node order should be counterclockwiseElement normals must all pointed at same direction ?Degree of freedom: translational DOFsElement properties: *SOLID SECTIONFormulation & Integration: Incompatible, Hybrid, Reduced ?Output variable: default in global coordinate system *ORIENTATION to define local coordinate systemDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications10Shell Elements (1)Used to model structure with one dimension small ?Element name begins with: S, SAX, SAXAThe first number indicates the number of nodesIf the last character is “5”, the element doesn’t use the rotational DOF around normal of middle plane ?Quadrilateral or triangular; Linear or quadratic elements ?Three different formulationsDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications 11Shell Elements (2)Degree of freedom: default 6 DOF, but also 5 DOF: S4R5 Axisymmetric shells have 3DOF:1Translation in the r-direction.2Translation in the z-direction.6Rotation in the r-z plane.*SHELL GENERAL SECTION: you define the properties*SHELL SECTION: ABAQUS calculates section properties ?Formulation & Integration: complicated, check before use ?Output variable: defined in the local material directions Lie on the surface of each shell elementAxes rotate with the element’s deformation in large-displacement simulationsDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications12Beam ElementsUsed to model structure with one dimension is quite large ?Element name begins with: B (e.g. B31)The first number indicates the dimensionalityThe third character indicates the interpolation order ?Degree of freedom: default 6 DOF, 2D beams have 3 DOF Open-section beams (e.g. B31OS) have DOF 7 for warping ?*BEAM GENERAL SECTION: you define the properties *BEAM SECTION: ABAQUS calculates section properties ?Formulation & Integration: Hybrid for very slender beams;B21, B31, B22, B32: shear deformable, and finite axial strain; B23 and B33 are not; Open section: B31OS, B32OS ?Output variable: axial stress (s 11), shear stress (s 12)Dr. Daming Zhang -Cal State Univ FresnoAdvanced Finite Element Analysis -And Applications13Truss ElementsModel rods that can carry only tensile or compressive loads ?Element name begins with: T (e.g. T2D3, T3D2) The next two characters indicates the dimensionality ?Thefinal character indicates the number of nodes ?Degree of freedom: has only translational DOFs ?*SOLID SECTION: specify the material properties The cross-sectional area is given on the data line ?Formulation & Integration: Hybrid for very rigid links ?Output variable: Axial stress and strainDr. Daming Zhang -Cal State Univ FresnoAdvanced Finite Element Analysis -And Applications14Rigid ElementsElement name begins with: R (e.g. R3D4, R3D3)The next two characters indicates the dimensionality ?The final character indicates the number of nodes ?The nodes have no independent degrees of freedomThe nodes defining rigid elements can have loads applied to them or can be connected to other elements but they cannot have any boundary conditions ?*RIGID BODY defines the rigid body reference node ?Pay attention to the ‘sides’ of the rigid body elements ?Formulation & Integration: none ?Output: motion onlyDr. Daming Zhang -Cal State Univ FresnoAdvanced Finite Element Analysis -And Applications15Continuum Elements OverviewThe biggest family with over 20 just for 3D models ?3D: Hexa, Penta, Tetra; 2D: triangles and quadrilaterals ?Linear and quadratic versions for each of these shapes ?Full-and reduced-integration elements for hexa and quad ?Standard or hybrid element formulationFor linear hexa or quad: incompatible mode formulation ?For quadratic tria or tetra: "modified" formulationThe accuracy of your simulation will depend strongly on the type of element you use in your modelDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications 16Full IntegrationThe accuracy of Gaussian Quadrature is (2n-1) for n=4?The Element Stiffness Matrix is calculated by:Fully integrated linear elements use two integration points in each directionFully integrated quadratic elements use three integration points in each directionDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications 17k []m e=B []V ∫TD []B []dVFull Integration ExampleUse a cantilever beam to show the accuracy of analysisDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications 18Shear LockingHappened on fully integrated, first-order, solid elements causes the elements to be too stiff in bending ?Deformation of material subjected to bending moment MDeformation of a fully integrated, linear element Dr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications 19Reduced IntegrationOnly quadrilateral and hexahedral elements can use a reduced-integration schemeuse one fewer integration point in each direction than the fully integrated elementsDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications 20Reduced Integration ResultsLinear reduced-integration elements tend to be too flexible But fine mesh will produce acceptable results ?Deformation of a linear element with reduced integrationDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications21Incompatible mode elementsAn attempt to overcome the problems of shear locking in fully integrated first-order elementsAdditional degrees of freedom enhance the element's deformation gradients as linear variationcan produce results in bending problems that are comparable to quadratic elements but at significantly lower computational costThe mesh distortion should be minimized as much as possible to improve the accuracy of the resultsDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications22Hybrid ElementsHybrid elements have the letter "H" in their names ?Hybrid elements are used when the material behavior is incompressible (Poisson's ratio = 0.5)The volume cannot change if thematerial is incompressibleThe pressure stress cannot becomputed from the displacementsof the nodesHybrid elements include an additional degree of freedom that determines the pressure stress in the element directly Dr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications23Selecting Continuum ElementsUse quadratic, reduced-integration elements (CAX8R, CPE8R, CPS8R, C3D20R, etc.) for general analysis workUse quadratic, fully integrated elements (CAX8, CPE8, CPS8, C3D20, etc.) locally where stress concentrations may exist Use a fine mesh of linear, reduced-integration elements (CAX4R, CPE4R, CPS4R, C3D8R, etc.) for large-strain analysis For contact problems use a fine mesh of linear, reduced-integration elements or incompatible elements (CAX4I, CPE4I, CPS4I, C3D8I, etc.)?Minimize the mesh distortion as much as possibleIn three dimensions use hexahedral (brick-shaped) elements wherever possible; Use C3D6 and C3D4 only when necessary ?modified quadratic tetrahedral element (C3D10M) is robust for large-deformation and contact problems and exhibits minimal shear and volumetric lockingDr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications24Thank You Dr. Daming Zhang -Cal State Univ Fresno Advanced Finite Element Analysis -And Applications25。
ABAQUS总结
ABAQUS总结1.在Help中搜索关键词时,⽤引号括起来表⽰把这些词当做整体搜索;不⽤引号,这些词将被分别搜索。
2.点击要删除的线时,按住Shift可以选择多条线同时操作。
要取消对某条线的选择,可以按住Ctrl再点击此线。
3.完成操作时,可以点击Done,或直接点击⿏标中键。
4.平⾯应⼒、平⾯应变问题的截⾯属性为Solid(实⼼体),不是Shell(壳)。
5.ABAQUS/CAE不把材料特性直接赋予单元或实体,⽽是先在Section(截⾯属性)中定义材料特性,再Assign Section(赋予截⾯属性),点取截⾯分别赋予。
6.Load(荷载类型)中的Pressure指单位⾯积上的⼒,压⼒为正,拉⼒为负。
7.Load⾥的BC指的是Boundary Condition(边界条件)。
8.遇到傻×死板的Assign Element Type时,调整显⽰器显⽰为纵向。
9.对话框底部经常出现Cancel和Dismiss,它们都是关闭当前对话框,区别在于:Dismiss出现在包含只读数据的对话框中;Cancel出现在允许修改的对话框中,但此时点Cancel表⽰不作修改⽽退出。
10.⼀个ABAQUS/CAE主窗⼝只能显⽰⼀个模型数据库。
如果想同时显⽰多个模型数据库,可以同时启动多个主窗⼝。
11.i nstance——实体。
所谓实体,是part在Assembly中的⼀种映射,⼀个部件可以对应多个实体。
材料和截⾯属性定义在part(部件)上,interaction(相互作⽤)、BCs(边界条件)Load(荷载)定义在instance(实体)上;mesh(⽹格)可以定义在部件或实体上。
12.A BAQUS中,Section指截⾯属性,包含⼴义的部件特性,⽽不是平时的梁板截⾯形状。
在ABAQUS中,梁板截⾯形状称为Profile。
13.S tep:默认的time period(分析步时间)是1,Nlgeom(⼏何⾮线性)是Off,如果模型中存在⼤的位移或转动,应设置Nlgeom为On。
abaqus系列教程-04应用实体单元
4. 应用实体单元在ABAQUS中,应力/位移单元的实体(continuum)单元族是包含最广泛的。
ABAQUS/Standard和ABAQUS/Explicit的实体单元库多少有所不同。
ABAQUS/Standard实体单元库ABAQUS/Standard的实体单元库包括二维和三维的一阶(线性)插值单元和二阶(二次)插值单元,它们应用或者完全积分或者减缩积分。
二维单元有三角形和四边形;在三维单元中提供了四面体、三角楔形体和六面体(砖型)。
也提供了修正的二阶三角形和四面体单元。
此外,在ABAQUS/Standard中还有杂交和非协调模式单元。
ABAQUS/Explicit实体单元库ABAQUS/Explicit的实体单元库包括二维和三维的减缩积分一阶(线性)插值单元,也有修正的二阶插值三角形和四面体单元。
在ABAQUS/Explicit中没有完全积分或者规则的二阶单元。
关于可选用的实体单元的详细信息,请参阅ABAQUS分析用户手册第14.1.1节“Solid (continumm) elements”。
当做出所有这些各种选项的排列,发现可供使用的实体单元的总数是相当大的,仅就三维模型而言就超过了20种。
模拟的精度将很大程度上依赖于在模型中采用的单元类型。
在这些单元中选择哪一个最适合于你的模型,可能是一件令人苦恼的事情。
特别是在初次使用时。
然而,你会逐渐认识到这种在20多件工具组中的选择,为你提供一种能力,对于一个特殊的模拟能够选择恰当正确的工具或单元。
本章讨论了不同的单元数学描述和积分水平对于一个特定分析的精度的影响,也给出了一些关于选择实体单元的一般性指导意见,这些为你积累ABAQUS的应用经验,并建立自己的知识库提供了的基础。
本章末尾的例子,当你建立和分析一个连接环构件模型时,将允许你应用这些知识。
4.1 单元的数学描述和积分通过考虑一个静态分析的悬臂梁,如图4-1所示,将演示单元阶数(线性或二次)、单元数学描述和积分水平对结构模拟的精度的影响。
abaqus 立方体 长方体 单元形状
abaqus中的有限元分析中,常用的单元形状包括立方体和长方体。
这两种单元形状在工程实践中有着广泛的应用,可以用来模拟不同类型的结构和材料行为。
下面我们将分别介绍abaqus中立方体和长方体单元的特点和应用。
1. 立方体单元立方体单元是abaqus中常用的一种三维单元形状,它具有以下特点:1)简单:立方体单元的几何形状和网格划分都比较简单,易于生成和处理。
2)稳定性好:立方体单元在受力情况下具有较好的稳定性,适用于不同类型的加载和边界条件。
3)适用范围广:立方体单元适用于各种材料和结构类型的有限元分析,包括钢结构、混凝土结构、土体等。
在工程实践中,立方体单元可用于模拟各种结构的应力、应变分布,分析结构的承载能力和变形情况。
立方体单元还可用于模拟材料的断裂、变形和破坏过程,对材料的性能进行评估和预测。
2. 长方体单元长方体单元是abaqus中另一种常用的三维单元形状,它具有以下特点:1)适用于不规则结构:长方体单元在处理不规则几何形状和复杂结构时具有较好的适用性,可以更精确地模拟实际工程中的结构。
2)网格划分灵活:长方体单元的网格划分相对灵活,可以根据实际需要进行细化或简化,适用于对结构细节要求较高的分析。
3)适用范围广:长方体单元适用于复杂结构和材料的有限元分析,包括汽车车身、航空航天结构、复合材料等。
在工程实践中,长方体单元可用于模拟复杂结构的应力、应变分布,分析结构的疲劳、强度和稳定性。
长方体单元还可用于模拟复合材料的层合结构和界面效应,对材料和结构的性能进行更精确的预测和评估。
abaqus中的立方体和长方体单元形状具有各自的特点和应用范围,在工程实践中都有着重要的作用。
合理选择和使用这两种单元形状,可以更准确地模拟和分析不同类型的结构和材料行为,为工程设计和科学研究提供可靠的分析基础。
立方体和长方体单元作为abaqus中常用的三维单元形状,在工程实践中有着广泛的应用。
它们不仅可以用来模拟结构的应力、应变分布,还可以用于分析结构的疲劳、强度和稳定性。
最新abaqus单元属性小结
最新abaqus单元属性小结
S4R 单元性能稳定,适用范围很广
对于复合材料,为模拟剪切变形的影响,应使用适于厚壳的单元(例如S4、S4R、S3、S3R、
S8R),并要注意检查截面是否保持平面。
对于几何非线性分析,在ABAQUS/Standard中的小应变壳单元(S4R5, S8R, S8R5, S8RT, S9R5, STRI3, 和STRI65)使用总体拉格朗日应变算法,应力应变可以相对于参考构型的材料方向改定。
垫片单
元是小
应变小位移单元,默认情况下其应力应变值也是以初始参考构型定义的行为方向输出。
对于有限膜应变单元(所有的膜单元以及S3/S3R, S4, S4R, SAX,和SAXA单元)和在ABAQUS/Explicit 中的小应变单元,其材料方向是随着曲面的平均刚性旋转运动而变以形成当前构型的材料方向。
此时这些单元的应力应变则是根据当前的参考构型中的材料方向给出的。
(更详细地说明可以参考ABAQUS相关手册)。
用户可以决定与*section print和*section file相关的局部坐标系统是固定不动还是随着曲面的平均刚性运动而旋转。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
在ABAQUS中,基于应力/位移的实体单元类型最为丰富:
(1)在ABAQUS/Sandard中,实体单元包括二维和三维的线性单元和二次单元,均可以采用完全积分或缩减积分,另外还有修正的二次Tri单元(三角形单元)和Tet单元(四面体单元),以及非协调模式单元和杂交单元。
(2)ABAQUS/Explicit中,实体单元包括二维和三维的线性缩减积分单元,以及修正的二次二次Tri单元(三角形单元)和Tet单元(四面体单元),没有二次完全积分实体单元。
------------------------------------------------------------------------------------------------------------
按照节点位移插值的阶数,ABAQUS里的实体单元可以分为以下三类:
线性单元(即一阶单元):仅在单元的角点处布置节点,在各个方向都采用线性插值。
二次单元(即二阶单元):在每条边上有中间节点,采用二次插值。
修正的二次单元(只有Tri 或Tet 才有此类型):在每条边上有中间节点,并采用修正的二次插值。
******************************************************************************* ***************
1、线性完全积分单元:当单元具有规则形状时,所用的高斯积分点的数目足以对单元刚度矩阵中的多项式进行精确积分。
缺点:承受弯曲载荷时,会出现剪切自锁,造成单元过于刚硬,即使划分很细的网格,计算精度仍然很差。
2、二次完全积分单元:
优点:
(1)应力计算结果很精确,适合模拟应力集中问题;
(2)一般情况下,没有剪切自锁问题(shear locking)。
但使用这种单元时要注意:
(1)不能用于接触分析;
(2)对于弹塑性分析,如果材料不可压缩(例如金属材料),则容易产生体积自锁(volumetric locking);
(3)当单元发生扭曲或弯曲应力有梯度时,有可能出现某种程度的自锁。
3、线性减缩积分单元:
减缩积分单元,比普通的完全积分单元在每个方向少用一个积分点;
线性缩减积分单元:
只在单元的中心有一个积分点,由于存在沙漏数值问题(hourglass)而过于柔软。
采用线性缩减积分单元模拟承受弯曲载荷的结构时,沿厚度方向上至少应划分四个单元。
优点:
(1)对位移的求解计算结果较精确;
(2)网格存在扭曲变形时(例如Quad 单元的角度远远大于或小于90º),分析精度不会受到明显的影响;
(3)在弯曲载荷下不易发生剪切自锁。
缺点:
(1)需要较细网格克服沙漏问题;
(2)如果希望以应力集中部位的节点应力作为分析目标,则不能选用此单元。
——因为线性缩减积分单元只在单元的中心有一个积分点,相当于常应力单元,在积分点上的应力结果实相对精确的,而在经过外插值和平均后得到的节点应力则不精确。
4、二次减缩积分单元
不但保持线性减缩积分单元的上述优点,还具有如下特点:
(1)即使不划分很细的网格也不会出现严重的沙漏问题;
(2)即使在复杂应力状态下,对自锁问题也不敏感。
使用这种单元要注意:
(1)不能用于接触分析;
(2)不能用于大应变问题;
(3)存在与线性减缩积分单元类似的问题,由于积分点少,得到的节点应力的精度往往低于二次完全积分单元。
5、非协调模式单元(imcompatible modes)
——仅在ABAQUS/Standard 有,可克服线性完全积分单元中的剪切自锁问题。
ABAQUS中的非协调模式单元和MSC.NASTRAN中的4节点四边形单元或8节点六面体单元很相似,所以在比较着两种有限元软件的计算结果时会发现,如果在ABAQUS中选择了非协调模式单元,得到的分析结果会和MSC.NASTRAN的结果一致。
优点:
(1)克服了剪切自锁问题,在单元扭曲比较小的情况下,得到的位移和应力结果很精确;(2)在弯曲问题中,在厚度方向上只需很少的单元,就可以得到与二次单元相当的结果,而计算成本却明显降低;
(3)使用了增强变形梯度的非协调模式,单元交界处不会重叠或开洞,因此很容易扩展到非线性、有限应变的位移。
但使用这种单元时要注意:如果所关心部位的单元扭曲比较大,尤其是出现交错扭曲时,分析精度会降低。
6、使用Tri 或Tet 单元要注意:
如果能用Quad 或Hex 单元,就尽量不要使用Tri或Tet 单元;
(1)线性Tri 或Tet 单元的精度很差,不要在模型中所关心的部位及其附近区域使用;
(2)二次Tri 或Tet 单元的精度较高,而且能模拟任意的几何形状,但计算代价比Quad 或Hex 单元大。
(3)二次Tet 单元(C3D10)适于ABAQUS/Standard 中的小位移无接触问题;
修正的二次Tet 单元(C3D10M)适于ABAQUS/Explicit 和ABAQUS/Standard 中的大变形和接触问题;
(4)使用自有网格不易通过布置种子来控制实体内部的单元大小。
7、杂交单元
在ABAQUS/Standard 中,每一种实体单元都有其对应的杂交单元,
用于不可压缩材料(泊松比为0.5,如橡胶)或近似不可压缩材料(泊松比大于0.475)。
除了平面应力问题之外,不能用普通单元来模拟不可压缩材料的响应,因为此时单元中的应力士不确定的。
ABAQUS/Explicit 中没有杂交单元。
******************************************************************************* ******************
混合使用不同类型的单元:
1、当三维实体几何形状复杂时,无法再整个实体上使用structure结构化网格或sweep扫略网格划分技术得到Hex单元网格,一种常用的做法是:
(1)对实体不重要的部分使用Free自由网格划分技术,生成Tet单元网格,而对于所关心的部分采用结构化网格或扫略网格划分技术,生成Hex单元网格。
(2)在生成这样的网格时,ABAQUS会给出提示信息,提示将生成非协调的网格,在不同单元类型的交界处将自动创建Tie绑定约束。
2、需要注意的是,在不同单元类型网格的交界处,即使单元角部节点是重合的,仍然有可能出现不连续的应力场,而且在交界处的应力可能大幅度的增大。
如果在同一实体中混合使用线性和二次单元,也会出现类似的问题。
因此在混合使用不同类型单元时,应确保其交界处远离所关心的区域,并仔细检查分析结果是否正确。