接触操作系统有哪些

LS-DYNA常见问题汇总2.
0资料整理:yuminhust2005职务:A2:LS-DYNA版主资料来源:网络和自己的总结更新时间:2008-10-12LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12目录1.
LS-DYNA简介.
12.
单位制度53.
质量缩放64.
长分析时间75.
准静态86.
计算不稳定97.
负体积108.
能量平衡119.
沙漏控制1310.
阻尼1411.
ASCIIoutputforMPPviabinout.
1512.
接触概述1613.
ContactSoft1接触Soft=11714.
LS-DYNA中的接触界面模拟.
1815.
LS-DYNA中夹层板(sandwich)的模拟2516.
怎样进行二次开发2617.
如何实现程序的连续计算2718.
如何得知LS-DYNA为一套完整的软件2819.
LS-DYNA与市面上其它的前处理软件兼容性如何2820.
LS-DYNA似乎很重视「ContactAlgorithm」,这是为什幺2821.
如果要利用LS-DYNA进行MPP(平行运算)的计算,硬件配备及操作系统有无特殊需求2822.
FEMB能够自动产生有限元素网格吗2823.
实际的产品CAD图文件,有许多复杂的几何造型,这些feature是否都该纳入分析考量2824.
LS-DYNA是否只能进行壳元素的分析2925.
在LS-DYNA里该如何强制刚体运动2926.
在LS-POST如何显示应变分布图2927.
强制结构体绕着特定轴旋转的时候,为何会出现奇怪的变形2928.
LS-DYNA该如何进行振动模态的分析设定3029.
究竟采用Implicit及Explicit的分析方式会有多大的差异性3030.
要如何设定才能让LS-DYNA自动进行Implicit/Explicit的切换3031.
LS-DYNA是唯一能够将Implicit与Explicitcode整合在一起的软件,这对实际的分析工作有何帮助3032.
请问以Explicit方式进行分析,影响运算时间的因素有哪些3033.
势流科技您好,请问弹簧/阻尼系统该采用什幺元素来仿真比较恰当3134.
请问安装LS-DYNA(MPP)版本时,操作系统需有哪些额外的设定3135.
在Linux系统下如何指定MPP的computingnode3136.
请问在执行MPP运算时,有没有类似PC版本Ctrl+C的功能以进行不同的switch3137.
如果我想要了解特定截面的力量分布情形,该设定哪些项目呢3238.
请问LS-DYNA在进行Implicit分析时,有没有什幺该注意的地方3239.
请教一下,LS-DYNA有没有什幺功能可以让上/下表面的网格移到中性平面来3240.
请问我该如何在后处理程序LS-Pre/Post里头看到壳元素的厚度3241.
请问在LS-Pre/Post里面该如何操作才能看到各相材料的变形及相关的力学资料3242.
请问LS-Pre/Post有没有办法改变SPH颗粒的大小3343.
请问LS-DYNA有没有办法仿真出「破裂」或「碎裂」的状况3344.
请问究竟是采用Automaticcontact好还是一一指定接触界面的方式比较好3345.
在完成前处理的设定后,如何知道每一个part的重量33LS-DYNA常见问题汇总2.
0yuminhust20052008-10-1246.
请问如果我自己有写程序产生出曲线资料,或是已将datapoint存成Excel档案,能不能直接读入FEMB3347.
请问贵公司所代理的LS-DYNA软件,其采用的单位制为何3348.
请问在loadcurve的输入方式除了一点一点慢慢key-in之外,有没有其它更快的方式3349.
如何处理LS-DYNA中的退化单元3550.
LS-DYNA中对于单元过度翘曲的情况有何处理方法.
3551.
在ANSYS计算过程中结果文件大于8GB时计算自动中断,如何解决这个问题3552.
关于梁、壳单元应力结果输出的说明3553.
LS-DYNA求解有时为什么有负的滑动能.
3654.
在DYNA中如何考虑材料失效.
3655.
在LS-DYNA中能否施加跟随力和跟随力矩3656.
如果在工程上遇到壳的厚度是坐标位置的函数时,这样的壳单元模型如何建立3657.
不同自由度的单元在merge时,ANSYS如何处理自由度的不匹配问题.
3758.
如何方便地建立压力容器的有限元模型3759.
ANSYS是否具有混合分网的功能3860.
在ANSYS中怎样给面施加一个非零的法向位移约束3861.
如何得到径向和周向的计算结果3862.
如何正确理解ANSYS的节点坐标系3863.
如何考虑结构分析中的重力3964.
对于具有高度不规则横截面的3D模型什么是最佳网格划分方法3965.
在交互方式下如何施加任意矢量方向的表面载荷4066.
LS-DYNA96版后在爆炸及流固耦合方面的功能增强.
4067.
ANSYS坐标系总结4068.
显式隐式分析转换的注意事项4169.
利用LS-DYNA进行接触分析应该注意的一些问题.
4270.
LS-DYNA求解中途退出的解决方案.
4271.
液面晃动4372.
复杂几何模型的系列网格划分技术5273.
LS-DYNA典型的射流模拟过程.
54LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第1页1.
LS-DYNA简介LS-DYNA是世界上最著名的通用显式动力分析程序,能够模拟真实世界的各种复杂问题,特别适合求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等非线性动力冲击问题,同时可以求解传热、流体及流固耦合问题.
在工程应用领域被广泛认可为最佳的分析软件包.
与实验的无数次对比证实了其计算的可靠性.
由J.
O.
Hallquist主持开发完成的DYNA程序系列被公认为是显式有限元程序的鼻祖和理论先导,是目前所有显式求解程序(包括显式板成型程序)的基础代码.
1988年J.
O.
Hallquist创建LSTC公司,推出LS-DYNA程序系列,并于1997年将LS-DYNA2D、LS-DYNA3D、LS-TOPAZ2D、LS-TOPAZ3D等程序合成一个软件包,称为LS-DYNA.
PC版的前后处理采用ETA公司的FEMB,新开发的后处理为LS-POST.
LS-DYNA的最新版本是2001年5月推出的960版.
(先前的说法,现在最新为971版)LS-DYNA功能特点LS-DYNA程序960版是功能齐全的几何非线性(大位移、大转动和大应变)、材料非线性(140多种材料动态模型)和接触非线性(50多种)程序.
它以Lagrange算法为主,兼有ALE和Euler算法;以显式求解为主,兼有隐式求解功能;以结构分析为主,兼有热分析、流体-结构耦合功能;以非线性动力分析为主,兼有静力分析功能(如动力分析前的预应力计算和薄板冲压成型后的回弹计算);军用和民用相结合的通用结构分析非线性有限元程序.
1.
分析能力:非线性动力学分析多刚体动力学分析准静态分析(钣金成型等)热分析结构-热耦合分析流体分析:欧拉方式任意拉格郎日-欧拉(ALE)流体-结构相互作用不可压缩流体CFD分析有限元-多刚体动力学耦合分析(MADYMO,CAL3D)水下冲击失效分析裂纹扩展分析实时声场分析设计优化隐式回弹多物理场耦合分析自适应网格重划并行处理(SMP和MPP)2.
材料模式库(140多种)金属塑料玻璃泡沫编制品橡胶(人造橡胶)蜂窝材料复合材料混凝土和土壤炸药推进剂粘性流体用户自定义材料3.
单元库体单元薄/厚壳单元梁单元焊接单元离散单元束和索单元安全带单元节点质量单元SPH单元4.
接触方式(50多种)柔体对柔体接触柔体对刚体接触刚体对刚体接触边-边接触侵蚀接触充气模型约束面刚墙面拉延筋5.
汽车行业的专门功能安全带滑环预紧器牵引器传感器加速计LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第2页气囊混合III型假人模型6.
初始条件、载荷和约束功能初始速度、初应力、初应变、初始动量(模拟脉冲载荷);高能炸药起爆;节点载荷、压力载荷、体力载荷、热载荷、重力载荷;循环约束、对称约束(带失效)、无反射边界;给定节点运动(速度、加速度或位移)、节点约束;铆接、焊接(点焊、对焊、角焊);二个刚性体之间的连接-球形连接、旋转连接、柱形连接、平面连接、万向连接、平移连接;位移/转动之间的线性约束、壳单元边与固体单元之间的固连;带失效的节点固连.
7.
自适应网格剖分功能自动剖分网格技术通常用于薄板冲压变形模拟、薄壁结构受压屈曲、三维锻压问题等大变形情况,使弯曲变形严重的区域皱纹更加清晰准确.
对于三维锻压问题,LS-DYNA主要有两种方法:自适应网格剖分和任意拉格朗日-欧拉网格(ALE)网格进行Rezoning),三维自适应网格剖分采用的是四面体单元.
8.
ALE和Euler列式ALE列式和Euler列式可以克服单元严重畸变引起的数值计算困难,并实现流体-固体耦合的动态分析.
在LS-DYNA程序中ALE和Euler列式有以下功能:·多物质的Euler单元,可达20种材料;·若干种Smoothing算法选项;·一阶和二阶精度的输运算法;·空白材料;·Euler边界条件:滑动或附着条件;·声学压力算法;·与Lagrange列式的薄壳单元、实体单元和梁单元的自动耦合.
9.
SPH算法SPH(SmoothedParticleHydrodynamics)光顺质点流体动力算法是一种无网格Lagrange算法,最早用于模拟天体物理问题,后来发现解决其它物理问题也是非常有用的工具,如连续体结构的解体、碎裂、固体的层裂、脆性断裂等.
SPH算法可以解决许多常用算法解决不了的问题,是一种非常简单方便的解决动力学问题的研究方法.
由于它是无网格的,它可以用于研究很大的不规则结构.
SPH算法适用于超高速碰撞、靶板贯穿等过程的计算模拟,下图是泰勒杆冲击试验模拟.
10.
边界元法LS-DYNA程序采用边界元法BEM(BoundaryElementMethod)求解流体绕刚体或变形体的稳态或瞬态流动,该算法限于非粘性和不可压缩的附着流动.
11.
隐式求解用于非线性结构静动力分析,包括结构固有频率和振型计算.
LS-DYNA中可以交替使用隐式求解和显式求解,进行薄板冲压成型的回弹计算、结构动力分析之前施加预应力等.
12.
热分析LS-DYNA程序有二维和三维热分析模块,可以独立运算,也可以与结构分析耦合,可进行稳态热分析,也可进行瞬态热分析,用于非线性热传导、静电场分析和渗流计算.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第3页热传导单元:8节点六面体单元(3D),4节点四边形单元(2D);材料类型:各向同性、正交异性热传导材料,可以与温度相关,以及各向同性热传导材料的相变;边界条件:给定热流flux边界,对流convection边界,辐射radiation边界,以及给定温度边界,它们可随时间变化;给定初始温度,可计算二个物体接触界面的热传导和热辐射,给定材料内部热生成(给定热源);热分析采用隐式求解方法,过程控制有:稳态分析还是瞬态分析;线性问题还是非线性问题;时间积分法:Crank-Nicholson法(a=0.
5)和向后差分法(a=1);求解器:直接法或迭代法;自动时步长控制.
13.
不可压缩流场分析LS-DYNA不可压缩流求解器是960版新增加的功能,用于模拟分析瞬态、不可压、粘性流体动力学现象.
求解器中采用了超级计算机的算法结构,在确保有限元算法优点的同时计算性能得到大幅度提高,从而在广泛的流体力学领域具有很强的适用性.
14.
多功能控制选项多种控制选项和用户子程序使得用户在定义和分析问题时有很大的灵活性.
输入文件可分成多个子文件;用户自定义子程序;二维问题可以人工控制交互式或自动重分网格(REZONE);重启动;数据库输出控制;交互式实时图形显示;开关控制-可监视计算过程的状态;对32位计算机可进行双精度分析.
15.
前后处理功能LS-DYNA利用ANSYS、LS-INGRID、ETA/FEMB及LS-POST强大的前后处理模块,具有多种自动网格划分选择,并可与大多数的CAD/CAE软件集成并有接口.
后处理:结果的彩色等值线显示、梯度显示、矢量显示、等值面、粒子流迹显示、立体切片、透明及半透明显示;变形显示及各种动画显示;图形的PS、TIFF及HPGL格式输出与转换等.
16.
支持的硬件平台LS-DYNA960版的SMP版本和MPP版本是同时发行的.
MPP版本使一项任务可同时在多台分布计算机上进行计算,从而最大限度地利用已有计算设备,大幅度减少计算时间.
计算效率随计算机数目增多而显著提高.
LS-DYNA960版的SMP版本和MPP版本可以在PC机(NT、LINUX环境)、UNIX工作站、超级计算机上运行.
LS-DYNA应用领域汽车工业:碰撞分析气囊设计乘客被动安全部件加工航空航天"鸟撞;叶片包容;飞机结构冲击动力分析;碰撞,坠毁;冲击爆炸及动态载荷;火箭级间分离模拟分析;宇宙垃圾碰撞;特种复合材料设计:制造业冲压锻造铸造切割,等建筑业:地震安全混凝土结构爆破拆除公路桥梁设计LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第4页国防:内弹道和终点弹道;装甲和反装甲系统;穿甲弹与破甲弹设计;战斗部结构设计;冲击波传播;侵彻与开坑;空气,水与土壤中爆炸;核废料容器设计,等电子领域:跌落分析包装设计热分析电子封装石油工业液体晃动;完井射孔;管道设计;爆炸切割;事故模拟;海上平台设计其它应用玻璃成型生物医学体育器材(高尔夫杆,高尔夫球,棒球杆,头盔)实例图片:http://forum.
simwe.
com/thread-826139-1-1.
htmlLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第5页2.
单位制度相信做仿真分析的人第一个需要明确的就是一致单位系统(ConsistentUnits).
计算机只认识0&1、只懂得玩数字,它才不管你用的数字的物理意义.
而工程师自己负责单位制的统一,否则计算出来的结果没有意义,不幸的是大多数老师在教有限元数值计算时似乎没有提到这一点.
见下面LS-DYNAFAQ中的定义:1力单位=1质量单位*1加速度单位1加速度单位=1长度单位/1时间单位的平方下表提供了一些单位换算关系:MASSLENGTHTIMEFORCESTRESSENERGYDENSITYYOUNG'sVelocity(56.
3KMPH)GRAVITYkgmsNPaJoule7.
83E+032.
07E+1115.
659.
806kgcms1.
e-02N7.
83E-032.
07E+091.
56E+039.
81E+02kgcmms1.
e+04N7.
83E-032.
07E+031.
569.
81E-04kgcmus1.
e+10N7.
83E-032.
07E-031.
56E-039.
81E-10kgmmmsKNGPaKN-mm7.
83E-062.
07E+0215.
659.
81E-03gmcmsdynedy/cm2erg7.
83E+002.
07E+121.
56E+039.
81E+02gmcmus1.
e+07NMbar1.
e7Ncm7.
83E+002.
07E+001.
56E-039.
81E-10gmmms1.
e-06NPa7.
83E-032.
07E+111.
56E+049.
81E+03gmmmmsNMPaN-mm7.
83E-032.
07E+0515.
659.
81E-03tonmmsNMPaN-mm7.
83E-092.
07E+051.
56E+049.
81E+03lbfs2/ininslbfpsilbf-in7.
33E-043.
00E+076.
16E+02386slugftslbfpsflbf-ft15.
24.
32E+0951.
3332.
17kgfs2/mmmmskgfkgf/mm2kgf-mm8.
02E-107.
00E+021.
56E+04(Japan)kgmmsmN1000Pa7.
83E-062.
07E+089.
81E+02gmcmms100000Pa7.
83E+002.
07E+06LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第6页3.
质量缩放质量缩放指的是通过增加非物理的质量到结构上从而获得大的显式时间步的技术.
在一个动态分析中,任何时候增加非物理的质量来增大时间步将会影响计算结果(因为F=ma).
有时候这种影响不明显,在这种情况下增加非物理的质量是无可非议的.
比如额外的质量只增加到不是关键区域的很少的小单元上或者准静态的分析(速度很小,动能相对峰值内能非常小).
总的来说,是由分析者来判断质量缩放的影响.
你可能有必要做另一个减小或消除了质量缩放的分析来估计质量增加对结果的灵敏度.
你可以通过人工有选择的增加一个部件的材料密度来实现质量缩放.
这种手动质量缩放的方法是独立于通过设置*Control_timestep卡DT2MS项来实现的自动质量缩放.
当DT2MS设置为一个负值时,质量只是增加到时间步小于TSSFAC*|DT2MS|的单元上.
通过增加这些单元的质量,它们的时间达到TSSFAC*|DT2MS|.
有无数种TSSFAC和DT2MS的组合可以得到同样的乘积,因而有相同的时间步,但是对于每一种组合增加的质量将是不一样的.
一般的趋势是TSSFAC越小,增加的质量越多.
作为回报,当TSSFAC减小时计算稳定性增加(就像在没有做质量缩放的求解中一样).
如果TSSFAC缺省的值0.
9会导致稳定性问题,可以试试0.
8或者0.
7.
如果你减小TSSFAC,你可以相应增加|DT2MS|,这样还是可以保证时间步乘积不变.
为了确定什么时候和位置质量自动增加了,可以输出GLSTAT和MATSUM文件.
这些文件允许你绘出完整的模型或者单独部件所增加的质量对时间的曲线.
为了得到由壳单元组成的部件增加的质量云图,将*database_extent_binary卡的STSSZ项设置为3.
这样你可以用ls-prepost绘出每个单元的质量增加量的云图,具体方法是通过选择Fcomp>Misc>timestepsize.
在*control_timestep中设置DT2MS正值和负值的不同之处如下:负值:初始时间步将不会小于TSSFAC*-DT2MS.
质量只是增加到时间步小于TSSFAC*|DT2MS|的单元上.
当质量缩放可接受时,推荐用这种方法.
用这种方法时质量增量是有限的.
过多的增加质量会导致计算任务终止.
正值:初始时间将不会小于DT2MS.
单元质量会增加或者减小以保证每一个单元的时间步都一样.
这种方法尽管不会因为过多增加质量而导致计算终止,但更难以作出合理的解释.
*control_timestep卡中的参数MS1ST控制是否只是在初始化时增加一次质量(MS1ST=1)还是任何需要维持由DT2MS所指定的时间步时都增加质量(MS1ST=0).
你可以通过在*control_termination卡片中设置参数ENDMAS来控制当质量增加到初始质量一定比率时终止计算(只对自动质量缩放有效)可变形点焊梁的质量缩放.
*mat_spotweld卡的质量缩放参数DT只影响点焊单元.
如果*control_timestep卡中没有指定质量缩放(DT2MS=0),而且时间由可变形点焊控制,可以用参数DT来在初始化时增加惯量到点焊单元上来提高时间步达到DT指定的值.
当DT不为0时,增加到可变形点焊梁元上的质量会输出到d3hsp文件里.
MATSUM中动量和动能不受增加到可变形点焊上的质量的影响.
GSLTAT中DOES和总的KE受增加的质量的影响.
考虑三种调用可变形点焊的质量缩放的情况:1.
当DT2MS为负值*mat_spotweld卡DT=0时,尽管在d3hsp文件中可变形点焊质量增量百分比不真实.
下面几个值是正确的:d3hsp中"addedspotweldmass";第一个时间步之后的"addedmass"&"percentageincrease";glstat和matsum中的"addedmass".
2.
当DT2MS为负值且*mat_spotweld卡DT≠0时,可变形点焊质量增加不会包含在d3hsp、glstat、matsum文件中的"addedmass"里.
这非常容易令人误解.
用户必须检查d3hsp文件的"addedspotweldmass".
建议不要同时使用两种质量缩放标准,推荐使用第一种方法(即负的DT2MS&DT=0).
3.
如果DT2MS=0且DT≠0,初始时间步将不考虑增加点焊的质量,但是之后每一个周期时间步都会增加10%,直到时间步达到正确的值(考虑点焊质量增加).
glstat&matsum不包含"addedmass"的行.
注意质量增加会引起能量比率增长.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第7页4.
长分析时间当用显式时间积分时,对于仿真非常小的部件而分析时间又要相当长时没有好的方法.
质量缩放(mass-scaling)增加了需要确认非物理质量的增加不会显著影响计算结果的负担.
当使用时间缩放(time-scaling)时也有同样的问题.
时间缩放(time-scaling)是指为了减小需要的时间步数,通过增加加载速率而缩短仿真时间.
要确认时间步不是仅由很少的小单元或者刚度大单元控制,可以通过在d3hsp文件中搜索"smallest"来显示100个最小的时间步单元.
如果只有很少的几个单元控制时间步,可以把那些单元及邻近区域重新remesh或者把它们变成刚体.
可是仅运行必要长的时间是很明显的.
这意味着在一个跌落分析的情况时,给跌落物体一个初速度,把它放在离地面一个非常小的距离.
冲击之后,仅运行足够得到需要的结果的时间.
值得注意的是对于一个长时间的仿真,如果时间步数超过了50万步,最好使用双精度版本的LS-DYNA求解器,使截断误差的影响最小化.
运行双精度版本要增加30%的时间.
对于长时间的分析,自动显式/隐式转换可能是一个选择.
使用这种方法,用户可以指定在一个时间段内使用隐式积分.
隐式积分的优点是时间步不由单元尺寸控制,所以可以得到大的时间步.
当然,隐式计算也非常点用cpu时间.
而且,目前并不是所有的LS-DYNA的功能和材料都在隐式分析中实现(大部分已经实现).
下面的FEAinformationnewsletter里讨论了显式/隐式转换(http://www.
feapublications.
com/pages/pdfnews/3feadec.
pdf).
Seealso:mass_scaling,quasistaticLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第8页5.
准静态动态松驰(Dynamicrelaxation)并不是有意为一般的准静态(quasi-static)分析设置的.
它适合于当预载只产生小的弹情况应变的施加预载,或者初始化系统到一个预定义的几何形状[1].
但对其它更多情况并不适合.
你可以通过做一个常规的显示仿真来模拟准静态分析,通过按需要调用时间/质量缩放(time-scaling,mass-scaling)来在可接受的时间内得到结果,但这种方法是需要技巧地.
你必须监测系统动能按希望的使惯性效应最小化.
基本上动能相对内能应该保持在一个较小的值.
时间缩放是指加载比在准静态实验里更快,以减少总的仿真时间.
关于质量缩放更多内容可以看"mass_scaling"一节.
或者你可以尝试用LS-DYNA运行一个隐式静力分析.
可以看用户手册里的卡片*control_implicit_…"和AppendixM.
Seealso:gravity.
txt,readme.
preload,mass_scaling,long_run_times,implicit.
general,quick_initialization.
Note[1]:初始化到预定义的几何1.
从第一次分析的最终状态输出一个节点位移文件.
(这一部分未按原文翻译)注意d3plot文件里不包含节点转动信息,因此转动输出为0.
这对初始化壳和梁单元会是个问题.
LS-Prepost有一个选项是输出节点位移,在Output->NodalDisplacements里.
但是这个输出是i8,3e16格式的,但需要的是i8,3e15,所以要注意修改一下.
如果你做了一个正常的动态松驰分析来得到初始状态,一个预定义位移和转动的drdisp.
sif文件在DR阶段结束时会自动创建.
2.
在第二次分析时,快速的初始化到第一步输出的预定义的几何.
你需要设置卡片*control_dynamic_relaxation里的参数IDRFLG=2,而且在命令行里指定"m=filename"(其中filename指第一步创建的文件).
这样在瞬态分析之前,LS-DYNA会自动做一个100步的预分析来使节点根据文件filename指定的数据移动到指定值.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第9页6.
计算不稳定一些表示计算不稳定的消息如:"out-of-rangevelocities"速度超出范围"negativevolumeinbrickelement"体单元负体积"terminationduetomassincrease"因质量增加而终止用来克服显式求解中的不稳定的方法如下:首先(也是最重要的)是使用可获得的最新的LS-DYNA版本.
最新的执行块可以从ftp://user@ftp.
lstc.
com上下载(注:前提是你有访问权限).
联系LSTC获得user帐号的密码.
最新的BETA版执行块可以在ftp://ftp.
lstc.
com/outgoing/ls971上找到(不需要密码,但lstc公司对ftp访问有IP限制).
其次是增加d3plot的输出频率到可以显示出不稳定的出现过程.
这可以提供导致不稳定性发生的线索.
其它的不些解决数值不稳定性的技巧:*试着用双精度LS-DYNA版本运行一次*试着减小时间步(timestep)缩放系数(即使使用了质量缩放mass-scaling)*单元类型和/或沙漏(hourglass)控制.
对出现不稳定的减缩体和壳单元,试着用沙漏控制type4和沙漏系数0.
05.
或者试着用类型16的壳单元,沙漏控制type8.
如果壳响应主要是弹性,设置BWC=1和PROJ=1(仅对B-T壳).
避免使用type=2体单元.
对体单元部件,在厚度方向最少用两个体单元.
*接触.
设置接触的bucketsorts之间周期数为0,这样会使用缺省的分类间隔.
如果参与接触的两个部件的相对速度异常的大,可能需要减小bucketsort的间隔(比如减小到5,2甚至1).
如果仿真过程中有明显的接触穿透出现,转换到使用*contact_automatic_surface_to_surface或者*contact_automatic_single_surface,并设置SOFT=1.
确保几何考虑了壳单元的厚度.
如果壳非常薄,比如小于1mm,放大或者设置接触厚度到一个更加合理的值.
*避免冗余的接触定义,也就是说不要对同样的两个部件定义多于一个的接触对.
*查找出现不稳定的部件的材料定义中的错误(比如误输入,不一致的单位系统等)*关掉所有的*damping这些技巧是一些通用的方法,可能并不适合于所有的情况.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第10页7.
负体积泡沫材料的负体积(或其它软的材料)对于承受很大变形的材料,比如说泡沫,一个单元可能变得非常扭曲以至于单元的体积计算得到一个负值.
这可能发生在材料还没有达到失效标准前.
对一个拉格朗日(Lagrangian)网格在没有采取网格光滑(meshsmoothing)或者重划分(remeshing)时能适应多大变形有个内在的限制.
LS-DYNA中计算得到负体积(negativevolume)会导致计算终止,除非在*control_timestep卡里面设置ERODE选项为1,而且在*control_termination里设置DTMIN项为任何非零的值,在这种情况下,出现负体积的单元会被删掉而且计算继续进行(大多数情况).
有时即使ERODE和DTMIN换上面说的设置了,负体积可能还是会导致因错误终止.
有助于克服负体积的一些方法如下:*简单的把材料应力-应变曲线在大应变时硬化.
这种方法会非常有效.
*有时候修改初始网格来适应特定的变形场将阻止负体积的形成.
此外,负体积通常只对非常严重的变形情况是个问题,而且特别是仅发生在像泡沫这样的软的材料上面.
*减小时间步缩放系数(timestepscalefactor).
缺省的0.
9可能不足以防止数值不稳定.
*避免用全积分的体单元(单元类型2和3),它们在包含大变形和扭曲的仿真中往往不是很稳定.
全积分单元在大变形的时候鲁棒性不如单点积分单元,因为单元的一个积分点可能出现负的Jacobian而整个单元还维持正的体积.
在计算中用全积分单元因计算出现负的Jacobian而终止会比单元积分单元来得快.
*用缺省的单元方程(单点积分体单元)和类型4或者5的沙漏(hourglass)控制(将会刚化响应).
对泡沫材料首先的沙漏方程是:如果低速冲击type6,系数1.
0;高速冲击type2或者3.
*对泡沫用四面体(tetrahedral)单元来建模,使用类型10体单元.
*增加DAMP参数(foammodel57)到最大的推荐值0.
5.
*对包含泡沫的接触,用*contact选项卡B来关掉shootingnodelogic.
*使用*contact_interior卡用partset来定义需要用contact_interior来处理的parts,在set_part卡1的第5项DA4来定义contact_interior类型.
缺省类型是1,推荐用于单一的压缩.
在版本970里,类型1的体单元可以设置type=2,这样可以处理压缩和减切混合的模式.
*如果用mat_126,尝试ELFORM=0*尝试用EFG方程(*section_solid_EFG).
因为这个方程非常费时,所以只用在变形严重的地方,而且只用于六面体单元.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第11页8.
能量平衡GLSTAT(参见*database_glstat)文件中报告的总能量是下面几种能量的和:内能internalenergy动能kineticenergy接触(滑移)能contact(sliding)energy沙漏能houglassenergy系统阻尼能systemdampingenergy刚性墙能量rigidwallenergyGLSTAT中报告的弹簧阻尼能"Springanddamperenergy"是离散单元(discreteelements)、安全带单元(seatbeltelements)内能及和铰链刚度相关的内能(*constrained_joint_stiffness…)之和.
而内能"InternalEnergy"包含弹簧阻尼能"Springanddamperenergy"和所有其它单元的内能.
因此弹簧阻尼能"Springanddamperenergy"是内能"Internalenergy"的子集.
由SMP5434a版输出到glstat文件中的铰链内能"jointinternalenergy"跟*constrained_joing_stiffness不相关.
它似乎与*constrained_joint_revolute(_spherical,etc)的罚值刚度相关连.
这是SMP5434a之前版本都存在的缺失的能量项,对MPP5434a也一样.
这种现象在用拉格朗日乘子(LagrangeMultiplier)方程时不会出现.
与*constrained_joint_stiffness相关的能量出现在jntforc文件中,也包含在glstat文件中的弹簧和阻尼能和内能中.
回想弹簧阻尼能"springanddamperenergy",不管是从铰链刚度还是从离散单元而来,总是包含在内能里面.
在MATSUM文件中能量值是按一个part一个part的输出的(参见*database_matsum).
沙漏能Hourglassenergy仅当在卡片*control_energy中设置HGEN项为2时才计算和输出.
同样,刚性墙能和阻尼能仅当上面的卡片中RWEN和RYLEN分别设置为2时才会计算和输出.
刚性阻尼能集中到内能里面.
质量阻尼能以单独的行"systemdampingenergy"出现.
由于壳的体积粘性(bulkviscosity)而产生的能量耗散(energydissipated)在版本970.
4748之前是不计算的.
在后续子版本中,设置TYPE=-2来在能量平衡中包含它.
最理想的情况下能量平衡:总能量totalenergy=初始总能量+外力功externalwork换句话说,如果能量比率energyratio(指的是glstat中的totalenergy/initialenergy,实际上是totalenergy/(initialenergy+externalwork))等于1.
0.
注意,质量缩放而增加质量可能会导致能量比率增加.
注意在LSprepost的History>Globalenergies中不包含删掉的单元(erodedelements)的能量贡献,然而GLSTAT文件中的能量包含了它们.
注意它们的贡献可以通过ASCII>glstat中的"ErodedKineticEnergy"&"ErodedInternalEnergy"来绘制.
侵蚀能量(Erodedenergy)是与删掉的单元相关的内能和删掉的节点相关的动能.
典型来说,如果没有单元删掉"energyratiow/oerodedenergy"等于1,如果有单元被删掉则小于1.
删掉的单元与"totalenergy/initialenergy"比率没有关系.
总能量比率增加要归于其它原因,比如增加质量.
重述一下,将一个单元删掉时,文件glstat中的内能和动能不会反映能量的丢失.
取而代之的是能量的丢失记录在glstat文件的"erodedinternalenergy"&"erodedkineticenergy"中.
如果用内能减去"erodedinternalenergy"将得到分析中还存在的单元的内能.
对动能也一样.
matsum文件中的内能和动能只包含余下(noneroded)的单元的贡献.
注意,如果在*control_contact卡中将ENMASS设置为2,则与删掉的单元的相关的节点不会删掉,"erodedkineticenergy"是0.
在LSprepost中History>Global只是动能和内能的简单相加,因此不包含接触能和沙漏能等的贡献.
壳的负内能:为了克服这种不真实的效应--关掉考虑壳的减薄(ISTUPDin*control_shell)--调用壳的体积粘性(setTYPE=-2在*control_bulk_viscosity卡中)--对在matsum文件中显示为负的内能的parts使用*damping_part_stiffness;先试着用一个小的值,比如0.
01.
如果在*control_energy中设置RYLEN=2,因为刚性阻尼而能会计算且包含在内能中.
正的接触能:当在接触定义中考虑了摩擦时将得到正的接触能.
摩擦将导致正的接触能.
如果没有设置接触LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第12页阻尼和接触摩擦系数,你将会看到净接触能为零或者一个很小的值(净接触能=从边和主边能量和).
所说的小是根据判断-在没有接触摩擦系数时,接触能为峰值内能的10%内可以被认为是可接受的.
负的接触能:突然增加的负接触能可能是由于未检测到的初始穿透造成的.
在定义初始几何时考虑壳的厚度偏置通常是最有效的减小负接触能的步骤.
查阅LS-DYNA理论手册的23.
8.
3&23.
8.
4节可得到更多接触能的信息.
负接触能有时候因为parts之间的相对滑动而产生.
这跟摩擦没有关系,这里说的负接触能从法向接触力和法向穿透产生.
当一个穿透的节点从它原来的主面滑动到临近的没有连接的主面时,如果穿透突然检测到,则产生负的接触能.
如果内能为负接触能的镜像,例如glstat文件中内能曲线梯度与负接触能曲线梯度值相等,问题可能是非常局部化的,对整体求解正确性冲击较小.
你可以在LS-prepost中分离出有问题的区域,通绘制壳单元部件内能云图(Fcomp>Misc>Internalenergy).
实际上,显示的是内能密度,比如内能/体积.
内能密度云图中的热点通常表示着负的接触能集中于那里.
如果有多于一个的接触定义,sleout文件(*database_sleout)将报告每一个接触对的接触能量,因此缩小了研究负接触能集中处的范围.
克服负接触能的一般的建议如下:-消除初始穿透(initialpenetration).
(在message文件中查找"warning")-检查和排除冗余的接触条件.
不应该在相同的两个parts之间定义多于一个的接触.
-减小时间步缩放系数-设置接触控制参数到缺省值,SOFT=1&IGNORE=1除外(接触定义选项卡C)-对带有尖的边的接触面,设置SOFT=2(仅用于segment-to-segment接触).
而且,在版本970中推荐设置SBOPT(之前的EDGE)为4对于部件之间有相对滑移的SOFT=2的接触.
为了改进edge-to-edgeSOFT=2接触行为,设置DEPTH=5.
请注意SOFT=2接触增加了额外的计算开消,尤其是当SBOPT或者DEPTH不是缺省值时,因此应该仅在其它接触选项(SOFT=0或者SOFT=1)不能解决问题时.
模型的细节可能会指示可用其它的一些方法.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第13页9.
沙漏控制沙漏(hourglass)模式是一种非物理的零能变形模式,产生零应变和应力.
沙漏模式仅发生在减缩积分(单积分点)体、壳和厚壳单元上.
LS-DYNA里面有多种算法用于抑制沙漏模式.
缺省的算法(type1)通常不是最有效的算法,但却是最经济的.
一种完全消除沙漏的方法是转换到全积分或者选择减缩积分(S/R)方程的单元.
但这种方法是一种下策.
例如,第一,类型2体单元比缺省的单点积分体单元计算开消大;其二,在大变形应用时更不稳定(更容易出现负体积);其三,类型2体单元当单元形状比较差时在一些应用中会趋向于剪切锁死(shear-lock),因而表现得过于刚硬.
三角形壳和四面体单元没有沙漏模式,但缺点是在许多应用中被认为过于刚硬.
减小沙漏的一个好的方法是细化网格,但这当然并不总是现实的.
加载方式会影响沙漏程度.
施加压力载荷优于在单点上加载,因为后者更容易激起沙漏模式.
为了评估沙漏能,在*control_energy卡片中设置HGEN=2,而且用*database_glstat和*database_matsum卡分别输出系统和每一个部件的沙漏能.
这一点是要确认非物理的沙漏能相对于每一个part的峰值内能要小(经验上来说Misc>hourglassenergy.
对于流体部件,缺省的沙漏系数通常是不合适的(太高).
因此对于流体,沙漏系数通常要缩小一到两个数量级.
对流体用基于粘性的沙漏控制.
缺省的沙漏方程(type1)对流体通常是可以的.
对于结构部件一般来说基于刚性的沙漏控制(type4,5)比粘性沙漏控制更有效.
通常,当使用刚性沙漏控制时,习惯于减小沙漏系数到0.
03~0.
05的范围,这样最小化非物理的硬化响应同时又有效抑制沙漏模式.
对于高速冲击,即使对于固体结构部件,推荐采用基于粘性的沙漏控制(type1,2,3).
粘性沙漏控制仅仅是抑制沙漏模式的进一步发展,刚性沙漏控制将使单元朝未变形的方向变形.
类型8沙漏控制仅用于单元类型16的壳.
这种沙漏类型激活了16号壳的翘曲刚度,因此单元的翘曲不会使解退化.
如果使用沙漏控制8,16号壳单元可以用于解被称为扭曲梁(TwistedBeam)问题.
对于单元类型1的体和减缩积分2D体(shelltypes13&15)类型6沙漏控制调用了一种假设应变协同转动方程.
使用沙漏控制类型6和系数1.
0,一个弹性部件在厚度方向仅仅需要划分一层类型1的体单元就可以获得正确的弯曲刚度.
在隐式计算里面,对于类型1的体单元应该总是使用类型6的沙漏控制(实际上,在V970里面这是自动设置的).
(Moreontype6HGcontrolfromLeeBindeman)类型6的沙漏控制与类型4,5不在于它用了一个假设应变场和材料属性来估算出假设应力场.
这个应力在单元封闭域内进行积分得到沙漏力,因此单元表现的像一个有同样假设应变场的全积分单元.
这种假设应变场设计成用来阻止纯弯曲中不真实的剪切变形和近似不可压材料中的体积锁死.
类型4和5的沙漏控制基于单元体积,波速和密度像在LS-DYNA理论手册中方程3.
21那样来计算沙漏刚度.
沙漏类型6主要的改进是应力场在单元域内积分.
这使得当使用大的长细比或者歪斜形状的体单元时沙漏控制非常鲁棒.
类型4和5的沙漏控制对大长细比和歪斜形状单元反应变不好,它趋向于对某些沙漏模式反应的过于刚硬而对其它模式反应得过弱.
沙漏控制类型6另一个理论上的优点是对在厚度方向只有一个单元的梁可以在弹性弯曲问题中得到准确的解.
要做到这一点,设置沙漏刚度参数为1.
0.
同样,对弹性材料方形截面杆的扭曲问题,当沙漏系数设为1.
0时可以用很少的单元来解.
然而,对于非线性材料,用粗糙的网格得到好的结果是不可能的,因为应力场不是像沙漏类型6假设的那样线性变化的.
在梁厚度方向上如果没有更多积分点的话,没有办法捕获应力场的非线性状态.
对于选择沙漏控制,下面几个问题要考虑.
对于单元有大的长细比或者明显歪斜(不管是初始还是变形过程中),推荐采用类型6的沙漏控制.
类型6的沙漏控制通常对软的材料更好,像泡沫或蜂窝材料在计算中会有非常明显的变形.
在材料不是特别软或者单元有合理的形状且网格不是太粗糙时,类型4,5和6沙漏控制似乎都能得到同样的结果.
这种情况推荐用类型4的沙漏控制,因为它比其它的更快.
类型6的沙漏控制在LS-DYNAUser'sManual中参考的Belytschko和Bindeman的论文中有更详细的描述.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第14页10.
阻尼在LS-DYNA中阻尼完全是可选的,通过使用一个*DAMPING卡片来调用.
应该知道能量可以通过其它的非*DAMPING的方式耗散,比如,因为沙漏力产生的能量,刚性墙的力产生的能量,接触摩擦力产生的能量,离散阻尼产生的内能等.
有时候,接触力可能将噪声引入到响应里.
在这种情况下,通过*CONTACT卡第二张卡的VDC参数来增加粘性阻尼,从而帮助减小噪声.
VDC以临界阻尼的百分比输入,典型的值是10到20.
*DAMPING卡片概览:LS-DYNA中的质量阻尼(Massdamping)包括*damping_global&*damping_part_mass,是用于抑止低频的结构振动模式,但此外它有抑制刚体模式的效应.
因此对经受明显刚体运动的部件,应该要么从质量阻尼中排除或者在部件经历大的刚体运动期间关掉质量阻尼;或者使用*damping_relative来替代.
通过使用*damping_relative,仅仅相对指定刚体的运动/振动被抑制.
在质量阻尼情况下临界阻尼系数是4*pi/T,其中T是要抑制的模态的周期(通常是最低阶(基频)模态).
周期可以通过特征值分析(eigenvalue)或者从一个无阻尼的瞬态分析结果来估计.
如果选择使用质量阻尼,建议使用小于临界阻尼系数的阻尼值.
取10%的临界阻尼的值,即输入0.
4*pi/T,是相当典型的值.
可以选择用同样的阻尼系数抑制所有的部件(*damping_global)或者对每一个部件指定不同的阻尼系数(*damping_part_mass).
在任何一种情况下,阻尼系数可能会随时间变化(在仿真中间关掉或打开阻尼时会有用).
*damping_part_stiffness是为了抑制高频振动和数值振荡,通常对结构振动没有明显的影响.
这种情况下阻尼系数COEF近似表示临界阻尼的一个系数.
典型的COEF值是0.
1.
如果使用刚性阻尼产生不稳定,消除阻尼或者减小COEF的值来使回复稳定(也许降低一个数量级或者在某些情况下更多).
质量和刚性阻尼在隐式瞬态分析中都实现了.
在版本970中另一个可选的阻尼是频率不相关的阻尼选项,它的目标是抑制一个范围的频率和一批部件(*damping_frequency_range).
Damping_frequency_range是由Arup的RichardSturt开发的,它的理论细节是私有的.
它开发的意图是帮助LS-DYNA来适当地处理振动预测问题中的阻尼--包括车辆NVH时间历程分析,某些地震问题和土木结构的振动问题.
*damping_frequency_range的关键点在于:-仅使用很小的阻尼,例如1%到2%-处理阻尼将轻微的减小了响应的刚度,那是因为阻尼力的应力滞后于理论上正确的阻尼力,由于需要估计频率内容.
-用户指定的频率范围最理想情况是不要超过最高值和最低值之间30%.
在这个范围之外同样可以获得阻尼,但阻尼值会减小.
这种阻尼是基于节点速度的,所以可能会由于结构模态或者刚体转动而有振荡.
在Rayleigh阻尼里,阻尼矩阵表达为质量和刚度矩阵的线性组合C=alpah*M+beta*KLS-DYNA为标准的非线性分析在单元级实现Rayleigh阻尼.
这是为了数值上的方便,因为在显式方法里不需要生成刚度矩阵K.
取而代之,通过简单的将应力在单元面积上积分得到内力.
Rayleigh阻尼作为这个应力的一个修正而实现.
版本960中的刚性阻尼(stiffnessdamping)完全重新实现,即使这样可以在960中提供COEF值,这个值与950中的BETA值是相符的,其中:COEF=BETA*(w/2)版本950和960中的刚性阻尼不完全相同.
在960中刚性阻尼方程在高频域提供一个近似的临界阻尼分数.
方程的详细信息是不公开的.
这个方程中的变化是因为使用旧版本方程时产生不稳定的频率而促使的.
在版本970的3510子版本(或者更高)中,旧的950版本的刚性阻尼方程作为一个选项提供,通过设置COEF值为负值来调用.
这个参数之后被解释成V950用户手册中所指的BETA值.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第15页11.
ASCIIoutputforMPPviabinout从970版本开始,有了二进制格式格式的所谓的ASCIIoutput文件(matsum,rcforc,etc.
).
--注单机版和smp版输出的为文本格式的ASCIIoutput文件.
MPP-DYNA不是直接写出ASCII输出文件,而是输出二进制格式数据.
这些数据可以有两种格式:dbout(参见V.
970用户手册AppendixL,P.
L3)和binout.
MPP缺省情况下会输出binout格式的数据.
卡片*database_matsum,*database_rcforc,etc的第二个参数控制输出哪一种类型的数据.
1=oldformat,例如,ASCII文件(SMP)或者dbout文件(MPP)2=newbinout格式(对MPP版这是缺省的)3=bothformats都输出LS-Prepost可以直接读取binout数据而不需要转换,但dbout文件需要用dumpbdb来转换成ASCII格式输出.
在LS-Prepost里,选项页"2″(启动界面右边的Group按钮下面),在CFD按钮下面有"Binout"按钮.
点击之后在左下面出现的子面板中点击"Load"按钮,选择求解器输出的binout文件.
之后可以进行相关曲线的绘制.
(此处修改原文部分)作为lsda软件包的一部分,有两个程序可以获得(以用户账号登陆lstcftpsite,在"lsda"目录下).
一个是"l2a",它可以用来从binout文件中提取不同的ascii文件.
另一个是"ioq",它是让你用来直接读取/浏览binout文件的小工具.
通常一个"l2a"执行块会包含在相应平台的MPP执行块的压缩包里,但你也可以用那个l2a来操作任何ls-dyna执行块(MPPorSMP)产生的binout文件.
要从binout文件中提取ASCII输出文件,在命令行执行l2a,后面包含binout文件的名字.
例如:".
/l2abinout.
0000″Binout文件是平台无关的,例如,你可以在与输出binout数据的平台同一平台或者任何其它平台下处理binout数据.
当MPPLS-DYNA执行块输出binout数据时,将会有多于一个的文件以"binout"作为根名字.
每一个这样的文件需要被独立的打开来得到所有需要的数据.
在d3hsp文件里,你可以看到像如下的部分所示的信息告许你每一个binout文件包含哪些数据.
>Thefollowingbinaryoutputfilearebeingcreated,>andcontaindataequivalenttotheindicatedasciioutputfiles>binout0000:(onprocessor0)>nodout>matsum>rcforc>abstat>rbdout>sleout>jntforc(type0)>binout0001:(onprocessor1)>jntforc>binout0003:(onprocessor3)>deforc对MPP可选的控制输出格式的方法(fromJason):如果你在pfile中输入如下行pfile:gen{dboutonly}toexecute:mpirun-np##mpp970i=…p=pfile程序将像之前一样输出dbout.
*文件(没有binout),之后你可以用dumpbdb来提取所有的ASCII文件.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第16页12.
接触概述这里提供LS-DYNA接触的一些简要信息.
讨论不包含所有的接触类型的接触选项.
更详细的LS-DYNA的接触回顾可以在www.
feapublicatoins.
com上找到FEAInformationNewsletters中关于接触模型的由四个部分组成的系列文章.
这一系列文章的四部分分别在2001年的8,9,10,12月份的新闻信里.
Automaticvs.
Non-automatic:对大多数显式分析推荐用自动接触类型.
非自动接触类型(接触方向是重要的情况)有时用于金属成型仿真,其中几何是非常直接的且接触面的方向在分析之间可以可靠的确定.
非自动接触通常推荐用于隐式分析.
类型:Type13接触(contact_automatic_single_surface)是一种单面接触(不需要定义主面),总是考虑壳的厚度且没有方向性.
因而参与接触的板壳面在建模的时候需要至少保持一个小的间隙.
为了避免初始穿透,间隙不能小于潜在会发生接触的两个壳之间厚度的平均值.
体单元之间不需要有间隙.
接触类型13的接触搜索算法比接触类型3(contact_surface_to_surface)或者接触类型a3(contact_automatic_surface_to_surface)更加复杂.
例如,type13可以处理例如壳边对面,梁对壳面的情况.
和任何其它单面接触类型一样,接触力不是直接从RCFORC文件中获取,你必需要定义contact_force_transducer_penalty来获取接触力.
文章"contact.
13vs26″提供了关于contact_automatic_single_surface和contact_automatic_general对比的一些额外信息.
(之后会贴出,敬请等待)Type3接触(contact_surface_to_surface)是双向的面对面接触,其中壳的厚度考虑选项可以通过*contact或者*control_contact(*contact优先)卡片打开或者关闭.
接触片的方向是很重要的,因为这种接触类型只在一个方向来检测潜在的接触.
在如面对面的一个双向接触中,先检测从面侧的节点对主面的穿透,然后再检测主面侧的节点对从面的穿透.
这种方法当用用设置SOFT=2来调用segment-based接触时会有例外.
接触类型a3没有方向性(从壳的中面检测任一侧的潜在接触)而且总是考虑壳厚度,从这一点考虑它非常类似于type13的接触.
接触参数的一些摘记:SOFTSOFT是*contact选项卡A的第一个参数.
SOFT的缺省值是0.
SOFT=1除了在确定接触刚度方面之外与SOFT=0差不多是一样的.
SOFT=2与SOFT=0根本上是不一样的,不但在确定接触刚度方法上,而且在搜索穿透的产生的方法上也不一样.
SOFT=2会调用所谓的"segment-basedcontact".
对于SOFT=1&2接触的摘记可以分别参考文章"contact.
soft1″&"contact.
soft2″.
(稍后推出,敬请期待)IGNORE在仿真分析中的任何一个时间点,如果一个节点突然检测到在面下面(比如说,节点运动的非常快,穿透之前没有检测到),原来的(IGNORE=0)算法仅仅是把节点移动到主面上,不施加任何力(称之为"shootingnodelogic").
如果shootingnodelogic选项被关掉(SNLOG=1),会突然出现一个相当大的力,而且产生负的接触能.
如果IGNORE设置为1,这样shootingnodelogic标志SNLOG没有作用.
相当于突然的穿透被留意到而且通过局部的调整接触厚度来补偿.
因此在任何分析的时间点,如果检测到突然的穿透,程序不会施加任何很大的力也不会移动任何节点.
但接触力会阻止进一步的穿透.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第17页13.
ContactSoft1接触Soft=1通过在*Contact卡的选项卡A中设置Soft=1所调用接触方程与缺省的罚函数接触方程的区别不如Soft=2接触算法大.
除了在确定接触刚度方面不一样外,Soft=1和Soft=0是差不多的.
Soft=1方法基于考虑稳定性的时间步来计算接触刚度.
换句话说,你可以把soft=1的接触看作是一组简单弹簧质量系统,其中每一个都有与实际仿真分析中用到的时间步相匹配的Courant时间步.
对软材料与刚硬材料的接触或两个接触面的网格密度不一致情况,Soft=1通常比Soft=0更有效.
当Soft=1时,采用通过soft=0&soft=1计算得到的最大的刚度值.
因此如果soft=0刚度更大,减小SOFSCL将没有效果.
K=max(SLSFAC*SFS*k0,SOFSCL*k1)其中:k是罚刚度SLSFAC是用户在*Control_contact中输入的选项SFS是用户在*Contact卡3中输入的选项SOFSCL是用户在*Contact选项卡A中输入的值k0是从材料体积模量和单元尺寸计算得到的刚度k1是从节点质量和计算时间步计算得到的刚度注意:对于双面接触,像*contact_automatic_surface_to_surface,在上面的方程中用SFM替换SFS(在*contact卡3中输入)LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第18页14.
LS-DYNA中的接触界面模拟1引言接触-碰撞问题属于最困难的非线性问题之一,因为在接触-碰撞问题中的响应是不平滑的.
当发生碰撞时,垂直于接触界面的速度是瞬时不连续的.
对于Coulcomb摩擦模型,当出现粘性滑移行为时,沿界面的切向速度也是不连续的.
接触-碰撞问题的这些特点给离散方程的时间积分带来明显的困难.
因此,方法和算法的适当选择对于数值分析的成功是至关重要的.
虽然通用商业程序LS-DYNA提供了大量的接触类型,可以对绝大多数接触界面进行合理的模拟,但用户在具体的工程问题中,面临接触类型的选择及棘手的接触参数控制等问题.
基于以上,本文对LS-DYNA中的接触-碰撞算法作了简要的阐述,对接触类型作了详尽的总结归纳,并对接触界面的模拟提出了一些建议.
2基本概念基本概念:"slave"、"master"、"segment".
在绝大多数的接触类型中,检查slavenodes是否与mastersegment产生相互作用(穿透或滑动,在TiedContacts中slave限定在主面上滑动).
因此从节点的连接方式(或从面的网格单元形式)一般并不太重要.
非对称接触算法中主、从定义的一般原则:粗网格表面定义为主面,细网格表面为从面;主、从面相关材料刚度相差悬殊,材料刚度大的一面为主面.
平直或凹面为主面,凸面为从面.
有一点值得注意的是,如有刚体包含在接触界面中,刚体的网格也必须适当,不可过粗.
3接触算法在LS-DYNA中有三种不同的算法处理碰撞、滑动接触界面,即:动态约束法(kinematicconstraintmethod)罚函数法(penaltymethod)分布参数法(distributedparametemethod)3.
1KinematicConstraintMethod采用碰撞和释放条件的节点约束法由Hughes等于1976年提出,同年被Hallquit首先应用在DYNA2|D中,后来扩展应用到DYNA3D中.
其基本原理是:在每一时间步Δt修正构形之前,搜索所有未与主面(mastersurface)接触的从节点(slavenode),看是否在此Δt内穿透了主面.
如是,则缩小Δt,使那些穿透主面的从节点都不贯穿主面,而使其正好到达主面.
在计算下一Δt之前,对所有已经与主面接触的从节点都施加约束条件,以保持从节点与主面接触而不贯穿.
此外还应检查那些和主面接触的从节点所属单元是否受到拉应力作用.
如受到拉应力,则施加释放条件,使从节点脱离主面.
这种算法存在的主要问题是:如果主面网格划分比从面细,某些主节点(masternode)可以豪无约束地穿过从面(slavesurface)(这是由于约束只施加于从节点上),形成所谓的"纽结"(Kink)现象.
当接触界面上的压力很大时,无论单元采用单点还是多点积分,这种现象都很容易发生.
当然,好的网格划分可能会减弱这种现象.
但是对于很多问题,初始构形上好的网格划分在迭代多次后可能会变得很糟糕,如爆炸气体在结构中的膨胀.
由于节点约束算法较为复杂,目前在LS-DYNA程序中仅用于固连与固连一断开类型的接触界面(统称固连界面),主要用来将结构网格的不协调两部分联结起来.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第19页3.
2DistributedParameterMethod分配参数法也是发展较早的一种接触界面算法,Wilkins在1964年将该算法成功地应用到HEMP程序中,Burton等在1982年将其应用于TENSOR分析程序中.
与节点约束法相比,这种算法具有较好的网格稳定性,因此被DYNA采用.
目前,在LS-DYNA程序中用来处理接触一滑动界面的问题.
该方法的基本原理是:将每一个正在接触的从单元(slaveelement)的一半质量分配到被接触的主面面积上,同时根据每个正在接触的从单元的内应力确定作用在接受质量分配的主面面积上的分布压力.
在完成质量和压力的分配后,修正主面的加速度.
然后对从节点的加速度和速度施加约束,以保证从节点在主面上滑动,不允许从节点穿透主表面,从而避免了反弹现象.
这种算法主要用来处理接触界面具有相对滑移而不可分开的问题.
因此,在结构计算中,该算法并没有太多的用处.
它最典型的应用是处理爆炸等问题,炸药爆炸产生的气体与被接触的结构之间只有相对滑动而没有分离.
3.
3PenaltyMethod该算法于1981年有Huag等人,1982年8月开始用于DYNA2D中[2].
现在,罚函数法已发展为一种非常用的接触界面算法,在数值计算中被广泛应用.
罚函数法的基本原理是:在每一个时间步首先检查各从节点是否穿透主面,如没有穿透不作任何处理.
如果穿透,则在该从节点与被穿透主面间引入一个较大的界面接触力,其大小与穿透深度、主面的刚度成正比.
这在物理上相当于在两者之间放置一法向弹簧,以限制从节点对主面的穿透.
接触力称为罚函数值.
"对称罚函数法"则是同时对每个主节点也作类似上述处理.
对称罚函数法由于具有对称性、动量守恒准确,不需要碰撞和释放条件,因此很少引起Hourglass效应,噪声小.
对称罚函数法在每一个时间步对从节点和主节点循环处理一遍,算法相同.
下面以从节点ns为例详细描述该算法的基本步骤:1搜索所有从动点,确定从动点是否穿透主面.
2如否,不作处理,搜索结束;如是,则在从节点与主面上的接触点间附加一法向接触力Fn.
3处理摩擦力.
4将接触力Fn和摩擦力投影到总体坐标,组集到总体载荷向量中.
4接触类型在具体介绍各种类型的接触前,先阐述几个基本的概念.
在壳单元中,自动接触通过法向投影中面的1/2"ContactThickness"来确定接触面.
这就是"shellthicknessoffsets".
接触厚度可以在接触的定义中明确指定.
如果接触厚度没有指定,则等于壳的厚度(在单面接触中,为壳厚度或单元边长的最小值).
相同的,在梁的接触中,接触面从梁的基线偏置梁截面等效半径距离.
因此,在有限元几何建模时,为考虑壳厚、梁截面尺寸必须在壳、梁的part间有适当的间隙,否则会有初始穿透现象发生(即发生不真实的接触现象).
虽然LS-DYNA可以通过移动穿透的从节点到主面上来消除初始穿透,但是并不是所有的初始穿透都能检查出.
DYNA中大多数的接触有一个"极限穿透深度",如侵彻超过这个深度则从节点被释放,接触力置为0.
这主要用在自动接触中,防止过大接触力的产生而引起数值不稳定性.
然而在有些情况下,因为这个阈值过早达到而使接触失效(常发生在非常薄的壳单元中).
此时应采取的措施是放大接触厚度因子或设置接触厚度为大于壳厚度的一个值,或者改变接触刚度的计算方法(如改为Soft=1).
LS-DYNA中的接触允许从节点与主段间压缩载荷的传递.
如接触摩擦激活,也允许切向载荷的传递.
Coulomb摩擦列式用来处理从静到动摩擦的转换,这种转换要求一个衰减系数、静摩擦系数大于动摩擦系数.
关于接触搜索方法,这里仅给出几个简单的要点,详细描述见TheoreticalManualofLS-DYNA.
DYNA中有两种搜索方法:IncrementalSearchTechnique与BucketSort.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第20页IncrementalSearchAlgorithmsGlobalBucketSort搜索方向仅在主段正方向从节点的穿透主面正、负方向检查穿透搜索步骤对每一个从节点的:找出最接近的主节点;搜索接近的主段(不止一个);搜索相邻的主段;局部利用IncrementalSearch确穿透检查;定最接近的主段;施加作用力.
穿透检查;施加作用力.
主面要求主面连续主面可以不连续特点简单、速度快非常有效,但耗时大所有的非自动LS-DYNA中的接触类型大体上可以分为五大类:One-WayContact(单向接触)Two-WayContact(双向接触)SingleContact(单面接触)EntityTiedContac(固-连接触)在以上接触类型中,前四种接触类型的接触算法均采用罚函数法.
固-连接触有的采用的罚函数法,有的采用动约束法,少部分采用分布参数法.
4.
1One-WayTreatmentofContactOne-Way、Two-Way是对接触搜索来讲的.
One-way仅检查从节点是否穿透主面,而不检查主节点.
在Two-WayContact中从节点与主节点是对称的,从节点与主节点都被检查是否穿透相应的主面或从面.
LS-DYNA中的_Node_To_Surface接触类型都属于单向接触,另外还有特别注明为单向接触的_Surface_To_Surface接触类型:z*Contact_Nodes_To_surfacez*Contact_Automatic_Nodes_To_Surfacez*Contact_Froming_Nodes_To_Surface(自动接触类型、主要用于金属拉压成形)z*Contact_Constraint_Nodes_To_Surface(现已很少用)z*Contact_Eroding_Nodes_To_Surfacez*Contact_One_Way_Surface_To_SurfacezContact_One_Way_Automatic_Surface_To_Surface由于在单向接触中,仅有从节点被检查是否穿透主面,而不考虑主节点,因此在使用时必须注意,应保证在接触过程中主节点不会穿过从面.
同样的原因,单向接触要比双向接触运行速度快得多,因此仍被广泛应用.
在以下情况中使用单向接触是合适的:主面是刚体相对细的网格(从)与相对平滑、粗的网格(主)接触beam_to_surface、Shelledge_to_surface接触.
beamnode、Shelledgenode作从点.
在接触分析中,由于问题的复杂性,判断接触发生的方向有时是很困难的,因此分析中应尽量使用自动接触(不需要人工干预接触方向).
但当面的方向在整个分析过程中都能确定的情况下,下面的非自动接触类型是非常有效的:*Contact_Nodes_To_Surface(5)*Contact_One_Way_Surface_To_Surface(10)*Contact_Constraint_Nodes_To_Surface(18)*Contact_Eroding_Nodes_To_Surface(16)LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第21页4.
2Two-WayTreatmentofContact主、从面的定义与算法处理上是完全对称的.
因此主面、从面可以随意定义.
计算资源大约是单向的2倍.
LS-DYNA中绝大多数_Surface_To_Surface接触都是双向接触类型.
双向接触除对主节点的搜索外,其它方面同单向接触是完全一样的.
与前述接触类型5、18、16相对应的双向接触为:*Contact_Surface_To_Surface(3)*Contact_Constraint_Surfaces_To_Surface(17)*Contact_Eroding_Surface_To_Surface(14)在CrashAnalysis中,*Contact_Automatic_Surface_To_Surface(a3)推荐使用.
在金属的拉压成形分析中推荐使用*Contact_Froming_Nodes_To_Surface.
4.
3SingleSurface单面接触是LS-Dyna中应用最为广泛的接触类型,尤其在Crashworthiness应用中.
在这中类型中,从面一般定义为Part或PartSetID.
各Part间及自身Part间的接触都考虑.
如果建模精确,该接触是可信、精确的.
在单面接触中,壳厚偏置总是考虑的,因此建模时不能有初始穿透存在.
单面接触有:*Contact_Single_Surface(4,不推荐使用)*Contact_Automatic_Single_Surface(推荐)*Contact_Automatic_General*Contact_General_Interior*Contact_Airbag_Single_Surface对于CrashAnalysis,推荐使用*Contact_Automatic_Single_Surface(13).
这个接触类型其性能随DYNA版本的提高不断改善.
4.
4TiedContact(TranslationalDOFonly,NoFailure,NoOffset)固-连接触用来将从节点约束、限定在主面上.
这种接触类型一般是非对称的,因此定义主、从是要符合§2中描述的一般规则.
在这种类型的接触中,主、从接触面最好不要以PartsID形式输入,应采用node/segment的形式.
固-连接触类型丰富,采用的接触算法也不一致,下面分别介绍.
TranslationalDOFonly,NoFailure,NoOffset这种接触仅约束从节点的平动自由度,且不考虑接触的失效,不允许从节点的偏置.
如果从节点与对应的主段间有微小的距离存在,则采用正交投影的方法将从节点移动到主面上.
因此,初始几何构形可能有微小的改变.
这种类型接触采用动态约束算法,因此不能将刚体约束到可变形体或刚体.
如下两个命令是常用的固连接触.
这两种接触在数值处理上是完全一样的,所不同的仅是输入数据格式.
*Contact_Tied_Nodes_To_Surface(6)*Contact_Tied_Surface_To_Surface(2)TranslationalDOFonly,NoFailure,WithOffset这种接触采用罚函数算法,允许从节点与主面间偏移(主、从面间存在微小的距离)存在,可以用于刚体相应的约束.
与上述接触类型2、6对应的为*Contact_Tied_Nodes_To_Surface_OFFSET(O6)*Contact_Tied_Surface_To_Surface_OFFSET(O2)由于从节点的偏置,可能会引起附加的动量矩.
但在这种类型的接触中,不考虑偏置引起的动量矩.
因此,主、从面必须相当的接近.
TranslationalDOF&RotationalDOF,WithFailure,NoOffset采用动态约束算法.
TranslationalDOF&RotationalDOF,WithFailure,WithOffset罚函数法.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第22页TranslationalDOFOnly,WithFailure,WithOffset动态约束算法.
5接触刚度的计算在基于罚函数算法的接触类型中,目前LS-DYNA有两种计算主、从面间接触刚度的方法.
5.
1Penalty-baseApproach(SOFT=0)该算法是LS-DYNA计算接触刚度的缺省方法.
它利用接触段的尺寸与其材料特性来确定接触刚度.
当两个接触面的材料刚度参数相差不大时,该方法是很有效的.
但当两个接触面的材料刚度相差很大时,由于接触刚度采用主、从面中较小的刚度,而使接触失效.
对于Crash分析,除非先验证明没有问题,否则一般不使用SOFT=0.
5.
2SoftConstraint-basedApproach(SOFT=1&2)计算接触刚度时,综合考虑了发生接触的节点之质量与整体时间步长,以保证接触的稳定性.
这样的处理,对于材料性质相差悬殊的接触问题是非常合适的.
Soft=1与Soft=0算法除刚度计算外,其它考虑是完全一致的.
当Soft=1时,采用下式计算接触刚度k(单向接触):k=max(SLSFAC*SFS*k0,SOFSCL*k1)其中:k0~根据材料弹性模量与单元尺寸确定的接触刚度;k1~根据节点之质量与整体时间步长来确定接触刚度.
对于Two-Way型接触,用SFM代替上式中的SFS.
5.
3Segment-basedContactVS.
StandardContact与Soft=0、1(以下简称"标准算法")不同,Soft=2是一种基于段(Segmentbased)的接触算法.
在标准算法中,检查从点穿透主段与否而施加罚力与从点及相应的主点;而在段接触算法中,直接检查段是否发生相互穿透而施加罚力与相应段的节点.
6接触控制参数LS-DYNA提供了多个与接触相关的控制参数.
根据不同接触问题的具体特点,设置不同的控制参数,对提高"接触模型"的精确性是非常必要的.
LS-DYNA中的接触控制参数可以在*Control_Contact、*Contact或*Part_Contact中设置,而有些参数也可以同时在多个命令中设定.
如一个参数在多个命令中设置,则这样的设置有一定的优先次序.
*Control_Contact对整个模型中的接触提供一种"全局性"的"缺省"参数设置;*Contact对每个具体的接触提供"局部"的参数设置,优先权较高;*Part_Contact则为某个具体的Part涉及的接触提供最高级别的参数控制.
6.
1Thicknessoffset:Automatic,SLTHK(Card1,*Control_Contact,OptionCardA)LS-DYNA中非自动接触类型:9*CONTACT_SURFACE_TO_SURFACE9*CONTACT_NODES_TO_SURFACE9*CONTACT_ONE_WAY_SURFACE_TO_SURFACE利用参数SHLTHK确定是否考虑"厚度偏置"(见下图),该参数可以在*CONTROL_CONTACT中全局定义,也可以在OptionalCardB中局部定义.
如果SHLTHK=0,不考虑厚度偏置,采用incrementalsearch方法来确定从节点最接近的主段;如SHLTHK=1,考虑变形体的厚度偏置,但不考虑刚体厚度偏置;如SHLTHK=2,变形体、刚体的厚度偏置都考虑.
如SHLTHK为1或2,程序采用globalbucketsearch来确定接触对.
接触建立以后,采用incrementalsearching来跟踪从节点在主面上的位置.
采用globalbucketsearching的优点是主、从面可以不连续(这对incrementalsearch是不可能的).
在非自动接触类型中,接触段的法向方向(符合右手法则,指向接触面)是非常重要的,必须保证所有接触段的法向一致指向接触面,这就是所谓的"orientedcontact".
一个简单的方向自动定位方法是激活LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第23页*CONTROL_CONTACT中的参数ORIEN(必须在两个中面间有一定的距离).
自动接触和单面接触总是考虑"壳厚偏置".
在这些类型的接触中,采用整体块搜索和局部增量搜索方法确定接触对.
在CrashAnalysis中,一般建议使用自动类型的接触,因为自动接触对于接触面的连续性、方向等基本没有限制,具有教强的适应能力.
6.
2ContactSlidingFriction:FS&FD(Card2)LS-DYNA中的摩擦采用Coulcomb摩擦列式与等效弹塑性弹簧模型.
摩擦通过设置*Contact或*Part_Contact中非零的静(FS)、动(FD)系数来激活.
如静、动摩擦系数不同,则FD应小于FS,同时必须指定非零的衰减系数DC.
对于伴有数值噪声的问题(如CrashAnalysis),FS、FD通常设为相同的值,以避免额外噪声产生.
为限制过大、不真实的摩擦力产生,通常设置3yVCσ=.
不同类型的问题对摩擦系数的敏感性是不同的,有时可能存在很大的差异.
在具体问题分析是,可以通过极限分析(设置FS和FD的上、下限)的方法确定摩擦的敏感性.
6.
3PenaltyScaleFactors:SFS&SFM(Card3)罚因子(SFS、SFM)用来增大或减小接触刚度.
在Soft=0、2时,真正的主、从罚因子是SFS、SFM分别与SLSFAC(*Control_Contact中定义)的乘积.
对于材料刚度相当、网格尺寸相差不大的两面间的接触问题,SFS、SFM、SLSFAC的缺省设置是可行的.
但相接触的两个面的材料、网格相差很远时,对于Soft=0的接触算法可能存在问题,此时一个简单的办法就是设置Soft=1,而不必考虑(或试算)罚因子的选择.
6.
4ContactThickness:SST&MST(Card3)SST、MST可以直接指定期望的"接触厚度".
如SST=MST=0(缺省值),则接触厚度等于*Section_Shell中定义的单元厚度.
有时通过设置非零的SST、MST值来消除"初始穿透"(尽量避免这样做).
SST和MST一般不应小于0.
6~0.
7.
ContactThicknessScaling(SFST&SFMT)同SST、MST作用相同.
6.
5ViscousDamping:VDC(Card2)粘性接触阻尼用来降低(高速)碰撞过程中接触力的高频振荡.
对于存在软材料(如泡沫材料)的接触问题,VDC设为40~60(临界阻尼的40~60%),通常能提高模型的稳定性.
对于金属间的碰撞接触问题,VDC一般可设为20.
6.
6Bucket-SortFrequency:BSORT(OptionalCardA,*Contact)&NSBCS(Card2,*Control_Contact)BucketSort是一种非常有效的接触搜索算法.
如果考虑"厚度偏置",则在所有的接触类型(自动、非自动接触)中,对于任一从节点均使用BucketSort方法搜索可能与之接触的主段.
由于接触搜索是"接触模拟"中非常耗时的一个步骤,因此应尽量减少搜索的次数.
BSORT用来指定两次搜索间的迭代时间步数,BucketSort的间隔一般为10~100(与具体的接触类型有关).
对于不连续面间的接触、高速碰撞等问题,应增加搜索的次数,即减小BSORT(或BSBCS),但一般不应小于10.
在这些问题中,如搜索间隔过大,一些从点就会在接触处理中被漏掉.
但对于相对平滑的面间接触问题,可以适当增加BSORT或NSBCS.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第24页6.
7MaximumPenetration:PENMAX(OptionalcardB,*ControlContact)&XPENE(Card2,*Control_Contact)为避免由于从节点穿透深度过大(罚力与穿透深度成正比)而引起的数值不稳定,当从节点穿透到一定的深度(MaximumPenetration),该节点从接触中自动释放(但依然参与其他的计算).
在对壳元的穿透中,为防止当从节点穿透壳的中面而引起的接触力方向的突然翻转,"壳厚偏置"的考虑也是非常必要的.
在非自动接触中,如SHLTHK=0,则缺省的最大穿透深度为1.
0e20,也就是说不考虑从节点的释放.
如SHLTHK=1or2,则参数XPENE确定节点释放准则:zMaxDistance(Solids)=XPENE(default=4.
0)*(thicknessofthesolidelement),SHLTHK=1zMaxDistance(Solids)=0.
05*(thicknessofthesolidelement),SHLTHK=2zMaxDistance(Shells)=XPENE(default=4.
0)*(thicknessoftheshellelement),SHLTHK=1zMaxDistance(Shells)=0.
05*(minimumdiagonallength),SHLTHK=2在自动接触、单面接触中(Automatic_General除外),最大穿透深度由PENMAX(缺省值为0.
4)确定:zMaxDistance=PENMAX*(thicknessofthesolid)zMaxDistance=PENMAX*(slavethickness+masterthickness)对于Automatic_General接触,PENMAX的缺省值为200(几乎不考虑节点的释放).
对于控制最大穿透深度的参数一般不要改动(使用缺省设置).
如果节点穿透过大而需要释放,可以采用增大接触刚度、改变罚函数算法(SOFT),或增加接触厚度等方法来实现.
7接触输出在LS-DYNA中,最常用的接触输出文件是RCFORC,它包含主、从面每一个节点接触力(GlobalCartesianCoordinateSystem)的ASCII文件.
为输出RCFORC必须在k文件中包含*Database_FCFORC,同时必须激活接触控制中的参数SPR、MPR(Card1).
注意:对于单面接触,RCFORC无效.
此时要输出接触节点力,必须通过*Contact_Force_Transducer_Penalty定义力传感器(forcetransducers).
力传感器仅用来输出接触力,对数值分析结果毫无影响.
接触面的能量通过*Database_Sleout输出到ASCII文件SLEOUT中.
该文件对于分析每个接触定义的可靠性是很有帮助的.
在某些情况下,有时需要接触界面的可视化(如应力云图等),这时必须通过以下控制输出二进制的接触界面文件:1)*Database_Binary_Intfor;2)设置接触面的输出标志SPR、MPR;3)在执行计算任务时,包含选项"s=filename".
8结束语LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第25页15.
LS-DYNA中夹层板(sandwich)的模拟在LS-DYNA中夹层板(sandwichplate)通常有三种方式来模拟1.
采用几层体单元来模拟,当板比较薄,沿厚度方向单元尺寸比较小,时间步很难提高,2.
中间的夹层材料用体单元模拟,外面的包层材料用shell单元模拟,在体单元和壳单元之间用Contact_tied(_offset)卡建立粘接关系,3.
用一层壳单元来模拟.
对材料*Mat_composite_damage和*mat_enhanced_composite_damage打开层板理论选项(*control_shell卡LAMSHT项设置为1)对最后一种方式,需要定义*intergration_shell卡片.
例如下面这个层合板:*INTEGRATION_SHELL1,8,0-.
9722,.
02778,1-.
9167,.
02778,1-.
6667,.
22222,2-.
2222,.
22222,2.
2222,.
22222,2.
6667,.
22222,2.
9167,.
02778,3.
9722,.
02778,3*PARTmaterial1111*PARTmaterial2112*PARTmaterial3111*SECTION_SHELL120.
00000008.
00000000.
00000-1.
000000018.
00000018.
00000018.
00000018.
0000000.
0000000*mat_composite_damage11,2.
7e-6,73.
4,73.
4,73.
4,0.
32,0.
32,0.
3227.
8,27.
8,27.
81.
e9,1.
e9,1.
e9,1.
e9*mat_composite_damage12,6.
3e-7,0.
286,0.
286,0.
286,.
3,.
3,.
30.
11,0.
11,0.
111.
e9,1.
e9,1.
e9,1.
e9其中*intergration_shell卡片中第一项S表示板厚度方向积分点的规一化坐标值,从-1到1.
积分点规一化坐标=(积分点全局Z坐标-中面Z坐标)/(thickness/2)第二项WF(Weightingfactors)为与积分点关联的厚度与整个板厚度的比值,WF和应该为1.
第三项指向各层对应的partID号,这个partID号仅用来指定这一层的材料.
完整的一个实例如下:shell8lam22.
rar14.
LS-DYNA隐式-显式自动转换例子在LS-DYNA中可以通过卡片*Control_implicit_general的选项imflag指定为一个负值(这个负值的绝对值为一条*define_curveID),来自动实现在指定的时间内作隐式或显式分析.
这里的例子为一个一端固定平板,另一端施加从0-0.
01秒从0增加到10.
0的力,之后释放.
曲线ID3指定从0-0.
01秒内用隐式计算,增量步为20步.
之后自动转为显式继续求解.
下载inputdeck-》implicit_explicit_switchk.
txt注意请用5434,5434aor971版本运行上面的k文件,不要用6763单精度或者双精度版本(这两个版本在implicit求解方面似乎存在bug,隐式分析部分无法完成)LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第26页16.
怎样进行二次开发1.
二次开发的过程和需要的文件在回答这个问题前,给大家一点背景知识(针对vwindows平台):1、DYNA只是一个求解器(一个固定的可执行的exe文件),它可以读取K文件,并进行计算,最终输出结果.
2、DYNA开发过程:编译连接DYNA源代码---->DYNA.
OBJ文件或DYNA.
lib————>DYNA.
exe如果进行二次开发,实际上就是要修改原有的exe文件,在其中加入自己的程序(子材料什么的).
而不管你是"一般用户"还是正版用户,都不太可能直接修改exe文件(除非你特别特别牛,能直接编写机器语言(比汇编语言还低一级别),开个玩笑,基本不可能实现).
所以目前DYNA开发商,给正版用户提供了下面的方式进行二次开发:编译连接自定义程序代码(F文件)————>OBJ文件+DYNA.
lib————>新的DYNA.
exe文件(其中包含自定义的部分)从上面可以看出二次开发实际上是重新生成了一个全新的可执行程序,这一过程需要F文件、DYNA.
lib和其它资源文件,而不需要原来的DYNA.
EXE文件,所以不管你是正版还是别的,如果没有DYNA.
lib文件等都是不可能进行二次开发的.
而"一般用户"和正版用户的区别在于:1、正版用户可以直接要求开发商提供F文件、DYNA.
lib和其它资源文件,而一般用户只有通过其它渠道.
2、将K文件提交DYNA.
EXE计算时,正版用户由于有license,可以直接通过检查并进行计算;而"一般用户"目前基本上是通过修改DYNA.
exe文件.
跳过检查,最终实现计算的.
一般用户"要实现二次开发就必需首先有F文件、DYNA.
lib和其它资源文件;通过编译连接后生成新的EXE文件;再修改EXE文件,跳过license检查(这个网上有说明具体过程,所以关键是要得到那几个二次开发文件).
"一般用户"完全有可能进行二次开发,关键是要通过其它渠道得到正版用户才有的F文件、DYNA.
lib和其它资源文件2.
子程序的编写见白金泽和赵海鸥的书.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第27页17.
如何实现程序的连续计算通过批处理设置连续计算可以省下很多时间.
首先创建一个记事本,将扩展名改成.
bat比如run.
bat.
编辑批处理文件,写下面的语句:f:cdF:\111\d:\lsdyna\program\ls971_s_7600_win32.
exei=*.
kmemory=内存大小cd.
.
第一步:f:这句很重要(假如你的k文件目录在f盘),因为批处理就是在命令提示符(DOS)下运行的.
在某个盘(除了c盘)建立一个目录,直接用cd命令是进入不到这个相应的目录的.
比如在f盘建立一个文件夹,命名为111.
用cdf:\111这句命令是进不到这个目录的,还是在c:.
必须先输入f:在输入cd111才能进入目录.
第二步:找到求解器所在的目录和命名,再设置好memory大小(也可以在k文件第一行*keyword后面设置),CPU个数即可.
D:\lsdyna\program\ls971_s_7600_win32.
exei=*.
kmemory=内存大小C:\ProgramFiles\LSDYNA\program\ls970_s_6763.
exei=*.
kmemory=内存大小上面是ls971求解器的,下面是ansys自带的970求解器的.
第三步:复制第二步的内容,有几个k文件就复制几次,改掉相应的路径和k文件名即可.
若想求解完毕后自动关机,还可以在最后一行写上"shutdown-s".
第四步:保存批处理文件,双击求解.
最后在提醒一点,每个k文件最好都提前试一下,运行无误后写入让它自动求解.
免得设置好高高兴兴走了,回来一看有问题,那就白高兴了,哈哈!
祝大家学业有成,工作顺利!
Timothy1985LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第28页18.
如何得知LS-DYNA为一套完整的软件解答:LS-DYNA一直以来都是独立软件,原因很简单:(1)美国原厂LivermoreSoftwareTechnologyCorp.
(简称LSTC)授权「势流科技」为台湾地区LS-DYNA软件的独家代理商.
(2)LS-DYNA拥有LSTC官方的前、后处理程序:FEMB及LS-Pre/Post,无需另外增购.
LS-DYNA强调「OneModel、OneCode、Multi-results」,单一核心程序即可进行Implicit、Explicit及混合Implicit/Explicit的求解方式,可以为您解决静力、准静力、掉落、冲击、模态分析、振动、预应力/预应变、流固耦合、流场、热应力.
.
.
.
等等的力学分析,使用者不必再额外添购模块程序,徒增预算的浪费与软件学习的负担.
19.
LS-DYNA与市面上其它的前处理软件兼容性如何解答:由于LS-DYNA是全球使用率最高的结构分析软件,因此其它的前处理软件与LS-DYNA是完全兼容的.
在此要强调的是:LS-DYNA的官方前处理程序为FEMB,因为FEMB是专门为LS-DYNA量身订作的前处理程序,有许多设定条件及架构逻辑是其它前处理软件所难望其项背的,为了避免在学习LS-DYNA的过程及操作上产生困扰,强烈建议使用者采用原厂出品的FEMB来做为LS-DYNA的前处理工具,使用者必定更能体会LS-DYNA直觉式的设定与强大的分析能力.
20.
LS-DYNA似乎很重视「ContactAlgorithm」,这是为什幺解答:是的,LS-DYNA很早以前就已经发展「接触算法」,这是因为基础力学所分析的对像均只考虑「力的受体」,故输入条件皆为外力量值.
然而在真实情况下,物体受力通常是因为与其它的物体发生「接触」(Contact)才受力,此时外力量值是无法预期的,应该输入的条件往往都是几何上的接触条件.
因为有完备的接触力学演算方式,LS-DYNA才得以忠实的仿真现实环境的复杂结构行为.
21.
如果要利用LS-DYNA进行MPP(平行运算)的计算,硬件配备及操作系统有无特殊需求解答:不论是PCcluster、工作站及一般的PC环境,都适合执行LS-DYNA的MPP平行运算功能,一般我们还是会建议要用来执行平行运算的计算机群组,彼此的等级宜尽量一致;操作系统方面并无特别需求,以一般的windows2000、LINUX或是UNIX皆可执行.
国外已有很多厂商利用非办公时间,将办公室内的计算机串连在一起,配合LS-DYNA来分析问题,宛如一部超级计算机,不仅可以有效提升研发的竞争力,同时亦可弹性地运用计算机资源,一举数得.
22.
FEMB能够自动产生有限元素网格吗解答:FEMB当然可以自动产生有限元素网格,使用者再也不必费心在每个几何边界上指定结点数量,仅需要输入元素尺寸的参考值,FEMB便会依此产生网格.
当然,如果能够先在CAD进行合理的简化,auto-mesh得到的元素品质会更好.
在国外CAE领域有两句名言:「WellDoneisQuicklyDone」、「QualityMeshleadstoQualityAnalysis」,因此势流科技建议使用者应该多花点心思在前处理上,这对后续的分析工作有着莫大的帮助.
23.
实际的产品CAD图文件,有许多复杂的几何造型,这些feature是否都该纳入分析考量解答:这样的问题事实上是属于「EngineerJudgement」的层次,而非软件层次的问题.
使用者当然可以将完全按照CAD几何来建立有限元素的模型,然而任何一位有经验的工程师一定会先将不必要的几何造型予以忽略,以提升分析效益,并将分析课题明确化.
以Block为例,有些工程师会以Solidelement来建构完整的模型,有些则采用Shellelement,甚至有些人会采用1Delement.
.
.
,这些方式都无伤大雅,重点是工程师如何从分析结果获得充份的信息.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第29页势流科技除了提供软件分析技术之外,同时更致力于协助业界将CAE导入于生产流程的正确位阶,培养正确的CAE分析观念,以提升产品研发品质,缩短生产流程.
24.
LS-DYNA是否只能进行壳元素的分析解答:LS-DYNA不只能分析实体元素(Tetra、Hexa)、薄/厚壳元素(Quad.
、Tri),同时还有梁元素、质点、Spring/Damper、spotweld等元素型态,另外还有SPHelement.
之所以建议台湾电子产业尽量以壳元素来仿真壳件产品,是由于电子产品的组件当中,板壳类占了90%以上,基于结构行为、准确性及分析效能的考量.
除非使用者有特殊的分析需求,否则并不建议采用实体元素来仿真板壳类型的产品组件.
我们藉由右边的简单案例,以了解不同元素型态在仿真壳件产品的结果表现:(1)实体元素在板壳厚度方向上至少需要分割三层以上,其位移变形量才能收敛至一定值;而壳元素及厚壳元素仅需一层,在位移量上即可达到理想精度.
(2)若是发生降伏产生大变形,在厚度方向上即使分割8层以上的实体元素,其掌握材料塑性应变的收敛速度仍然相当缓慢,而付出的时间代价为shellelement的25倍以上.
厚壳元素的运算速度及精度与三层实体元素相当,然而就塑性应变的精度而言,厚壳元素的表现比实体元素要好太多.
25.
在LS-DYNA里该如何强制刚体运动解答:CAE分析常需要藉由刚体运动,以简化复杂物理行为.
应用技巧有二:给定初始速度或是强制运动轨迹(prescribedmotion).
在LS-DYNA中,预设的刚体旋转轴位置为刚体的质心,若要强制刚体绕特定轴旋转,只要利用FEMB在PART选单里设定INERTIA相关的数值即可.
其中的XC、YC、ZC为旋转中心所在位置,即旋转轴通过的地方(内定值为刚体质心).
TM为TranslationalMass,一般状况输入质量即可.
IXX~IZZ为此刚体质心位置相对于旋转中心的惯性矩,一般CAD软件应已具备此项计算功能,因此使用者无须再额外进行人工运算.
以Pro/E的操作为例:1.
确认图文件的单位(公制、英制).
2.
于旋转中心建立坐标系.
3.
选择「Analysis」「ModelAnalysis」「ModelMassProperties」.
4.
选择步骤2.
产生之坐标系.
Pro/E计算结果会有二组InertiaTensor,其中「InertiaatCenterofGravitywithrespecttoXXXcoordinateFrame」才是我们要的,再将相对应的数值填入IXX~IZZ即可.
26.
在LS-POST如何显示应变分布图解答:LS-DYNA内定的输出控制并不包含应变值.
如果您希望能有应变值的输出的话,请在FEMB里依照下列步骤完成设定即可.
步骤1:点选「DynaMisc」,进入「DATABASE」选项.
步骤2:请在「Extent」选项中新增一个「Binary」的输出控制.
步骤3:设定STRFLG=1即可.
步骤4:分析完成后,使用者可直接在LS-Pre/Post中,点选「Fcomp」「Strain」观察计算后的应变分布图.
LS-Pre/Post还有其它关于应变的显示按键,如:「Infin」、「Green」及「Almans」等,基本上这些应变量是根据d3plot档案中所记录的位置资料,经过微分计算而得到.
如果点位取得够密的话,理论上会与板壳中性平面的应变值相当接近.
27.
强制结构体绕着特定轴旋转的时候,为何会出现奇怪的变形解答:在定义物体旋转自由度的时候,常会以*Initial_velocity_generation来设定角速度,然而这样的方式会将离心的效果纳入考量,特别是当角速度值很大的时候,,旋转物体会由旋转中心沿半径方向膨胀变形.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第30页因此,如果不希望有离心膨胀的效果的话,可以利用*Load_body_rz(orrx、ry),并配合*Define_curve给定角速度对应时间的关系曲线.
几个要注意的地方如下:步骤1:起始速度的设定选项,请设定Phase=1,以进行dynamicrelaxation.
步骤2:定义角速度的曲线中,必须开启stressinitialization的功能(SIDR=1).
步骤3:在controlcard的部份须开启accuracy及dynamicrelaxation的功能.
步骤4:参考范例请至下载区download.
28.
LS-DYNA该如何进行振动模态的分析设定解答:模态分析即为特征值分析,须采取Implicit的方式求解,您只需要在前处理程序FEMB里依照下列方骤进行设定即可.
步骤1:点选「DynaMisc」,进入「Control」选项.
步骤2:将一般Implicit功能选项开启,其中一项即为「IMP.
EIGENVALUE」.
步骤3:最基本的设定,您只需要设定模态个数(NEIG)的值即可.
步骤4:LS-DYNA针对特征值,还有许多非常实用的输出控制,详细用法请参考LS-DYNAKeywordUser'sManualVol.
I,或是直接与势流科技的工程师联络.
29.
究竟采用Implicit及Explicit的分析方式会有多大的差异性解答:依照静力的观念,其破坏点必发生于试片中点;然而若是以瞬间拉扯试片的话,因为应力传递速度与路径的关系,其破坏点未必会发生在中点,一般的静力软件必定不能满足真实的结构行为.
事实上,ASTM对于材料性质试验,都有相关的规范以定义外力加载的速率,也就是因为这个原因.
30.
要如何设定才能让LS-DYNA自动进行Implicit/Explicit的切换解答:在LS-DYNA960的版本即有Implicit/Explicit切换的功能,使用者可在FEMBDynaMiscIMP.
General的第一个字段输入「-100」,表示Implicit/Explicit的切换时间是依据#100的loadcurve来定义的.
而在LS-DYNA970最新版本更提供了自动切换的功能,做法如下:步骤1:进入FEMBDynaMiscCONTROL步骤2:进入IMP.
AUTO选项,给定DTEXP的值,此值即为explicit分析所需的计算步长.
步骤3:进入IMP.
General选项,设定IMFLAG=4or5即可.
31.
LS-DYNA是唯一能够将Implicit与Explicitcode整合在一起的软件,这对实际的分析工作有何帮助解答:LS-DYNA将Implicit及Explicit两种不同的求解方式整合在一起,对于初学者而言可以降低学习的困扰、缩短软件学习的时程、省去面临不同软件的适应问题;对进阶的使用者而言,可以从事更详实的力学分析.
最实际的效益在于:(1)预变形及预应力的仿真:真实的结构行为通常都是以连续作动的方式在进行,因为第一阶段的力学行为而衍生出后续阶段的变形.
传统的作法,使用者必须凭借个人经验或是以简化的方式在结构体上施予一假设的预应力(预变形)量,然而所施加的量值却未必能符合真实的数值.
LS-DYNA整合Implicit及Explicit解法,可以忠实的仿真各种结构反应.
以卡勾卡合的动作机制而言,在第一阶段可利用Implicit解法指定卡勾运动方向,接着可马上切换至Explicit解法,仿真卡合瞬间的反应.
(2)相信以往从事Implicit分析的工程师都有调整收敛参数的不愉快经验,而LS-DYNA970可以视运算过程收敛性的难易程度,自动将solver由Implicit切换到Explicit,以克服收敛性的问题,这个功能绝对是非常吸引人的!
!
32.
请问以Explicit方式进行分析,影响运算时间的因素有哪些解答:采用Explicit方式最直接影响运算时间的因素是「计算步长」(timestep)的大小,LS-DYNA决定计LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第31页算步长的方式为:另外还有许多因素会影响运算的时间:(1)分析历时(TerminationTime),分析历时较短者,可缩减运算时间.
(2)接触算法:LS-DYNA采用最有效率的接触演算理论.
(3)机器等级:计算机等级高者,运算时间较短.
(4)元素理论:采用较多高斯积分点的元素,其运算时间会增加.
(5)材料模式之复简:材料性质直接影响计算步长.
(6)模型的复简:适当的简化CADmodel可大幅提升运算品质.
所以当您在评估软件时,若有任何关于LS-DYNA计算速度的疑问,都欢迎与势流科技联络,我们将提供最正确的CAE信息与您分享.
33.
势流科技您好,请问弹簧/阻尼系统该采用什幺元素来仿真比较恰当解答:若您是古早版本(LS-DYNA940)的使用者,可以直接选用discreteelement来仿真弹簧或是阻尼元素.
不过LSTC原厂建议采用最新版本LS-DYNA970所提供的BeamElement(#6)来仿真Spring/Damper,配合*MAT(#66、#67)可一次定义多轴的劲度与阻尼值,并且不用额外定义LumpMass.
例如右图即是采用单一BeamElement(#6)来仿真单自由度(SDOF)系统在NoDamped、UnderDamped、CriticleDamped及OverDamped状况下的振动行为.
34.
请问安装LS-DYNA(MPP)版本时,操作系统需有哪些额外的设定解答:若你所采用的操作系统是Window2000,可以请贵单位的网管人员协助下列设定:(1)确定参与平行运算的各机器皆处于同一网域内.
(2)各机器须提供公用帐号及密码给予MPP的使用者登入.
(3)于每一台机器上安装最新版本的MPICH程序(Mpich-1.
2.
5),利用所附的mpiregister程序设定MPP的使用者帐号及密码.
若您所用的系统是UNIX/Linux,亦请贵单位的网管人员协助以下的设定:(1)确定所有参与MPP计算的机器皆位于同一网域内.
(2)建立NFSserver让cluster内的所有机器将LS-DYNA的所在目录挂载进来.
(3)进行相关安全设定(/etc/hosts、/etc/hosts.
equiv、.
rhosts.
.
.
.
etc),让MPP的使用帐号可以利用rsh或是ssh的方式在Cluster的机器间互相登入而无须输入密码.
(4)安装MPICH或是LAM最新版本(LAM-6.
5.
6)35.
在Linux系统下如何指定MPP的computingnode解答:若你是采用MPICH的话,请修改/usr/local/mpich/share/machines.
LINUX的档案即可,参考格式如下:node1:2node2:2若你是采用LAM的话,请先编辑hostfile,参考格式如下:node1cpu=2node2cpu=2执行运算前请别忘了要先激活LAM,请执行lamboot–vhostfile;若是激活成功,屏幕上会显示参与MPP计算的Nodenumber.
36.
请问在执行MPP运算时,有没有类似PC版本Ctrl+C的功能以进行不同的switch解答:当然是有的.
你可以利用另一个terminal窗口在工作目录产生名为「d3kil」的档案,将正确的switch(如sw2)填写在第一行,LS-DYNA程序会周期性的读取d3kil档案,并进行相关的switchreport.
若是Linux/UNIX系统的话,还可以有另一种方式:请另外开启新的terminal,执行sourcebg_switch,LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第32页此时正在运算的程序会先暂停,并且等待你的switch讯号.
37.
如果我想要了解特定截面的力量分布情形,该设定哪些项目呢解答:您所想获得的信息应该是属于Sectionforce的输出结果.
其做法如下:(1)进入Set选单,将特定截面上的所有结点设成同一个群组(NodeSet).
(2)将这些结点所在的元素(在截面的同一侧)也设成同一个群组(ElementSet).
(3)进入DynaMisc选单点选DatabaseCrosssection建立Crosssectionset,将步骤(1)(2)所设定的结点及元素群组点选进来即可.
(4)最后再到ACSII选单定义secforc的记录输出时距即可至于资料的读取方式与一般利用LS-Pre/Post读取ACSIIdata的作法一样.
一般常见的应用在于获得支承的扭转性质(如LCD、铰接托架等结构).
38.
请问LS-DYNA在进行Implicit分析时,有没有什幺该注意的地方解答:您这个问题好象有点「蛮广泛」.
.
.
在这边先跟您说明一下大原则:Implicitmethod分析主要就是在求解矩阵式,给定的力学条件不外乎是「ForcePrescribed」与「DisplacementPrescribed」,一般最常犯的毛病就是边界条件给定不足,造成model发生「刚体运动」的情形,您可以将这种情形想象成物体浮在空中一样,宜尽量避免这种情形产生.
在此提供您一个判断的方式:您可以先利用LS-DYNAeigenvalueanalysis的功能先进行模态分析,藉后处理程序检查有无刚体模态,以判断Implicitmodel的束制条件是否充足.
39.
请教一下,LS-DYNA有没有什幺功能可以让上/下表面的网格移到中性平面来解答:有的,请依照下列步骤进行:(1)利用Check的功能,将元素法向量调整一致.
(2)进入Element选单,点选transform.
(3)接着会出现向右图一样的窗口,请选择NORMALOFFSET.
(4)输入半个壳厚即可.
Normaloffset与直接平移的不同是在于Normaloffset会依照壳厚修正曲率,能够让网格正确的移至中性平面的位置.
所以你大可放心在CAD的外表面产生igessurface,再利用Normaloffset的功能让网格产生在中性平面上,以加快建构网格的时间.
40.
请问我该如何在后处理程序LS-Pre/Post里头看到壳元素的厚度解答:你好,请依照下列步骤进行:(1)进入Appear选单.
(2)勾选Thickness的方框.
(3)选择Allvis(Allvisible)或是直接于画面上点取你要显示壳厚的Part亦可.
41.
请问在LS-Pre/Post里面该如何操作才能看到各相材料的变形及相关的力学资料解答:(1)进入Fcomp选单.
(2)将下方的Fringe钮切换成Ffds.
(3)于CFD选项内选择欲观测的物理量.
(4)进入Selpar选单后,您会发现多了一些F1,F2之类的part,此即为参与ALE分析的各相材料您可以LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第33页直接点选欲观察的part.
42.
请问LS-Pre/Post有没有办法改变SPH颗粒的大小解答:在Setting的选项里面有提供SPH颗粒大小的控制方式,请依照下列步骤进行:(1)点选SPHnodes.
(2)下方字段即可选择SPH的半径比例及显示方式.
43.
请问LS-DYNA有没有办法仿真出「破裂」或「碎裂」的状况解答:您可以参考LS-DYNAUser'sManualVol.
II选择具备failurecriterion的材料模式.
以最简单的*MAT_KINEMATIC_PLASTIC模式而言,您只要给定FS(plasticstrain)的值,当结构体变形达到此值即会进行elementelimination的的动作,若mesh够密即可看到「破裂」的样子.
话说回来,「破裂」现像还有许多判别的方式,利用小技巧来产生碎裂的行为给大家看,通常是「作秀」的成分居多,对于电子产品结构系统的改善并无直接的帮助.
要进行真正的大变形仿真,其实可以考虑采用EFG或是ALE的方式来进行,右图所示为LSTC原厂所提供「metalcutting」的案例.
此例即是利用ALE的方式仿真金属切削的行为,突破一般FEMLagrangemesh的变形限制及建构网格的困扰.
44.
请问究竟是采用Automaticcontact好还是一一指定接触界面的方式比较好解答:基本上这个问题也是属于工程师判断的问题.
早期的接触算法是需要一一指定接触界面,甚至工程师得非常小心各个接触界面的法向量必须要两两相对.
这样的方式可想而知,很容易发生人为的疏忽,进而产生结构体发生穿透的情形,面对复杂的结构系统而言几乎是不太可行的方式.
个人的建议是除非有特殊理由,必须明显区分不同的接触界面,如打印机进纸机构等产品,否则可先以*contact_automatic_single_surface的方式来进行分析.
以附图为例,这是很单纯的积木堆栈结构,然而却让人很难预测在球击的过程中,哪些界面会发生接触行为.
若是采用一一指定接触界面的方式来进行设定其难度可想而知.
45.
在完成前处理的设定后,如何知道每一个part的重量解答:你可以在FEMB的Part选项里找到MassCalculation的功能.
程序允许你将各个part的质量信息存成文字文件,同时也会立即将这些信息显示在屏幕上.
除了FEMB自动计算质量信息的功能外,当LS-DYNA开始进行分析的时候,会自动产生d3hsp的档案,里头也会记录质量、质心位置等相关信息.
46.
如何将Excel档案直接读入FEMB中解答:FEMB建立数值曲线的方式有两种,一种是在FEMB里用key-in的方式来输入点位,另一种方式就是读入Excel的格式,其方式如下:(1)请先将Excel档另外储存成附档名为*.
csv的格式,再将附档名改成「*.
cvr」.
(2)进入FEMB「DynaMisc」的功能选项,选择「LoadCurve」CreateOKFileImport即可.
47.
请问贵公司所代理的LS-DYNA软件,其采用的单位制为何解答:事实上LS-DYNA并不强制你采用何种单位制,只要model的尺度及相关力学常数能彼此搭配一致即可.
在此提供单位换算的简表,方便您进行查核.
个人的习惯是采用Ton-mm-sec制,因为这样所得到力量、应力及能量单位皆是一般惯用的单位制.
48.
请问在loadcurve的输入方式除了一点一点慢慢key-in之外,有没有其它更快的方式解答:如果你的线型很单纯,用key-in的方式会比较快速;如果你的曲线是属于函数图形之类的话,可LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第34页依下列方式,直接将Excel资料读入FEMB当中:(1)利用Excel产生两栏(x,y)的数值,将其存成*.
csv的格式.
(2)利用档案总管将*.
csv的扩展名改成*.
cvr(3)在FEMB的DynaMiscLoadcurveCreate一直按OK进入直角坐标图FILEImport将*.
cvr的曲线点位謮入即可.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第35页49.
如何处理LS-DYNA中的退化单元在网格划分过程中,我们常遇到退化单元,如果不对它进行一定的处理,可能会对求解产生不稳定的影响.
在LS-DYNA中,同一PartID下既有四面体,五面体和六面体,则四面体,五面体既为退化单元,节点排列分别为N1,N2,N3,N4,N4,N4,N4,N4和N1,N2,N3,N4,N5,N5,N6,N6.
这样退化四面体单元中节点4有5倍于节点1-3的质量,而引起求解的困难.
其实在LS-DYNA的单元公式中,类型10和15分别为四面体和五面体单元,比退化单元更稳定.
所以为网格划分的方便起见,我们还是在同一PartID下划分网格,通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元.
50.
LS-DYNA中对于单元过度翘曲的情况有何处理方法有两种方法:1.
采用默认B-T算法,同时利用*control_shell控制字设置参数BWC=1,激活翘曲刚度选项;2.
采用含有翘曲刚度控制的单元算法,第10号算法.
该算法是针对单元翘曲而开发的算法,处理这种情况能够很好的保证求解的精度.
除了上述方法外,在计算时要注意控制沙漏,确保求解稳定.
翘曲翘曲是由板厚内的应力——应变分布不均匀产生的.
成形的材料通过rd进入侧壁部时,在被成形的侧壁部产生的弯曲,称为翘曲.
如图1所示.
消除方法①不用压板弯曲带凸缘制品时,凸缘在行程最后阶段压出.
如果rd大,就可能成形.
如果rd小,首先用大rd进行粗弯曲.
然后增加一道精整工序.
②由于用厚板产生的翘曲很大,所以先弯曲凸缘部,后进行内侧壁弯曲.
③使用薄板时,可用拉深筋或多次拉深增加张力,同时成形凸缘和侧壁.
51.
在ANSYS计算过程中结果文件大于8GB时计算自动中断,如何解决这个问题解决超大结果文件的方案:1.
将不同时间段内的结果分别写入一序列的结果记录文件;2.
使用/assign命令和重启动技术;3.
ANSYS采用向指定结果记录文件追加当前计算结果数据方式使用/assign指定的文件,所以要求指定的结果记录文件都是新创建的文件,否则造成结果文件记录内容重复或混乱.
特别是,反复运行相同分析命令流时,在重复运行命令流文件之前一定要删除以前生成的结果文件序列.
具体操作方法和过程参见下列命令流文件的演示.
52.
关于梁、壳单元应力结果输出的说明问题:怎样显示梁单元径向和轴向的应力分布图(我作的梁单元结果只有变形图DOFSOLUTIN–Translation,但是没有stress等值线图,只有一种颜色)和壳单元厚度方向的应力、变形图(我们只能显示一层应力、变形,不知道是上下表层或中间层的结果).
解答:如果想显示梁单元的应力等值线图,请打开实际形状显示功能(PLotCtrl->Style->SizeandLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第36页Shape->/ESHAPE选为ON),然后即可绘制.
注意梁单元(如BEAM188,BEAM189)的应力结果是在单元坐标系中显示的,即SXX为轴向正应力,SXY,SXZ为截面剪应力,没有其他应力分量.
另外,缺省情况下,只输出SXX,如果想观察SXY,SXZ,请将BEAM188或189的KEYOPT(4)选为Includeboth(以这两个单元为例,其他单元可能不同,请看帮助文件,推荐使用BEAM188,BEAM189,这是功能最强的梁单元).
至于壳的应力显示也类似,请打开实际形状显示功能,即可如同在实体上一样显示结果,您可以很清楚地看出不同位置、高度的应力值.
当然如果你只想画出顶部、中部或底部的应力图也可以,以shell63为例,首先需关闭powergraphics(Toolbar上点POWRGRPH,选择OFF),然后进入GeneralPostProc->Optionforoutp->SHELL中选择位置即可.
53.
LS-DYNA求解有时为什么有负的滑动能这是由于在建立模型时PART与PART之间有初始穿透,尤其是壳单元模型时很容易发生,应当避免这种情况的出现,否则容易在有初始穿透的地方产生塑性铰,原因是程序在求解的开始阶段给与穿透相应的接触力消除穿透,使材料发生局部塑性变形.
解决方法见2002年11月的应用技巧.
54.
在DYNA中如何考虑材料失效问题:在LS-DYNA的材料库中,能考虑失效的材料其失效模式往往比较单一,或者是应力失效,或是应变失效,如果材料本身较为复杂,在破坏过程中可能涉及多种失效模式,能否在一种材料中同时定义多种实效模式答:可以.
LS-DYNA材料库中提供了专门定义失效准则的命令,即*mat_add_erosion,利用该命令,可以同时定义压力、主应力、等效应力、主应变、临界应力以及应力脉冲六种失效准则,在加载过程中满足任何一种失效准则都会使材料发生破坏.
55.
在LS-DYNA中能否施加跟随力和跟随力矩答:能,对于一些应用,施加的载荷相对与坐标系不仅大小变化,而且方向变化,此时按照通常的施加力方向(X、Y、Z)不能满足要求,在LS-DYNA中,可以方便的施加跟随力和跟随力矩,在关键字*LOAD_NODE_OPTION中,对DOF选择4和8就可以施加跟随力和跟随力矩.
56.
如果在工程上遇到壳的厚度是坐标位置的函数时,这样的壳单元模型如何建立我们常用到等厚度的壳单元,如果在工程上遇到壳的厚度是坐标位置的函数时,这样的壳单元模型如何建立要用到RTHICK命令.
/PREP7ET,1,63RECT,,10,,10ESHAPE,2ESIZE,,20AMESH,1EPLOMXNODE=NDINQR(0,14)*DIM,THICK,,MXNODE*DO,NODE,1,MXNODE*IF,NSEL(NODE),EQ,1,THENTHICK(node)=0.
5+0.
2*NX(NODE)+0.
02*NY(NODE)**2*ENDIF*ENDDONODE=MXNODE=RTHICK,THICK(1),1,2,3,4/ESHAPE,1.
0/USER,1/DIST,1,7/VIEW,1,-0.
75,-0.
28,0.
6/ANG,1,-1LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第37页/FOC,1,5.
3,5.
3,0.
27EPLO57.
不同自由度的单元在merge时,ANSYS如何处理自由度的不匹配问题ANSYS可以将在给定容差(tolerance)内的节点通过merge可以合并到一起,它可以是不同的单元类型和不同的自由度(如:壳或梁(6自由度)与块单元(3自由度),ANSYS在处理共节点的自由度关系使,类似于藕合自由度(CoupleDOFs)而非约束方程,即只把相同的自由度联系起来,这样就可以方便的处理如柱销、铰链之类的约束问题.
58.
如何方便地建立压力容器的有限元模型在建立压力容器模型时,充分考虑模型的对称性,灵活利用ANSYS软件的工作平面和坐标系,利用合理的分网工具,可以得到漂亮的有限元模型,如下面的命令流所示:/prep7et,1,45mp,ex,1,2e11mp,nuxy,1,0.
3!
cylind,0.
89,0.
8,0,-1.
7125,90,270,wpoff,0,0,-1.
7125wprot,0,90,cylind,0.
47/2,0.
37/2,-1.
5,0,90,180,vovlap,allvsel,s,loc,x,-0.
11,0vdel,all,,,1vsel,s,loc,y,0.
3,0.
5vdel,all,,,1vsel,s,loc,y,-0.
3,0.
3vsbw,allafillt,21,11,0.
1al,33,50,5al,37,53,7alls,allva,15,13,25,24,11kgen,2,35,,,,-0.
2,,,0lstr,35,15adrag,54,,,,,,21vsba,2,16lsel,s,loc,y,1.
5lsel,u,length,,0,0.
06lesize,all,,,16,alls,allvsel,s,loc,y,0.
9,1.
5vsweep,allvsel,s,loc,y,0.
89,0.
99vsweep,allvsel,s,loc,y,0.
8,0.
89vsweep,allasel,s,loc,z,asel,r,loc,y,0,1asel,a,loc,y,0accat,allalls,allvsel,s,loc,y,0,0.
6vmesh,allvsel,s,loc,y,0,-0.
6vmesh,allalls,all!
WPCSYS,-1,0wpoff,0,0,-1.
7125CSWPLA,11,0,1,1,VSYMM,z,all,nummrg,allnumcmp,allwpoff,0,0,-1.
7125CSWPLA,12,0,1,1,VSYMM,z,all,nummrg,allnumcmp,all!
wpoff,0,0,-3.
425wprot,0,-90sphere,0.
89,0.
8,90,180,csys,0vsel,s,loc,z,-8.
0,-6.
85vsbw,allalls,allnummrg,allnumcmp,allvsel,s,loc,z,-8.
0,-6.
85vsweep,allLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第38页59.
ANSYS是否具有混合分网的功能ANSYS具有混合网格剖分的功能.
例如两个粘在一起的面,可以对一个面进行三角形划分,再对另一个面进行四边形划分.
过程见下列命令:/prep7et,1,42rect,,1,,1rect,1,2,,1aglue,allmshape,0,2damesh,1mshape,1,2damesh,360.
在ANSYS中怎样给面施加一个非零的法向位移约束给面施加非零的法向位移约束的过程如下:1)在面上施加一个对称约束条件(DA,2,SYMM)2)将实体模型上的载荷传递到有限元模型(SBCTRAN)3)选择需要施加约束条件的面(ASEL,S,,,2)4)选择附在面上的节点(NSLA,S,1)5)创建节点组元(CM,AREA2_N,NODE)6)删除面上的对称约束条件(DADELE,2,SYMM)7)删除节点上的零位移约束(DDELE,AREA2_N,UY)8)在节点组元上施加一个非零的法向位移约束(D,AREA2_N,UY,.
05)9)图示节点验证约束是否正确(NPLOT)61.
如何得到径向和周向的计算结果在圆周对称结构中,如圆环结构承受圆周均布压力.
要得到周向及径向位移,可在后处理/POST1中,通过菜单GeneralPostproc>OptionsforOutp>Rsys>Globalcylindric或命令Rsys,1将结果坐标系转为极坐标,则X方向位移即为径向位移,Y向位移即为周向位移.
62.
如何正确理解ANSYS的节点坐标系节点坐标系用以确定节点的每个自由度的方向,每个节点都有其自己的坐标系,在缺省状态下,不管用户在什么坐标系下建立的有限元模型,节点坐标系都是与总体笛卡尔坐标系平行.
有限元分析中的很多相关量都是在节点坐标系下解释的,这些量包括:输入数据:1自由度常数2力3主自由度4耦合节点5约束方程等输出数据:1节点自由度结果2节点载荷3反作用载荷等但实际情况是,在很多分析中,自由度的方向并不总是与总体笛卡尔坐标系平行,比如有时需要用柱坐标系、有时需要用球坐标系等等,这些情况下,可以利用ANSYS的"旋转节点坐标系"的功能来实现节点坐标系的变化,使其变换到我们需要的坐标系下.
具体操作可参见ANSYS联机帮助手册中的"分析过程指导手册->建模与分网指南->坐标系->节点坐标系"中说明的步骤实现.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第39页63.
如何考虑结构分析中的重力在结构分析中,如何模拟结构自重和设备重量是一个经常遇到的问题,对于结构自重有两点要注意:1.
在材料性质中输入密度,如果不输入密度,则将不会产生重力效果.
2.
因为ANSYS将重力以惯性力的方式施加,所以在输入加速度时,其方向应与实际的方向相反.
对于结构上的设备重量可以用MASS21单元来模拟,该单元为一个空间"点"单元.
设备重量可通过单元实常数来输入.
下面附上一个小例子(设重力方向向下).
/prep7et,1,42et,2,21r,2,10,10,10mp,ex,1,2e5mp,nuxy,1,0.
3mp,dens,1,1rect,,10,,1esize,.
5amesh,alltype,2real,2e,node(5,1,0)fini/soludk,1,alldk,2,uy,acel,,10solvefini/post1plnsol,u,sum,2/SOLUANTYPE,MODALMODOPT,subspa,2,,,2,ONMXPAND,2,,,YESSOLVEFINISH/post1set,1,1etabl,kene,kenessum*get,keneval1,ssum,,item,kene*get,freqval1,mode,1,freqeigen1=(2*3.
14159*freqval1)**2pmass1=2*keneval1/eigen1set,1,2etabl,kene,kenessum*get,keneval2,ssum,,item,kene*get,freqval2,mode,2,freqeigen2=(2*3.
14159*freqval2)**2pmass2=2*keneval2/eigen2finish64.
对于具有高度不规则横截面的3D模型什么是最佳网格划分方法答:在横截面上自由划分四边形网格,然后在体内扫掠成六面体单元.
在扫掠前可对四边形网格加密(如需要).
确认加密后生产的单元保持四边形以保证扫掠成六面体单元.
(ANSYS5.
6/FLOTRAN支持锲形单元,所以无此要求.
)set,1,1etabl,kene,kenessum*get,keneval1,ssum,,item,kene*get,freqval1,mode,1,freqeigen1=(2*3.
14159*freqval1)**2LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第40页pmass1=2*keneval1/eigen1set,1,2etabl,kene,kenessum*get,keneval2,ssum,,item,kene*get,freqval2,mode,2,freqeigen2=(2*3.
14159*freqval2)**2pmass2=2*keneval2/eigen265.
在交互方式下如何施加任意矢量方向的表面载荷答:若需在实体表面上施加任意方向的表面载荷,可通过在实体表面生成表面效应单元(比如SURF154单元)的方法来完成.
施加面载荷时,可施加在表面效应单元上,这样可以任意控制面力的方向.
加载过程中,选定表面效应单元,对话框中LKEY取值不同,则所加表面载荷的方向不同.
(请仔细看一看surf154的单元手册).
比如:LKEY=1(缺省),载荷垂直于表面;LKEY=2,载荷为+X切向;LKEY=3,载荷为+Y切向;LKEY=4,载荷垂直于表面;LKEY=5,则可输入任意矢量方向的载荷.
特别地:LKEY=5,VALUE项为均布压力值VAL2、VAL3、VAL4三项的值确定矢量的方向.
66.
LS-DYNA96版后在爆炸及流固耦合方面的功能增强在LS-DYNA中,处理爆炸和流固耦合单元一般采用ALE列式和Euler列式(也可采用Lagrange),从而克服单元严重畸变引起的数值计算困难,并实现流体-固体耦合的动态分析.
ALE列式先执行一个或几个Lagrange时步计算,此时单元网格随材料流动而产生变形,然后执行ALE时步计算:(1)保持变形后的物体边界条件,对内部单元进行重分网格,网格的拓扑关系保持不变,称为SmoothStep;(2)将变形网格中的单元变量(密度、能量、应力张量等)和节点速度矢量输运到重分后的新网格中,称为AdvectionStep.
用户可以选择ALE时步的开始和终止时间,以及其频率.
Euler列式则是材料在一个固定的网格中流动,在LS-DYNA中只要将有关实体单元标志Euler算法,并选择输运(advection)算法.
LS-DYNA还可将Euler网格与全Lagrange有限元网格方便地耦合,以处理流体与结构在各种复杂载荷条件下的相互作用问题,并在95和96版中得到了极大的增强.
67.
ANSYS坐标系总结工作平面(WorkingPlane)工作平面是创建几何模型的参考(X,Y)平面,在前处理器中用来建模(几何和网格)总体坐标系在每开始进行一个新的ANSYS分析时,已经有三个坐标系预先定义了.
它们位于模型的总体原点.
三种类型为:CS,0:总体笛卡尔坐标系CS,1:总体柱坐标系CS,2:总体球坐标系数据库中节点坐标总是以总体笛卡尔坐标系,无论节点是在什么坐标系中创建的.
局部坐标系局部坐标系是用户定义的坐标系.
局部坐标系可以通过菜单路径Workplane>LocalCS>CreateLC来创建.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第41页激活的坐标系是分析中特定时间的参考系.
缺省为总体笛卡尔坐标系.
当创建了一个新的坐标系时,新坐标系变为激活坐标系.
这表明后面的激活坐标系的命令.
菜单中激活坐标系的路径Workplane>ChangeactiveCSto>.
节点坐标系每一个节点都有一个附着的坐标系.
节点坐标系缺省总是笛卡尔坐标系并与总体笛卡尔坐标系平行.
节点力和节点边界条件(约束)指的是节点坐标系的方向.
时间历程后处理器/POST26中的结果数据是在节点坐标系下表达的.
而通用后处理器/POST1中的结果是按结果坐标系进行表达的.
例如:模型中任意位置的一个圆,要施加径向约束.
首先需要在圆的中心创建一个柱坐标系并分配一个坐标系号码(例如CS,11).
这个局部坐标系现在成为激活的坐标系.
然后选择圆上的所有节点.
通过使用"Prep7>Move/Modify>RotateNodalCStoactiveCS",选择节点的节点坐标系的朝向将沿着激活坐标系的方向.
未选择节点保持不变.
节点坐标系的显示通过菜单路径Pltctrls>Symbols>NodalCS.
这些节点坐标系的X方向现在沿径向.
约束这些选择节点的X方向,就是施加的径向约束.
注意:节点坐标系总是笛卡尔坐标系.
可以将节点坐标系旋转到一个局部柱坐标下.
这种情况下,节点坐标系的X方向指向径向,Y方向是周向(theta).
可是当施加theta方向非零位移时,ANSYS总是定义它为一个笛卡尔Y位移而不是一个转动(Y位移不是theta位移).
单元坐标系单元坐标系确定材料属性的方向(例如,复合材料的铺层方向).
对后处理也是很有用的,诸如提取梁和壳单元的膜力.
单元坐标系的朝向在单元类型的描述中可以找到.
结果坐标系/Post1通用后处理器中(位移,应力,支座反力)在结果坐标系中报告,缺省平行于总体笛卡尔坐标系.
这意味着缺省情况位移,应力和支座反力按照总体笛卡尔在坐标系表达.
无论节点和单元坐标系如何设定.
要恢复径向和环向应力,结果坐标系必须旋转到适当的坐标系下.
这可以通过菜单路径Post1>Optionsforoutput实现.
/POST26时间历程后处理器中的结果总是以节点坐标系表达.
显示坐标系显示坐标系对列表圆柱和球节点坐标非常有用(例如,径向,周向坐标).
建议不要激活这个坐标系进行显示.
屏幕上的坐标系是笛卡尔坐标系.
显示坐标系为柱坐标系,圆弧将显示为直线.
这可能引起混乱.
因此在以非笛卡尔坐标系列表节点坐标之后将显示坐标系恢复到总体笛卡尔坐标系.
68.
显式隐式分析转换的注意事项运用ANSYS/LS-DYNA进行隐、显式分析时,由于隐、显式分析过程中所涉及的对象一般会有所不同,ANSYS/LS-DYNA使用手册中只介绍了一种方法,即下面所述的第一种.
实际上,共有两种方法可以选择使用:1、将隐、显式分析过程中涉及到的所有对象都在隐式分析前建好模型,把隐式分析不需要的对象的所有节点自由度都约束住,进行隐式求解,转换单元类型,进入显式求解阶段,将显式part的约束去除,执行动力松弛求解以便对相应part进行应力初始化,并按照需要施加新的边界和载荷条件,进行显式分析.
2、在隐式分析时只对隐式分析涉及的对象建模,而不考虑显式分析需要的part,完成隐式分析后,单元类型转换完成后,通过定义新的单元类型和材料,创建显式分析所需的模型,生成新的part列表,选择所有节点,读入隐式求解结果文件进行动力松弛求解,对相应part进行应力初始化,施加必要的约束和载荷条件,执行显式求解.
实际上,动力松弛过程是执行一次稳态或是准静态分析,目的就是将隐式分析的结果中的位移、温度结果作为体载荷施加到相关节点上,实现相应部件的应力初始化,作为后续分析的初始条件.
需要注意的是,LS-DYNA中无高阶单元,所以在进行隐式求解时要选择缩减积分的低阶单元.
如果隐式分析使用高阶单元,则程序无法自动转换单元类型,需要手动转换.
上面所述的是利用ANSYS作为隐式求解器时的操作方法.
我们知道,近几年来,LSTC公司不断加强LS-DYNA程序本身的隐式分析能力,所以我们也可以利用LS-DYNA本身的隐式求解器来完成隐式分析,也基本有两种方法:1、进行隐式分析时,涉及的关键字主要有:*control_implicit_solverLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第42页*control_implicit_general*control_implicit_solution*control_implicit_auto*control_implicit_dynamics等.
在这些命令中,设置隐式求解的求解方法(波前、迭代)、时间步长等控制参数.
在dyna的输入文件中加入下列命令,*interface_springback_nike3D在该关键字中,声明需要进行应力初始化的part号,完成隐式求解后,生成一个nikin文件,包含了相关part的应力应变信息.
在后续的显式分析中,在inputdeck中加入下列命令,*includenikin程序就会自动将存在应力、应变的相关part导入,进行显式分析.
2、另外,可以LS-DYNA的动力松弛方式来对某一构件进行应力初始化.
相关的关键字为:*DEFINE_CURVE将此卡片的SIDR参数设置为1即可启动动力松弛分析.
*CONTROL_DYNAMIC_RELAXATION此卡片在随后的显式分析中用来进行应力初始化操作.
*LOAD_BODY_RX(RY、RZ)等运行后收敛的结果即为初始化应力,同时生成动力松弛文件drdisp.
sif,该文件与drelax文件结构、用法完全一致,只是精度上较差.
建议:使用ANSYS作为隐式求解器,因为它的隐式功能和计算精度都优于LS-DYNA.
69.
利用LS-DYNA进行接触分析应该注意的一些问题在定义材料特性时确保使用了协调单位.
不正确的单位将不仅决定材料的响应,而且影响材料的接触刚度.
确保模型中使用的材料数据是精确的.
大多数非线性动力学问题的精度取决于输入材料数据的质量.
多花点时间以得到精确的材料数据.
对所给模型选择最合适的材料模型.
如果不能确定某个part的物理响应是否应该包含某个特殊特性(例如:应变率效应),定义一种包含所有可能特点的材料模型总是最好的.
在两个接触面之间不允许有初始接触,确保在定义接触的地方模型没有任何重叠.
总是使用真实的材料特性和壳厚度值,接触面的材料特性和几何形状被用来决定罚刚度.
在相同的part之间不要定义多重接触.
对壳单元,除非需要接触力否则使用自动接触.
无论何时尽可能使用自动单面接触(ASSC),此接触是最容易定义的接触类型而不花费过多的CPU时间.
在求解之前列示所定义的接触面以保证定义了合适的接触.
避免单点载荷,它们容易引起沙漏模式.
既然沙漏单元会将沙漏模式传给相邻的单元,应尽可能避免使用点载荷.
在定义载荷曲线之后,使用EDLDPLOT命令进行图形显示以确保其精确性.
因为LS-DYNA可能会多算几个微秒,将载荷扩展至超过最后的求解时间(终止时间)常常是有用的.
对准静态问题,施加一个高于真实情况的速度常常是有利的,这能极大的缩减问题的求解时间.
不允许约束刚体上的节点,所有的约束必须加在刚体的质心(通过EDMP,RIGID命令).
70.
LS-DYNA求解中途退出的解决方案LS-DYNA在求解过程中由于模型的各种问题常发生中途退出的问题,归纳起来一般有三种现象:一是单元负体积,二是节点速度无限大,三是程序崩溃.
LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第43页1.
单元负体积:这主要是由于人工时间步长设置的不合理,调小人工时间步长可解决该问题.
还有就是材料参数和单元公式的选择合理问题.
2.
节点速度无限大:一般是由于材料等参数的单位不一致引起,在建立模型时应注意单位的统一,另外还有接触问题,若本该发生接触的地方没有定义接触,在计算过程中可能会产生节点速度无限大.
3.
程序崩溃:该现象不常发生,若发生,首先检查硬盘空间是否已满,二是检查求解的规模是否超过程序的规模.
最后就是对于特定的问题程序本身的问题.
当然对于程序中途退出问题原因是比较复杂的,不过对于其他一些刚开始就中断的现象LS-DYNA都会提示用户怎样改正,如格式的不对,符号的缺少等等.
71.
液面晃动液面晃动(sloshing)问题的研究在实际工程中有重要的意义,比如在石化工业中广泛应用的大型储罐,一般直径在几十米,甚至上百米.
在地震或其他意外条件下液面的波动情况如何,是否存在安全隐患,都需要进行数值模拟研究.
下面就ANSYS/LS-DYNA软件在这方面的应用.
众所周知,ANSYS/LS-DYNA在显式计算领域占据主导地位,随着各种新的算法的不断采用,在求解的广度、精度以及效率上,ANSYS/LS-DYNA具有同类软件所无法比拟的优势.
针对液面晃动问题,ANSYS/LS-DYNA提供以下三种方法:1、流固耦合流固耦合是ANSYS/LS-DYNA计算流体和结构间相互作用的最常用的方法,包括单物质+空材料和多物质耦合两大类,流体单元有Euler和ALE两种.
其涉及的主要命令如下:*control_ale算法选择有两种2、3,分别为Euler和ale实质上此处二者没有区别,只是因为兼容性进行的设置;两种精度供选择-单精度、双精度.
*section_solid_ale对单物质+空材料为12号算法,对多物质耦合为11号算法.
*ale_multi-material_group进行多物质的定义,最多可以定义20种材料.
可以根据物质间能否混合将各种材料定义在不同的材料组ID中.
*ale_multi-material_system_group该命令决定流体物质的算法(Euler或Ale),或是在运算过程中切换使用两种算法,并可对流体物质进行自由度约束.
该命令多与下列三个命令结合使用:*ale_multi-material_system_curve定义ale系统的运动曲线.
*ale_multi-material_system_node通过一系列节点定义ale的运动参考坐标系统.
*ale_multi-material_system_switch定义euler和ale参考系统的切换.
上述命令是流体物质涉及的关键字,而我们知道,结构采用Lagrange单元来离散,二者之间的耦合通过下列命令来实现:*constrained_lagrange_in_solid耦合算法分为两种:罚耦合和运动约束.
前者遵循能量守恒,后者遵循动量守恒.
一般令结构网格较流体网格密以保证界面不出现渗透,否则可以增大NQUAD参数值来增加耦合点,如设置该值为4或5.
在970中,此命令第三行又增加了一个控制字ILEAK-0,1或2,一般可设置为1.
最后给出一个典型算例-水箱跌落的部分关键字:*KEYWORD*TITLEboxwater2.
k:droppingawaterboxontoarigidplatformLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第44页$[1]EXECUTIONCONTROLS*CONTROL_TERMINATION$ENDTIMENDCYCDTMINENDENGENDMAS0.
050000000.
000000000.
0000000*CONTROL_TIMESTEP$DTINITTSSFACISDOTSLIMTDT2MSLCTMERODEMS1ST0.
00000000.
200000000.
00000000.
0000000000*CONTROL_ENERGY$HGENRWENSLNTENRYLEN2222$[3]OUTPUTCONTROLS*DATABASE_BINARY_D3PLOT$DTCYCLLCDTBEAM0.
00050000*DATABASE_GLSTAT0.
0001000$[5]|SECTIONS|PARTS|DEFs*PARTwaterinthebox$PIDSECIDMIDEOSIDHGIDGRAVADPOPTTMID11110000*SECTION_SOLID_ALE$SECIDELFORMAET112$AFACBFACCFACDFACSTARTENDAAFAC0.
00000000.
00000000.
00000000.
00000000.
00000000.
0000000*MAT_NULL$MIDRHOPCMUTERODCERODYMPR11000.
0000-1.
000+100.
00000000.
00000000.
00000000.
00000000.
0000000*EOS_LINEAR_POLYNOMIAL$EOSIDC0C1C2C3C4C5C610.
00000001.
50000+90.
00000000.
00000000.
00000000.
00000000.
0000000$E0V00.
00000001.
0000000*PARTvoidportioninthebox21110000*INITIAL_VOID_PART2*PARTLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第45页rigidboxcontainingwater$PIDSECIDMIDEOSIDHGIDGRAVADPOPTTMID33300000*SECTION_SOLID$SECIDELFORMAET30*MAT_RIGID32000.
00001.
00000+80.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
00000000.
0000000*PARTrigidsuper-heavyplatform$PIDSECIDMIDEOSIDHGIDGRAVADPOPTTMID444*SECTION_SHELL$SIDELFORMSHRFNIPPROPTQR/IRIDICOMP40$T1T2T3T4NLOC0.
0110.
0110.
0110.
011*MAT_ELASTIC$MIDRHOEPRDADBK41000000.
01.
0000+14$[8]BC's+IC's+BODYLOADS+FORCEFIELDS*INITIAL_VELOCITY$NSIDNSIDEXBOXID0$VXVYVZVXRVYRVZR0.
0-20.
00.
0*LOAD_BODY_Y$LCIDSFLCIDDRXCYCZC11.
00*DEFINE_CURVE$LCIDSIDRSFOOFFAOFFODATTYP1$X=abcissaY=ordinate0.
0981.
01.
0981.
0$[9]LAGRANGIANCONTACTSCONSTRAINTS,.
.
.
$SFS=scalefactondfltSLAVEpenalstifns(seeCONTROLL_CONTACT)$SFM=scalefactondfltMASTERpenalstifns(seeCONTROLL_CONTACT)*CONTACT_AUTOMATIC_NODES_TO_SURFACE$SSIDMSIDSSTYPMSTYPSBOXIDMBOXIDSPRMPRLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第46页3433$FSFDDCVCVDCPENCHKBTDT$SFSSFMSSTMSTSFSTSFMTFSFVSF100.
100.
$[10]EULERIAN&ALECONTACTSCONSTRAINTS,.
.
.
*CONTROL_ALE$DCTNADVMETHAFACBFACCFACDFACEFAC214-1.
00000000.
00000000.
00000000.
0000000$STARTENDAAFACVFACTVLIMITEBC0.
00000000.
00000000.
00000000.
0*ALE_REFERENCE_SYSTEM_GROUP$SIDSTYPEPRTYPPRIDBCTRANBCEXPBCROTICOORD1051$XCYCZCEXPLIM*SET_PART_LIST$SIDDA1DA2DA3DA41$PID1PID2PID3PID4PID5PID6PID7PID812*ALE_REFERENCE_SYSTEM_NODE$NSID1$N1N2N3N4N5N6N7N85672、SPH算法SPH算法作为DYNA中第一种无网格(meshfree)算法,在连续体的破碎或分离分析中得到了广泛的关注和应用.
在解决极度变形和破坏类型的问题上SPH有着其他方法无法比拟的优势,可以说无网格算法正在成为数值分析领域的研究热点,具有很好的发展前景.
我们知道传统的有限单元法中,单元的形状对结果的精度影响很大,如果单元因为变形过大可能造成矩阵奇异,使得精度降低甚至无法计算下去.
而SPH算法则是把每个粒子作为一个物质的插值点,各个粒子间通过规则的内插函数计算全部质点即可得到整个问题的解.
主要的关键字如下:*section_sph提供算法选择,以及sph粒子的滑顺长度的定义;*control_sph提供sph算法的控制,如粒子排序后的循环次数、计算空间、中止时间以及维数;处理sph粒子与其它结构的相互作用采用接触算法.
下面给出某一算例的部分命令流:*KEYWORD*TITLEsphtest$*DATABASE_FORMATLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第47页0$units:cm,gm,us$CONTROLOPTIONS$*CONTROL_PARALLEL1*CONTROL_ENERGY2222*CONTROL_SHELL20.
01-11221*CONTROL_TIMESTEP0.
00000.
900000.
000.
00*CONTROL_TERMINATION$1000.
000000.
000000.
000000.
000000.
800E+0500.
000000.
000000.
00000*CONTROL_SPH20$$TIMEHISTORY$$*DATABASE_BINARY_D3PLOT200.
0E+00$0.
500E+00*DATABASE_BINARY_D3THDT0.
8000E+02*DATABASE_EXTENT_BINARY00310000004000$$$SECTIONDEFINITIONS$$*SECTION_SOLID21*SECTION_SPH1$$$MATERIALDEFINITIONS$$*MAT_RIGIDLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第48页27.
802.
100.
3000000.
00.
00.
01.
007.
007.
00*MAT_RIGID37.
802.
100.
3000000.
00.
00.
01.
006.
007.
00*MAT_RIGID47.
802.
100.
3000000.
00.
00.
01.
006.
007.
00*MAT_NULL11.
00*EOS_GRUNEISEN1.
14840001.
9790000.
0000000.
0000000.
11000003.
0000000.
0000000.
0000000$$$$PARTSDEFINITIONS$$$*PARTPart1forMat4andElemType11111000$*PARTPart2forMat2andElemType12220000$*PARTPart3forMat3andElemType13230000*PARTPart3forMat3andElemType14240000$$RIGIDBOUNDRIES$*LOAD_BODY_Y2,1.
0*DEFINE_CURVE20.
0,9.
8E-101.
0,9.
8E-10LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第49页$*DEFINE_CURVE101.
0001.
0000.
0000.
0000.
000000000000E+001.
000000000000E-041.
000000000000E+051.
000000000000E-04*BOUNDARY_PRESCRIBED_MOTION_RIGID3201-1.
0000.
0000.
000$$NODEDEFINITIONS$$*CONTACT_AUTOMATIC_NODES_TO_SURFACE113200000.
0000.
0000.
0000.
0000.
00000.
0000.
0000E+080.
0000.
0000.
1000.
0000.
0000.
0000.
0000.
00010.
10000003*SET_PART_LIST12,3,4*INCLUDEmesh.
k*END3、ALE(接触算法)采用接触算法分析流固耦合问题也是一种选择,在液面波动幅度较小时可以采用此种方法进行分析,流体用ALE算法描述,结构采用Lagrange算法;需要注意的一点:对ALE网格要进行滑顺处理,以控制网格形态,保证求解精度.
下面是某算例的部分命令流:*KEYWORD*TITLEALE$*DATABASE_FORMAT0$$CONTROLOPTIONS$$*CONTROL_PARALLEL1*CONTROL_ENERGY2222*CONTROL_ALE3121.
00000001.
00000000.
0000001.
00000001.
0000e+90.
0000000.
0000002*CONTROL_TIMESTEPLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第50页0.
00000.
900000.
000.
00*CONTROL_CONTACT0.
00000000.
0000000102000000*CONTROL_TERMINATION0.
100E+0500.
000000.
000000.
00000$$TIMEHISTORY$$*DATABASE_BINARY_D3PLOT0.
1000E+03*DATABASE_BINARY_D3THDT0.
1000E+02*DATABASE_EXTENT_BINARY00310000004000$$$SECTIONDEFINITIONS$$*SECTION_SOLID21*SECTION_SOLID_ALE15$$$MATERIALDEFINITIONS$$*MAT_RIGID27.
802.
100.
2800000.
00.
00.
01.
007.
007.
00*MAT_RIGID37.
802.
100.
2800000.
00.
00.
01.
006.
007.
00*MAT_RIGID47.
802.
100.
2800000.
00.
00.
01.
006.
007.
00*MAT_NULL11.
00000000.
00000001.
00000-80.
00000000.
00000000.
00000000.
0000000LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第51页*EOS_LINEAR_POLYNOMIAL11.
00000-61.
92100-30.
00000000.
00000000.
40000000.
40000000.
00000000.
00000000.
0000000$$$$PARTSDEFINITIONS$$$*PARTPart1forMat4andElemType11111000$*PARTPart2forMat2andElemType12220000$*PARTPart3forMat3andElemType13230000*PARTPart3forMat3andElemType14240000$$RIGIDBOUNDRIES$$*CONTACT_SURFACE_TO_SURFACE123300000.
00000.
00000.
00000.
00000.
000000.
00000.
1000E+081.
00000.
00010.
00000.
00001.
00001.
00001.
00001.
0000*CONTACT_SURFACE_TO_SURFACE133300000.
00000.
00000.
00000.
00000.
000000.
00000.
1000E+081.
00000.
00010.
00000.
00001.
00001.
00001.
00001.
0000*CONTACT_SURFACE_TO_SURFACE143300000.
00000.
00000.
00000.
00000.
000000.
00000.
1000E+081.
00000.
00010.
00000.
00001.
00001.
00001.
00001.
0000$*DEFINE_CURVE101.
0001.
0000.
0000.
0000.
000000000000E+001.
000000000000E-041.
000000000000E+051.
000000000000E-04*BOUNDARY_PRESCRIBED_MOTION_RIGIDLS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第52页3201-1.
0000.
0005.
00e3*BOUNDARY_PRESCRIBED_MOTION_RIGID4201-1.
0005.
00e30.
000$通常我们处理液面晃动采用上述方法,LS-DYNA求解器在下一个版本LS970中在MESHFREE功能上增加了更为稳定和高效的EFG技术,在解决此类问题上将给予我们更大的灵活性和更多的选择.
72.
复杂几何模型的系列网格划分技术众所周知,对于有限元分析来说,网格划分是其中最关键的一个步骤,网格划分的好坏直接影响到解算的精度和速度.
在ANSYS中,大家知道,网格划分有三个步骤:定义单元属性(包括实常数)、在几何模型上定义网格属性、划分网格.
在这里,我们仅对网格划分这个步骤所涉及到的一些问题,尤其是与复杂模型相关的一些问题作简要阐述.
1.
自由网格划分自由网格划分是自动化程度最高的网格划分技术之一,它在面上(平面、曲面)可以自动生成三角形或四边形网格,在体上自动生成四面体网格.
通常情况下,可利用ANSYS的智能尺寸控制技术(SMARTSIZE命令)来自动控制网格的大小和疏密分布,也可进行人工设置网格的大小(AESIZE、LESIZE、KESIZE、ESIZE等系列命令)并控制疏密分布以及选择分网算法等(MOPT命令).
对于复杂几何模型而言,这种分网方法省时省力,但缺点是单元数量通常会很大,计算效率降低.
同时,由于这种方法对于三维复杂模型只能生成四面体单元,为了获得较好的计算精度,建议采用二次四面体单元(92号单元).
如果选用的是六面体单元,则此方法自动将六面体单元退化为阶次一致的四面体单元,因此,最好不要选用线性的六面体单元(没有中间节点,比如45号单元),因为该单元退化后为线性的四面体单元,具有过刚的刚度,计算精度较差;如果选用二次的六面体单元(比如95号单元),由于其是退化形式,节点数与其六面体原型单元一致,只是有多个节点在同一位置而已,因此,可以利用TCHG命令将模型中的退化形式的四面体单元变化为非退化的四面体单元,减少每个单元的节点数量,提高求解效率.
在有些情况下,必须要用六面体单元的退化形式来进行自由网格划分,比如,在进行混合网格划分(后面详述)时,只有用六面体单元才能形成金字塔过渡单元.
对于计算流体力学和考虑集肤效应的电磁场分析而言,自由网格划分中的层网格功能(由LESIZE命令的LAYER1和LAYER2域控制)是非常有用的.
2.
映射网格划分映射网格划分是对规整模型的一种规整网格划分方法,其原始概念是:对于面,只能是四边形面,网格划分数需在对边上保持一致,形成的单元全部为四边形;对于体,只能是六面体,对应线和面的网格划分数保持一致;形成的单元全部为六面体.
在ANSYS中,这些条件有了很大的放宽,包括:1面可以是三角形、四边形、或其它任意多边形.
对于四边以上的多边形,必须用LCCAT命令将某些边联成一条边,以使得对于网格划分而言,仍然是三角形或四边形;或者用AMAP命令定义3到4个顶点(程序自动将两个顶点之间的所有线段联成一条)来进行映射划分.
2面上对边的网格划分数可以不同,但有一些限制条件.
3面上可以形成全三角形的映射网格.
4体可以是四面体、五面体、六面体或其它任意多面体.
对于六面以上的多面体,必须用ACCAT命令将某些面联成一个面,以使得对于网格划分而言,仍然是四、五或六面体.
5体上对应线和面的网格划分数可以不同,但有一些限制条件.
对于三维复杂几何模型而言,通常的做法是利用ANSYS布尔运算功能,将其切割成一系列四、五或六面体,然后对这些切割好的体进行映射网格划分.
当然,这种纯粹的映射划分方式比较烦琐,需要的时间和精力较多.
面的三角形映射网格划分往往可以为体的自由网格划分服务,以使体的自由网格划分满足一些特定的要求,比如:体的某个狭长面的短边方向上要求一定要有一定层数的单元、某些位置的节点必须在一条直线上、等等.
这种在进行体网格划分前在其面上先划分网格的方式对很多复杂模型可以进行良好的控制,但别忘了在体网格划分完毕后清除面网格(也可用专门用于辅助网格划分的虚拟单元类型-MESH200-来划分面网LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第53页格,之后不用清除).
3.
拖拉、扫略网格划分对于由面经过拖拉、旋转、偏移(VDRAG、VROTAT、VOFFST、VEXT等系列命令)等方式生成的复杂三维实体而言,可先在原始面上生成壳(或MESH200)单元形式的面网格,然后在生成体的同时自动形成三维实体网格;对于已经形成好了的三维复杂实体,如果其在某个方向上的拓扑形式始终保持一致,则可用(人工或全自动)扫略网格划分(VSWEEP命令)功能来划分网格;这两种方式形成的单元几乎都是六面体单元.
通常,采用扫略方式形成网格是一种非常好的方式,对于复杂几何实体,经过一些简单的切分处理,就可以自动形成规整的六面体网格,它比映射网格划分方式具有更大的优势和灵活性.
4.
混合网格划分混合网格划分即在几何模型上,根据各部位的特点,分别采用自由、映射、扫略等多种网格划分方式,以形成综合效果尽量好的有限元模型.
混合网格划分方式要在计算精度、计算时间、建模工作量等方面进行综合考虑.
通常,为了提高计算精度和减少计算时间,应首先考虑对适合于扫略和映射网格划分的区域先划分六面体网格,这种网格既可以是线性的(无中节点)、也可以是二次的(有中节点),如果无合适的区域,应尽量通过切分等多种布尔运算手段来创建合适的区域(尤其是对所关心的区域或部位);其次,对实在无法再切分而必须用四面体自由网格划分的区域,采用带中节点的六面体单元进行自由分网(自动退化成适合于自由划分形式的单元),此时,在该区域与已进行扫略或映射网格划分的区域的交界面上,会自动形成金字塔过渡单元(无中节点的六面体单元没有金字塔退化形式).
ANSYS中的这种金字塔过渡单元具有很大的灵活性:如果其邻接的六面体单元无中节点,则在金字塔单元四边形面的四条单元边上,自动取消中间节点,以保证网格的协调性.
同时,应采用前面描述的TCHG命令来将退化形式的四面体单元自动转换成非退化的四面体单元,提高求解效率.
如果对整个分析模型的计算精度要求不高、或对进行自由网格划分区域的计算精度要求不高,则可在自由网格划分区采用无中节点的六面体单元来分网(自动退化成无中节点的四面体单元),此时,虽然在六面体单元划分区和四面体单元划分区之间无金字塔过渡单元,但如果六面体单元区的单元也无中节点,则由于都是线性单元,亦可保证单元的协调性.
5.
利用自由度耦合和约束方程对于某些形式的复杂几何模型,可以利用ANSYS的约束方程和自由度耦合功能来促成划分出优良的网格并降低计算规模.
比如,利用CEINTF命令可以将相邻的体在进行独立的网格划分(通常是采用映射或扫略方式)后再"粘结"起来,由于各个体之间在几何上没有联系,因此不用费劲地考虑相互之间网格的影响,所以可以自由地采用多种手段划分出良好的网格,而体之间的网格"粘结"是通过形函数差值来进行自由度耦合的,因此连接位置处的位移连续性可以得到绝对保证,如果非常关注连接处的应力,可以如下面所述再在该局部位置建立子区模型予以分析.
再如,对于循环对称模型(如旋转机械等),可仅建立一个扇区作为分析模型,利用CPCYC命令可自动对扇区的两个切面上的所有对应节点建立自由度耦合条件(用MSHCOPY命令可非常方便地在两个切面上生成对应网格).
6.
利用子区模型等其它手段子区模型是一种先总体、后局部的分析技术(也称为切割边界条件方法),对于只关心局部区域准确结果的复杂几何模型,可采用此手段,以尽量小的工作量来获得想要的结果.
其过程是:先建立总体分析模型,并忽略模型中的一系列细小的特征,如导角、开孔、开槽等(因为根据圣维南原理,模型的局部细小改动并不特别影响模型总的分析结果),同时在该大模型上划分较粗的网格(计算和建模的工作量都很小),施加载荷并完成分析;其次,(在与总体模型相同的坐标系下)建立局部模型,此时将前面忽略的细小特征加上,并划分精细网格(模型的切割边界应离关心的区域尽量远),用CBDOF等系列命令自动将前面总体模型的计算结果插值作为该细模型的边界条件,进行求解计算.
该方法的另外好处是:可以在小模型的基础上优化(或任意改变)所关心的细小特征,如改变圆角半径、缝的宽度等;总体模型和局部模型可以采用不同的单元类型,比如,总体模型采用板壳单元,局部模型采用实体单元等.
子结构(也称超单元)也是一种解决大型问题的有效手段,并且在ANSYS中,超单元可以用于诸如各种非线性以及装配件之间的接触分析等,有效地降低大型模型的求解规模.
巧妙地利用结构的对称性对实际工作也大有帮助,对于常规的结构和载荷都是轴对称或平面对称的问题,毫无疑问应该利用其对称性,对于一些特殊情况,也可以加以利用,比如:如果结构轴对称而载荷非轴对称,则可用ANSYS专门用于处理此类问题的25、83和61号单元;对于由多个部件构成装配件,如果其每个零件都满足平面对称性,但各对称平面又不是同一个的情况下,则可用多个对称面来处理模型(或至少可用此LS-DYNA常见问题汇总2.
0yuminhust20052008-10-12第54页方法来减少建模工作量:各零件只需处理一半的模型然后拷贝或映射即可生成总体模型).
总之,对于复杂几何模型,综合运用多种手段建立起高质量、高计算效率的有限元模型是极其重要的一个步骤,这里介绍的注意事项仅仅是很少一部分,用户自己通过许多工程问题的不断摸索、总结和验证才是最能保证有效而高效地处理复杂模型的手段.
73.
LS-DYNA典型的射流模拟过程问题对于正在设计过程中的破甲弹,在试验之前了解其起爆、射流的形成及穿甲的全过程是非常重要的.
计算模型一个完整的弹体,从起爆到射流穿甲全过程都可以在ANSYS/LS-DYNA完全模拟.
在ANSYS/LS-DYNA中计算时,如果采用三维实体进行计算时,炸药、药型罩、空气、泡沫等材料都最好考虑有状态方程,弹壳可采用一般的弹塑性材料.
Keyword文件中会涉及到如下材料相关的关键字:*MAT_HIGH_EXPLOSIVE_BURN*EOS_JWL*MAT_NULL*EOS_LINEAR_POLYNOMIAL*MAT_STEINBERG*EOS_GRUNEISEN*MAT_ELASTIC_PLASTIC_HYDRO*MAT_PLASTIC_KINEMATIC计算区域包括破甲弹弹体部分以及射流在空中行程路径周围的空气.
其中炸药、空气和药形罩采用ALE算法,泡沫、弹壳等可以采用Lagrange算法.
也可全采用Euler算法.
计算到一定程度,逐渐删除对计算影响越来越小的part.
涉及到的关键字:*SECTION_SOLID*SECTION_SOLID_ALE*CONSTRAINED_LAGRANGE_IN_SOLID*ALE_MULTI-MATERIAL_GROUP_PART*INITIAL_DETONATION为了降低计算量,要充分考虑到弹体的对称性,计算模型只考虑取几分之一即可.
如图2所示.
对各个部分划分非常细的均匀网格(特别是药形罩部分和弹体中心线射流拉长流动的区域),网格要求都是六面体的,选用SOLID164单元.
载荷及边界条件由于计算是采用几分之一的弹体进行的,在切开的对称面上需要施加滑移边界条件.
涉及到的关键字:*BOUNDARY_SLIDING_PLANE*SET_NODE_LIST*BOUNDARY_SPC_SET计算结果处理计算过程中可逐渐删除泡沫、炸药、弹体外壳材料等影响已不大的PART.
在LSPOST中可以非常方便的绘制系统动能、内能和总能量随时间的变化曲线.
下图为系统动能随时间变化的曲线图:计算完成后,应和理论分析结果进行比较,并让专业设计工程师判断,如果有条件,可与试验结果进行比较.
对实际工程设计应有较大的指导作用,特别适合于试验之前设计初期的方案选择,设计后期的方案验证等等.
应该选择ansys分析结果里面的哪些数值来判别材料会不会失效判断材料是否失效的准则与材料特性有关,一般来讲,对于各向同性塑性材料可以用VonMises应力与拉伸极限来比较;对于各向同性脆性材料,可以用第一和第三主应力分别与抗拉和抗压强度比较.
对于复合材料需要比较拉压、剪切等不同情况.
对于各向异性则需对不同方向的强度分别比较.