地震波输入的问题

F_watcher

请问各位大侠，在输入地震加速度时程时，加速度是作用在节点上、单元上、面上还是作用在实体上？另外，在ansys中，如何指定加速度的作用方向，例如，我要沿X方向输入地震加速度时程，该怎么操作，请给予指点，先行谢过！！！

shijin0301

首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：

   -0.227109E-02   -0.209046E+00    0.467072E+01
   -0.413893E-02   -0.168195E+00    0.261523E+01
   -0.574753E-02   -0.157890E+00    0.809014E-01
   -0.731227E-02   -0.152996E+00    0.119975E+01
   -0.876865E-02   -0.138102E+00    0.130902E+01
   -0.101067E-01   -0.131582E+00    0.143611E+00
.......................

然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：

    0.100000E-01
    0.200000E-01
    0.300000E-01
    0.400000E-01
    0.500000E-01
    0.600000E-01
.......................

编写如下的命令流文件，并命名为acce.inp

1. *dim,ACCEXYZ,TABLE,2000,3                     !01行
   
2. *vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000   !02行(3e16.6)                                      !03行
   
3. *vread,ACCEXYZ(1,0),time,txt                  !04行
   
4. (e16.6)                                       !05行
   
5. ACCEXYZ(0,1)=1                                !06行
   
6. ACCEXYZ(0,2)=2                                !07行，同上
   
7. ACCEXYZ(0,3)=3                                !08行，同上
   
8. finish
   
9. /SOLU
   
10. ANTYPE,trans
    
11. btime=0.01               !定义计算起始时间
    
12. etime=15.00              !定义计算结束时间
    
13. dtime=0.01               !定义计算时间步长
    
14. *DO,itime,btime,etime,dtime
    
15. time,itime
    
16. AUTOTS,0
    
17. NSUBST,1, , ,1
    
18. KBC,1
    
19. acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3)     !施加三个方向的地震加速度
    
20. SOLVE
    
21. *ENDDO

复制代码

最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。

说明和讨论：
1、命令流中各行说明：
01行：定义2000行，3列的数组，(行数根据数据文件里加速度点数来定)ACCEXYZ用来存放三个方向的加速度值。
02行：从数据文件accexyz.txt里读加速度值到数组ACCEXYZ,2000为行数，可根据情况修改。其中的JIK，3，2000非常重要，它决定着将从加速度数据文件中数据输入到ACCEXYZ数组时的存放格式，这里用到按JIK方式变化,即读进来的数据依次放入ACCEXYZ(I,J)（J从1到3，I从1到2000，K默认从1到1）。根据这行命令下面(3e16.6)的格式，每次从数据文件里读一行三个数据，分别放入ACCEXYZ(1,1)、ACCEXYZ(1,2)、ACCEXYZ(1,3)，接着再读下一行的三个数据，分别放入ACCEXYZ(2,1)、ACCEXYZ(2,2)、ACCEXYZ(2,3)，依次类推。当然依据数据文件的格式，也可以选用IJK, IKJ, JIK, JKI, KIJ, KJI等格式，其中IJK为默认。

03行：读数据的格式，每行三个数值。由于数据文件中的数据是用科学记数年法表示的，因此，这里也用相应FORTRAN的科学记数法的格式。如果数据文件里的数值是如“0.2876   1.2333   2.9938”这样的格式，此行的数据格式也就修改为“(3f10.4)”这样的格式。
04行：从数据文件time.txt里读时间值到数组ACCEXYZ的第零列。
06行：将数组ACCEXYZ的第零行赋值，如果不对行插值的话也可以不赋值

F_watcher

十分感谢，请问您有qq吗？有的话留一个吧，日后有问题好再次请教。

笑话大哥

太好了，雪中送资料呀

笑话大哥

十分感谢，请问您有qq吗？有的话留一个吧，日后有问题好再次请教。

风中的誓言

shijin0301 发表于 2009-2-13 16:16
首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为 ...

你好，能加你QQ吗我也是研究地震作用的想和你讨论一些问题，我的qq是460645499