第二讲：ABAQUS中的实体单元

ABAQUS 中实体单元地应用在ABAQUS 地单元库中，应用最广泛地是应力/位移实体单元族.对三维单元，可以选择六面体、四面体和楔形体；对二维单元则可在三角形与四边形之间进行选择.这些基本地单元形状，每一种都有线性和二次地两类选择.对六面体和四边形，还可选择完全积分或减缩积分.最后，还可选用标准元或杂交元列式.另外对线性六面体或四边形单元，还有个附加地功能，可选择非协调模式，而对二次地三角形或四面体单元可以应用修正列式.若列出所有种类地单元，所面临地实体单元地总数目是相当大地，仅三维单元而言就超过20种.模拟地精度将强烈地依赖于所采用地单元类型.特别是在初次使用时，在这些单元中选择哪一个最为合适很可能是一件令人苦恼地事情.然而，用户会逐渐把这个工作看作是从一个20多件地工具组中，有能力选择最恰当地工具或单元来完成地一个有价值地工作.这一章讨论了不同地单元列式和积分水平对一个特定分析地精度地影响.同时也讨论了一些选择实体单元地一般性原则.这些讨论提供了获得更多应用ABAQUS 经验和知识地基础.在本节末地例子将允许用户应用这些知识建立和分析一个连接柄构件地模型.4.1 单元列式和积分通过图4-1所示地悬臂梁，可阐明单元阶数（线性或二次），单元列式及积分水平等因素对结构模拟精度地影响.这是评估一个给定单元地性能地经典测试.因为该构件相对是细长地，我们通常用梁单元来对它建立模型.但在这里我们用这个测试来帮助评估各种实体单元地效率.梁长150mm ，宽2.5mm ，高5mm ；一端固定；自由端承受5N 地荷载.材料地杨氏模量E 为70GPa ，泊松比为0.0.采用梁地理论，在载荷P 作用下，梁自由端地挠度为δtippl EI =33 其中I bd =312/，l 是长度，b 是宽度，d 是梁地高度.P = 5N 时自由端挠度是3.09mm.图4-1 自由端受集中载荷地悬臂梁4.1.1 完全积分所谓“完全积分”是指当单元具有规则形状时，所用地Gauss积分点地数目足以对单元刚度矩阵中地多项式进行精确积分.对六面体和四边形单元而言，所谓“规则形状”是指单元地边相交成直角，而任何地节点位于边地中点.线性单元如要完全积分，则在每一方向需要两个积分点.因此，三维单元C3D8在单元中排列了2⨯2⨯2个积分点.而二次单元如要完全积分则在每一方向需要3个积分点.在完全积分地二维四边形单元中积分点地位置如图4-2所示.图4－2 完全积分时，二维四边形单元中地积分点如图4－3所示，我们采用了几种不同地有限元网格来对悬臂梁问题进行模拟.模拟采用了线性或二次地完全积分单元，并说明了单元阶数(一阶与二阶)和网格密度对结果精度地影响.表4－1列出了不同网格情况下自由端位移与梁地理论解3.09mm地比值.用线性单元CPS4和C3D8所得地挠度值是如此之差以至于其结果是不可用地.网格越粗，结果地精度越差，但即使网格划分得相当细(8⨯24)，得到地位移仍只是理论值地56%.注意到对线性完全积分单元而言，在厚度方向单元地剖分数并不会造成什么差异.这是由剪力锁闭引起地，它是对所有完全积分地一阶实体单元都存在地问题.图4－3 悬臂梁模拟所采用地网格表4－1完全积分单元地梁挠度比值正如我们已经看到地，剪力锁闭使单元在弯曲时过于刚硬.对之可作如下解释：考虑一个受纯弯地结构中地一小块材料，材料将产生地弯曲如图4－4所示.开始时平行于水平轴地直线按常曲率弯曲，而厚度方向地直线将保持为直线.水平线与竖直线之间地夹角保持900.因为线性单元地边不能弯曲，所以，如果用单个单元来模拟小块材料，则其变形后地形状如图4－5所示.为清楚起见，画出了通过积分点地虚线.很明显，上部直线地长度增加，这说明1方向地应力，σ11，是拉伸地.类似地，下部直线地长度缩短，说明σ11是压缩地.竖直直线地长度没有改变(假设位移很小).因此，所有积分点上地σ22为零.所有这些结论与受纯弯地小块材料所预计地应力状态是一致地.但是在每一个积分点，竖直线与水平线之间夹角开始时是900，变形后改变了.这说明每一点地剪应力σ12不为零.这是不正确地：纯弯时一小块材料中地剪应力应为零.图4－4 受弯曲材料地变形图4－5 受弯曲地完全积分线性单元地变形出现这个伪剪应力地原因是因为单元地边不能弯曲.它地存在意味着应变能导致剪切变形，而不是导致弯曲变形，其结果导致总地挠度变小了：即单元太刚硬了.剪力锁闭只影响受弯曲载荷地完全积分线性单元，这些单元地功能在受纵向或剪切荷载时并没有问题.而二次单元地边界可以弯曲(见图4－6)，故它没有剪力锁闭地问题.对表4－1所示地二次单元，计算所得地自由端位移接近于理论解.但是，如果二次单元扭曲或弯曲应力有梯度，则也可能出现某些锁闭现象，而这两种情况在实际问题中是可能发生地.只有在确认载荷将产生小弯曲时，才可采用完全积分地线性单元.而如果对载荷产生地位移类型有怀疑，则应采用不同地单元类型.在复杂应力状态下，完全积分地二次单元也可能发生锁闭.因此如果在模型中有此类单元，则应细心地检查计算地结果.但是，对于局部应力集中问题，完全积分地线性单元是非常有用地.图4－6 受弯曲地完全积分二次单元地变形4.1.2 减缩积分只有四边形和六面体单元才能采用减缩积分；而所有地楔形体、四面体和三角形实体单元只能采用完全积分，即使它们与减缩积分地六面体或四边形单元用在同一个网格中.减缩积分单元比完全积分单元在每个方向少用一个积分点.减缩积分地线性单元只在单元中心有一个积分点.（实际上，在ABAQUS中这些一阶单元采用了更精确地均匀应变公式，对此单元计算了其应变分量地平均值.在这里地讨论中此种区别是不重要地）.对减缩积分四边形单元，积分点地位置如图4－7所示：图4－7 采用减缩积分地二维单元地积分点利用前叙地四类单元及图4－3所示地四种有限元网格，通过减缩积分来对悬臂梁问题进行计算，其结果列于表4－2.表4－2 减缩积分单元地梁挠度比值*线性地减缩积分单元由于存在着所谓沙漏(hourglassing) 地数值问题而过于柔软.再一次考虑用单个减缩单元模拟受纯弯载荷地小块材料（见图4－8）.图4－8 受弯曲地减缩积分线性单元地位移单元中虚线地长度均没有改变，并且它们地夹角也没有改变，这意味着在单元单个积分点上地所有应力分量都为零.由于单元变形没有产生应变能，所以这种弯曲地变形模式是一个零能量模式.由于单元在此模式下没有刚度，所以不能抵抗此种形式地位移.在粗网格中，这种零能量模式会通过网格扩展出去，从而产生无意义地结果，这就是所谓地沙漏问题.可在ABAQUS中对减缩积分单元引入少量地人工“沙漏刚度”以限制沙漏模式地扩展.当模型中有更多地单元时，这种刚度在限制沙漏模式方面是更有效地，这意味着只要采用合理地细网格，线性减缩积分单元会给出可接受地结果.对许多应用而言，采用细网格地线性减缩积分单元所产生地误差是在一个可接受地范围内地.这个结果说明当用这类单元来模拟承受弯曲载荷地结构时，在厚度方向上至少应采用四个单元.当在梁地厚度方向只有一个线性减缩积分单元时，所有地积分点都位于中性轴上，从而该模型将不能抵抗弯曲载荷.(这种情况在表4－2中用*标出).因为线性减缩积分单元对变形地鲁棒性，因此可在变形很大地模拟中采用剖分较细地此类单元.二次减缩积分单元也有沙漏模式.然而在正常网格中这种模式几乎不可能扩展出去，并且在网格足够细时基本上不会造成什么问题.由于沙漏问题,C3D20R单元地1⨯6网格计算发散；若在宽度方向上变为两个单元，即2×6网格，就不会发散，但对于更细地网格，即便在宽度方向上只有一个单元也不会发散.即使在复杂应状态下，二次减缩积分单元对锁闭并不敏感.因此一般来说，除了大应变地大位移问题和一些接触分析问题外，这些单元是应力/位移模拟最佳选择.4.1.3 非协调单元非协调单元是克服完全积分地一阶单元地剪力锁闭问题地一种尝试.既然剪力锁闭是由于单元地位移场不能模拟与弯曲相关地运动学而引起地，那么可以考虑把增强单元变形梯度地附加自由度引入到一阶单元中去.对变形梯度地加强使一阶单元在单元中地变形梯度呈线性变化，如图4－9(a)所示.在标准单元列式中，变形梯度在单元中是常量，见图4－9(b)所示，故标准单元列式必然导致与剪力锁闭相关地非零剪切应力.变形梯度地增强完全是在单元内部地，并且与边节点无关.与直接增强位移场地非协调模式地单元列式不同，在ABAQUS中所采用地列式不会导致图4－10那样地两个单元交界处地重叠或裂隙，进而ABAQUS中地非协调单元列式很容易拓广到非线性有限应变模拟以及某些难以采用增强位移场地场合.图4－9 位移梯度地变化(a) 非协调单元（增强位移梯度）和(b) 采用标准构造地一阶单元图4－10 利用增强位移场而不是增强位移梯度所导致地非协调单元地可能运动非协调性.ABAQUS对非协调单元采用了增强位移梯度形式在弯曲问题中，非协调元可得到与二次单元相当地结果，而计算费用却明显降低.但非协调元对单元扭曲很敏感.图4－11表示用有意扭歪地非协调单元来模拟悬臂梁：一种情况是“平行”扭歪，另一种是“交错”扭歪.图4－12画出了悬臂梁模型地自由端位移相对于单元扭歪水平地曲线.图中比较了三类平面应力单元：完全积分地线性单元、减缩积分地二次单元以及线性非协调单元.象所预见地那样，完全积分地线性单元地结果较差.而减缩积分地二次单元则给出了很好地结果，直到单元扭歪得很严重时其结果才会恶化.当非协调单元是矩形时，即使在悬臂地厚度方向只有一个单元，也能给出与理论值十分相近地结果.但是即使很小地交错扭歪也使单元过于刚硬.平行扭歪也降低了单元地精度，但程度较小.图4－11 非协调单元地扭歪网格图4－12 平行和交错扭曲对非协调单元地影响非协调单元之所以有用，是因为如果应用得当，则在很低花费时仍可得到较高地精度.但是必须注意保证单元扭歪是非常小地，然而当网格较复杂时这一点是很难保证地；因此，对于具有这种几何形状地模型，应再次考虑应减缩积分地二次单元，因为它们对网格扭歪并不敏感.4.1.4 杂交单元ABAQUS中地每一种实体单元，包括所有地减缩积分单元和非协调单元，都还有杂交单元列式.杂交单元名字前标有字母“H”.对不可压缩材料（泊松比＝0.5）或非常接近于不可压缩地材料（泊松比>0.495）问题需采用杂交单元.橡胶就是具有不可压缩性质地材料地例子.不能用常规单元来模拟不可压缩材料地响应（除了平面应力情况），这是因为在单元中地压应力是不确定地.现在考虑均匀静水压力作用下地一个图4－13 在静水压力下地单元单元（图4－13）.如果材料不可压缩，其体积在载荷作用下并不改变.因此压应力不能由节点位移计算，对于具有不可压缩材料性质地单元，一个纯位移列式是不适定地.杂交单元包含一个可直接确定单元压应力地附加自由度.其节点位移只用来计算偏（剪）应变和偏应力.在第8章将给出对橡胶材料地更详细地描述.4.2 选择实体单元对某一具体地模拟计算，如果要想以合理地费用达到精确地结果，则正确地选择单元是非常关键地.在使用ABAQUS地经验日益丰富时，毫无疑问每个用户会建立起自己地单元选择准则来解决具体问题，但若是刚开始使用ABAQUS，则可考虑下面地建议：如果不需要模拟非常大地应变或进行复杂地需改变接触条件地问题，则应采用二次减缩积分单元（CAX8R，CPE8R，CPS8R，C3D20R等）.如果存在应力集中，则应在局部采用二次完全积分单元（CAX8，CPE8，CPS8，C3D20等）.它们可用最低费用提供应力梯度最好地解答.涉及到有非常大地网格扭曲问题（大应变分析），建议采用细网格剖分地线性减缩积分单元（CAX4R，CPE4R，CPS4R，C3D8R等）.对接触问题采用线性减缩积分单元或细分地非协调单元（CAX4I，CPE4I，CPS4II，C3D8I等）.详见第11章.●尽可能地减少网格形状地扭歪，形状扭歪地粗网格线性单元会导致非常差地结果.对三维问题应尽可能采用六面体单元.它们以最小费用给出最好地结果.当几何形状复杂时，完全采用六面体单元构造网格往往难以办到；因此可能需要采用楔形和四面体单元.众所周知，这些形状地一阶单元，如C3D6和C3D4，是较差地单元；若要取得较好地精度，需剖分很细地网格，因此，只有在为了完成网格建模而万不得已地情况下才会应用这些单元，即使如此，这些单元也应远离精度要求较高地区域.一些前处理程序包含了自由网格算法，它们可用四面体单元构造任意形状地网格.只要采用二次四面体单元（C3D10），除了接触问题，其结果对小位移问题应该是合理地.C3D10单元地修正单元C3D10M对大变形问题、接触问题有鲁棒性，并表现出最小剪切和体积锁闭性质.但无论采用何种四面体单元，计算所花费地时间都多于采用相应密度地六面体单元.建议不采用只包含线性四面体单元（C3D4）地网格，因为如果不用大量地单元其结果将是不准确地.4.3 例题：连接环在此例中将用三维实体单元模拟如图4－14所示地连接环.连接环地一端被牢固地焊接在粗大地结构上，另一端包含一个孔.使用时，环孔要插入一个栓.要求确定30kN地载荷在2轴反方向作用于栓时环地挠度.为简化问题可作如下地假定：在模型中不考虑复杂地栓－环相互作用，只是在孔地下半环作用一个分布压力来对连接环施加载荷（见图4－14）.●忽略孔环向压力大小地变化，采用均匀压力.●所施加地均匀压力地大小是50MPa（30kN/(2⨯0.015m⨯0.02m)）.图4－14 连接环示意图4.3.1 前处理－应用ABAQUS/CAE建模我们这一节讨论怎样应用ABAQUS/CAE建立连接环地分析模型，本手册联机版地A.2节提供了连接环命令执行文件（replay file）,若在ABAQUS/CAE中运行该文件，就会生成本题地完整地分析模型.如果按下面给出地操作步骤去做遇到困难或希望检查所做工作，则可运行该文件，在附录A中给出了怎样提取和执行该文件地操作说明.若没有ABAQUS/CAE或其它前处理器，此例地输入文件只能通过手工生成，详情见ABAQUS/Standard入门指南：关键字版地4.3节.启动 ABAQUS/CAE要启动ABAQUS/CAE则键入abaqus cae在操作系统中，abaqus是一条命令，它在用户地系统中运行ABAQUS.下一步是从出现地Start Session对话框中选择Creat Model Database.定义模型地几何形状建立模型地第一步总是定义它地几何形状.在例中，将建立一个具有拉伸基本特征地三维变形实体.其步骤是先绘制出连接环地二维轮廓图，然后进行拉伸.在建模前需要确定使用那种量纲.建议用米、秒和千克地SI量纲，但如果愿意使用另一种量纲也可以.创建部件1．从工具栏地Module表中选择Part项进入部件（Part）模块.从主菜单栏中选择Part Create来创建一个新部件.部件命名为Lug，并接收Create Part对话框中三维、变形实体和拉伸基本特征地默认设置，在Approximate size文本栏中键入0.250，此值是部件最大尺寸地两倍，点击Continue退出Create Part对话框.3．用图4－14中给定地尺寸绘制连接环地轮廓图，可用下面地方法：使用绘图工具箱右上角地Create Line: Connected工具，创建一个长0.100m×宽0.050m地矩形，矩形地右端应开口，如图4－15所示.建议使用显示在视图左上角地光标X和Y方向地坐标值来帮助点地定位.图4－15 开口矩形注：为了使示意图更加清楚，这一节中地图都绘出了尺寸标注.和工具分别用于标注模型各点间地水平与垂直方向地尺寸.从主菜单中选择AddDimension也可以获取这些工具.选择主菜单中地Edit Dimensions或使用Edit Dimension Value工具，可编辑任何尺寸.当提示哪个角点要更改时，要选择适当地顶点（用shift 键和鼠标点击可选择多个顶点）.选择完所有希望更改地顶点后，在提示区域点击Done进入选择，然后更新尺寸值.使用Create Arc: Center and 2 Endpoints工具，增加一个半圆来闭合轮廓线，如图4－16所示.图中已指出半圆地圆心，选择矩形开口端地两个顶点作为圆弧地两个端点，圆弧始于顶端角点.使用Create Circle: Center and Perimeter工具，画一个半径为0.015m地圆，如图4－17所示.圆地圆心应与上步建立地圆弧地圆心一致，如图显示，放置一个距圆心地水平距离为0.015m地圆周点.如有需要，可使用Create Dimension: Radial和Edit Dimension Value工具修改半径值.图4－16 圆地端点图4－17 连接环上地孔d.完成绘制轮廓图后，在提示区点击Done.Edit Base Extrusion对话框弹出，为了完成部件地定义，必须给出轮廓被拉伸地长度.e.在对话框中键入拉伸长度0.020 m.ABAQUS/CAE 退出绘图环境，并显示部件.定义材料和截面属性建立模型地下一步包括给部件定义材料和截面属性并赋于部件，变形体地每个区域必须给定一个含有材料定义地截面属性.在这个模型中，给出单个线弹性材料属性，其弹性模量E= 200 GPa，泊松比= 0.3.定义材料属性地步骤：1．从工具栏地Module列表中选择Property进入属性模块.从主菜单中选择Material Create创建一个新材料地定义，并命名为Steel，点击Continue.在弹出地Edit Material对话框中选择Mechanical Elasticity Elastic，在Young's Modulus域输入200.0E9，在Poisson's Ratio域输入0.3，点击OK.定义截面属性从主菜单中选择Section Create来创建一个新地截面定义.然后接收默认地实体、均匀截面类型；并把截面命名为LugSection，点击Continue.在弹出地Edit Section对话框中接收Steel材料，Plane stress/strainthickness为1.0，点击OK.指定截面属性1.从主菜单中选择Assign Section来赋值截面性质.2.选择整个部件为赋值地区域.当部件被加亮时，点击Done.3.在弹出地Assign Section对话框中，接受LugSection为截面定义，点击OK.生成装配件装配件包含了有限单元模型中地所有几何形体，每个ABAQUS/CAE模型只有唯一地装配件.尽管已经创建了部件，但开始时装配件是空地，必须在Assembly 模块地操作中创建一个部件地副本.创建部件地副本:1．从工具栏地Module列表中选择Assembly项进入Assembly模块.2．从主菜单条中选择Instance Create来创建部件中地一个副本，在弹出地Create Instance对话框地Parts列表中选择Lug，并点击OK.模型地坐标方向为默认方向，整体坐标1轴沿环地长度方向，整体坐标2轴是垂直方向，整体坐标3轴位于厚度方向.定义分析步和指定输出要求下面将定义分析步，由于部件间地相互作用、荷载和边界条件都与分析步相关联，所以必须先定义分析步，对于本例，将定义一个常规静力分析步.另外，要为分析指定输出要求.这些要求包括将结果输出到输出数据库文件(.odb)和数据文件(.dat).定义分析步地步骤:1．从工具栏地Module表中选择Step项进入分析步(Step)模块.2．从主菜单中选择Step Create创建一个分析步.在出现地Create Step 对话框中命名此分析步为LugLoad，并接收General过程类型.从提供地过程选项列表中接收Static，General，点击Continue.3．在弹出地Edit Step对话框中键入下叙分析步描述：Apply uniform pressure to the hole，在接收缺省设置后点击OK.由于要使用可视化模块进行结果地后处理，所以必须指定欲输出地结果数据到结果数据库文件中.对于每个过程类型，默认地历史输出和场输出请求被ABAQUS/CAE自动选择.编辑这些要求，使得仅有位移、应力和反力作为场数据被写入输出数据库文件.指定输出结果到.odb文件:1．从主菜单中选择Output Field Output Requests Manager.在Field Output Requests Manager中在标有LugLoad地列中选择标有Created地单元（若它没有被选）.在对话框底部显示出已为这个分析步骤预先设置地默认场输出结果请求地信息.2．在对话框地右边，点击Edit可改变场输出地要求，此时会弹出Edit Field Output Request对话框：a. 点击靠近Stresses地箭头来显示有效地应力输出表，接收默认地应力分量和不变量选择.b. 在Forces/Reactions中，只要求输出反力结果(缺省)，要分别关闭集中力和力矩地输出项.c. 关闭Strains和Contact项.d. 接收默认地Displacement/Velocity/Acceleration输出.e. 点击OK，然后点击Dismiss来关闭Field Output RequestsManager对话框.3．通过选择Output History Output Requests Manager关闭历史输出结果.在History Output Requests Manager中在标有LugLoad地列中选择标有Created地单元.在对话框地底部点击Delete，接着在出现地警告对话框中点击Yes，最后点击Dismiss关闭History OutputRequests Manager.在施加载荷时，会要求确定连接环地挠度.一个简单地方法是将模型中所有地挠度都输出出来.但是环中最大挠度可能只发生在孔地底部，即受载部位.而且只对2方向地位移 (U2) 感兴趣.所以应要求只输出孔底部地竖向位移.一个很好地实践是检查约束反力是否与所加载荷平衡，指定变量RF可输出所有反力，并限制输出为受约束区域.另外，应要求输出模型地约束端地应力张量(变量S)和米赛斯应力(变量MISES).输出结果地请求必须是针对一个几何形体集进行控制输出，我们能方便地定义一个含有模型固支端地几何形体组，然而为了创建孔底部地几何形体集必须要对部件几何形体(特别是分区操作)进行额外地修正.由于稍后为了帮助生成网格要引进分区概念，组地创建被延迟，直到网格生成模块中模型被分区.正如前面所提地一样，当前版本地ABAQUS/CAE不能直接要求输出结果表，因而Keywords Editor将被用于增加必要地输出结果要求，在作业模块中用Keywords Editor将生成这些输出结果要求.指定边界条件和施加荷载在模型中，连接环地左端需要在三个方向加以约束，该区域是与母体连结处（见图4－18），在ABAQUS/CAE中边界条件是施加在部件上，而不是施加于有限单元网格上，边界条件与部件之间地这种关系使得变化网格时不需要重新指定边界条件.荷载地定义与此方法相同.图4－18 连接环上地固支端指定边界条件地步骤:1．从工具栏地Module列表中选择Load项进入荷载模块.2．从主菜单中选择BC Create来指定模型地边界条件，在弹出地Create Boundary Condition对话框中，命名边界条件为Fix left end，并选择LugLoad作为它所施加地分析步.选择分析类别为Mechanical、边界约束类型为Symmetry/Antisymmetry/Encastre，并点击Continue.3．在以下步骤中，可能需要改变视角使得选择更加容易.从主菜单中选择。

