找回密码
 注册
Simdroid-非首页
查看: 395|回复: 25

[接触碰撞] 求助泡沫填充管冲击模拟

[复制链接]
发表于 2012-12-29 14:53:37 | 显示全部楼层 |阅读模式 来自 四川乐山
悬赏30仿真币已解决
小弟最近模拟了一个泡沫铝填充薄壁管的冲击,发现结构的变形不对称,并且在模拟过程中由于负体积问题而导致终止。急!!!求解决,麻烦各位帮忙看下!!!
附件: 您需要 登录 才可以下载或查看,没有账号?注册

最佳答案

查看完整内容

哦, 我这里CPU不空, 大概说个思路,你来试试看 1 trigger 不要, 那个肯定破坏压溃模式 2 foam 的部分possion比设0 , 让它不要膨胀开来 3 铝外壳的部分,设置好single surface的接触厚度, 5mm 4 给铝壳材料曲线的最后加上一个强硬化, 要不然肯定压溃成纸一样,弯曲面出不来 加一个网格细分的 control adaptive. ...

点评

变模式不对,建议增加一些诱导槽  发表于 2012-12-30 11:47
发表于 2012-12-29 14:53:38 | 显示全部楼层 来自 上海
Simdroid开发平台
哦, 我这里CPU不空, 大概说个思路,你来试试看

1 trigger 不要, 那个肯定破坏压溃模式
2 foam 的部分possion比设0 , 让它不要膨胀开来
3 铝外壳的部分,设置好single surface的接触厚度, 5mm
4 给铝壳材料曲线的最后加上一个强硬化, 要不然肯定压溃成纸一样,弯曲面出不来
  加一个网格细分的 control adaptive.
回复

使用道具 举报

 楼主| 发表于 2012-12-30 13:54:20 | 显示全部楼层 来自 四川乐山
版主,求教关于泡沫材料的负体积解决方法,这个问题已经困扰近一个多星期了,也看了论坛上关于负体积的解决办法,但是都没能够解决
回复

使用道具 举报

发表于 2013-1-3 08:24:36 | 显示全部楼层 来自 湖南长沙
Negative volume in soft materials

In materials that undergo extremely large deformations, such as soft foams, an element may become so distorted that the volume of the element is calculated as negative. This may occur without the material reaching a failure criterion. There is an inherent limit to how much deformation a Lagrangian mesh can accommodate without some sort of mesh smoothing or remeshing taking place. A negative volume calculation in LS-DYNA will cause the calculation to terminate unless ERODE in *CONTROL_TIMESTEP is set to 1 and DTMIN in *CONTROL_TERMINATION is set to any nonzero value in which case the offending element is deleted and the calculation continues (in most cases). Even with ERODE and DTMIN set as described, a negative volume may cause an error termination. Some approaches that can help to overcome negative volumes include the following:

1. Simply stiffen up the material stress-strain curve at large strains. This approach can be quite effective.

2. Sometimes tailoring the initial mesh to accomodate a particular deformation field will prevent formation of negative volumes. Again, negative volumes are generally only an issue for very severe deformation problems and typically occur only in soft materials like foam.

3. Reduce the timestep scale factor. The default of 0.9 may not be sufficient to prevent numerical instabilities.

4. Avoid fully-integrated solids (formulations 2 and 3) which tend to be less stable in situations involving large deformation or distortion. (The fully integrated element is less robust than a 1-point element when deformation is large because a negative Jacobian can occur at one of the integration points while the element as a whole maintains a positive volume. The calculation with fully integrated element will therefore terminate with a negative Jacobian much sooner than will a 1-point element. (lpb))

5. Use the default element formulation (1 point solid) with type 4 or 5 hourglass control (will stiffen response).Preferred hourglass formulations for foams are type 6 with hourglass coef = 0.1 if low velocity impact, and type 2 or 3 with hourglass coef. = 0.1 if high velocity impact. The hourglass type and hourglass coefficient may warrant modification based on observed hourglass modes and on reported hourglass energy in matsum

6. Model the foam with tetrahedral elements using solid element formulation 10 although this approach may give an overly stiff response.

7. Increase the DAMP parameter (foam model 57) to the maximum recommended value of 0.5.

8. Use optional card B of *CONTACT to turn shooting node logic off for contacts involving foam.

9. Use *CONTACT_INTERIOR. A part set defines the parts to be treated by *CONTACT_INTERIOR. Attribute 4 (DA4 = 5th field of Card 1) of the part set defines the TYPE of *CONTACT_INTERIOR used. The default TYPE is 1 which is recommended for uniform compression. In version 970, solid formulation 1 elements can be assigned TYPE=2 which treats combined modes of shear and compression.

10. If *MAT_126 is used, try ELFORM = 0.

11. Try EFG formulation (*SECTION_SOLID_EFG). Use only where deformations are severe as this formulation is very expensive. Use only with hex elements.

--------From LS-DYNA Support "http://www.dynasupport.com/howtos/material/negative-volume-in-soft-materials"
回复

使用道具 举报

 楼主| 发表于 2013-1-3 09:30:00 | 显示全部楼层 来自 四川成都
本帖最后由 风龙 于 2013-1-3 09:42 编辑
goalyhfeng 发表于 2013-1-3 08:24
Negative volume in soft materials

In materials that undergo extremely large deformations, such as  ...

