11_LS-DYNA碰撞解释
LS-DYNA软件对某半潜式海洋平台发生碰撞的分析

t he c o l l i s i on pr oc e s s,t he me c ha ni c a l c ha r a c t e r i s t i c s o f s hi p a nd o c e a n pl a t f or m we r e a na l y z e d.And t he n a c —
11_LS-DYNA碰撞解释

三节点仅拉伸单元,第三个节点定义单元初始方向,用于索绳建模MASS166 是一个有9个自由度的单点质量单元:在x,y,z方向的平动、速度、加速度这个单元还有附加的选项用来定义无质量的转动惯量:KEYOPT(1)=0 无惯量的3-D质量: 输入质量KEYOPT(1)=1 3-D转动惯量(无质量):输入6 个惯量值这个单元用来调整例如汽车碰撞这样复杂模型的质量,其中许多组件:(如座位,车灯,控制工具和假人等)未被建模(以质量单元替代)•使用两个节点和离散的材料模式来定义•能与其它所有显式单元连接•具有平动和转动自由度•能定义复杂的力-位移关系•而COMBIN14, 弹簧和阻尼必须是不同的单元..由于只能同时定义一个弹簧或阻尼选项,所以定义弹簧-阻尼集合体时需要重叠定义两个单元有两种实体单元算法:–单点积分实体(整个单元中常应力)•缺省形式•对于单元大变形单元非常快和有效•通常需要沙漏控制来阻止沙漏模式–全积分实体(2x2x2 积分)•比较慢,但无沙漏•对于高的泊松比时会同时出现剪切锁定和体积锁定,得到比较差的结果•精度比缺省算法对单元形状更敏感•在特定区域被选用来降低病态效应极力反对用退化的四面体网格////对显式动力学单元使用映射网格SHELL163 有12 种不同的单元算法,重要的包括:•Belytschko-Tsay ( BT, KEYOPT(1)=0 or 2, 缺省):–简单壳单元–非常快(相对速度= 1.0)–翘曲时易出错•Belytschko-Wong-Chiang ( BWC, KEYOPT(1)=10 ):–相对速度= 1.28 * BT–设用于翘曲分析–推荐使用•Belytschko-Leviathan ( BL , KEYOPT(1)=8 ):–相对速度= 1.25 * BT–较新,仍在开发中–第一个有物理沙漏控制的单元–(对于EDMP,HGLS,Mat,Val1无参数)•S/R co-rotational Hughes-Liu (S/R CHL, KEYOPT(1)=7):–没有沙漏控制的壳–相对速度= 8.84 * BT•单元算法BT, BWC, BL 仅适用于平面内单点积分,而S/R CHL 用于平面内4点积分。
LS-DYNA 理论及功能(简介)

LS-DYNA 理论及功能LS-DYNA 的理论及功能LS-DYNA 发展概况 (LS-DYNA Introduction)LS-DYNA是以显式为主、隐式为辅的通用非线性动力分析有限元程序,特别适合求解 各种二维、三维非线性结构的高速碰撞、爆炸和金属成形等非线性动力冲击问题,同时可以 求解传热、流体及流固耦合问题。
DYNA 程序系列最初是 1976 年在美国 Lawrence Livermore National Lab. 由 J.O.Hallquist 博士主持开发完成的,主要目的是为武器设计提供分析工具,后经 1979、1981、1982、1986、 1987、1988 年版的功能扩充和改进,成为国际著名的非线性动力分析软件,在武器结构设 计、内弹道和终点弹道、军用材料研制等方面得到了广泛的应用。
1988 年 J.O.Hallquist 创建 LSTC 公司,推出 LS-DYNA 程序系列,主要包括显式 LS-DYNA2D、LS-DYNA3D、隐式 LS-NIKE2D、LS-NIKE3D、热分析 LS-TOPAZ2D、 LS-TOPAZ3D、前后处理 LS-MAZE、LS-ORION、LS-INGRID、LS-TAURUS 等商用程序, 进一步规范和完善 DYNA 的研究成果,陆续推出 930 版(1993 年)、936 版(1994 年)、940 版(1997 年),950 版(1998 年)增加了汽车安全性分析(汽车碰撞、气囊、安全带、假人)、 薄板冲压成形过程模拟以及流体与固体耦合(ALE 和 Euler 算法)等新功能,使得 LS-DYNA 程序系统在国防和民用领域的应用范围进一步扩大,并建立了完备的质量保证体系。
1997 年LSTC公司将LS-DYNA2D、LS-DYNA3D、LS-TOPAZ2D、LS-TOPAZ3D等程序 合成一个软件包,称为LS-DYNA,PC版的前后处理采用ETA公司的FEMB,新开发的后处 理器为LS-POST。
LS_DYNA在汽车碰撞模拟过程中的应用_陈海树