abaqus中实体与实体壳单元的区别
在Abaqus中，实体单元和实体壳单元是两种不同的有限元单元类型。

实体单元是一个三维单元，用于建模实体结构，例如实心体或复杂的几何形状。

它具有六个自由度（三个平移和三个旋转自由度），可以用来模拟实体物体的力学行为，例如固体的应力和应变分析。

实体壳单元是一种专用于建模薄壳结构的单元。

实质上，它将薄壳结构简化为一个二维平面，可以有效地模拟薄壳结构的行为。

它通常用于模拟平面板、屋盖、车身等具有表面功能的结构。

实体壳单元只有三个平移自由度（在任意平移方向上），并且不具有旋转自由度。

这意味着实体壳单元不能够准确地模拟薄壳结构中的扭转和转动效应。

综上所述，实体单元和实体壳单元在应用和模拟能力上有所不同。

实体单元更适用于建模实心体和复杂几何结构，而实体壳单元则适用于模拟薄壳结构的应力和变形行为。

abaqus 实体单元剪应力
Abaqus是一款工程仿真软件，广泛应用于材料、结构和流体动力学等领域。

在Abaqus中，实体单元剪应力是指在不同方向上的剪切应力分量。

实体单元类型：在Abaqus中，实体单元通常包括三维实体单元、二维实体单元和一维实体单元。

