最全的负体积讨论方法.docVIP

最近被负体积整的郁闷,所以找到了很多关于负体积的资料,.现分享.顺便请教; |$ q) g% g; _% r1 k! ]关于用ANSYS做为前处理工具的??我想对模型中的实体进行包壳处理.在原来体单元上9 c) M* Y5 F0 c我在没有另外建立面的情况下对体表面进行壳划分.生成K文件编辑.但是还是会有负体积??h( M8 T8 A% L+ J! m1 O问题.请教各位大神你们关于包壳处理方法.欢迎各位加QQ群176071936探讨(整车碰撞技术群).??y. f??a( g( P+ k# I- ~% I d9 N# b1 Y, v Z) X2 T+ p6 t2 d) O. I$ T; p% M2 E O- v+ \! e7 L* U负体积定义？ Negitive volume7 f; b# r- B2 q) F, i* t; w负体积是由于element本身产生大变形造成自我体积的内面跑到外面接着被判断为负体积。8 s! {3 F% C3 z, J关于负体积的解决办法？0 c1 M/ z( I. @! f3 c负体积多是网格畸变造成的，和网格质量以及材料、载荷条件都有关系。有可能的原因和解决的方法大概有几种：0 w9 j _6 n0 j m??m D3 h$ Z（1）材料参数设置有问题，选择合适的材料模式) s v# N: T5 K P, n3 q（2）沙漏模式的变形积累，尝试改为全积分单元( c! O8 I e* m7 E+ E! N! O. h) f（3）太高的局部接触力（不要将force施在单一node上，最好分散到几个node上以pressure的方式等效施加），尝试调整间隙，降低接触刚度或降低时间步。/ H; x. Z! P. Y! F; k（4）在容易出现大变形的地方将网格refine。6 i- b8 a; w5 {（5）材料换的太软,是不是也会出现负体积!. s3 _0 Y- U7 w6 }; ~2 I% M6 E（6）另外也可以采用ALE或是euler单元算法，用流固耦合功能代替接触，控制网格质量。例如在承受压力的单元在受压方向比其他方向尺寸长。. X m; S- W h j8 u（7）尝试减小时间步长从0.9减小到0.6或更小。 ~. |2 Y! G e5 w??B经验总结：* Q3 J ? S+ A时间步长急剧变小，可能是因为单元产生了严重的畸变而导致的负体积现象，如果采用的是四面体单元，你可以用网格重划分的方法来解决。如果你采用的是六面体单元，那目前就没有很有效的方法，可以试一下*ELEMENT_SOLID_EFG，那对机器的要求相对就会比较高了。# e# l Z$ o6 B: z* {5 K) L7 w2 X+ tQ1:材料负体积解决方法（全面、有效） \9 M4 {) H _2 l, B- G% g材料负体积解决方法8 c e# e q3 j5 T2 k( S4 e* [2 e( h8 N在仿真中，通常有材料的大变形问题，如泡沫材料，由于单元大扭曲而出现了单元负体积，这种情况一般出来在材料失效之前。在没有网格光滑和网格从划分的情况下，ls-dyna有一个内部的限制来调节lagrange单元的变形。负体积一般都会导致计算中止，除非你设置时间步长控制中的erode=1和设置终止控制中的dtmin为一非零数，这种情况下，出现负体积的单元将被自动删除，计算也不会中止。不过就算你如上设置了erode与dtmin，负体积有时候也会导致计算出错停止。5 { \ L3 O* {( a- p; |8 Z; n( k: X* l1 `% O2 b1 K一些常用的解决负体积的方法如下：8 d7 I, t$ ~) y! `! P0 E在材料出现大应变的情况下增强材料的应力-应变曲线中材料应力。这种方法往往非常有效果。 Z! i j( E4 ?6 l% u??C9 a2、重新划分网格，在出现大变形的地方把网格加密。.) d4 d4 b. g2 T1 d7 h0 e, b9 i# u3、减小时间步长系数。默认的0.9系数可能不足以避免数值的不稳定。3 Z- { N3 ^, A4、避免采用全积分体单元（算法2和3），这会导致大变形和大扭曲的情况下计算相对不稳定。5 g: {# L4 O J6 C; `# m2 y5、采用默认的单元算法（单点体单元），采用沙漏控制type4和5。泡沫的沙漏控制算法为：在低速冲击问题中采用type6，系数为1；在高速冲击问题中采用type2和3。3 o1 C6 R! q) X9 u6、泡沫材料网格划分采用四面体网格，单元算法为10，虽然这样会导致材料相对比较刚性。增大材料（泡沫材料57号材料）的阻尼系数，推荐采用系数为0.5。6