![ANSYS中单元类型与形函数及计算精度问题[1]](https://img.doczj.com/imgf1/078sm3wp3p79ti7wznm0-11.webp)
![ANSYS中单元类型与形函数及计算精度问题[1]](https://img.doczj.com/imgf1/078sm3wp3p79ti7wznm0-62.webp)
ANSYS中单元类型与形函数及计算精度问题有限元计算中，单元类型的选择，最终归结为形函数（插值函数）的选择直接从理论上影响着分析的精度。
在ANSYS中，单元类型与形函数之间存在什么关系呢？通过参考其理论手册发现，单元类型与形函数之间并不是一一对应关系，在某些情况下，一种单元类型在采用不同离散形式（例如体单元分别采用四面体，五面体和六面体）时对应不同的形函数（插值函数）。
以下针对常用的几种体单元对单元类型与形函数之间的关系作一简单说明，便于以后在分析问题时根据工程需要采用不同的单元类型及形函数。
1．SOLID45
对于SOLID45单元采用六面体、五面体和四面体形状时与形函数的对应关系为：
采用六面体时：
五面体时：
四面体时：
从六面体到四面体，虽然都采用的是SOLID45，但其对应的形函数阶次依次降低，因而在同等计算规模下计算精度也依次降低。当采用四面体离散时，相应形函数对应的单元退化为常应变四面体单元，要准确模拟则需要相当数量的单元方可。
同理对于SOLID95在采用不同形状离散时也分别对应不同的形函数。当退化为四面体时，SOLID95和SOLID92（只具有四面体形式）采用相同的形函数（可参考理论手册）：
因而从理论上来看，在采用相同数目的四面体单元进行离散时其计算结果应该是一致的。以下为在ANSYS中对其进行的验证：
从图1、2中可以明显看出，采用相同形函数而不同单元类型离散后，计算结果是相同的。
图1 使用SOLID92计算结果
以上仅为本人一点见解，欢迎一起探讨。Swifthorse1979@https://www.doczj.com/doc/f81751845.html,
1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中； 3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元； 4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性； 3. 关于下降段的问题 1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。 2）下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用rush模型。 3）利用最大压应变准则来判断混凝土是否破坏。 4. Solid65单元中的破坏准则 1）采用Willam&Warnke五参数破坏准则 2）需要参数： 单轴抗拉强度，单轴，双轴抗压强度，围压压力，在围压作用下双轴，单轴抗压强度
ANSYS中单元类型介绍和单元的选择原则 ANSYS中单元类型的选择 初学ANSYS的人，通常会被ANSYS所提供的众多纷繁复杂的单元类型弄花了眼，如何选择正确的单元类型，也是新手学习时很头疼的问题。 类型的选择，跟你要解决的问题本身密切相关。在选择单元类型前，首先你要对问题本身有非常明确的认识，然后，对于每一种单元类型，每个节点有多少个自由度，它包含哪些特性，能够在哪些条件下使用，在ANSYS的帮助文档中都有非常详细的描述，要结合自己的问题，对照帮助文档里面的单元描述来选择恰当的单元类型。 1.该选杆单元（Link）还是梁单元(Beam)？ 这个比较容易理解。杆单元只能承受沿着杆件方向的拉力或者压力，杆单元不能承受弯矩，这是杆单元的基本特点。 梁单元则既可以承受拉，压，还可以承受弯矩。如果你的结构中要承受弯矩，肯定不能选杆单元。 对于梁单元，常用的有beam3,beam4,beam188这三种，他们的区别在于： 1)、beam3是2D的梁单元，只能解决2维的问题。 2)、beam4是3D的梁单元，可以解决3维的空间梁问题。 3)、beam188是3D梁单元，可以根据需要自定义梁的截面形状。（常规是6个自由度，比如是用于桁架等框架结构，如鸟巢，飞机场的架构） 2.对于薄壁结构，是选实体单元还是壳单元？ 对于薄壁结构，最好是选用shell单元，shell单元可以减少计算量，如果你非要用实体单元，也是可以的，但是这样计算量就大大增加了。而且，如果选实体单元，薄壁结构承受弯矩的时候，如果在厚度方向的单元层数太少，有时候计算结果误差比较大，反而不如shell单元计算准确。 实际工程中常用的shell单元有shell63,shell93。shell63是四节点的shell单元(可以退化为三角形)，shell93是带中间节点的四边形shell单元(可以退化为三角形),shell93单元由于带有中间节点，计算精度比shell63更高，
Ansys的热载荷及热单元类型 Ansys的6种热载荷 ANSYS共提供了6种载荷，可以施加在实体模型或单元模型上，包括：温度、热流率、对流、热流密度、生热率和热辐射率。 1. 温度 作为第一类边界条件，温度可以施加在有限元模型的节点上，也可以施加在实体模型的关键点、线段及面上。 2. 热流率 热流率（Heal Flow)—种节点集中载荷，只能施加在节点或关键点上，主要用于线单元模型。提示：如果温度与热流率同时施加在某一节点上，則ANSYS读取温度值进行计算。 3.对流 对流（Convection)是一种面载荷，用于计算流体与实体的热交换。它可以施加在有限元模型的节点及单元上，也可以施加在实体模型的线段和面上。 4.热流密度 热流密度，又称热通量（Heat Flux)，单位为W/m2。热流密度是一种面载荷，表示通过单位面积的热流率。当通过单位面积的热流率己知时，可在模型相应的外表面施加热流密度。若输入值为正，则表示热流流入单元：反之，则表示热流流出单元。它可以施加在有限元模型的节点及单元上，也可以施加在实体模型的线段和面上。 提示：热流密度与对流可以施加在同一外表面，但ANSYS将读取最后施加的面载荷进行计算。 5. 生热率 如前所述，生热率既可看成是材料的一种基本属性，又可作为载荷施加在单元上，它可以施加在有限元模型的节点及单元上，也可以施加在实体模型的关键点、线段、面及体上。 6. 热辐射率 热辐射率也是一种面载荷，通常施加于实体的外表面。它可以施加在有限元模型的节点及单元上，也可以施加在实体模型的线段和面上。
Ansys的热单元类型 ANSYS 10.0热分析共提供了 40余种单元，其中包括辐射单元、对流单元、特殊单元以及前面所介绍的耦合场中-元等。其中常见的用于热分析的单元有16种： 下面一次对各单元进行介绍●MASS71 维度：1D、2D、3D 节点数：1 自由度：温度 性质：质量单元 几何形状 ●LINK31 维度：2D、3D 节点数：2 自由度：温度 性质：热辐射单元 几何形状
初学ANSYS的人，通常会被ANSYS所提供的众多纷繁复杂的单元类型弄花了眼，如何选择正确的单元类型，也是新手学习时很头疼的问题。 单元类型的选择，跟你要解决的问题本身密切相关。在选择单元类型前，首先你要对问题本身有非常明确的认识，然后，对于每一种单元类型，每个节点有多少个自由度，它包含哪些特性，能够在哪些条件下使用，在ANSYS的帮助文档中都有非常详细的描述，要结合自己的问题，对照帮助文档里面的单元描述来选择恰当的单元类型。 1.该选杆单元（Link）还是梁单元(Beam)？ 这个比较容易理解。杆单元只能承受沿着杆件方向的拉力或者压力，杆单元不能承受弯矩，这是杆单元的基本特点。梁单元则既可以承受拉，压，还可以承受弯矩。如果你的结构中要承受弯矩，肯定不能选杆单元。对于梁单元，常用的有beam3,beam4,beam188这三种，他们的区别在于： 1)beam3是2D的梁单元，只能解决2维的问题。 2)beam4是3D的梁单元，可以解决3维的空间梁问题。 3)beam188是3D梁单元，可以根据需要自定义梁的截面形状。 2.对于薄壁结构，是选实体单元还是壳单元？ 对于薄壁结构，最好是选用shell单元，shell单元可以减少计算量，如果你非要用实体单元，也是可以的，但是这样计算量就大大增加了。而且，如果选实体单元，薄壁结构承受弯矩的时候，如果在厚度方向的单元层数太少，有时候计算结果误差比较大，反而不如shell单元计算准确。 实际工程中常用的shell单元有shell63,shell93。shell63是四节点的shell 单元(可以退化为三角形)，shell93是带中间节点的四边形shell单元(可以退化为三角形),shell93单元由于带有中间节点，计算精度比shell63更高，但是由于节点数目比shell63多，计算量会增大。对于一般的问题，选用shell63 就足够了。除了shell63,shell93之外，还有很多其他的shell单元，譬如shell91,shell131，shell163等等，这些单元有的是用于多层铺层材料的，有的是用于结构显示动力学分析的，一般新手很少涉及到。通常情况下，shell63单元就够用了。 3.实体单元的选择。 实体单元类型也比较多，实体单元也是实际工程中使用最多的单元类型。常用的实体单元类型有solid45, solid92,solid185,solid187这几种。其中把 solid45,solid185可以归为第一类，他们都是六面体单元，都可以退化为四面体和棱柱体，单元的主要功能基本相同,(SOLID185还可以用于不可压缩超弹性材料)。Solid92, solid187可以归为第二类，他们都是带中间节点的四面体单元，单元的主要功能基本相同。 实际选用单元类型的时候，到底是选择第一类还是选择第二类呢？也就是到底是选用六面体还是带中间节点的四面体呢？如果所分析的结构比较简单，可以很方便的全部划分为六面体单元，或者绝大部分是六面体，只含有少量四面体和棱柱体，此时，应该选用第一类单元，也就是选用六面体单元；如果所分析的结构比较复杂，难以划分出六面体，应该选用第二类单元，也就是带中间节点的四面体单元。
摘要 本文介绍ANSYS 模拟钢筋混凝土梁的过程，讨论了有限元模型的建立以及在 ANSYS 中的实现，给出了用分离式配筋方法对混凝土梁的分析的一般过程。并给出了详细的命令流过程。并在此基础上对混凝土梁进行了分析，讨论了在力的作用下混凝土梁的塑形变形和裂缝的发展过程。 关键词 Ansys 混凝土梁 分离式配筋 The analysis of mechanics of a reinforced concrete based on ANSYS Abstract This paper introduces ANSYS simulation of the reinforced concrete beam process, discusses the establishment of the finite element model and the realization, and gives the ANSYS reinforcement method with separate the analysis of concrete beams of the general process. And gives the detailed command flow process. Based on the analysis of concrete beams, and discussed the concrete beam under the action of forces of the body deformation and fracture process. Keywords Ansys concrete beams reinforced separated 1 引言 由于钢筋混凝上材料性质复杂，使其表现出明显的非线性行为[1]。长期以来采用线弹性理论的设计方法来研究钢筋混凝上结构的应力或内力，显然不太合理，尽管有此理论是基于人量试验数据上的经验公式，还是不能准确反映混凝上的力学性能，特别是受力复杂的重要结构，必须采用三维钢筋混凝上非线性有限元方法才能很好地掌握其力学性能。利用ANSYS 对钢筋混凝上结构弹塑性的仿真分析，可以对结构自开始受荷载直到破坏的全过程进行分析，获得不同阶段的受力性能。本文将以混凝土梁的弹塑性分析为例，介绍在Ansys 中分析材料非线性问题的具体实现方法。 2 问题介绍 如图所示的钢筋混凝土梁[2]，横截面尺寸为200400b h mm mm ?=?,梁的跨度为3.0L m =，支座宽度为250mm 采用C20混凝土，梁内受拉纵筋3φ20，架立筋采用2φ12， 箍筋采用φ6@150，钢筋保护层厚度为25mm 。如图一。 图一 对于梁中所采用的所有钢筋，弹性模量为5 2.110MPa ?，抗拉强度设计值210MPa ， 密度33 7.810/kg m ?，泊松比为0.3。
ANSYS软件中常用的单元类型 一、单元 （1）link(杆)系列： link1(2D)和link8(3D)用来模拟珩架，注意一根杆划一个单元。 link10用来模拟拉索，注意要加初应变，一根索可多分单元。 link180是link10的加强版，一般用来模拟拉索。 （2）beam(梁)系列： beam3(2D)和beam4(3D)是经典欧拉梁单元，用来模拟框架中的梁柱，画弯据图用etab 读入smisc数据然后用plls命令。注意：虽然一根梁只划一个单元在单元两端也能得到正确的弯矩图，但是要得到和结构力学书上的弯据图差不多的结果还需多分几段。该单元需要手工在实常数中输入Iyy和Izz，注意方向。 beam44适合模拟薄壁的钢结构构件或者变截面的构件，可用"/eshape,1"显示单元形状。 beam188和beam189号称超级梁单元，基于铁木辛科梁理论，有诸多优点：考虑剪切变形的影响，截面可设置多种材料，可用"/eshape,1"显示形状，截面惯性矩不用自己计算而只需输入截面特征，可以考虑扭转效应，可以变截面（8.0以后），可以方便地把两个单元连接处变成铰接（8.0以后，用ENDRELEASE命令）。缺点是：8.0版本之前beam188用的是一次形函数，其精度远低于beam4等单元，一根梁必须多分几个单元。8.0之后可设置“KEYOPT(3)=2”变成二次形函数，解决了这个问题。可见188单元已经很完善，建议使用。beam189与beam188的区别是有3个结点，8.0版之前比beam188精度高，但因此建模较麻烦，8.0版之后已无优势。 (3)shell(板壳)系列 shell41一般用来模拟膜。 shell63可针对一般的板壳，注意仅限弹性分析。它的塑性版本是shell43。加强版是shell181（注意18*系列单元都是ansys后开发的单元，考虑了以前单元的优点和缺陷，因而更完善），优点是：能实现shell41、shell63、shell43...的所有功能并比它们做的更好，偏置中点很方便（比如模拟梁板结构时常要把板中面望上偏置），可以分层，等等。 (4)solid(体)系列 土木中常用的就solid45、solid46、solid65、solid95等。 solid45就不用多说了，solid95是它的带中结点版本。
一、关于模型 钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。 就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。 二、关于本构关系 混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。 就ANSYS而言，其问题比较复杂些。 1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？ 采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。而工程上又常将二者等同，其原因是工程结构不容许有很大的塑性变形，且混凝土等材料的屈服点不够明确，但破坏点非常明确。 定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系，即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系，而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系，如tb,MKIN,则定义的屈服准则是Von Mises,流动法则、硬化法则也就确定了)。 2 定义tb,concr后可否定义其它的应力应变关系 当然是可以的，并且只有在定义tb,concr后，有些问题才好解决。例如可以定义tb,miso,输入混凝土的应力应变关系曲线(多折线实现)，这样也就将屈服准则、流动法则、硬化法则等确定了。 这里可能存在一点疑问，即ANSYS中的应力应变关系是拉压相等的，而混凝土材料显然不是这样的。是的，因为混凝土受拉段非常短，认为拉压相同影响很小，且由于定义的tb,concr 中确定了开裂强度，所以尽管定义的是一条大曲线，但应用于受拉部分的很小。 三、具体的系数及公式 1 定义tb,concr时候的两个系数如何确定？ 一般的参考书中，其值建议先取为0.3~0.5(江见鲸)，原话是“在没有更仔细的数据时，不妨先取0.3~0.5进行计算”，足见此0.3~0.5值的可用程度。根据我的经验和理由，建议此值取大些，即开裂的剪力传递系数取0.5,（定要>0.2）闭合的剪力传递系数取1.0。支持此说法的还有 现行铁路桥规的抗剪计算理论，以及原公路桥规的容许应力法的抗计剪计算。
a n s y s单元类型
在结构分析中，“结构”一般指结构分析的力学模型。 按几何特征和单元种类，结构可分为杆系结构、板 壳结构和实体结构。 杆系结构：其杆件特征是一个方向的尺度远大于其它两个方向的尺度，例如长度远大于截面高度和宽度的 梁。单元类型有杆、梁和管单元（一般称为线单元）板壳结构：是一个方向的尺度远小于其它两个方向尺度的结构，如平板结构和壳结构。单元为壳单元 实体结构：则是指三个方向的尺度约为同量级的结构，例如挡土墙、堤坝、基础等。单元为3D实体单元和2D 实体单元 杆系结构： ①当构件15>L/h≥4时，采用考虑剪切变形的梁单元。 ②当构件L/h≥15时, 采用不考虑剪切变形的梁单元。 ③BEAM18X系列可不必考虑的上限，但在使用时必须 达到一定程度的网格密度。 对于薄壁杆件结构，由于剪切变形影响很大，所以必 须考虑剪切变形的影响。 板壳结构： 当L/h<5～8时为厚板，应采用实体单元。 当5～8<L/h<80～100时为薄板，选2D体元或壳元 当L/h>80～100时,采用薄膜单元。 对于壳类结构，一般R/h≥20为薄壳结构，可选择薄 壳单元，否则选择中厚壳单元。 对于既非梁亦非板壳结构，可选择3D实体单元。 杆单元适用于模拟桁架、缆索、链杆、弹簧等构件。该类单元只承受杆轴向的拉压，不承受弯矩，节点只有 平动自由度。不同的单元具有弹性、塑性、蠕变、膨胀、 大转动、大挠度（也称大变形）、大应变（也称有限应变）、应力刚化（也称几何刚度、初始应力刚度等）等 功能 ⑴杆单元均为均质直杆，面积和长度不能为零（LINK11 无面积参数）。仅承受杆端荷载，温度沿杆元长线性变 化。杆元中的应力相同，可考虑初应变。 ⑵LINK10属非线性单元，需迭代求解。LINK11可作用线 荷载；仅有集中质量方式。 ⑶LINK180无实常数型初应变，但可输入初应力文件， 可考虑附加质量；大变形分析时，横截面面积可以是变 化的，即可为轴向伸长的函数或刚性的。 ⑷通常用LINK1和LINK8模拟桁架结构，如屋架、网架、 网壳、桁架桥、桅杆、塔架等结构，以及吊桥的吊杆、 拱桥的系杆等构件，必须注意线性静力分析时，结构
ANSYS中单元类型介绍和单元的选择原则ANSYS中单元类型的选择 初学ANSYS的人，通常会被ANSYS所提供的众多纷繁复杂的单元类型弄花了眼，如何选择正确的单元类型，也是新手学习时很头疼的问题。 类型的选择，跟你要解决的问题本身密切相关。在选择单元类型前，首先你要对问题本身有非常明确的认识，然后，对于每一种单元类型，每个节点有多少个自由度，它包含哪些特性，能够在哪些条件下使用，在ANSYS的帮助文档中都有非常详细的描述，要结合自己的问题，对照帮助文档里面的单元描述来选择恰当的单元类型。 1.该选杆单元（Link）还是梁单元(Beam)？ 这个比较容易理解。杆单元只能承受沿着杆件方向的拉力或者压力，杆单元不能承受弯矩，这是杆单元的基本特点。 梁单元则既可以承受拉，压，还可以承受弯矩。如果你的结构中要承受弯矩，肯定不能选杆单元。 对于梁单元，常用的有beam3,beam4,beam188这三种，他们的区别在于： 1)、beam3是2D的梁单元，只能解决2维的问题。 2)、beam4是3D的梁单元，可以解决3维的空间梁问题。 3)、beam188是3D梁单元，可以根据需要自定义梁的截面形状。（常规是6个自由度，比如是用于桁架等框架结构，如鸟巢，飞机场的架构） 2.对于薄壁结构，是选实体单元还是壳单元？ 对于薄壁结构，最好是选用shell单元，shell单元可以减少计算量，如果你非要用实体单元，也是可以的，但是这样计算量就大大增加了。而且，如果选实体单元，薄壁结构承受弯矩的时候，如果在厚度方向的单元层数太少，有时候计算结果误差比较大，反而不如shell单元计算准确。 实际工程中常用的shell单元有shell63,shell93。shell63是四节点的shell单元(可以退化为三角形)，shell93是带中间节点的四边形shell单元(可以退化为三角形),shell93单元由于带有中间节点，计算精度比shell63更高，但是由于节点数目比shell63多，计算量会增大。对于一般的问题，选用shell63就足够了。
ANSYS 钢筋混凝土建模 一、简介 钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。希望大家一起讨论、批评指正(wang.jian@https://www.doczj.com/doc/f81751845.html,)。 程序：ANSYS 单元：SOLID65、BEAM188 建模方式：分离 暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。 二、单元选择 以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。只有计算的开裂荷载与实验还算是比较接近，但这个手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？ 所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。建模时发现，只要充分、灵活地运用APDL的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。 暗支撑剪力墙数值模型 看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如
Ansys单元类型设置 一、单元类型选择概述: ANSYS的单元库提供了100多种单元类型，单元类型选择的工作就是将单元的选择范围缩小到少数几个单元上； 单元类型选择方法： 1.设定物理场过滤菜单，将单元全集缩小到该物理场涉及的单元； 二、单元类型选择方法（续一） 2.根据模型的几何形状选定单元的大类，如线性结构则只能用“Plane、Shell”这种单元去模拟； 3.根据模型结构的空间维数细化单元的类别，如确定为“Beam”单元大类之后，在对话框的右栏中，有2D和3D的单元分类，则根据结构的维数继续缩小单元类型选择的范围； 三、单元类型选择方法（续二） 4.确定单元的大类之后，又是也可以根据单元的阶次来细分单元的小类，如确定为“Solid-Quad”，此时有四种单元类型：Quad 4node 42 Quad 4node 183 Quad 8node 82 Quad 8node 183 前两组即为低阶单元，后两组为高阶单元； 四、单元类型选择方法（续三） 5.根据单元的形状细分单元的小类，如对三维实体，此时则可以根据单元形状是“六面体”还是“四面体”，确定单元类型为“Brick”还是“Tet”； 五、单元类型选择方法（续四） 6.根据分析问题的性质选择单元类型，如确定为2D的Beam单元后，此时有三种单元类型可供选择，如下：2D elastic 3 2Dplastic 23 2D tapered 54，根据分析问题是弹性还是塑性确定为“Beam3”或“Beam4”，若是变截面的非对称的问题则用“Beam54”。 六、单元类型选择方法（续五） 7.进行完前面的选择工作，单元类型就基本上已经定位在2-3种单元类型上了，接下来打开这几种单元的帮助手册，进行以下工作： 仔细阅读其单元描述，检查是否与分析问题的背景吻合、
ANSYS中混凝土的计算问题【精华】 最近做了点计算分析，结合各论坛关于这方面的讨论，就一些问题探讨如下，不当之处敬请指正。 一、关于模型 钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。 就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。 二、关于本构关系 混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。 就ANSYS而言，其问题比较复杂些。 1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？ 采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。而工程上又常将二者等同，其原因是工程结构不容许有很大的塑性变形，且混凝土等材料的屈服点不够明确，但破坏点非常明确。 定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系，即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系，而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系，如tb,MKIN,则定义的屈服准则是V on Mises,流动法则、硬化法则也就确定了)。
一、用钢筋混凝土简支梁的数值模拟为实例，对ANSYS的使用方法进行说明 钢筋混凝土简支梁，尺寸为长2000mm，宽150mm，高300mm。混凝土采用C30，钢筋全部采用HRB335，跨中集中荷载P作用于一刚性垫板上，垫板尺寸为长150mm，宽100mm。 建立分离式有限元模型，混凝土采用SOLID65单元，钢筋采用LINK8单元，不考虑钢筋和混凝土之间的粘结滑移。创建分离式模型时，将几何实体以钢筋位置切开，划分网格时将实体的边线定义为钢筋即可。加载点以均布荷载近似代替钢垫板，支座处则采用线约束和点约束相结合。单元尺寸以50mm左右为宜。 二、命令流 !钢筋混凝土简支梁数值分析 !分离式模型 FINISH /CLEAR /PREP7 !1.定义单元与材料属性 ET,1,SOLID65,,,,,,,1 ET,2,LINK8 MP,EX,1,13585 !混凝土材料的初始弹模以及泊松比 MP,PRXY,1,0.2 FC=14.3 !混凝土单轴抗压强度和单轴抗拉强度 FT=1.43 TB,CONCR,1 TBDA TA,,0.5,0.95,FT,-1 !定义混凝土材料及相关参数，关闭压碎 TB,MISO,1,,11 !定义混凝土应力应变曲线，用MISO模型 TBPT,,0.0002,FC*0.19 TBPT,,0.0004,FC*0.36 TBPT,,0.0006,FC*0.51 TBPT,,0.0008,FC*0.64 TBPT,,0.0010,FC*0.75 TBPT,,0.0012,FC*0.84 TBPT,,0.0014,FC*0.91 TBPT,,0.0016,FC*0.96 TBPT,,0.0018,FC*0.99 TBPT,,0.002,FC TBPT,,0.0033,FC*0.85 MP,EX,2,2.0E5 !钢筋材料的初始弹模以及泊松比 MP,PRXY,2,0.3 TB,BISO,2 TBDA TA,,300,0 !钢筋的应力应变关系，用BISO模型
在结构分析中，“结构”一般指结构分析的力学模型。按几何特征和单元种类，结构可分为杆系结构、板壳结构和实体结构。杆系结构：其杆件特征是一个方向的尺度远大于其它两个方向的尺度，例如长度远大于截面高度和宽度的梁。元类型有杆、梁和管单元（一般单称为线单元）。板壳结构：是一个方向的尺度远小于其它两个方向尺度的结构，如平板结构和壳结构。单元为壳单元。实体结构：则是指三个方向的尺度约为同量级的结构，例如挡土墙、堤坝、基础等。单元为3D实体单元和2D 实体单元。 杆系结构： ①当构件15>L/h≥4时，采用考虑剪切变形的梁单元。（h为杆系的高度） ②当构件L/h≥15时, 采用不考虑剪切变形的梁单元。 ③BEAM18X系列可不必考虑L/h的值，但在使用时必须达到一定程度的网格密度。对于薄壁杆件结构，由于剪切变形影响很大，所以必须考虑剪切变形的影响。 板壳结构： 当L/h<5～8时为厚板，应采用实体单元。（h为板壳的厚度）当5～8
对于壳类结构，一般R/h≥20为薄壳结构，可选择薄壳单元，否则选择中厚壳单元。 对于既非梁亦非板壳结构，可选择3D实体单元。 杆单元适用于模拟桁架、缆索、链杆、弹簧等构件。该类单元只承受杆轴向的拉压，不承受弯矩，节点只有平动自由度。不同的单元具有弹性、塑性、蠕变、膨胀、大转动、大挠度（也称大变形）、大应变（也称有限应变）、应刚化（也称几何刚度、初始应力刚度等）等功能 ⑴杆单元均为均质直杆，面积和长度不能为零（LINK11无面积参数）。仅承受杆端荷载，温度沿杆元长线性变化。杆元中的应力相同，可考虑初应变。 ⑵LINK10属非线性单元，需迭代求解。LINK11可作用线荷载；仅有集中质量方式。 ⑶LINK180无实常数型初应变，但可输入初应力文件，可考虑附加质量；大变形分析时，横截面面积可以是变化的，即可为轴向伸长的函数或刚性的。 ⑷通常用LINK1和LINK8模拟桁架结构，如屋架、网架、网壳、桁架桥、桅杆、塔架等结构，以及吊桥的吊杆、拱桥的系杆等构件，必须注意线性静力分析时，结构不能是几何可
【原创】钢筋混凝土分离式建模方法（含ANSYS命令流） 钢筋混凝土, 分离式, 建模, ANSYS, 命令钢筋混凝土, 分离式, 建模, ANSYS, 命令 一、简介 钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。希望大家一起讨论、批评指正(wang.jian@https://www.doczj.com/doc/f81751845.html,)。 程序：ANSYS 单元：SOLID65、BEAM188 建模方式：分离 暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。 二、单元选择 以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。只有计算的开裂荷载与实验还算是比较接近，但这个 手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？ 所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。建模时发现，只要充分、灵活地运用APDL 的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。 [center] 暗支撑剪力墙数值模型[/center] 看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如下图箭头方向），从而导致总刚矩阵小主元地出现影响计算精度，或者干脆形成瞬变体系导致计 算提前发散。 [center] LINK8+SOLID65的问题[/center] 如果采用梁单元模拟暗钢筋，就算包裹钢筋的混凝土破坏了，钢筋单元本身仍可对连接点提供一定的侧向刚度（其实钢筋本身就是有一定抗弯刚度的），保证计算进行下去。ANSYS中的梁单元比较多，建议选取beam188单元。beam188支持弹塑性分析、自定义截面。可以用内力计算结果按截面插值得出应力结果，这样，SOLID65+beam188不仅解决了SOLID65+beam188的小主元问题，而且可以方便地控 制钢筋单元的划分密度，也扩充了钢筋单元输出信息。 三、单元组合方式 将剪力墙中所有钢筋单元（包括暗柱、梁的纵、箍筋、暗支撑钢筋、暗支撑箍筋、暗分布筋）单独建模，为了能够与混凝土单元节点共享，将混凝土单元细化，单元高度设为暗柱箍筋间距与墙片分布筋间距的最大公约数。 钢筋与混凝土单元节点共享。不考虑粘接-滑移影响。其实由于混凝土单元已经细化过了，钢筋周围的混凝土由于钢筋作用而开裂之
一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可，不需要约束方程，否则需要耦合自由度与约事方程。例如： （1）杆与梁、壳、体单元有公共节点即可，不需要约束方程。 （2）梁与壳有公共节点怒可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。 （3）梁与体则要在相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。 （4）壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。 上面所述的不同单元之间的接连方法主要是用耦合自由度和约束方程来实现的，有一定的局限性，只适用于小位移，下面介绍一种支持大位移算法的方法，MPC法。 MPC即Multipoint Constraint,多点约束方程，其原理与前面所说的方程的技术几乎一致，将不连续、自由度不协调的单元网格连接起来，不需要连接边界上的节点完全一一对应。 MPC能够连接的模型一般有以下几种。 solid 模型-solid 模型 shell模型-shell模型 solid 模型-shell 模型 solid 模型-beam 模型 shell 模型-beam模型 在 ANSYS中，实现上述MPC技术有三种途径。 （1）通过MPC184单元定义模型的刚性或者二力杆连接关系。定义MPC184单元模型与定义杆的操作完全一致，而MPC单元的作用可以是刚性杆（三个自由度的连接关系）或者刚性梁（六个自由度的连接关系）。 （2）利用约束方程菜单路径Main Menu>preprocessor>Coupling/Ceqn>shell/solid Interface创建壳与实体模型之间的装配关系。 （3）利用ANSYS接触向导功能定义模型之间的装配关系。选择菜单路径Main
ansys单元类型种类统计 单元名称种类单元号 LINK (共12种) 1,8,10,11,31,32,33,34,68,160,167,180 PLANE (共20种)2,13,25,35,42,53,55,67,75,77,78,82,83,121,145,146,162,182,183,223 BEAM (共09种)3,4,23,24,44,54,161,188,189 SOLID (共30 种)5,45,46,62,64,65,69,70,87,90,92,95,96,97,98,117,122,123,127,128,147,148,164,168, 185,186,187,191,226,227 COMBIN (共05种)7,14,37,39,40 INFIN (共04种)9,47,110,111 CONTAC (共05种)12,26,48,49,52 PIPE (共06种)16,17,18,20,59,60 MASS (共03种)21,71,166 MATRIX (共02种)27,50 SHELL (共19种)28,41,43,51,57,61,63,91,93,99,131,132,143,150,157,163,181,208,209 FLUID (共14种)29,30,38,79,80,81,116,129,130,136,138,139,141,142 SOURC (共01种)36 HYPER (共06种)56,58,74,84,86,158 VISCO (共05种)88,89,106,107,108 CIRCU (共03种)94,124,125 TRANS (共02种)109,126 INTER (共05种)115,192,193,194,195 HF (共03种)118,119,120 ROM (共01种)144 SURF (共04种)151,152,153,154 COMBI (共01种)165 TARGE (共02种)169,170 CONTA (共06种)171,172,173,174,175,178 PRETS (共01种)179 MPC (共01种)184 MESH (共01种)20
ANSYS土工结构计算案例 ANSYS－CHINA广州办事处 2019年6月24日
目录 计算题目及计算要求说明 (1) 题目一 (4) 一、计算说明 (4) 二、计算所用ANSYS邓肯－张的E－B模型说明 (5) 三、计算有限元模型及计算结果 (6) 题目二 (7) 一、用三维有限元模型计算 (7) 二、用三维有限元模型计算 (8) 题目三 (10) 一、计算说明 (10) 二、计算有限元模型及边界条件 (10) 三、强夯地基固结计算 (10) 题目四 (17) 一、计算说明 (17) 二、计算几何模型和有限元模型 (17) 三、计算结果 (18) 1、计算边界条件 (18) 2、计算结果 (19) 3、结论 (20)
计算题目及计算要求说明 题目一：高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析 题目二：大圆筒结构、波浪与地基的相互作用分析（大圆筒作为重力式码头结构，波浪为水平动荷载，门吊为竖向动荷载，地基为三层以上地基包括自抛碎石垫层、粘土层、粉细砂层和岩层，粉细砂层可能在波浪动荷载作用下液化造成圆筒倾覆） 题目三：（冲击荷载下）强夯地基固结有限元分析（提供固结方程或固结方程处理方案，孔隙水压力消散计算方案、沉降计算方案及其他一些处理技巧） 题目四：在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析（由上至下为公路面层，垫层，挡墙，挡墙面板采用预制混凝土块0.6?0.6?0.6m3，混凝土后方为钩挂式土工格栅，边坡比较陡，边坡有一定排水特性）。 具体处理方案包括： 1、提供计算输入界面 2、计算模型或采用本构情况 3、前处理方案及网格划分技巧 4、特殊材料或模型嵌入技术 5、计算技巧及解决方案 6、后处理提供内容
如果模型中加入（或删除）材料，模型中相应的单元就“存在”（或消亡）。单元生死选项就用于在这种情况下杀死或重新激活选择的单元。（可用的单元类型在表6－1中列出。）本选项主要用于钻孔（如开矿和挖通道等），建筑物施工过程（如桥的建筑过程），顺序组装（如分层的计算机芯片组装）和另外一些用户可以根据单元位置来方便的激活和不激活它们的一些应用中。单元生死功能只适用于ANSYS/Multiphysics，ANSYS/Mechanical和ANSYS/Structural产品。 Table 6-1 Elements with birth and death capability LINK1 SURF19 SHELL41 SOLID64 LANE83 SHELL143 PLANE2 IPE20 LANE42 SOLID65 SOLID87 SURF151 BEAM3 MASS21 SHELL43 LANE67 SOLID90 SURF152 BEAM4 SURF22 BEAM44 LINK68 SOLID92 SURF153 SOLID5 BEAM23 SOLID45 SOLID69 SHELL93 SURF154 LINK8 BEAM24 LANE53 SOLID70 SOLID95 SHELL157 LINK10 PLANE25 BEAM54 MASS71 SOLID96 TARGE169 LINK11 MA TRIX27 PLANE55 SOLID72 SOLID97 TARGE170 PLANE13 LINK31 SHELL57 SOLID73 SOLID98 CONTA171 COMBIN14 LINK32 PIPE59 PLANE75 SHELL99 CONTA172 PIPE16 LINK33 PIPE60 PLANE77 PLANE121 CONTA173 PIPE17 LINK34 SOLID62 PLANE78 SOLID122 CONTA174 PIPE18 PLANE35 SHELL63 PLANE82 SOLID123 在一些情况下，单元的生死状态可以根据ANSYS的计算数值决定，如温度，应力，应变等。可以用ETABLE命令（Main Menu>General Postproc>Element Table>Define Table）和ESEL命令（Utility Menu>Select>Entities）来确定选择的单元的相关数据，也可以改变单元的状态（溶和，固结，俘获等）。本过程对于由相变引起的模型效应（如焊接过程中原不生效的熔融材料变为生效的模型体的一部分），失效扩展和另外一些分析过程中的单元变化是有效的。 单元生死是如何工作的？ 要激活“单元死”的效果，ANSYS程序并不是将“杀死”的单元从模型中删除，而是将其刚度（或传导，或其他分析特性）矩阵乘以一个很小的因子[ESTIF]。因子缺省值为1.0E-6，可以赋为其他数值（详见“施加载荷并求解”一章）。死单元的单元载荷将为0，从而不对载荷向量生效（但仍然在单元载荷的列表中出现）。同样，死单元的质量，阻尼，比热和其他类似效果也设为0值。死单元的质量和能量将不包括在模型求解结果中。单元的应变在“杀死”的同时也将设为0。 与上面的过程相似，如果单元“出生”，并不是将其加到模型中，而是重新激活它们。用户必须在PREP7 中生成所有单元，包括后面要被激活的单元。在求解器中不能生成新的单元。要“加入”一个单元，先杀死它，然后在合适的载荷步中重新激活它。 当一个单元被重新激活时，其刚度，质量，单元载荷等将恢复其原始的数值。重新激活的单元没有应变记录（也无热量存储等）。但是，初应变以实参形式输入（如LINK1 单元）的不为单元生死选项所影响。而且，除非是打开了大变形选项[NLGEOM,ON],一些单元类型将以它们以前的几何特性恢复（大变形效果有时用来得到合理的结果）。单元在被激活后第一个求解过程中同样可以有热应变（等于a*(T-TREF)),如果其承受热量体载荷。