移动载荷作用下桥梁的动态响应

Chelsea

1. 
   
2. /filname,mobileforce,1
   
3. keyw,pr_struc,1
   
4. /prep7
   
5. /config,nrst,2000  !定义最大的计算结果子步为2000步
   
6. !定义单元类型
   
7. et,1,beam3
   
8. !定义实常数
   
9. r,1,0.1,0.0001/12,0.1
   
10. !定义材料属性
    
11. mp,ex,1,207e9
    
12. mp,prxy,1,0.3
    
13. mp,dens,1,2000
    
14. !定义相关标量参数，路径parameter>scalar parameter
    
15. *SET,deltl,32/100   
    
16. *SET,v,120  
    
17. *SET,deltt,deltl/v*3.6  
    
18. *SET,f,1000
    
19. *SET,w,10
    
20. !建立网格化模型,将整个简支梁划分为100个单元，则车子的前后轮之间将包含2.56/（32/100）=8个单元
    
21. n,1
    
22. n,101,32
    
23. fill,1,101
    
24. *do,i,1,100,1
    
25. e,i,i+1
    
26. *enddo
    
27. !加载与求解
    
28. !匀速常量力 的情况
    
29. !(1)指定分析类型
    
30. /solu
    
31. antype,4
    
32. trnopt,full
    
33. LUMPM,0
    
34. 
    
35. !(2)位移边界条件
    
36. d,1,all
    
37. d,101,all
    
38. !(3)施加常量力并求解
    
39. *do,ii,1,101+8,1
    
40. outres,all,all          !输出所有加载子步的所有结果
    
41. time,ii*deltt           !每个载荷步的最终求解时间
    
42. nsubst,5                !每个载荷步包含五个载荷子步
    
43. *if,ii,lt,8,then
    
44. fdele,all,all           !删除所有集中力载荷
    
45. nsel,s,loc,x,(ii-1)*deltl    !选择集中力作用节点
    
46. f,all,fy,-f                  !施加集中载荷
    
47. allsel
    
48. solve
    
49. *elseif,ii,lt,101
    
50. fdele,all,all
    
51. nsel,s,loc,x,(ii-1)*deltl
    
52. nsel,a,loc,x,(ii-7-1)*deltl
    
53. f,all,fy,-f
    
54. allsel
    
55. solve
    
56. *else
    
57. fdele,all,all
    
58. nsel,s,loc,x,(ii-7-1)*deltl
    
59. f,all,fy,-f
    
60. allsel
    
61. solve
    
62. *endif
    
63. *enddo
    
64. 
    
65. 
    
66. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
67. !匀速简谐力的情况
    
68. ！简谐力的施加与常量力的施加过程差不多，只需将上面的第三步路径命令流改为如下形式
    
69. !注意：以下命令流不是续接第（3）步，而是代替第（3）步的命令流
    
70. *do,ii,1,101+8,1
    
71. outres,all,all          !输出所有加载子步的所有结果
    
72. time,ii*deltt           !每个载荷步的最终求解时间
    
73. nsubst,5                !每个载荷步包含五个载荷子步
    
74. *if,ii,lt,8,then
    
75. fdele,all,all           !删除所有集中力载荷
    
76. nsel,s,loc,x,(ii-1)*deltl    !选择集中力作用节点
    
77. f,all,fy,-f*cos(w*(ii-1)*deltt)                  !施加集中载荷
    
78. allsel
    
79. solve
    
80. *elseif,ii,lt,101
    
81. fdele,all,all
    
82. nsel,s,loc,x,(ii-1)*deltl
    
83. nsel,a,loc,x,(ii-7-1)*deltl
    
84. f,all,fy,-f*cos(w*(ii-1)*deltt)
    
85. allsel
    
86. solve
    
87. *else
    
88. fdele,all,all
    
89. nsel,s,loc,x,(ii-7-1)*deltl
    
90. f,all,fy,-f*cos(w*(ii-1)*deltt)
    
91. allsel
    
92. solve
    
93. *endif
    
94. *enddo
    

复制代码