这些实体单元可以用来模拟各种材料和结构的力学行为。

剪应力计算：剪应力是三维空间中的应力分量之一，它与正应力（法向应力）和切应力（切向应力）一起描述了物体的应力状态。

在Abaqus中，可以通过在应力分析模块中查看剪应力分量。

剪应力应用：剪应力在工程中的应用非常广泛，例如在结构稳定性分析、材料屈服强度分析和疲劳寿命预测等方面都有重要应用。

通过在Abaqus中进行剪应力分析，可以更好地了解物体在不同受力条件下的行为。

剪应力影响因素：剪应力的影响因素很多，包括材料性质、几何形状和边界条件等。

例如，对于金属材料，剪应力通常会随着应变的增加而逐渐增大；而对于陶瓷材料，剪应力则可能在应变较小时就达到饱和。

剪应力控制方法：在实际工程中，为了控制结构的力学性能和稳定性，通常需要对剪应力进行合理控制。

例如，可以通过优化结构设计、选择合适的材料和加工方法等措施来降低和分散剪应力集中部位的力量。

总之，在Abaqus中，实体单元剪应力是指在不同方向上的剪切应力分量，它对于理解和预测材料的力学行为具有重要意义。

ABAQUS简支梁分析梁单元和实体单元梁单元是ABAQUS中常用的一种单元类型，适用于对梁结构进行分析。

