隐式-显示 初始应力岩石爆破模拟

惆怅依旧

在ansys给岩石两端加上压力进行隐式静力分析，导出隐式分析生成的drelax文件，显式分析时导入到ls-dyna中计算，爆破过程可以计算，但是计算完成后出现了不正常的情况，如下 是1/2模型                     

           0时刻，岩石初始应力场正确                                                                爆破5us时应力云图，问题就出在这，怎么两端有应力波往里传呢，按理说不应该啊   
         显示分析的时候除了炸药没有再加其他的外载荷，大家帮忙看看问题出在哪？考虑挺长时间了，还是不得其解。

        附件太大，上传不了，只能把主要的关键字给贴出来了
     *CONTROL_DYNAMIC_RELAXATION
       250  0.001000  0.995000     0.000     0.000         0  0.040000         2

*SECTION_SOLID
         1         1         0
*SECTION_SOLID
         2         1         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_ADD_EROSION
         1 12350.000     0.000     0.000     0.000     0.000     0.000     0.000
-2.000E-4 3.9000E-5 12350.000 12350.000 12350.000 12350.000 12350.000     0.000
*MAT_PLASTIC_KINEMATIC
         1  0.002600  0.051800  0.3300+00 7.5000E-5  0.004000     0.000
     0.000     0.000  1.250000     0.000
*MAT_HIGH_EXPLOSIVE_BURN
         2  0.001000  3.000000  0.00000     0.000     0.000     0.000     0.000
*EOS_JWL
         2   0.15000 2.2300E-3  4.150000  0.950000  0.300000 2.5000E-3  1.000000
*INITIAL_DETONATION
         2     0.000     0.000   130.000     0.000
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         1 and Elem Type         1
         1         1         1         0         0         0         0         0
*PART
Part          2 for Mat         2 and Elem Type         1
         2         2         2         2         0         0         0         0

*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_TIMESTEP
    0.0000    0.9000         0  0.00      0.00   
*CONTROL_TERMINATION
  100.             0   0.00000   0.00000   0.00000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_BINARY_D3PLOT
1.000   
*DATABASE_BINARY_D3THDT
1.000   
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         4         0         0         0
其他的地方都是节点单元 无反射信息了


还请大家帮忙看看

我的橙子

楼主，显式分析读取drelax文件后，初始应力与静力分析时的不同，位移相同，这是哪里错了吗？

zhuxiangy

楼主可否做个教程啊？用这个实例讲一讲ansys隐显式顺序分析过程，谢谢

惆怅依旧

经过几天的尝试，前几天终于把初始压力和无反射边界条件加载一块了，用完全重启动可以实现，先把外载荷线性加到所需的值，在保持一段时间，是的整个模型的应力场保持稳定，然后完全重启动分析，模型继承上个分析中的应力场（外载荷还是得加上哦），现在加上无反射边界条件，就不会出现隐式-显式方法出现的初始应力场和无反射边界条件的扰动问题了。

惆怅依旧

有做过类似的吗? 求帮忙啊

zhenghonghao88

这个是因为边界没有设计投射边界，导致应力波发生反射造成的，加上透射边界就好了！能不能把k文件给我发一份过来啊，邮箱jbn111@126.com，谢谢了！

惆怅依旧

zhenghonghao88 发表于 2014-1-7 22:44
这个是因为边界没有设计投射边界，导致应力波发生反射造成的，加上透射边界就好了！能不能把k文件给我发一 ...

无反射边界条件我是加了的，那个不是应力波的反射吧，应力波还没传到那呢

惆怅依旧

好像是因为隐式的初始应力场在显式分析中没有一直保持下去 然后造成了扰动  谁有比较好的解决办法呢  真的是不知道怎么办了

惆怅依旧

经过尝试，好像有个解决办法了，在显式分析中，隐式分析加的外载荷不能被撤掉，还要继续在显式中加上，而且define curve中有个代表瞬态还是动力松弛的那个关键字要设置成1.

惆怅依旧