沈阳大学 学报 JOURNAL OF SHENYANG U NIVERSIT Y
文章编号: 1008-9225( 2006) 04- 0010- 03
Vol118, No. 4 Aug. 2 0 0 6
LS-DYNA 在汽车碰撞模拟过程中的应用
陈海树, 赖征海, 邸建卫
18 0
10 000
10 000
20 0
A 立柱后移 量/ mm
35 50
B 立柱下端最 是汽车产品开发设计链中的常 规流程, 没有有限元分析的设计不能进入下一个 技术流程#使用 Dyna 软件在设计初 期对产品的 安全性能进行验证, 及时发现新产品的问题, 这为 设计工程师提供了更大的创造空间, 使设计质量 大幅度提高#
使用实际材 料的试验 结果值; 对于钢, 使用 24 号材料模式; 可恢复的泡沫材料 使用 57 号材 料模式, 若要考虑应变率的影响, 使用 83 号材料; 发动机等在碰撞过程中不变形物体采用 20 号刚 体材料#
收稿日期: 2006- 05- 12 作者简介: 陈海树( 1977- ) , 男, 辽 宁朝阳人, 沈阳华晨金杯汽车有限公司研发中心工程 师, 硕士#
第4期
陈海树等: LS- DYNA 在汽车碰撞模拟过程中的应用
11
2 乘员约束系统有限元模型
2. 1 假人有限元模型 仿真用的假人有限元模型是完全基于试验的
50 百分位的 Hybrid Ó型假人#模型细分为 98 个 部件#主要部件有: 头部、颈部、胸部、腹部、臀部和 四肢等共约 6 400 个节点, 3 900 个单元#假人各 部分质量和转动惯量都符合试验用假人#人体各 部分的关节用球形铰链( Spherical Joint ) 和旋转铰 链( Revolut e Joint ) 单元定义, 假人有限元模型铰 链总计 47 个#用非线 性的扭矩弹簧 Spring 和阻 尼 Damper 单元连接来模拟人体各部分之间的弹 性#胸部和颈部这些关键部位都设置成柔性, 能对 发生碰撞后假人的受伤进行评价#使 用 Oasys 软 件调整假人的 H 点及坐姿, 并保证假人与地板、 坐椅、防火墙等不发生穿透干涉# 2. 2 安全带模型
LS-Dyna碰撞分析资料

LS-DYNA 碰撞分析调试LS-DYNA碰撞计算模型的主要检查、调试项目有:a、质量增加百分比小于5%;b、总沙漏能小于5%;c、滑移界面能;d、检查各部件之间的连接、接触关系是否定义正确,检查模型的完整性;e、检查数值输出的稳定性。
一、质量缩放Mass scale的检查:质量缩放——对于时间步长小于控制卡片中设置的最小时间步长的单元,我们通常采取增加单元材料密度的方法来增大其时间步长,以减短模型的计算时间。
关于LS-DYNA中单元时间步长的计算方法请参见附录一。
1、初步检查。
让模型在dyna中运行2个时间步,在Hyper view中调出glstat 文件并检查mass scaling项(质量增加应该小于5%);调出matsum文件并检查各部件的质量增加情况,对于质量增加过大以及有快速增长趋势的部件应检查此部件的网格质量和材料参数设置(质量增加一般是由于单元的特征长度太小或者是材料参数E、ρ设置错误,导致该单元的时间步长低于控制卡片中设置的最小时间步长,从而引起质量缩放)。
2、全过程检查。
调整模型使其符合初步检查的标准,计算模型至其正常结束。
再按[初步检查]的要求检查调试整个模型直至达到要求。
一个计算收敛的模型在其整个计算过程中,最大质量缩放应小于总质量的5% 。
二、沙漏能Hourglass energy的检查:沙漏能的出现是因为模型中采用了缩减积分引起的,我们常用的B-T单元采用的是面内单点积分,这种算法会引起沙漏效应(零能模式)。
具体介绍参见附录二。
检查:在dyna中计算模型至其正常结束。
在Hyper view中调出glstat文件并检查energy的total energy 、Hourglass energy两项,整个计算过程中沙漏能应小于总能量的5% 。
三、滑移界面能sliding interface energy的检查:滑移界面能是由摩擦和阻尼所引起的。
剧烈的滑动摩擦会引起大的正值的滑移界面能;未能检测到的穿透(undetected penetrations)常常会引起大的负值的滑移截面能。
基于LS-DYNA 的移动式压力容器侧翻碰撞分析

