[分享]地震波的输入和求解

cxb305

!对于地震波的输入，可以把荷载记录做成文件，利用apdl的读取功能读入倒数据库中。下面的例子是自己编的一个小文件。修改一下可以更简洁。有用到的朋友自己作一下把。

1. fini
   
2. /config,nres,1000
   
3. *dim,aceX,TABLE,3000,1
   
4. *dim,aceY,TABLE,3000,1
   
5. *dim,aceZ,TABLE,3000,1
   
6. 
   
7. *creat,ff
   
8. *vread,aceX(1,1),acex,txt,,1
   
9. (e16.6)
   
10. *vread,aceX(1,0),ACETT,,,1
    
11. (e17.6)
    
12. ACEX(0,1)=1
    
13. *end
    
14. /input,ff
    
15. 
    
16. *creat,ff
    
17. *vread,aceY(1,1),acey,txt,,1
    
18. (e16.6)
    
19. *vread,aceY(1,0),ACETT,,,1
    
20. (e17.6)
    
21. ACEY(0,1)=1
    
22. *end
    
23. /input,ff
    
24. 
    
25. *creat,ff
    
26. *vread,aceZ(1,1),acez,txt,,1
    
27. (e16.6)
    
28. *vread,aceZ(1,0),ACETT,,,1
    
29. (e17.6)
    
30. ACEZ(0,1)=1
    
31. *end
    
32. /input,ff  

复制代码

！地震波时程记录分成了3个文件，每个文件是一列。分别记录x,y,z方向的加速度。acett是时间记录。
这样就可以把加速度记录读取倒ansys数据库中作为数组。
也可以把加速度记录做成一个文件，这样程序就简单多了。大家可以试看看修改一下。
下面是计算部分语句：

1. /SOLU
   
2. 
   
3. ANTYPE,trans
   
4. ！求解其自己选了
   
5. TM_START=0.01
   
6. TM_END=15.00
   
7. TM_INCR=0.01
   
8. *DO,TM,TM_START,TM_END,TM_INCR
   
9. TIME,tm
   
10. alpha,
    
11. BETAD,
    
12. ACEL,acex(tm),acey(tm),acez(tm)
    
13. SOLVE
    
14. *ENDDO
    
15. fini
    

复制代码

vigorao

在地震输入中怎样考虑行波效应呢？

sysh320

我正在用ANSYS做抗震分析，有用！

AaronSpark

1）是正弦波的地震波反应

1. /solu
   
2. *dim,baseace1,array,600 ! 定义正弦波的“基础加速度”, 500时间点
   
3. *do,i,1,600
   
4. baseace1(i)=0.1*sin(i*2) !
   
5. *enddo
   
6. NSUBST,1, , ,1 !1个子步
   
7. OUTRES,ALL,1 !输出每个子步的结果
   
8. ANTYPE,TRANS
   
9. 
   
10. CSYS,1
    
11. NSEL,S,LOC,X,13000,13000
    
12. CSYS,0
    
13. *do,i,1,600
    
14. f,all,fx,3e4*baseace1(i)
    
15. f,all,fz,3e4*baseace1(i)
    
16. TIME,i*0.02
    
17. solve
    
18. *enddo

复制代码

2）做z向行波分析

1. /nerr,,99999999
   
2. fini
   
3. /solu
   
4. *dim,baseace1,array,300 ! 定义正弦波的“基础加速度”, 500时间点
   
5. *dim,baseace2,array,300
   
6. *do,i,1,300
   
7. baseace1(i)=0.1*sin(i*0.2)
   
8. baseace2(i)=0.1*sin(i*0.2-31)
   
9. *enddo
   
10. 
    
11. ANTYPE,TRANS !时程分析
    
12. CSYS,1
    
13. NSEL,S,LOC,X,13000,13000
    
14. CSYS,0
    
15. *do,i,1,300
    
16. f,all,fx,3e4*baseace1(i)
    
17. f,all,fz,3e4*baseace1(i) ! F=m*a
    
18. TIME,i*0.02
    
19. *enddo
    
20. 
    
21. NSUBST,300, , ,1 !500个子步
    
22. OUTRES,ALL,1 !输出每个子步的结果
    
23. LSWRITE
    
24. 
    
25. !在Z方向加行波效应
    
26. ANTYPE,TRANS
    
27. CSYS,1
    
28. NSEL,S,LOC,X,13000,13000
    
29. CSYS,0
    
30. *do,j,1,300
    
31. f,all,fz,3e4*baseace2(j)
    
32. f,all,fx,0
    
33. TIME,(j+301)*0.02
    
34. *enddo
    
35. 
    
36. NSUBST,300, , ,1 !500个子步
    
37. OUTRES,ALL,1 !输出每个子步的结果
    
38. LSWRITE
    
39. 
    
40. LSSOLVE,1,2

复制代码

sysh320

感谢阿！

shenliying

我也要做地震分析 谢谢拉！

wawa

谢谢！

sunjianpeng2001

如何获得位移地震波

sunjianpeng2001

有没有位移地震波啊

yangmingyu

我也正在考虑行波效应，谢谢了！

dreamskyweb

谢谢大侠的代码！！！！！！！！！！！！！