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

利用Ansys/Ls-Dyna模拟的若干问题

[复制链接]
发表于 2003-11-8 16:08:41 | 显示全部楼层 |阅读模式 来自 陕西西安
本帖最后由 gdyu_yu 于 2010-7-12 21:54 编辑

眼看毕业日期临近,好坏出不了结果!
一、如何在Ansys/Ls-Dyna中加对称边界?
    我只在Solution->constrains->Apply->Symmetry B.C.-Bndry Plane 下看到施加滑移对称(Sliding symmetry)和
循环对称(cyclic symmetry)的菜单选项,其中的Xc,Yc,Zc分别代表什么意思呢?
    比如我的模型是外圆环土体,内圆炸药(很小一块),现只取其1/4做模型,
其中左边沿y轴对称,下边沿x轴对称,那Xc,Yc,Zc该怎样写?
    我用的是cyclic symmetry,
*BOUNDARY_CYCLIC
-1.0000    -1.0000    0.0000         1         2
(分别是左边和下边的Node set ID)
这样对不对?
  
二、负体积的问题
    看了以前的帖子,说解决办法有1、减小时间步长; 2,增大罚因子; 3,控制网格质量,例如在承受压力的单元在受压方向比其他方向尺寸长。
    现在我的时间步长因子是0.5,更小的也试过,还是不行!网格质量
我也试了,因为我的模型如下图,无论网格疏密,还是受压单元的受压
方向的尺寸,都还出现负体积。
  
三、材料参数
    最终我怀疑到这方面来,我用的都是国际单位制,几乎所有爆炸的
例子都是用克-厘米-微妙-兆巴为单位制的,对炸药的eos_jwl卡,
参数例如A、B,我这里的都是GPa,换算成兆巴后,跟其它能算的k文件
相比,其值非常小,我老觉得这不正常。(曾换算成克-厘米……的单位
只是算过,不成功,所以又换了回来)。
     
    还有土材料的参数,我用5号材料,老感觉着参数有点问题,请高手看看,是不是真的有些不对劲?
  
还有什么问题,望各位老大提出!

本帖子中包含更多资源

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

×
 楼主| 发表于 2003-11-8 16:14:07 | 显示全部楼层 来自 陕西西安

回复: 【注意】弄了一个多月了,还是问题不断,麻烦高手们看看吧,受不了了!

Simdroid开发平台
k文件如下
*KEYWORD
*TITLE
                                                                                
$
*DATABASE_FORMAT
         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               NODE DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*NODE
       1 2.500000000E+00 0.000000000E+00 0.000000000E+00       0       0
       2 1.530757942E-16 2.500000000E+00 0.000000000E+00       0       0
       3 2.478611299E+00 3.263219696E-01 0.000000000E+00       0       0
       4 2.414814566E+00 6.470476128E-01 0.000000000E+00       0       0
       5 2.309698831E+00 9.567085811E-01 0.000000000E+00       0       0
       6 2.165063509E+00 1.250000000E+00 0.000000000E+00       0      
……
  
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             SECTION DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SHELL
         1        14    1.0000       1.0       0.0       0.0         0         1
  0.00      0.00      0.00      0.00      0.00     
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*MAT_HIGH_EXPLOSIVE_BURN
         1      1000      4300  3.85E+09      0.00      0.00      0.00      0.00
*EOS_JWL
         1  0.42E+11  0.44E+09      3.55      0.16      0.41  3.15E+09      1.00
*INITIAL_DETONATION
         1         0         0         0         0
*MAT_SOIL_AND_FOAM
         2      1800 1.601E+07 1.328E+13 5.479E+08 9.035E+04    0.4104       0.0
       0.0       0.0
       0.0      0.05      0.09      0.11      0.15      0.19      0.21      0.22
      0.25      0.30
       0.0  3.42E+09  4.53E+09  6.76E+09  1.27E+10  2.08E+10  2.71E+10  3.92E+10
  5.66E+10  1.23E+11   
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         2 and Elem Type         1
         1         1         1         1         0         0         0
$
*PART
Part          2 for Mat         1 and Elem Type         1
         2         1         2         0         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DISCRETE MATERIALS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             ELEMENT DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*ELEMENT_SHELL
       1       1     123       1       3     402
       2       1     400       4       5     401
       3       1     398       6       7     399
……
  
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              COORDINATE SYSTEMS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*BOUNDARY_NON_REFLECTING_2D
         3
*SET_NODE_LIST
         1     0.000     0.000     0.000     0.000
         2         8         9        10        20        26        27        28
        29        30        31        32        33        34        35        36
        37        38       118       125       127       129       131       133
       135       137       139       141       143       145       147       149
       152