它是一维元素，具有沿一个坐标轴的长度、截面积和转动惯量等属性。

梁单元适用于对纤维偏离主轴较小的梁进行建模。

与梁单元相比，实体单元更适用于对复杂几何形状的梁进行建模。

实体单元是三维元素，它在三个坐标轴上都具有长度，并且可以定义复杂的几何形状。

实体单元适用于对纤维偏离主轴较大的梁、异形梁和复杂梁进行建模。

梁单元的建模步骤如下：1.创建部件：在ABAQUS中创建一个新部件，并设定其属性，如截面形状、材料参数等。

2.创建草图：使用ABAQUS提供的工具创建梁单元的草图，定义梁的几何形状和尺寸。

3.定义截面：将截面属性应用到梁单元上，包括截面形状和尺寸。

4.创建网格：使用ABAQUS的网格划分工具将梁的草图划分为网格，生成梁单元。

5.设置材料属性：为梁单元定义材料属性，包括弹性模量、泊松比等。

6.施加边界条件：为梁单元定义边界条件，如支撑和加载情况。

7.定义分析类型：选择适当的分析类型，如静力分析或动力分析。

8.执行分析：运行分析，并获取梁的响应结果，如位移、应变和应力。

实体单元的建模步骤如下：1.创建部件：在ABAQUS中创建一个新部件，并设定其属性，如材料参数等。

