共享：联系梁的影响线的例子

djaity

和大家分享一下，希望对大家有所帮助

如果用sap2000来画的画，实在是非常简单。但如果在ansys里面画，就实在是两个世界。以下是一个联系梁的影响线的例子（在别人中拿来，忘记是谁了，不好意思）：

1. !!连续梁桥影响线的计算示例
   
2. !!日期：2005-11-12
   
3. !!计算单位：kN,m
   
4. !!!!********************************
   
5. !!!!初始化分析
   
6. FINI
   
7. /CLEAR
   
8. SPAN=3 !定义连续梁跨数
   
9. L=20 !定义单跨跨径
   
10. /FILNAME,MODEL
    
11. /TITLE,The %SPAN%*%L%M Continuous Bridge Analysis
    
12. /NOPR
    
13. /PREP7
    
14. *AFUN,DEG
    
15. /VIEW,1,0,0,1
    
16. !!!! 初始化结束
    
17. !!!!********************************
    
18. !!!!********************************
    
19. !!!!定义相关参数
    
20. SEGMENT=5 !每米单元个数
    
21. NMAX=SEGMENT*L*SPAN+1 !计算最大节点号
    
22. X0=0.0
    
23. Y0=0.0 !定义原点
    
24. SEC_H=L/20
    
25. SEC_W=L/10 !计算截面高、宽
    
26. !!!!定义相关参数定义结束
    
27. !!!!********************************
    
28. !!!!********************************
    
29. !!!!建立有限元模型
    
30. ET,1,BEAM3
    
31. MP,EX,1,3.5E7
    
32. MP,PRXY,1,0.167 !定义单元和材料属性
    
33. R,1,SEC_W*SEC_H,SEC_H**3*SEC_W/12,SEC_H !定义梁截面
    
34. N,1,X0,Y0
    
35. N,NMAX,X0+L*SPAN,Y0
    
36. FILL,1,NMAX !定义节点
    
37. 
    
38. TYPE,1
    
39. MAT,1
    
40. REAL,1
    
41. ESYS,0 !定义单元属性
    
42. *DO,I,1,NMAX-1
    
43. E,I,I+1
    
44. *ENDDO !生成单元
    
45. *DO,J,1,NMAX,SEGMENT*L
    
46. D,J,UY
    
47. *IF,J,EQ,1,THEN
    
48. D,J,UX
    
49. *ENDIF
    
50. *ENDDO !对节点施加约束
    
51. ALLS
    
52. /PBC,ALL,,1
    
53. EPLOT
    
54. SAVE !存储模型
    
55. !!!!模型建立结束
    
56. !!!!********************************
    
57. !!!!********************************
    
58. !!!!模型求解
    
59. /SOLU
    
60. *DO,I,1,NMAX
    
61. F,I,FY,-1
    
62. *IF,I,GT,1,THEN
    
63. FDELE,I-1,FY
    
64. *ENDIF !对节点施加单元荷载
    
65. LSWRITE,I !存储载荷步
    
66. *ENDDO
    
67. ANTYPE,0
    
68. OUTRES,ALL,LAST,
    
69. OUTPR,BASIC,ALL,
    
70. TIME,1
    
71. AUTOTS,-1
    
72. NSUBST,1, , ,1
    
73. KBC,0 !定义有关求解参数
    
74. ALLS,ALL
    
75. LSSOLVE,1,NMAX !求解所有载荷步
    
76. !!!!模型求解结束
    
77. !!!!********************************
    
78. !!!!********************************
    
79. !!!!进入后处理器计算影响线
    
80. /POST1
    
81. *DIM,NODE_X,ARRAY,NMAX !定义节点X坐标存储数组
    
82. 
    
83. *DIM,N_QY,TABLE,NMAX,NMAX !定义QY 影响线存储表
    
84. *DIM,N_MZ,TABLE,NMAX,NMAX !定义MZ 影响线存储表
    
85. *DIM,QY_EXTR,ARRAY,4,NMAX !定义QY 影响线极值存储数组
    
86. *DIM,MZ_EXTR,ARRAY,4,NMAX !定义MZ 影响线极值存储数组
    
87. *DIM,TRANS_V,ARRAY,NMAX-1 !定义中间传递数据向量
    
88. ESEL,ALL
    
89. *DO,I,1,NMAX
    
90. N_QY(0,I)=I
    
91. N_MZ(0,I)=I
    
92. *ENDDO !给N_QY和N_MZ表的0行赋值（节点号）
    
93. *VGET,NODE_X,NODE,ALL,LOC,X !给NODE_X赋值（节点X坐标）
    
94. *ABSET,'Creating Efficacy Data ......',BAR !初始化状态条
    
95. *DO,I,1,NMAX
    
