负体积(L-DYNA).docVIP

负体积定义？ Negative volume8 i6 g, W# P7 _ z) e负体积是由于element本身产生大变形造成自我体积的内面跑到外面接着被判断为负体积。% L6 H P. l1 P关于负体积的解决办法？# i+ c$ P1 K$ F: R负体积多是网格畸变造成的，和网格质量以及材料、载荷条件都有关系。有可能的原因和解决的方法大概有几种：, m$ B8 E h% [ K1 j$ V4 @（1）材料参数设置有问题，选择合适的材料模式), C! {# M; ?% y: H0 U（2）沙漏模式的变形积累，尝试改为全积分单元 G: d, s7 F% `（3）太高的局部接触力（不要将force施在单一node上，最好分散到几个node上以pressure的方式等效施加），尝试调整间隙，降低接触刚度或降低时间步。6 u( o2 U8 a9 C/ ^) R2 N9 O（4）在容易出现大变形的地方将网格refine。8 M L* a; U- _ q. P4 M3 A# ]（5）材料换的太软,是不是也会出现负体积!: I( R% k+ n??e) u, `（6）另外也可以采用ALE或是euler单元算法，用流固耦合功能代替接触，控制网格质量。例如在承受压力的单元在受压方向比其他方向尺寸长。??x. s4 c* n7 G) [! H. k??h（7）尝试减小时间步长从0.9减小到0.6或更小。! J5 t0 U4 U: K0 M经验总结：, d `% Z |9 d+ b, U$ L时间步长急剧变小，可能是因为单元产生了严重的畸变而导致的负体积现象，如果采用的是四面体单元，你可以用网格重划分的方法来解决。如果你采用的是六面体单元，那目前就没有很有效的方法，可以试一下*ELEMENT_SOLID_EFG，那对机器的要求相对就会比较高了。- l) f- U+ q5 G4 J. ]$ C. O0 k) y# t??]3 D+ w, @7 r; OQ1:材料负体积解决方法（全面、有效）. S# J# \$ ]# Q材料负体积解决方法$ G+ D- @4 C- }( P6 F1 `* \ 在仿真中，通常有材料的大变形问题，如泡沫材料，由于单元大扭曲而出现了单元负体积，这种情况一般出来在材料失效之前。在没有网格光滑和网格从划分的情况下，ls-dyna有一个内部的限制来调节lagrange单元的变形。负体积一般都会导致计算中止，除非你设置时间步长控制中的erode=1和设置终止控制中的dtmin为一非零数，这种情况下，出现负体积的单元将被自动删除，计算也不会中止。不过就算你如上设置了erode与dtmin，负体积有时候也会导致计算出错停止。! w- Z4 n g+ s9 I h0 \ T 一些常用的解决负体积的方法如下：! M L) a2 _6 E8 E V在材料出现大应变的情况下增强材料的应力-应变曲线中材料应力。这种方法往往非常有效果。 B# @% W, C T$ f3 |5 I2、重新划分网格，在出现大变形的地方把网格加密。.% G$ m4 \1 L* a c, D5 r3、减小时间步长系数。默认的0.9系数可能不足以避免数值的不稳定。??{$ y/ G0 i C- U N, m$ H1 j4、避免采用全积分体单元（算法2和3），这会导致大变形和大扭曲的情况下计算相对不稳定。: P: w8 R% t5 P- h5、采用默认的单元算法（单点体单元），采用沙漏控制type4和5。泡沫的沙漏控制算法为：在低速冲击问题中采用type6，系数为1；在高速冲击问题中采用type2和3。??b$ x, e/ R: q??J, T, Q6、泡沫材料网格划分采用四面体网格，单元算法为10，虽然这样会导致材料相对比较刚性。增大材料（泡沫材料57号材料）的阻尼系数，推荐采用系数为0.5。) d, X2 m4 z0 [( h8、在泡沫接触计算中，采用*contact中的选项卡B，关闭shooting node logic。- C3 Q0 }/ p k, ~! m9、如果你采用的是126号材料，设置elform=0。- J+ l t$ ] m$ z6 r10、尝试使用EFG算法（*SECTION_SOLID_EFG）。3 n d, H5 q- ]9 C k ( l4 K K A$ ~: d LQ2:各位高手:我做分析时将terminate time设为0.006s时没有负体积出现，但一旦延长求解终止时间就会出现负体积，比如将terminate time设为0.01s时就会出现负体积，请问是什么原因呢，怎么解决好呢？+ a4 \. D6 C9 g1 j1 把ADMAP的参数值设置为0.1，在材料属性里面设置。