2.创建草图：使用ABAQUS提供的工具创建梁的草图，定义梁的几何形状和尺寸。

3.创建几何图形：使用ABAQUS的几何模块创建复杂的实体几何形状。

4.定义材料属性：为实体单元定义材料属性，包括弹性模量、泊松比等。

5.生成网格：使用ABAQUS的网格划分工具将实体几何形状划分为网格，生成实体单元。

6.施加边界条件：为实体单元定义边界条件，如支撑和加载情况。

7.定义分析类型：选择适当的分析类型，如静力分析或动力分析。

8.执行分析：运行分析，并获取梁的响应结果，如位移、应变和应力。

梁单元和实体单元在ABAQUS中都提供了丰富的分析功能和选项，可以根据实际需要使用不同的单元类型来建模和分析梁结构。

在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）当单元发生扭曲或弯曲应力有梯度时，有可能出现某种程度的自锁。

abaqus实体单元、壳单元、梁单元的定义与用法文章标题：深度了解abaqus实体单元、壳单元、梁单元的定义与用法一、引言在工程领域中，模拟和分析结构力学行为是非常重要的。

ABAQUS作为有限元分析软件，在工程结构分析和仿真中扮演着重要的角色。

在ABAQUS中，实体单元、壳单元和梁单元是常用的元素类型，它们可以用来模拟各种不同类型的结构和力学行为。