谢谢goalyfeng的热心,你贴的这些方法我都试过了,但是没有从根本上解决负体积的问题,可能我还是有些地方设置不对吧

网络卡了一下,多回复了几次,怎么删掉多余的回复?

点评

点击 编辑,里面右侧有删除帖子。  发表于 2013-1-3 12:30
回复

使用道具 举报

 楼主| 发表于 2013-1-4 15:08:43 | 显示全部楼层 来自 四川乐山
xiaopiujy版主,我加了诱导槽之后,计算的时候还是会发生负体积,外面的壳单元变形模式和单独的壳单元时候不一样,这是怎么会事啊?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

发表于 2013-1-4 16:35:10 | 显示全部楼层 来自 上海
大概是你泡沫表面铝材料的定义有问题

变形模式是应该这样的吗?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

 楼主| 发表于 2013-1-4 18:39:41 | 显示全部楼层 来自 四川乐山
本帖最后由 风龙 于 2013-1-4 18:40 编辑
necrosiscai 发表于 2013-1-4 16:35
大概是你泡沫表面铝材料的定义有问题

变形模式是应该这样的吗?

恩,你计算的过程有没有出现那负体积?
回复

使用道具 举报

发表于 2013-1-4 23:14:32 | 显示全部楼层 来自 上海
风龙 发表于 2013-1-4 18:39
恩,你计算的过程有没有出现那负体积?

没有,把你的铝材料定义修改一下

LCID 2 被用了两次, 当中有个是好像应变速度有关的,把那个曲线去掉
回复

使用道具 举报

 楼主| 发表于 2013-1-5 10:20:57 | 显示全部楼层 来自 四川乐山
necrosiscai 发表于 2013-1-4 23:14
没有,把你的铝材料定义修改一下

LCID 2 被用了两次, 当中有个是好像应变速度有关的,把那个曲线去掉 ...

忘记这个是先前传的K文件,你把那个termination time改为0.02看下,这样好像会出现负体积,0.01的时候刚体的动能只减少了一半左右
回复

使用道具 举报

发表于 2013-1-5 13:15:01 | 显示全部楼层 来自 上海
没问题,

附上K文件

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

 楼主| 发表于 2013-1-5 16:17:27 | 显示全部楼层 来自 四川乐山
necrosiscai 发表于 2013-1-5 13:15
没问题,

附上K文件

control timestep 的TSSFAC设置为0.3是不是有点太小了?这样会增加很多的计算时间啊,SNLOG设置为1有没有什么道理?
回复

使用道具 举报

 楼主| 发表于 2013-1-8 10:19:38 | 显示全部楼层 来自 四川乐山
necrosiscai 发表于 2013-1-5 13:15
没问题,

附上K文件

necrosiscai ,我在模型的前端加了诱导槽之后,计算的时候会出现负体积,怎么解决啊,附上K文件,做诱导槽的时候,泡沫材料要不要因为外面的壳单元而做成一个诱导槽?急求

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

发表于 2013-1-8 12:39:00 | 显示全部楼层 来自 上海
载荷均匀分布的话, 负体积不容易出现, 你这样应力集中以后, 不出负体积我到感觉奇怪了

非得做诱导槽吗? 等我看看K文件
回复

使用道具 举报

发表于 2013-1-8 12:43:05 | 显示全部楼层 来自 上海
我终于明白你说加trigger以后算得慢是为什么了
你trigger的地方元素很小, 最小timestep 1e-7 .你又没在contrl里面设定, 那是要很慢的
回复

使用道具 举报

发表于 2013-1-8 12:57:31 | 显示全部楼层 来自 上海
风龙 发表于 2013-1-8 10:19
necrosiscai ,我在模型的前端加了诱导槽之后,计算的时候会出现负体积,怎么解决啊,附上K文件,做诱导 ...

你画网格熟练吗? 个人建议你在trigger的地方用三角形和正方形来代替现在的梯形
就是5面体和6面体混合网格, 你现在的网格质量不太好

你要不熟练的话, 把foam里面材料曲线的最后拉一个大阳线,让它硬起来
同时设定timestep大概在1e-6的样子便于计算 control time step 里面大概第四项 -1e-6

评分

1

查看全部评分

回复

使用道具 举报

 楼主| 发表于 2013-1-9 10:49:40 | 显示全部楼层 来自 四川乐山
necrosiscai 发表于 2013-1-8 12:57
你画网格熟练吗? 个人建议你在trigger的地方用三角形和正方形来代替现在的梯形
就是5面体和6面体混合网 ...

重新划分了网格及设定了质量缩放系数都会发生因负体积而终止计算,

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

发表于 2013-1-9 12:25:57 | 显示全部楼层 来自 上海
风龙 发表于 2013-1-9 10:49
重新划分了网格及设定了质量缩放系数都会发生因负体积而终止计算,

材料曲线改了吗?
回复

使用道具 举报

 楼主| 发表于 2013-1-9 21:13:11 | 显示全部楼层 来自 四川乐山
本帖最后由 风龙 于 2013-1-9 21:13 编辑
necrosiscai 发表于 2013-1-9 12:25
材料曲线改了吗?

修改了,材料曲线修改了的话,会不会对结果有很大的影响?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

发表于 2013-1-11 15:38:14 | 显示全部楼层 来自 上海
挖坟, 要实现什么样的变性模式? 有个照片什么的吗?
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-5-15 10:36 , Processed in 0.054992 second(s), 16 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表