*SET_NODE_LIST
         2     0.000     0.000     0.000     0.000
         1         8        11        12        19        39        40        41
        42        43        44        45        46        47        48        49
        50        51       123       126       128       130       132       134
       136       138       140       142       144       146       148       150
       157
*BOUNDARY_CYCLIC
   -1.0000   -1.0000    0.0000         1         2
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               LOAD DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                RIGID BOUNDRIES                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SET_part_list
        1
1
*SET_part_list
        2
2
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             BOUNDARY DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SET_NODE_LIST
         3     0.000     0.000     0.000     0.000
        19        20        21        22        23        24        25       151
       153       154       155       156       158
*BOUNDARY_SPC_SET
         3         0         1         1         0         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTACT DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTACT_2D_AUTOMATIC_SURFACE_TO_SURFACE
         2         1
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTROL OPTIONS                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_TIMESTEP
    0.0000    0.5000         0  0.00      0.00     
*CONTROL_TERMINATION
  10.0             0   0.00000   0.00000   0.00000
*DATABASE_HISTORY_NODE
         1                                                  
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_GLSTAT   
0.1000E-01
*DATABASE_BINARY_D3PLOT
0.1000     
*DATABASE_BINARY_D3THDT
0.1000E-01
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         4         0         0         0
*END

本帖子中包含更多资源

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

×
发表于 2003-11-8 16:40:14 | 显示全部楼层 来自 新加坡

回复: 【注意】弄了一个多月了,还是问题不断,麻烦高手们看看吧,受不了了!

你的土模型有点问题。你换成空气的看看?可能就可以了
发表于 2003-11-8 18:18:55 | 显示全部楼层 来自 台湾

回复: 【注意】弄了一个多月了,还是问题不断,麻烦高手们看看吧,受不了了!

MAT_HIGH_EXPLOSIVE_BURN  
MAT_SOIL_AND_FOAM  
這兩個材模可用在實體元素不能用在殼元素呀
這個應該對你有幫助
是3d

本帖子中包含更多资源

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

×
 楼主| 发表于 2003-11-9 09:14:37 | 显示全部楼层 来自 陕西西安

回复: 【注意】弄了一个多月了,还是问题不断,麻烦高手们看看吧,受不了了!

to  kerlit
不好意思,能说清楚点吗?我搞得就是土中爆炸,怎么换成空气呀?是不是要弄成不耦合装药?
  
to ye01.chen
  
我用的是solid162(plant)单元,这个也不行吗?那我再弄个3D的试试!
  
还有一个Femb的问题,在Femb中,网格划分怎样画成上述图形所示的呢,我画出的单元尺寸总是平均的,有没有象ansys中的那个比例因子的?
发表于 2003-11-9 18:12:31 | 显示全部楼层 来自 新加坡

回复: 【注意】弄了一个多月了,还是问题不断,麻烦高手们看看吧,受不了了!

哦,你要用3维的才有点意思。2维的太简单了点。用6面体单元,比较稳定。
不过你的爆炸和土的材料都好像没有。什么毕业?要大学 毕业了?
  
土的模型我也没有搞到,不过快了,我上次拿到的以后,正在测试它们,如果可以我就给你一份。
  
祝你快乐
 楼主| 发表于 2003-11-10 11:26:26 | 显示全部楼层 来自 陕西西安

回复: 【注意】弄了一个多月了,还是问题不断,麻烦高手们看看吧,受不了了!

哦,我想两个都算一下,2D的简单一些,就先弄它了
3D的我做做看
谢谢!:)
发表于 2003-11-11 16:28:14 | 显示全部楼层 来自 北京

回复: 【注意】弄了一个多月了,还是问题不断,麻烦高手们看看吧,受不了了!

(1)BOUNDARY_CYCLIC  
   -1.0000 -1.0000 0.0000 1 2  
是循环对称边界条件,我想你在命令流中输入DSYM,SYMM,X,和DSYM,SYMM,Y,就可以加对称条件了;
(2)负体积问题可能由你的接触方式引起CONTACT_2D_AUTOMATIC_SURFACE_TO_SURFACE  
         2 1  
个人观点,我觉得你所选的这个接触方式对你这种情况下不太理想,可以用节点对节点的接触方式,就像3D模型种的SEGMENT一样,用SET_NODE_COLUMN试试;
(3)你用2D-PLANE162把问题简化成平面应变问题是完全行得通的,不管你是2D还是3D,只要跟你的试验结合起来,就相当不错了。
发表于 2004-2-7 11:51:33 | 显示全部楼层 来自 北京理工大学

回复: 【注意】弄了一个多月了,还是问题不断,麻烦高手们看看吧,受不了了!

你做得是爆炸问题,变形较大不适合用L算法的,你用ALE算法应该就不会出现负体积了。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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

GMT+8, 2024-3-29 18:56 , Processed in 0.070047 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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