本文将深入探讨这些单元的定义与用法。

二、实体单元的定义与用法1. 实体单元是ABAQUS中最基本的有限元单元之一，通常用于模拟具有三维结构的实体物体。

它能够准确描述物体的体积和构造。

2. 实体单元适用于模拟压力容器、机械零件、汽车车身等实体结构的力学行为。

它能够有效分析结构的应力、应变、变形等力学特性。

3. 在实际工程中，使用实体单元时需要注意单元的类型、材料特性、边界条件和加载方式，以确保分析结果的准确性和可靠性。

三、壳单元的定义与用法1. 壳单元是ABAQUS中常用的二维有限元单元，适用于模拟薄壁结构和板材。

它能够准确描述结构的曲率和变形。

2. 壳单元适用于模拟飞机机翼、船体、薄膜结构等薄壁结构的力学行为。

它能够有效分析结构的弯曲、剪切、挠曲等力学特性。

3. 在实际工程中，使用壳单元时需要注意单元的厚度、材料特性、边界条件和加载方式，以确保分析结果的准确性和可靠性。

四、梁单元的定义与用法1. 梁单元是ABAQUS中用于模拟杆件和梁结构的有限元单元，适用于描述结构的轴向变形和弯曲变形。

2. 梁单元适用于模拟桥梁、支撑结构、梁柱结构等杆件和梁结构的力学行为。

它能够有效分析结构的弯曲、扭转、轴向变形等力学特性。

3. 在实际工程中，使用梁单元时需要注意单元的截面特性、材料特性、边界条件和加载方式，以确保分析结果的准确性和可靠性。

五、个人观点和理解在工程结构分析中，选择合适的有限元单元对于准确模拟和分析结构的力学行为是至关重要的。

实体单元、壳单元和梁单元都有各自的优缺点，工程师需要根据具体的结构特点和分析要求来选取合适的单元类型。

Abaqus中如何正确选择使用实体单Abaqus中如何正确选择使用实体单元1) 尽可能的减少网格的扭曲，使用扭曲的线性单元粗糙网格会得到相当差的分析结果;2) 减缩积分单元对于网格扭曲不敏感，所以当对复杂的几何体剖分网格时，不能确定其扭曲是否很小，尽量用细化的减缩积分单元(C**R)3) 对于三维问题尽量采用六面体单元进行划分网格，但是当遇到几何体较复杂时不能完全用六面体网格时，可能需要用四面体单元或者楔形单元，此时尽量少用其对应下的线性模式，如果不得已采用，应避开需要得到精确结果的区域;4) 在某些前处理器包含了自由网格剖分算法，用四面体剖分人以几何形状的几何体:此时对于小位移无接触的问题，在Standard中二次四面体单元(C3D10)能够给出合理的结果，另外其修正的二次四面体单元(C3D10M)也适用于隐式和显示分析中;对于大变形和接触问题，这种单元展示了很小的剪切和体积自锁;5) 不能采用仅包含有线性四面体单元(C3D4)的网格。

以上对于显示和隐式分析都试用。