96. *ABCHECK,NINT(100*I/NMAX) !更新状态条
    
97. SET,I
    
98. ETABLE,QYI,SMISC,2
    
99. ETABLE,QYJ,SMISC,8
    
100. ETABLE,MZI,SMISC,6
     
101. ETABLE,MZJ,SMISC,12 !建立Beam3的单元表数据
     
102. N_QY(I,0)=NODE_X(I,1)
     
103. N_MZ(I,0)=NODE_X(I,1) !给N_QY和N_MZ的0列赋值（节点X坐
     
104. 标）
     
105. *VGET,TRANS_V,ELEM,ALL,ETAB,QYI !从读取各单元I端剪力QYI赋予TRANS_V
     
106. *MFUN,N_QY(I,1),TRAN,TRANS_V !将TRANS_V中的值赋予N_QY的第一列
     
107. *GET,N_QY(I,NMAX),ELEM,NMAX-1,SMISC,8 !读取最后一个单元的J端剪力QY
     
108. *VGET,TRANS_V,ELEM,ALL,ETAB,MZI !从读取各单元I端弯矩MZI赋予TRANS_V
     
109. *MFUN,N_MZ(I,1),TRAN,TRANS_V !将TRANS_V中的值赋予N_MZ的第一列
     
110. *GET,N_MZ(I,NMAX),ELEM,NMAX-1,SMISC,12 !读取最后一个单元的J端弯矩MZ
     
111. *ENDDO
     
112. *ABFINISH !终止状态条（*ABSET *ABCHECK *ABFINISH须联用）
     
113. TRANS_V(0)=
     
114. !!!!QY和MZ的影响线建立结束
     
115. !!!!********************************
     
116. !!!!********************************
     
117. !!!!计算影响线的最大、最小值并决定其位置
     
118. *DIM,COLUMN_V,ARRAY,NMAX !定义中间临时数组
     
119. *ABSET,'Geting Extremum ......',BAR !初始化状态条
     
120. *DO,I,1,NMAX
     
121. *ABCHECK,NINT(100*I/NMAX) !更新状态条
     
122. *VCOL,1
     
123. *MFUN,COLUMN_V(1),COPY,N_QY(1,I) !从N_QY中读取第1列数据给COLUMN_V
     
124. *VSCFUN,INDEX,LMAX,COLUMN_V !从COLUMN_V中提取最大值所处位置序号
     
125. QY_EXTR(1,I)=NODE_X(INDEX) !根据序号从NODE_X中提取坐标值
     
126. *VSCFUN,QY_EXTR(2,I),MAX,COLUMN_V !从COLUMN_V表中提取QY最大值
     
127. 
     
128. *VSCFUN,INDEX,LMIN,COLUMN_V !从COLUMN_V中提取最小值所处位置序号
     
129. QY_EXTR(3,I)=NODE_X(INDEX) !根据序号从NODE_X中提取坐标值
     
130. *VSCFUN,QY_EXTR(4,I),MIN,COLUMN_V !从COLUMN_V表中提取QY最小值
     
131. *VCOL,1
     
132. *MFUN,COLUMN_V(1),COPY,N_QY(1,I) !从N_MZ中读取第1列数据给COLUMN_V
     
133. *VSCFUN,INDEX,LMAX,COLUMN_V !从COLUMN_V中提取最大值所处位置序号
     
134. MZ_EXTR(1,I)=NODE_X(INDEX) !根据序号从NODE_X中提取坐标值
     
135. *VSCFUN,MZ_EXTR(2,I),MAX,COLUMN_V !从COLUMN_V表中提取MZ最大值
     
136. *VSCFUN,INDEX,LMIN,COLUMN_V !从COLUMN_V中提取最小值所处位置序号
     
137. MZ_EXTR(3,I)=NODE_X(INDEX) !根据序号从NODE_X中提取坐标值
     
138. *VSCFUN,MZ_EXTR(4,I),MIN,COLUMN_V !从COLUMN_V表中提取MZ最小值
     
139. *ENDDO
     
140. *ABFINISH !终止状态条
     
141. COLUMN_V(0)=
     
142. !!!!计算影响线的最大、最小值结束
     
143. !!!!********************************
     
144. /eof !结束该宏的执行，以下可以选择执行
     
145. /AXLAB,X,Span_X
     
146. /XRANGE,0,SPAN*L
     
147. /GROPT,DIVX,12
     
148. /GROPT,CURL,1
     
149. /AXLAB,Y,VALUE
     
150. *VPLOT,N_QY(1,0),N_QY(1,1), 151
     
151. *VPLOT,N_MZ(1,0),N_MZ(1,1),51,101,151
     
152. !!!!利用图表画出影响线命令示例，
     
153. !!!!********************************

复制代码