桥梁的地震响应分析【转】

Chelsea

1. 
   
2. /filname,earthquake,1
   
3. keyw,pr_struc,1
   
4. /prep7
   
5. !定义单元类型
   
6. et,1,beam3
   
7. !定义实常数
   
8. r,1,0.0256,0.546e-4,0.16
   
9. !定义材料属性
   
10. mp,ex,1,2.1e11
    
11. mp,prxy,1,0.3
    
12. mp,dens,1,7800
    
13. !建立几何模型
    
14. k,1
    
15. k,2,100
    
16. k,3,50,25
    
17. l,1,2
    
18. larc,1,2,3
    
19. ldiv,1,,,10
    
20. ksel,s,kp,,4,12,1
    
21. kgen,2,all,,,,30
    
22. kdele,3
    
23. 
    
24. *do,i,4,12,1
    
25. l,i,i+9
    
26. *enddo
    
27. lsel,s,line,,12,20,1
    
28. lsel,a,line,,2
    
29. lsbl,2,all
    
30. l,4,22
    
31. *do,i,5,11,1
    
32. l,i,i+19
    
33. *enddo
    
34. l,12,23
    
35. allsel
    
36. numcmp,all
    
37. !网格划分、施加边界条件
    
38. lesize,all,,,1
    
39. type,1
    
40. real,1
    
41. mat,1
    
42. lmesh,all
    
43. !边界条件
    
44. nsel,s,loc,x,0
    
45. nsel,a,loc,x,100
    
46. d,all,uy
    
47. nsel,s,loc,x,0
    
48. d,all,ux
    
49. allsel
    
50. !自重载荷
    
51. acel,,9.81
    
52. allsel
    
53. !地震谱响应分析
    
54. /solu
    
55. !模态分析
    
56. antype,2
    
57. MODOPT,SUBSP,28  !子空间法，模态提取数为28
    
58. solve
    
59. finish
    
60. !谱分析
    
61. /solu
    
62. antype,8     !spectrum
    
63. SPOPT,SPRS,28,1   !单点谱分析，参与计算的模态数为28，打开盈利响应开关
    
64. SVTYP,2,1,        !设置反应谱类型为加速度谱，
    
65. SED,1,1,0,        !激励方向为x和y
    
66. FREQ,0.333,0.357,0.416,0.5,0.625,0.833,1.25,2.5,2.778   
    
67. FREQ,3.333,3.571,4.167,5,6.25,8.333,10,12.5,16.667  
    
68. FREQ,25,50  
    
69. 
    
70. SV,0,0.013,0.013884,0.01595,0.018794,0.022974,0.029763,0.042871,0.08,0.08,  
    
71. SV,0,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.0712,0.0624,  
    
72. SV,0,0.0536,0.0448,
    
73. solve
    
74. finish
    
75. !模态扩展
    
76. /solu   
    
77. antype,2
    
78. MODOPT,SUBSP,28
    
79. mxpand,28,0,0,1,0.005  !扩展模态数设为28，计算单元的计算结果设置为yes，临界点设为0.005
    
80. solve
    
81. finish
    
82. 
    
83. !模态合并
    
84. /solu
    
85. antype,8
    
86. SRSS,0.02,DISP  !组合方式为srss，significant threshold为0.02
    
87. solve
    
88. finish
    
89. 
    
90. 
    
91. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
92. !地震波瞬态分析
    
93. 
    
94. *DIM,tjx,ARRAY,2,50,1
    
95. 
    
96. *DIM,tjy,ARRAY,2,50,1  
    
97.   
    
98. *CREATE,ansuitmp
    
99. *VREAD,tjx(1,1),'tjx','txt',' ',50
    
100. (e9.3,e11.3)
     
101. *END
     
102. /INPUT,ansuitmp
     
103. 
     
104. *CREATE,ansuitmp
     
105. *VREAD,tjy(1,1),'tjy','txt',' ',50
     
106. (e9.3,e11.3)
     
107. *END
     
108. /INPUT,ansuitmp
     
109. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下是一个矩阵形成的例子
     
110. !*dim,aa,,2,3
     
111. !*CREATE,ansuitmp
     
112. !*vread,aa(1,1),'data','dat',,jik,3,2
     
113. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     
114. antype,4
     
115. TRNOPT,FULL
     
116. LUMPM,0
     
117. *do,t,1,50,1    !地震波时间5秒，时间间隔0.1
     
118. time,0.1*t      !每一个加载步的求解时间
     
119. kbc,0           !每个加载步的加载形式为rump
     
120. nsub,10         !每一个加载步分为10个子步
     
121. acel,tjx(2,t),tjy(2,t)   !施加横向、竖向加速度
     
122. nsel,all
     
123. solve            !求解
     
124. *enddo
     

复制代码

李子

先谢谢了，能否请教一下，怎样在瞬态分析中施加正弦荷载？