lsdyna-问题解决.docVIP

1、hypermesh找CAE模型的重心的方法Y( ]! G, c9 V3 {post-summary，点击load工具条，找到Altair目录hw7.0-templates-summary选择你所应用的模块，hm,dyna,abaques,.etc.ctr_of_gravity,选择完之后点一下summary，就会弹出质量分布窗口，所有的信息都在这里了，可以显示整个模型的，也可以显示当前显示单元的。最后几行就是整个模型的重心坐标还有一种方法，就是你把模型调试好了，计算一下，在写出的d3hsp文件中，有重心的坐标?负体积定义？??^+ [8 h0 l, h7 \8 V+ z# z负体积是由于element本身产生大变形造成自我体积的内面跑到外面接着被判断为负体积。2 q M5 } e5 f: { ]5 h# k; w$ G I2 ~关于负体积的解决办法？# j6 x0 c: |8 y7 N4 z8 v3 ?. K2 g负体积多是网格畸变造成的，和网格质量以及材料、载荷条件都有关系可能的原因和解决的方法大概有几种：/ t??j7 p3 J* Z# k( y（1）材料参数设置有问题，选择合适的材料模式（2）沙漏模式的变形积累，尝试改为全积分单元3）太高的局部接触力（不要将force施在单一node上，最好分散到几个node上以pressure的方式等效施加），尝试调整间隙，降低接触刚度或降低时间步。9 [, ~7 n8 Q. ^* \（4）在容易出现大变形的地方将网格refine3 `6 E3 k+ G???8 ~) V! q0 J（5）材料换的太软,是不是也会出现负体积! G, ^1 v- q7 U（6）另外也可以采用ALE或是euler单元算法，用流固耦合功能代替接触，控制网格质量。例如在承受压力的单元在受压方向比其他方向尺寸长9 @# Y$ v @0 v9 \9 \! n/ A+ @7 Y N( n（7）尝试减小时间步长从0.9减小到0.6或更小。J! J% I: U% \3 j4 H8 `2 K* g( e4 T- |# l, { T! y经验总结：$ c) b( E7 [) f0 e L2 L( r; i时间步长急剧变小，可能是因为单元产生了严重的畸变而导致的负体积现象，如果采用的是四面体单元，你可以用网格重划分的方法来解决。如果你采用的是六面体单元，那目前就没有很有效的方法，可以试一下*ELEMENT_SOLID_EFG，那对机器的要求相对就会比较高了。out-of-range velocities错误??s J4 |0 E3 q（1）、划分的网格质量差，如长宽比超过约定；- j k( q2 l W4 J7 _, Y（2）、材料参数输入格式出错，如多出一个空格或一个逗号等（3）、输入的参数单位不协调，如在长度--mm；质量--kg；时间-s单位制下，压力的单位以Mpa输（4）、自定义的子程序存在问题，如岩石爆破中只考虑压，不考虑拉等（5）、介质太软，计算网格严重变形，如在淤泥中爆炸（6）、接触定义上出错。# @3 R% J8 {$ I( i f% S$ f( T我补充一个现象，不过说不出原因。5 F2 x: ?: L; Y8 r曾经做的算例中，网格划分与单位制绝对没有问题，计算时出现out-of-range velocities错误提示。检查出错节点时，发现都分布在边界上。想不出原因，干脆删掉了无反射边条，错误消失，可以计算。但不知原因. (More on type 6 HG control from Lee Bindeman) 类型6的沙漏控制与类型4,5不在于它用了一个假设应变场和材料属性来估算出假设应力场。这个应力在单元封闭域内进行积分得到沙漏力，因此单元表现的像一个有同样假设应变场的全积分单元。这种假设应变场设计成用来阻止纯弯曲中不真实的剪切变形和近似不可压材料中的体积锁死。 类型4和5的沙漏控制基于单元体积，波速和密度像在LS-DYNA理论手册中方程3.21那样来计算沙漏刚度。 沙漏类型6主要的改进是应力场在单元域内积分。这使得当使用大的长细比或者歪斜形状的体单元时沙漏控制非常鲁棒。类型4和5的沙漏控制对大长细比和歪斜形状单元反应变不好，它趋向于对某些沙漏模式反应的过于刚硬而对其它模式反应得过弱。 沙漏控制类型6另一个理论上的优点是对在厚度方向只有一个单元的梁可以在弹性弯曲问题中得到准确的解。要做到这一点，设置沙漏刚度参数为1.0。同样，对弹性材料方形截面杆的扭曲问题，当沙漏系数设为1.0时可以用很少的单元来解。然而，对于非线性材料，用粗糙的网格得到好的结果是不可能的，因为应力场不是像沙漏类型6假设的那样线性变化的。在梁厚度方向上如果没有更多积分点的话，没有办法捕获应力场的非线性状态。 对