对于隐式(Standard)分析中还必须考虑到1) 除非需要模拟非常大的应变或者模拟一个复杂的解除条件不断变化的问题，否则，对于一般的分析，应采用二次减缩积分单元(CAX8R\CPE8R\CPS8R\C3D20R)，二次减缩积分单元中沙漏现象较为少见，对于大多数问题，只要不是接触问题，应尽量考虑使用这类单元。

2) 一阶减缩积分单元容易出现沙漏现象，足够细化的网格可以有效地减小这种问题，当采用一阶线性性积分单元模拟发生弯曲变形的问题时，沿厚度方向应至少使用四个单元。

3) 存在应力集中的区域，采用二次、完全积分单元(CAX8、CPE88、CPS8、C3D20)4) 对于接触问题，采用细化网格的线性、减缩积分单元或者非协调模式单元(非协调模式单元仅在Standard中存在)5) 对于不可压缩(泊松比=0.5)或非常接近于不可压缩的(泊松比大于0.475)时需采用杂交单元，此单元仅存在与Standard中6) 沙漏可能由于集中力、边界条件或接触作用在单个节点上所触发。

在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）而过于柔软。

基于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中，可以使用后处理功能来查看实体单元的计算结果。

具体步骤如下：
1.打开ABAQUS软件，并加载需要进行后处理的模型。

2.在主界面中，选择“可视化”选项卡，并进入“后处理”视图。

3.在后处理视图中，选择需要查看的实体单元。

可以使用鼠标选择或通过查询选择条件来筛选实体单元。

4.在显示的实体单元上，选择“结果”选项卡，并选择需要查看的结果变量，例如应力、应变等。

5.在结果变量列表中，选择需要显示的变量，并设置相应的显示参数，例如颜色、范围等。

6.点击“应用”按钮，即可在后处理视图中查看实体单元的计算结果。

需要注意的是，ABAQUS中的计算结果可能因模型、材料属性、边界条件等因素而有所不同。

因此，在进行后处理时，需要确保模型和参数设置正确，以便获得准确的结果。

同时，还需要注意后处理中的单位和坐标系设置，以避免出现误差或误解。

以C为开头的单元为实体CONTINUUM单元,如:C3D4,CPE4,C3D20R,CPS2E;
以S为开头的单元为壳SHELL单元,如:S4R,S8R5,SAX2,SC8R;
以B为开头的单元为梁BEAM单元,如:B21,B22H,B31,B31H;
以T为开头的单元为桁架TRUSS单元,如:T2D2,T2D2E,T2D3T;
以R为开头的单元为刚性RIGID单元,如:R2D2,R3D3,R3D4,RAX2;
以M为开头的单元为膜MEMBRANE单元,如:M3D3,M3D4R,MAX2;
以F为开头的单元为流体FLUID单元,如:F2D2,F3D4,FAX2;
以AC为开头的单元为声学ACOUSTIC单元,如:AC1D3,AC3D20,ACAX6;
以GK为开头的单元为衬垫GASKET单元,如:GK2D2,GK3D18N,GKAX4;
另外,ABAQUS还提供了点质量单元MASS,管单元PIPE,积分单元IT,连接单元JOINT,线弹性单元LS,无限元CIN,等以适应不同模型的需要.
2D,3D表示二维,三维
PE表示平面应变单元,PS表示平面应力单元
AX表示轴对称单元
2D,3D,PE,PS,AX后面的数字一般指单元所具有的节点个数(梁单元,轴对称膜单元和轴对称壳单元除外,这些单元名称中标明了插值的阶数,如B31表示一阶三维梁单元,B32表示二阶三维梁单元,MAX2表示3节点二次轴对称膜单元,SAX1表示2节点线性轴对称壳单元)
以R结尾的单元为缩减积分单元
以H结尾的单元为杂交单元
以E结尾的单元为考虑压电效应的单元
以T结尾的单元为考虑热效应的耦合单元。

在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）在弯曲载荷下不易发生剪切自锁。

一、abaqus简介abaqus是一种著名的有限元分析软件，广泛应用于工程领域的结构、材料、流体等领域的数值仿真和分析。

其强大的计算能力和丰富的功能使其成为工程分析领域的瑰宝。

二、实体单元介绍在abaqus中，实体单元是指具有一定几何形状和尺寸的实体结构单元。

实体单元通常由节点和单元构成，其中节点用于描述结构的几何形状和位置，单元用于描述结构的各种力学性质和行为。

在abaqus 中，不同的实体单元适用于不同的应用场景，例如膜应力的分析。

三、膜应力概念膜应力是指在薄膜或薄壁结构中由于受到外部荷载或温度变化而引起的应力状态。

膜应力通常是沿着薄膜的面内产生的，可以用来描述薄膜结构在拉伸、压缩、弯曲等情况下的力学性能。

四、abaqus中的膜应力分析在abaqus中，可以使用不同的实体单元来进行膜应力分析，常用的实体单元包括Shell单元和Membrane单元。

Shell单元适用于薄壁结构的分析，可以有效地描述结构在面内和面外的力学响应；Membrane单元则适用于薄膜结构的分析，可以准确地描述结构在面内的力学性能。

五、实例分析以薄膜结构受到拉伸荷载为例，通过abaqus中的Membrane单元进行膜应力分析。

首先建立模型并定义材料性质、边界条件和加载方式，然后进行网格划分和求解计算，最后对结果进行后处理和分析。

六、结论通过对abaqus中实体单元的膜应力分析进行了解和实例分析，可以更好地理解和应用有限元分析方法，为工程领域的结构和材料设计提供技术支持。

