(求助)看看我的流固耦合的效果.

我就是大头

我做的是炸药在岩石介质中爆炸,用的是多物质流固耦合算法.空气单元与炸药共节点.空气和岩石单元搭接.计算结果好象有问题,请高手指教.

我就是大头

顺序发错了,按照图片的名称看吧.
我在计算结果中发现,地表的震动速度不是随距离爆源的距离增大而减小,而是一个震荡的过程.



$
$
*SECTION_SOLID
         1         1
*SECTION_SOLID_ALE
$      SID    ELFORM       AET
         2        11
$     AFAC      BFAC      CFAC      DFAC     START       END     AAFAC
0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
*SECTION_SOLID_ALE
$      SID    ELFORM       AET
         3        11
$     AFAC      BFAC      CFAC      DFAC     START       END     AAFAC
0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_PLASTIC_KINEMATIC
         1      2700 0.600E+11  0.300000 5.000E+07 0.500E+09      0.05   
      0.00      0.00      0.10   
*MAT_HIGH_EXPLOSIVE_BURN
$      MID        RO         D       PCJ      BETA         K         G      SIGY
         2 1.0000E+3 0.3600E+4 0.361E+10 0.0000000 0.0000000 0.0000000 0.0000000
*EOS_JWL
$    EOSID         A         B        R1        R2      OMEG        E0        V0
         1 0.476E+11 0.592E+09 3.5000000 0.9000000 1.3000000 0.045E+11 0.0000000
$         
$air from landmine (chen 1.225)      
*MAT_NULL                              
         3   1.29000      0.00      0.00      0.00      0.00      0.00      0.00   
*EOS_LINEAR_POLYNOMIAL
         2   -1.0E-6      0.00      0.00      0.00       0.4       0.4      0.00   
   0.25000    1.0E+0
*INITIAL_DETONATION
$      PID         X         Y         Z        LT
         2       0.0     -0.50      0.02       0.0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              ALE CONTROL                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*CONTROL_ALE
$      DCT      NADV      METH      AFAC      BFAC      CFAC      DFAC      EFAC
         2         1         2-1.0000000 0.0000000 0.0000000 0.0000000
$    START       END     AAFAC     VFACT    VLIMIT       EBC
0.0000000 0.0000000 0.0000000
$
*CONSTRAINED_LAGRANGE_IN_SOLID
$    SLAVE    MASTER     SSTYP     MSTYP     NQUAD     CTYPE     DIREC     MCOUP
         1         2         0         0         3         5         1         0
$    START       END      PFAC      FRIC
         0    1.0e10       0.1      0.00
                                    
*SET_PART_LIST
         1
         1
*SET_PART_LIST
         2
         2         3         
*ALE_MULTI-MATERIAL_GROUP
         2         1
         3         1                 
*HOURGLASS
         1         0     0.100         0      1.50 0.600E-01      0.00      0.00           
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         1 and Elem Type         1
         1         1         1         0         1         0         0
$
*PART
Part          2 for Mat         2 and Elem Type         2
         2         2         2         1         1         0         0
$
*PART
Part          3 for Mat         3 and Elem Type         3
         3         3         3         2         1         0         0
$
$

liushuangq05

你能不能把k文件上传一下.
地震波本身就是一种震荡波,波本身就是一种周期运动.当你计算的时间段只有周期的1/4时,那应该没有什么波动效果.因为它不是抛物线或者椭圆方程.

danlex

大头兄，我现在也要接手一个地震局的一个项目，作地震对建筑物倒塌影响的模拟。希望能和你多多交流，我的qq是419258365，我的电子邮箱是danlex@163.com