现在又出现了问题，只要一加上无反射边界条件，就还会出现边界的扰动，看来无反射边界跟初始应力两个还是存在问题

n3ch

发K文件到n3ch@163.com

zhuxiangy

你好，看你的做法也是在Ansys中进行隐式分析后导出K文件，再修改K文件进行解算的吧？我也是这样做的，但进行显式分析时却没有初始应力存在，不知道什么原因，我的步骤是这样的：
1 隐式分析
2 隐式分析结束后
Main Menu>Solution>Solve
Main Menu>Finish
存储隐式分析数据到Jobname.DB文件
3 将作业名改为Jobname2，以防止覆盖隐式结果文件。
4 重新进入前处理： /PREP7
5 将隐式单元转换到对应的显式单元类型
6 重新对显式单元定义关键选项、实常数、材料性质、边界条件和载荷值
7 从隐式计算结果文件读入节点位移（.rst文件）、转角和温度，并将此数据记入ANSYS/LS-DYNA文件“drelax”。
Main Menu>Solution>Cosraints>Read Disp
8 根据存储在“drelax”文件中的位移、转角将结构初始化为给定几何构形
Main Menu>Solution>Analysis Options>Dynamic Relax
9 设置显式分析选项导出K文件
10 修改K文件进行求解

不知道其中哪个地方与你的不同，急切盼望你的回复，谢谢

惆怅依旧

zhuxiangy 发表于 2014-2-21 22:13
你好，看你的做法也是在Ansys中进行隐式分析后导出K文件，再修改K文件进行解算的吧？我也是这样做的，但进 ...

你可以搜下隐式-显示相关帖子，有的帖子里有详细的过程，照着教程一步步来应该就没问题了

zhuxiangy

惆怅依旧 发表于 2014-2-23 10:19
你可以搜下隐式-显示相关帖子，有的帖子里有详细的过程，照着教程一步步来应该就没问题了 ...

你在做静力学分析的时候用的是solid185单元吧？施加的力是多大的？有没有检查过静力学分析后的结果与施加的力之间的量级相差多大？

惆怅依旧

zhuxiangy 发表于 2014-2-24 09:38
你在做静力学分析的时候用的是solid185单元吧？施加的力是多大的？有没有检查过静力学分析后的结果与施加 ...

是，静力学分析后的结果和施加的力基本上是一致的啊

zhuxiangy

惆怅依旧 发表于 2014-2-24 18:48
是，静力学分析后的结果和施加的力基本上是一致的啊

好的，谢谢，我知道问题所在了，是约束设置的事，已经解决了

zhuxiangy

惆怅依旧 发表于 2014-1-15 21:54
经过几天的尝试，前几天终于把初始压力和无反射边界条件加载一块了，用完全重启动可以实现，先把外载荷线性 ...

你好，”先把外载荷线性加到所需的值，在保持一段时间“这个操作是指的隐式分析么？还是隐式分析结束后显式分析开始前单独进行的一个步骤？“保持一段时间”是怎么实现的啊？谢谢

惆怅依旧

zhuxiangy 发表于 2014-2-27 16:11
你好，”先把外载荷线性加到所需的值，在保持一段时间“这个操作是指的隐式分析么？还是隐式分析结束后显 ...

是显式分析，就是定义一个上线性上升再保持不变的载荷曲线。

geguangzhou

惆怅依旧 发表于 2014-1-15 21:54
经过几天的尝试，前几天终于把初始压力和无反射边界条件加载一块了，用完全重启动可以实现，先把外载荷线性 ...

亲 我也做到你这一步了  但是我调用的处理器 输入i=jobname2.k m=drelax  在计算过程中炸药单元出现畸形  小启动和完全重启动能解决问题 但是怎么输入呢   你显示分析的过程中 用的什么求解的呢？？

ansys/ls-dyna

惆怅依旧 发表于 2014-2-27 19:02
是显式分析，就是定义一个上线性上升再保持不变的载荷曲线。

那为什么还要做前面隐式分析的一步啊，这样直接在显示分析中做不就好了，有啥区别？