混凝土抗冲击数值模拟问题

mulin

ANSYS/LS-DYNA7.0,冲击模拟简要描述：装置如图1所示，冲击杆、入射杆和透射杆为钢杆，采用线弹性模型；图1中的试件为混凝土试件，混凝土采用Nonlinear--Inelastic--Pressure Dependment--Hydrodynamic模型(即Elastic-PlasticHydrodynamic Model)。混凝土抗压强度为50MPa。更详细的信息请见附件（命令流）。
   求解后，有以下5个问题：
1、输入杆与透射杆严重变形，如图2所示。即使将钢的弹性模量扩大10倍(即2000GPa)，钢杆还是变形。
2、进行时间历程处理，无弹性应变选择项，而塑性应变总是“invalid”,如图3所示,导致无法显示求解结果混凝土试件的应力--应变曲线。
3、应力波形状奇怪，透射波居然与入射波波幅相等，而且是三角波，如图4所示。（通常来说，大多数试验结果都是矩形波）。
4、混凝土破坏是中心先破坏，而试验一般是边缘破坏。如图5所示。这是不是与扫掠网格剖分有关？
5、如何使高速冲击模拟效果出现混凝土碎块溅飞的情况？
  
    命令流(for ANSYS/LS-DYNA)与K文件(for LS-DYNA)请见附件。
    除此之外，执行求解命令后，出现警告提示：“Warning！ SOLID164 wedeges are recommended only inregions of relatively low stress gradients”。是不是应力梯度太大？这会影响求解结果吗？如何将应力
梯度变小一点？
    请高手们不吝赐教。
    我的信箱：xxiangsun@sina.com.cn
    非常感谢。
  
    图1如下所示，图2至图5请见后续贴子。

mulin

图2：输入杆与透射杆严重变形，中间网格密集的是混凝土试件。

mulin

图3：进行时间历程处理，无弹性应变选择项，而塑性应变总是“invalid”。

mulin

图4：应力波形状奇怪，透射波居然与入射波波幅相等，而且是三角波。

mulin

如图5：混凝土破坏是中心先破坏，而试验一般是边缘破坏。这是不是与扫掠网格剖分有关？

mulin

     附件：命令流(for ANSYS/LS-DYNA)
    

mulin

    附件：K文件(for LS-DYNA)请见附件。  
    如何使高速冲击模拟效果出现混凝土碎块溅飞的情况？
    除此之外，执行求解命令后，出现警告提示：“Warning！ SOLID164 wedeges are recommended only inregions of relatively low stress gradients”。是不是应力梯度太大？这会影响求解结果吗？如何将应力  
梯度变小一点？  
    图2和图3中，取的是1/4的造型，图5中是全型。

极地猎人

命令流文件在哪里，俺看看

mulin

    命令流文件就是圆圆的混凝土试件图形下面的附件啊。
  
    第一个问题与第四个问题，即输入杆与透射杆严重变形、混凝土破坏是中心先破坏的问题，已解决，其实是我粗心忘了在对称面上加约束。惭愧惭愧。
    其它问题还搞不清楚。
    另外，将混凝土试件沿半径方向剖分密集一点，比如超过输入杆与透射杆的径向划分数目，就出现多外负体积而停止计算。这该如何解决？
    谢谢！

haiyang7948

老兄能不能把你的软件和大家分享一下！！！

mulin

我所用的就是ANSYS7.0中的LS-DYNA模块，很容易得到的啊。