52 2. 5 2. 5. 1 条件加载 仿真条件的设定
石家庄铁道大学学报( 自然科学版)
第 27 卷
本文仿真的移动式压力容器运输车侧翻试验特指其在正常运行状态下因过度转弯而导致的侧翻 , 运 输车从开始侧倾到触地有一个较长的运动过程 , 这不是研究的重点, 对仿真结果也没有影响, 故无需进行 仿真分析。从运输车侧翻触地的最后一刻开始模拟 , 直到运输车基本不再变形为止。只需算出运输车触 [9 ] 地一瞬间的角速度和直线速度, 作为仿真初始条件加载在模型即可 。 罐体内部已充满液氨, 在本仿真 液氨在模型中以点质量的形式均匀地分布在罐体内部 。整个模型的仿 中忽略其振荡对侧翻碰撞的影响 , , 。 真过程中 考虑重力的影响 2. 5. 2 侧翻角速度的计算 7] 本模型可通过能量守恒定律计算出触地时运输车的角速度 , 并将之加载。 文献[ 中提到, 运输车接 触地面时的动能为 E = 0 . 75 Mgh 也可表示为 E = 0 . 5 J ω2 ( 2) M 为运输车满载整车质量, 45 990 kg, 式中, 仿真假设液氨完全充满, 并忽略其对侧翻造成的影响; h 为运 890 mm; g 为重力加速度, 9 . 8 m / s2 ; J 为运输车 输开始侧倾时重心的高度和即将触地时的重心高度之差, 3 . 0 × 10 5 kg·m2 ; ω 为运输车绕转轴角速度。 绕转轴的转动惯量, ( 2 ) 式求解得 联立( 1 ) 、 ω = 计算得 2. 5. 3 ( 1)
槡
1 . 5 Mgh J
( 3)
ω = 1 . 42 rad / s
( 4)
其他关键字的设置 DYNA EXPORTE 模块输出 K 文件, 对模型定义速度 36 km / h, 和整体施加重力速度 g 后, 由 LS打开 K 文件设置内存空间、 合理的时间步长控制、 沙漏控制、 能量控制和输出控制等相关参数 , 以及为提高计算
LS-Dyna碰撞分析调试指南
LS-DYNA 碰撞分析调试LS-DYNA碰撞计算模型的主要检查、调试项目有:a、质量增加百分比小于5%;b、总沙漏能小于5%;c、滑移界面能;d、检查各部件之间的连接、接触关系是否定义正确,检查模型的完整性;e、检查数值输出的稳定性。
一、质量缩放Mass scale的检查:质量缩放——对于时间步长小于控制卡片中设置的最小时间步长的单元,我们通常采取增加单元材料密度的方法来增大其时间步长,以减短模型的计算时间。
关于LS-DYNA中单元时间步长的计算方法请参见附录一。
1、初步检查。
让模型在dyna中运行2个时间步,在Hyper view中调出glstat 文件并检查mass scaling项(质量增加应该小于5%);调出matsum文件并检查各部件的质量增加情况,对于质量增加过大以及有快速增长趋势的部件应检查此部件的网格质量和材料参数设置(质量增加一般是由于单元的特征长度太小或者是材料参数E、ρ设置错误,导致该单元的时间步长低于控制卡片中设置的最小时间步长,从而引起质量缩放)。
2、全过程检查。
调整模型使其符合初步检查的标准,计算模型至其正常结束。
再按[初步检查]的要求检查调试整个模型直至达到要求。
一个计算收敛的模型在其整个计算过程中,最大质量缩放应小于总质量的5% 。
二、沙漏能Hourglass energy的检查:沙漏能的出现是因为模型中采用了缩减积分引起的,我们常用的B-T单元采用的是面内单点积分,这种算法会引起沙漏效应(零能模式)。
具体介绍参见附录二。
检查:在dyna中计算模型至其正常结束。
在Hyper view中调出glstat文件并检查energy的total energy 、Hourglass energy两项,整个计算过程中沙漏能应小三、滑移界面能sliding interface energy的检查:滑移界面能是由摩擦和阻尼所引起的。
剧烈的滑动摩擦会引起大的正值的滑移界面能;未能检测到的穿透(undetected penetrations)常常会引起大的负值的滑移截面能。
汽车碰撞分析LS_DYNA控制卡片设置
控制卡片参数说明
*CONTROL_TIMESTEP(时间步长控制卡片) $ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST 0.0 0.9 2 0.0 -0.001 0 1 1 $ DT2MSF DT2MSLC 计算所需时间步长时,要检查所有的单元。出于稳定性原因,用0.9(缺省)来 减小时间步:Δt = 0.9 l/c ,特征长度l,和波的传播速度c,都与单元的类型有关。 DTINIT:初始时间步长,如为0.0,由DYNA自行决定初始步长; TSSFAC:时间步长缩放系数,用于确定新的时间步长。默认为0.9,当计算不稳定时,可以减小该值,但同时 增加计算时间; ISDO:计算4节点壳单元时间步长的(不同的值对应特征长度的不同算法,推荐使用2,因为此选项可以获得 最大的时间步长,但有三角形单元存在时会导致计算不稳定); TSLIMT:壳单元最小时间步分配 ,使单元的时间步长控制在最小时间步长之上;只适用于使用 *mat_plastic_kinematic,*mat_power_law_plasticity*mat_strain_rate_dependent_plasticity,*mat_piecewise_linear_pla sticity等材料模型的壳单元,不建议使用该选项,因为使用DT2MS选项更好。 DT2MS:因质量缩放计算得到的时间步长。当设置为一个负值时,初始时间将不会小于TSSFAC*|DT2MS|。质 量只是增加到时间步小于TSSAFC*|DT2MS|的单元上。当质量缩放可接受时,推荐用这种方法。用这种方法时 质量增加是有限的,过多的增加质量会导致计算终止。当设置为正值时,初始时间步长不会小于DT2MS。单 元质量会增件或者减小以保证每一个单元的时间步都一样。这种方法尽管不会因为过多增加质量而导致计算终 止,但更难以作出合理的解释。默认为0.0,不进行质量缩放; LCTM:限制最大时间步长的Load-curve,该曲线定义最大允许时间步长和时间的关系(可选择) ; ERODE:当计算时间步长小于TSMIN(最小时间步长)时体单元和t-shell被自动删除。
基于LSDYNA的汽车正碰分析
中图分类号:U467.1+4
0前言
在汽车被动安全性研究中,汽车碰撞是一个 十分复杂的力学问题,它是一个动态的大位移和 大变形的瞬态接触过程。接触和高速冲击载荷影 响碰撞的全过程,碰撞系统具有大位移、大转动 和大应变的非线性特征,以材料弹塑性变形为典 型特征的材料非线性和以接触摩擦为特征的接 触非线性的多重非线性特点。这些非线性物理陛 能的综合,使得利用传统的手段对汽车碰撞过程 的精确描述和求解变得非常复杂和困难。文中进 行的是汽车100%ODB正面碰撞的仿真分析。采 用动态显式非线性有限元技术,利用NX、ANSA、 LS.DYNA等有限元软件对某汽车进行几何建 模、网格建模和计算求解,依据计算结果分析了 汽车正面主要结构在碰撞中的耐撞性,并将计算 结果与试验结果进行对比,以验证仿真模型的准 确性。
1碰撞模型的建立
1.1汽车几何模型 建立汽车整车几何模型是进行有限元分析
的第一步,模型质量的好环对后续的有限元分析
·34·
文献标识码:A
起着至关重要的作用,直接决定有限元模型网格 的质量和分析时间的长短以及分析进展是否顺 利,尤其在碰撞分析中尤为突出。另外,在保证模 型正确及几何元素相关联的基础上,可以进行简 化,以达到事半功倍的效果。文中利用NX软件建 立某汽车整车的几何模型,该汽车三维几何模型 如图1所示。
[2]陈海树,赖征海,邸建卫.Ls—DYNA在汽车碰撞模拟过程 中的应用[J].沈阳大学学报,2006.
[3]何涛,杨竞.ANSYS/LS—DYNA非线性有限元分析实例指 导教程[M].机械工业出版社.
[4]张金换,杜汇良,马春生.汽车碰撞安全性设计[M].清华 大学出版社.
[5]Oaryl L.Logan.有限元方法基础教程[M].北京电子工业 出版社.
汽车碰撞精确分析LSDYNA控制卡片设置
THKCHG:在单面接触时考虑壳厚度的改变(默认时不考虑)。
ORIEN:在初始化时可选择性的对接触面部分自动再定位。
控制卡片参数说明
ENMASS:接触单元被腐蚀的质量处理。0-节点被移除,1-体单元节点被保留,2-体单元壳单元节点被保留。 USRSTR:每个接触面分配的存储空间,针对用户提供的接触控制子程序。 USRFRC:每个接触面分配的存储空间,针对用户提供的接触摩擦子程序。 NSBCS:接触搜寻的循环数(使用三维Bucket分类搜索),推荐使用默认项。 INTERM:间歇搜寻主面和从面接触次数。 XPENE:接触面穿透检查最大乘数,默认4.0。 SSTHK:在单面接触中是否使用真实壳单元厚度,默认0,不使用真实厚度。 ECDT:时间步长内忽略腐蚀接触。
16. DATABASE_BINARY_RUNRSF 设置如下:
控制卡片参数设置
17. DATABASE_BINARY_RUNRSF 设置如下:
控制卡片参数说明
*CONTROL_TERMINATION
$ ENDTIM ENDCYC DTMIN ENDENG ENDMAS
150
0
0.0
0.0
0.0
SLSFAC:滑动接触惩罚系数 ,默认为0.1。当发现穿透量过大时,可以调整该参数;
RWPNAL: 刚体作用于固定刚性墙时,刚性墙罚函数因子系数,为0.0时,不考虑刚体与刚性墙的作用,>0时, 刚体作用于固定的刚性墙,建议选择1.0;
ISLCHK:接触面初始穿透检查,为0或1(默认)时,不检查。为2时,检查。
后面将逐一介绍碰撞分析中经常用到的控制卡片,并对每个卡片的作 用进行说明。
控制卡片使用规则
卡片相应的使用规则如下:
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
三节点仅拉伸单元,第三个节点定义单元初始方向,用于索绳建模
MASS166 是一个有9个自由度的单点质量单元:在x,y,z方向的平动、速度、加速度
这个单元还有附加的选项用来定义无质量的转动惯量:
KEYOPT(1)=0 无惯量的3-D质量: 输入质量
KEYOPT(1)=1 3-D转动惯量(无质量):输入6 个惯量值
这个单元用来调整例如汽车碰撞这样复杂模型的质量,其中许多组件:(如座位,车灯,控制工具和假人等)未被建模(以质量单元替代)
•使用两个节点和离散的材料模式来定义
•能与其它所有显式单元连接
•具有平动和转动自由度
•能定义复杂的力-位移关系
•而COMBIN14, 弹簧和阻尼必须是不同的单元
..由于只能同时定义一个弹簧或阻尼选项,所以定义弹簧-阻尼集合体时需要重叠定义两个单元
有两种实体单元算法:
–单点积分实体(整个单元中常应力)
•缺省形式
•对于单元大变形单元非常快和有效
•通常需要沙漏控制来阻止沙漏模式
–全积分实体(2x2x2 积分)
•比较慢,但无沙漏
•对于高的泊松比时会同时出现剪切锁定和体积锁定,得到比较差的结果•精度比缺省算法对单元形状更敏感
•在特定区域被选用来降低病态效应
极力反对用退化的四面体网格////对显式动力学单元使用映射网格
SHELL163 有12 种不同的单元算法,重要的包括:
•Belytschko-Tsay ( BT, KEYOPT(1)=0 or 2, 缺省):
–简单壳单元
–非常快(相对速度= 1.0)
–翘曲时易出错
•Belytschko-Wong-Chiang ( BWC, KEYOPT(1)=10 ):
–相对速度= 1.28 * BT
–设用于翘曲分析
–推荐使用
•Belytschko-Leviathan ( BL , KEYOPT(1)=8 ):
–相对速度= 1.25 * BT
–较新,仍在开发中
–第一个有物理沙漏控制的单元
–(对于EDMP,HGLS,Mat,Val1无参数)
•S/R co-rotational Hughes-Liu (S/R CHL, KEYOPT(1)=7):
–没有沙漏控制的壳
–相对速度= 8.84 * BT
•单元算法BT, BWC, BL 仅适用于平面内单点积分,而S/R CHL 用于平面内4点积分。
•所有的壳单元沿着厚度方向有任意多数目的积分点(NIP) 。
对于弹性行为NIP = 2 (缺省)
对于塑性行为,3 < NIP < 5 (推荐NIP=5)
实常数用来定义积分点的数目
•R, NSET, SHRF,NIP,T1
–NSET = 实常数组参考号
–SHRF = 剪切因子(对于薄壳推荐为5/6)
–NIP = 积分点数
–T1 = 单元厚度
用EDINT命令定义结果输出的沿厚度方向的积分点数目
Solution > Output Controls > Integ Pt Storage …
–EDINT, SHELLIP, BEAMIP
•SHELLIP 是输出中壳的积分点数目|| SHELLIP 3
•每一个积分点与一个LAYER 相关|| 缺省值是 3 (顶层,中层和底层)
•BEAMIP是输出的梁积分点数目
对于三角形壳单元有两种算法:
C0 三角形壳( KEYOPT(1)=4 ),基于Mindlin-Reissner 平板理论。
此算法刚度偏大,不推荐用于整个壳体网格中。
BCIZ 三角形壳( KEYOPT(1)=3 ),基于Kirchhoff 平板理论,较慢,在混合网格中, C0 三角形单元通常比退化的4节点单元算法更好。
所以当混合划分(自由划分)通常使用下面的命令:•EDSHELL, , ,ITRST
–ITRST = 1 : 退化的四边形单元被当作三角形单元(缺省)
–ITRST = 2 : 退化的四边形单元保持不变
•Preprocessor > Shell Elem Ctrls > Triangular Shell Sorting > Full Sorting > OK
有两种膜单元算法:
单点积分的膜单元:Belytschko-Tsay-膜( KEYOPT(1)=5 );全积分Belytschko-Tsay-膜
(KEYOPT(1)=9):具有4点积分的膜单元
•全积分Belytschko-Tsay 壳( KEYOPT(1)=12 ): –不需要沙漏控制
–对横向剪切,假设的小应变弥补了剪切锁定
–平面内4点积分(2 X 2 积分) ,但速度仍然很快–比缩减积分的Belytschko-Tsay 壳慢2.5 倍
–当沙漏模式难以控制时推荐使用
–在每一层的单元中心平均各层应力结果
•PLANE162 – 2-D, 4-节点体
–3节点三角形单元(不推荐)
–仅支持Lagrangian 算法
–UX, UY, VX, VY, AX, AY 自由度
–对于轴对称模式,Y 轴= 对称轴
–不允许混用2D 和3D 单元类型
–不允许全积分选项
•PLANE162 KEYOPT 设置:
–Keyopt(2) –面积加权或体积加权( AXISYM )
–Keyopt(3) –平面应力,轴对称或平面应变
|__在给定的分析中仅仅可以使用一种2-D类型(如在一个模型中不能同时轴对称和平面应力单元)
•在X-Y 平面建立PLANE162 单元
•PLANE162单元不要定义实常数
•支持许多材料模式(如塑性,复合材料,Mooney-Rivlin橡胶材料)
•RSYS 支持位移和应力(不包括应变)
•Lagrangian 算法基于大应变理论,根据此理论实体被离散化,并且当网格随时间物理变形时几何体不断更新。
该算法同样使用于隐式ANSYS中。
•3-D 梁单元适用于刚体旋转,因为它不产生应变
–用三个节点定义单元
–第三个节点用来定义梁的方向
–可以定义许多标准的梁横截面
•3-D 杆单元只能承受轴向载荷
–用三个节点定义单元
–第三个节点用来定义初始杆方向
只要可能尽量避免小单元,因为它将大大减小时间步,从而增加求解时间。
如果小单元不可避免,使用质量缩放(见第2-2章)。
减少使用三角形、四面体和棱柱单元,尽管程序支持,但不推荐使用。
避免尖角单元和翘曲的壳,因为它们将降低结果精度。
在需要沙漏控制的地方使用全积分单元,但是全积分六面体单元会导致体积锁定(由于泊松比接近于0.5)和剪切锁定(如剪支梁的弯曲)。