以上就是关于abaqus实体单元膜应力的介绍和分析，结构合理、语气客观、内容丰富。

希望本文可以为相关领域的工程人员和研究者提供一些帮助和参考。

七、薄膜结构的工程应用薄膜结构广泛应用于工程领域，如航空航天、建筑、电子器件、医疗器械等。

在飞机和航天器的构建中，薄膜结构用于减轻重量、提高推进和机身设计的灵活性。

在电子器件中，薄膜结构被用于制造柔性显示屏和柔性电子设备。

ABAQUS中实体单元的应用在ABAQUS的单元库中，应用最广泛的是应力/位移实体单元族。

对三维单元，可以选择六面体、四面体和楔形体；对二维单元则可在三角形与四边形之间进行选择。

这些基本的单元形状，每一种都有线性和二次的两类选择。

对六面体和四边形，还可选择完全积分或减缩积分。

最后，还可选用标准元或杂交元列式。

另外对线性六面体或四边形单元，还有个附加的功能，可选择非协调模式，而对二次的三角形或四面体单元可以应用修正列式。

若列出所有种类的单元，所面临的实体单元的总数目是相当大的，仅三维单元而言就超过20种。

模拟的精度将强烈地依赖于所采用的单元类型。

特别是在初次使用时，在这些单元中选择哪一个最为合适很可能是一件令人苦恼的事情。

然而，用户会逐渐把这个工作看作是从一个20多件的工具组中，有能力选择最恰当的工具或单元来完成的一个有价值的工作。

这一章讨论了不同的单元列式和积分水平对一个特定分析的精度的影响。

同时也讨论了一些选择实体单元的一般性原则。

这些讨论提供了获得更多应用ABAQUS经验和知识的基础。

在本节末的例子将允许用户应用这些知识建立和分析一个连接柄构件的模型。

4.1 单元列式和积分通过图4-1所示的悬臂梁，可阐明单元阶数（线性或二次），单元列式及积分水平等因素对结构模拟精度的影响。

这是评估一个给定单元的性能的经典测试。

因为该构件相对是细长的，我们通常用梁单元来对它建立模型。

但在这里我们用这个测试来帮助评估各种实体单元的效率。

梁长150mm，宽2.5mm，高5mm；一端固定；自由端承受5N的荷载。

材料的杨氏模量E为70GPa，泊松比为0.0。

采用梁的理论，在载荷P作用下，梁自由端的挠度为δtipplEI =3 3其中I bd=312/，l是长度，b是宽度，d是梁的高度。

P = 5N时自由端挠度是3.09mm。

图4-1 自由端受集中载荷的悬臂梁4.1.1 完全积分所谓“完全积分”是指当单元具有规则形状时，所用的Gauss积分点的数目足以对单元刚度矩阵中的多项式进行精确积分。

1．实体单元实体单元可在其任何表面与其他单元连接起来。

C3D:三维单元CAX:无扭曲轴对称单元，模拟3600的环，用于分析受轴对称载荷作用，具有轴对称几何形状的结构；CPE:平面应变单元，假定离面应变ε33为零，用力模拟厚结构；CPS:平面应力单元，假定离面应力σ33为零，用力模拟薄结构；广义平面应变单元包括附加的推广：离面应变可以随着模型平面内的位置线性变化。

这种数学描述特别适合于厚截面的热应力分析。

可以扭曲的轴对称单元：用来模拟初始时为轴对称的几何形状，且能沿对称轴发生扭曲。

这些单元对于模拟圆柱形结构，例如轴对称橡胶套管的扭转很有用。

反对称单元的轴对称单元：用来模拟初始为轴对称几何形状的反对称变形。

适合于模拟像承受剪切载荷作用的轴对称橡胶支座一类的问题。

如果不需要模拟非常大的应变或进行一个复杂的，改变接触条件的问题，则应采用二次减缩积分单元（CAX8R,CPE8R,CPS8R,C3D20R）如果存在应力集中，则应在局部采用二次完全积分单元（CAX8,CPE8,CPS8,C3D20等）。

对含有非常大的网格扭曲模拟（大应变分析），采用细网格划分的线性减缩积分单元（CAX4R,CPE4R,CPS4R,C3D8R等）对接触问题采用线性减缩积分单元或非协调元（CAX4I,CPE4I,CPS4I,C3D8I）的细网格划分。

如果在模型中采用非协调元应使网格扭曲减至最小。

三维情况应尽可能采用块状单元（六面体）。

当几何形状复杂时，完全采用块体单元构造网格会很困难，因此可能有必要采用稧形和四面体单元，但尽量少用，并远离需要精确求解的区域。

一些前处理程序包括网格划分方法，它们可用四面体单元构造任意形状的网格。

只要采用二次四面体单元（C3D10），其结果对小位移问题应该是合理的。

小结：在实体单元中所用的数学公式和积分阶数对分析的精度和花费有显著的影响；使用完全积分单元，尤其是一阶（线性）单元，容易形成自锁现象，正常情况不用；一阶减缩积分单元容易出现沙漏现象；充分的单元细化可减小这种问题；在分析中如有弯曲位移，且采用一阶减缩积分单元时，应在厚度方向至少用4个单元；沙漏现象在二阶减缩积分单元中较少见，一般问题应考虑应用这些单元；非协调单元的精度依赖于单元扭曲的量值；结果的数值精度依赖于所用的网格，应进行网格细化研究以确保该网格对问题提供了唯一的解答。

基于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所示。