Fortran不想做程序员了如下，我想对NG,NE,MC等等赋值该怎么做
点击联系发帖人
时间：2018-05-19 05:11
豆丁微信公众号
君，已阅读到文档的结尾了呢~~
2016年度精品--FORTRAN 90 程序编程规范
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
2016年度精品--FORTRAN 90 程序编程规范
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='http://www.docin.com/DocinViewer--144.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口fortran程序警告 nuused dummy argument NE，见图，怎么解决。_百度知道
fortran程序警告 nuused dummy argument NE，见图，怎么解决。
我有更好的答案
一般来说，Fortran的数组在定义时需要指定其大小。如果你要用real::h(i)那么i必须是确定的值，例如常数，虚参，COMMON中的变量。如果定义时尚不能确定大小，你需要使用可变的数组。这需要书写allocatable属性和allocate语句。
为您推荐：
其他类似问题
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。以下试题来自：
单项选择题阅读下列FORTRAN程序：
INTEGER X(7)
DATA X/4，5，7，8，10,12,11/
DO 10K=1，7
IF(MOD(X(K)，2).NE.0)THEN
END IF 10 CONTINUE
WRITE(*．20)S，T 20 FORMAT(2X，3F6.1)
程序运行后的输出结果是：A． 23.0 7.7B． 47.0 6.7C． 34.0
8.5D． 57.0 11.4
为您推荐的考试题库
你可能感兴趣的试题
1A． K4=28B． K4=27C． K4=30D． K4=292A． SIN(30°)B． SIN(π/6)C． SIN(0.523)D． SIN(π+30°)3A． Internet网B． 广域网、互连网和城域网C． 局域网、互连网和Internet网D． 广域网，局域网和城域网4A． (125)10B． (1101011)2C． (347)8D． (FF)165A． 硬盘中的一块区域B． 光盘中的一块区域C． 内存中的一块区域D． 优盘中的一块区域
热门相关试卷
最新相关试卷当前位置： >>
Fortran 讲解-平面有限元源程序
Fortran 讲解FORTRAN 是科学计算语言。很多年不碰它了（87 年用）。 。。你怎么不去 FORTRAN 论坛？ 试着说说： SUBROUTINE KE(IO,NE,NWE,T,A1,A2,V,EK,BCA) DIMENSION B(7),BCA(7,NE),EK(6,6) DO 10 I=1,7 B(I)=BCA(I,IO) 10 CONTINUE A=A1/B(7)*T DO 20 I=1,3 DO 20 J=I,3 I1=2*I J1=2*J EK(I1-1,J1-1)=A*(B(I)*B(J)+A2*B(I+3)*B(J+3)) *2 维数组通过表达式计算后赋值 EK(I1-1,J1)=A*(V*B(I)*B(J+3)+A2*B(I+3)*B(J)) EK(I1,J1-1)=A*(V*B(I+3)*B(J)+A2*B(I)*B(J+3)) EK(I1,J1)=A*(B(I+3)*B(J+3)+A2*B(I)*B(J)) 20 CONTINUE DO 30 I=3,6 DO 30 J=1,I EK(I,J)=EK(J,I) 30 CONTINUE IF(NWE.EQ.0) GOTO 60 WRITE(6,40) IO 制输出格式 40 FORMAT(/1X,'EK NE='I5) WRITE(6,50) EK 50 60 FORMAT(1X,6E11.4) RETURN END *格式说明，似乎 1X 前多个/ *输出 EK *同上 *返回操作系统 *程序结束 *循环未完成继续 *如果 NWE=0 转标号 60 处 *输出 IO，6=显示器或打印机，40 是表控格式，就是由 40 标号语句控 *循环未完成继续 *有是一个双重循环 *子例行程序 KE(可以 CALL 调用) *定义 1 个一维数组和 2 个两维数组*循环，10 是下面的标号 *给一维数组 B(I)赋值。但 BCA 函数我没见过，外部函数？ *未完成 7 次循环继续 *表达式计算结果赋值给 A *下面是双重循环――外循环 *内循环EK(I1-1,J1-1)=A*(B(I)*B(J)+A2*B(I+3)*B(J+3)) 假如 I=2,J=3 I1=2*I=4,J1=2*J=6，这样上面的语句就相当于： EK(3,5)=A*(B(2)*B(3)+A2*B(5)*B(6)) 呵呵，就这意思。一维、二维数组你总该知道吧？B(2)就是元素在数组中第二个，EK（3，5）第三行第五 列那个元素(行列式)。有限元离不开行列式，甚至需要用到张量（N 维数组） 。 FORTRAN 在计算上得天独厚。 平面弹性力学有限元源程序（FORTRAN） 1. $DEBUG 2. PROGRAM PLANE 3. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 4. ALLOCATABLE::IJK(:,:),XY(:,:),BCA(:,:),SK(:,:),STR( :,:),MB(:,:),ZB(:),B(:) 5. ALLOCATABLE::DELD(:,:,:),TOD(:,:),DELST(:,:,:),TOST (:,:),DELSUP(:,:),TOTSUP(:) 6. DIMENSION EK(6,6) 7. CHARACTER PN*40,FN*12 8. 9. WRITE(*,'(A)') ' 本程序为计算平面问题的有限元程序' 10. WRITE(*,'(A)') ' 特点：(1)采用三结点三角形单元；' 11. WRITE(*,'(A)') ' (2)采用等带宽存贮技术；' 12. WRITE(*,'(A)') ' (3)采用高斯消元法解线性方程组。' 13. WRITE(*,'(/A)') ' 输入计算问题名(PN):' 14. READ(*,'(A)') PN 15. CALL FNAME(PN,'.DAT',FN) 16. WRITE(*,'(2A)') ' 输入数据文件名为：',FN 17. OPEN(5,FILE=FN,STATUS='OLD') 18. CALL FNAME(PN,'.OUT',FN) 19. WRITE(*,'(/2A)') ' 结果输出数据文件名为： ',FN 20. OPEN(6,FILE=FN,STATUS='UNKNOWN') 21. CALL FNAME(PN,'.OU1',FN) 22. WRITE(*,'(/2A)') ' 参数输出数据文件名为： ',FN 23. OPEN(7,FILE=FN,STATUS='UNKNOWN') 24. 25. READ(5,*) NG,NE,MC,NX,NB,EO,VO,DENSITY,T 26. WRITE(6,120) NG,NE,MC,NX,NB 27. WRITE(6,130) EO,VO,DENSITY,T 28. READ(5,*) NWA,NWE,NWK,NWP 29. NT=2*NG 30. ALLOCATE (IJK(3,NE),XY(2,NG),BCA(7,NE),STR(3,NE),MB(2,NB),ZB(NB),B(NT)) 31. ALLOCATE (DELD(2,NG,NX),TOD(2,NG),DELST(3,NE,NX),TOST(3,NE),DELSUP(NB,NX),TO TSUP(NB)) 32. CALL CLEAR(2,NG,TOD) 33. CALL CLEAR(3,NE,TOST) 34. CALL CLEAR1(NB,TOTSUP) 35. IF (NG.EQ.0) THEN 36. STOP 000 37. ENDIF 38. 39. 40. 41. 42. 43. 44. 45. 46. 47. 48. 49. 50. 51. 52. 53. 54. 55. 56. 57. 100 58. 59. 110 60. 61. 62. 63. 64. 65. 66. 67. 68. 69. 70. 71. 72. 73. 74. 75. 76. 77. 78. 79. 80. 81.CALL INPUT(NG,NE,NB,IJK,XY,MB,ZB) CALL CALND(NG,NE,IJK,ND) ALLOCATE (SK(NT,ND)) IF(MC.EQ.0) THEN E=EO V=VO ELSE E=EO/(1-VO*VO) V=VO/(1-VO) ENDIF IP=0 NX1=NX A1=E/(1-V*V)/4.0 A2=0.5*(1-V) CALL ABC(NG,NE,IJK,XY,BCA,NWA) CALL CLEAR(NT,ND,SK) DO 100 K=1,NE CALL KE(K,NE,T,V,A1,A2,IJK,BCA,EK,NWE) CALL SUMK(K,NE,NT,ND,IJK,EK,SK) CONTINUE CALL CHECK(NT,ND,SK) CONTINUE IP=IP+1 WRITE(6,'(///1X,A,I4)') &荷载工况=&,IP READ(5,*) NF,NP,NM DO I=1,NT B(I)=0.0 ENDDO IF(NF.GT.0) THEN CALL PF(NF,NT,B) ENDIF IF(NP.GT.0) THEN CALL PP(NP,NG,NT,T,XY,B) ENDIF IF(NM.GT.0) THEN CALL PG(NM,NE,NT,DENSITY,T,IJK,BCA,B) ENDIF DO I=1,NB I1=2*(MB(1,I)-1)+2-MB(2,I) DELSUP(I,IP)=-B(I1) ENDDO IF(NF.NE.0.OR.NP.NE.0.OR.NM.NE.0) THEN CALL DBC(NB,ND,NT,NX,NX1,MB,ZB,SK,B) CALL GAUSS(NT,ND,NX1,NX,SK,B,NWK,NWP) 82. 83. 84.CALL STRESS(NE,NT,V,A1,A2,IJK,BCA,B,STR) CALL TRES(IP,NG,NE,NX,NT,B,STR,DELD,TOD,DELST,TOST) CALL SUPT(IP,NG,NE,NB,NX,T,V,A1,A2,IJK,MB,BCA,DELD,DELSUP,TOTSUP) 85. CALL OUTPUT(IP,NG,NE,NB,NX,MB,DELD,TOD,DELST,TOST,DELSUP,TOTSUP) 86. ELSE 87. WRITE(*,'(2x,A,I4,A)') &荷载工况=&,IP,&没有荷载!& 88. WRITE(6,'(2x,A,I4,A)') &荷载工况=&,IP,&没有荷载!& 89. ENDIF 90. NX1=NX1-1 91. IF (NX1.GT.0) GOTO 110 92. 120 FORMAT(1X, '结点总数=', I3, 1X, '单元总数=', I3, 1X, '问题类型=', & 93. & I1, 1X, '荷载工况数=', I2, 1X, '支承位移数=', I2) 94. 130 FORMAT(/1X, '弹性模量=', E10.3, 1X, '泊松比=', F5.3, 1X, '密度=', E10.3, & 95. & 1X, '厚度=', F6.3) 96. END 97. 98. 99. SUBROUTINE GAUSS(NT,ND,NX1,NX,SK,B,NWK,NWP) 100. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 101. DIMENSION SK(NT,ND),A(NT,NT),B(NT) 102. CALL CLEAR(NT,NT,A) 103. DO I=1,NT 104. DO J=1,ND 105. IF((I+J-1).LE.NT) THEN 106. A(I,I+J-1)=SK(I,J) 107. ENDIF 108. ENDDO 109. ENDDO 110. IF(NWK.EQ.1.AND.NX1.EQ.NX) THEN 111. WRITE(7,'(A)') &结构总刚(列出右上三角元素)& 112. DO I=1,NT 113. WRITE(7,100) (A(I,J),J=1,NT) 114. ENDDO 115. ENDIF 116. IF (NWP.EQ.1) THEN 117. WRITE(7,'(1X,A,I3,A)') &第&,NX-NX1+1,&荷载工况 的荷载列阵& 118. DO I=1,NT 119. WRITE(7,'(1x,E11.4)') B(I) 120. ENDDO 121. ENDIF 122. DO 50 K=1,NT-1 123. DO 60 I=K+1,NT 124. C1=A(K,I)/A(K,K) 125. DO 70 J=I,NT 126. A(I,J)=A(I,J)-C1*A(K,J) 127. 70 CONTINUE 128. B(I)=B(I)-C1*B(K) 129. 60 CONTINUE 130. 50 CONTINUE 131. B(NT)=B(NT)/A(NT,NT) 132. DO 80 I=NT-1,1,-1 133. DO 90 J=I+1,NT 134. B(I)=B(I)-A(I,J)*B(J) 135. 90 CONTINUE 136. B(I)=B(I)/A(I,I) 137. 80 CONTINUE 138. 100 FORMAT(1x, ) 139. END 140. 141. 142. SUBROUTINE CHECK(NT,ND,SK) 143. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 144. DIMENSION SK(NT,ND) 145. M=0 146. DO I=1,NT 147. IF(SK(I,1).LE.0.0) THEN 148. M=M+1 149. ENDIF 150. ENDDO 151. IF (M.GT.0) THEN 152. WRITE(*,*) &总刚矩阵的对角元素为负值!!& 153. STOP . ENDIF 155. END 156. 157. SUBROUTINE STRESS(NE,NT,V,A1,A2,IJK,BCA,B,STR) 158. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 159. DIMENSION D1(6),D2(3),B(NT),S(3,6),IJK(3,NE),BCA(7,NE),STR(3,NE) 160. CALL CLEAR(3,NE,STR) 161. DO 10 K=1,NE 162. DO I=1,3 163. II=IJK(I,K) 164. D1(2*(I-1)+1)=B(2*(II-1)+1) 165. D1(2*(I-1)+2)=B(2*(II-1)+2) 166. ENDDO 167. CALL CLEAR(3,6,S) 168. C1=2*A1/BCA(7,K) 169. DO 20 I=1,3 170. S(1,2*(I-1)+1)=C1*BCA(I,K) 171. S(1,2*(I-1)+2)=C1*V*BCA(I+3,K) 172. S(2,2*(I-1)+1)=C1*V*BCA(I,K) 173. S(2,2*(I-1)+2)=C1*BCA(I+3,K) 174. S(3,2*(I-1)+1)=C1*A2*BCA(I+3,K) 175. S(3,2*(I-1)+2)=C1*A2*BCA(I,K) 176. 20 CONTINUE 177. CALL MATMUL(3,6,1,S,D1,D2) 178. DO I=1,3 179. STR(I,K)=D2(I) 180. ENDDO 181. 10 CONTINUE 182. END 183. 184. SUBROUTINE OUTPUT(IP,NG,NE,NB,NX,MB,DELD,TOD,DELST,TOST,DELSUP,TOTSUP) 185. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 186. DIMENSION MB(2,NB),DELD(2,NG,NX),TOD(2,NG) 187. DIMENSION DELST(3,NE,NX),TOST(3,NE),DELSUP(NB,NX),TOTSUP(NB) 188. CHARACTER VE*6 189. WRITE(6,20) IP 190. WRITE(6,30) 191. DO I=1,NG 192. WRITE(6,40) I,DELD(1,I,IP),TOD(1,I),DELD(2,I,IP),TOD(2,I) 193. ENDDO 194. WRITE(6,50) 195. DO J=1,NE 196. WRITE(6,60) J,(DELST(L,J,IP),TOST(L,J),L=1,3) 197. ENDDO 198. WRITE(6,70) 199. DO J=1,NB 200. IF (MB(2,J).EQ.1) THEN 201. VE='x 方向' 202. ELSE 203. VE='Y 方向' 204. ENDIF 205. WRITE(6,80) MB(1,J),VE,DELSUP(J,IP),TOTSUP(J) 206. ENDDO 207. 20 FORMAT(/1X, '荷载工况=', I3, '的计算结果') 208. 30 FORMAT(/, 1X, '结点号', 5X, 'X 位移增量', 5X, 'X 累计位移', 5x, & 209. & 'Y 位移增量', 5X, 'Y 累计位移') 210. 40 FORMAT(1X, I3, 6X, F10.4, 4x, F10.4, 4X, F10.4, 4x, F10.4) 211. 50 FORMAT(/, 1X, '单元号', 6x, 'X 应力增量', 6x, 'X 累计应力', 6x, & 212. & 'Y 应力增量', 6x, 'Y 累计应力', 6x, '剪应力增量', 6x, & 213. & '累计剪应力') 214. 60 FORMAT(2X, I3, 6X, E10.3, 5X, E10.3, 5X, E10.3, 5X, E10.3, 6X, E10.3, & 215. & 6X, E10.3) 216. 70 FORMAT(/, 1X, '支座结点号', 6x, '支反力方向', 6x, '支反力增量', & 217. & 6x, '支反力累计量') 218. 80 FORMAT(5X, I3, 12X, A, 6x, E10.3, 8X, E10.3) 219. END 220. 221. 222. 223. 224. 225. 226. 227. 228. 229. 230. 231. 232. 233. 234. 235. 236. 237. 238. 239. 240. 241. 242.100SUBROUTINE MATMUL(M,N,L,A,B,C) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION A(M,N),B(N,L),C(M,L) DO 100 I=1,M DO 100 J=1,L C(I,J)=0.0 DO 100 K=1,N C(I,J)=C(I,J)+A(I,K)*B(K,J) END SUBROUTINE CLEAR(M,N,A) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION A(M,N) DO I=1,M DO J=1,N A(I,J)=0.0 ENDDO ENDDO END SUBROUTINE ABC(NG,NE,IJK,XY,BCA,NWA) 243. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 244. DIMENSION X(2,5),XY(2,NG),IJK(3,NE),B(7),BCA(7,NE) 245. DO 10 I=1,NE 246. DO 20 K=1,3 247. DO 20 J=1,2 248. X(J,K)=XY(J,IJK(K,I)) 249. 20 CONTINUE 250. DO 30 J=1,2 251. X(J,4)=X(J,1) 252. X(J,5)=X(J,2) 253. 30 CONTINUE 254. DO 40 K=1,3 255. B(K)=X(2,K+1)-X(2,K+2) 256. B(K+3)=X(1,K+2)-X(1,K+1) 257. 40 CONTINUE 258. B(7)=(B(1)*B(5)-B(4)*B(2))*0.5 259. IF (B(7).LE.0.0) THEN 260. WRITE(*,*) &单元号=&,I 261. WRITE(*,*) &单元的结点 I,J AND K 编号： &,(IJK(K,I),K=1,3) 262. STOP . ELSE 264. DO J=1,7 265. BCA(J,I)=B(J) 266. ENDDO 267. ENDIF 268. 10 CONTINUE 269. IF(NWA.EQ.1) THEN 270. WRITE(7,100) 271. DO I=1,NE 272. WRITE(7,110) I,(BCA(J,I),J=1,7) 273. ENDDO 274. ENDIF 275. 100 FORMAT(1X, '单元号', 3x, 'bi', 7x, 'bj', 7x, 'bk', 7x, 'ci', 7x, & 276. & 'cj', 7x, 'ck', 7x, '面积 A') 277. 110 FORMAT(2X, I3, 3X, F6.3, 3X, F6.3, 3X, F6.3, 3X, F6.3, 3X, F6.3, & 278. & 3X, F6.3, 3X, F6.3) 279. END 280. 281. SUBROUTINE SUMK(K,NE,NT,ND,IJK,EK,SK) 282. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 283. ,ND) 284. 285. 286. 287. 288. 289. 290. 291. 292. 293. 294. 295. 296. 297. 298. 299. 20 300. 10 301. 302. 303. 304. 305. 306. 307. 308. 309. 310. 311. 312. 313. 314. BCA(J+3,K) 315. ,K)*BCA(J,K) 316. )*BCA(J,K) 317. 318. )*BCA(J+3,K) 319. 20 320. 10 321.DIMENSIONIJ(3),EK(6,6),IJK(3,NE),SK(NTDO I=1,3 IJ(I)=IJK(I,K) ENDDO DO 10 I=1,3 DO 20 J=1,3 IF (IJ(I).GT.IJ(J)) GOTO 20 M=2*(IJ(I)-1)+1 N=2*(IJ(J)-1)+1-(2*(IJ(I)-1)+1)+1 MO=2*(I-1)+1 NO=2*(J-1)+1 SK(M,N)=SK(M,N)+EK(MO,NO) SK(M,N+1)=SK(M,N+1)+EK(MO,NO+1) SK(M+1,N)=SK(M+1,N)+EK(MO+1,NO+1) IF (IJ(I).EQ.IJ(J)) GOTO 20 SK(M+1,N-1)=SK(M+1,N-1)+EK(MO+1,NO) CONTINUE CONTINUE END SUBROUTINE KE(K,NE,T,V,A1,A2,IJK,BCA,EK,NWE) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION IJ(3),IJK(3,NE),BCA(7,NE),EK(6,6) CALL CLEAR(6,6,EK) DO I=1,3 IJ(I)=IJK(I,K) ENDDO DO 10 I=1,3 DO 20 J=I,3 IO=2*(I-1) JO=2*(J-1) EK(IO+1,JO+1)=BCA(I,K)*BCA(J,K)+A2*BCA(I+3,K)* EK(IO+1,JO+2)=V*BCA(I,K)*BCA(J+3,K)+A2*BCA(I+3 EK(IO+2,JO+2)=BCA(I+3,K)*BCA(J+3,K)+A2*BCA(I,K IF (I.EQ.J) GOTO 20 EK(IO+2,JO+1)=V*BCA(I+3,K)*BCA(J,K)+A2*BCA(I,K CONTINUE CONTINUE DO I=2,6 322. 323. 324. 325. 326. 327. 328. 329. 330. 331. 332. 333. 334. 335. 336. 337. 338. 339. 340. 341. 342. 343. 344. 345. 346. 347. 348. 349. 350. 351. 352. 353. 354. 355. 356. 357. 358. 359. 360. 361. 362. 363. 364. 365.40 50DO J=1,I-1 EK(I,J)=EK(J,I) ENDDO ENDDO DO I=1,6 DO J=1,6 EK(I,J)=EK(I,J)*A1*T/BCA(7,K) ENDDO ENDDO IF(NWE.EQ.1) THEN WRITE(7,40) K, &号单元的单刚& DO I=1,6 WRITE(7,50) (EK(I,J),J=1,6) ENDDO ENDIF FORMAT(1X, I3, A) FORMAT(1X, 6(E10.4, 4X)) END SUBROUTINE PF(NF,NT,B) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION MF(2,NF),ZF(NF),B(NT) DO I=1,NF READ(5,*) K,MF(1,I),MF(2,I),ZF(I) ENDDO WRITE(6,'(/A,I3)') &集中荷载个数 NF=&,NF WRITE(6,20) DO I=1,NF WRITE(6,30) MF(1,I),MF(2,I),ZF(I) ENDDO DO 10 I=1,NF II=2*(MF(1,I)-1)-MF(2,I)+2 B(II)=B(II)+ZF(I) CONTINUE FORMAT(1X, '结点号', 4x, '方向', 4x, '荷载值') FORMAT(2X, I3, 5X, I3, 4X, E10.4) END SUBROUTINE FNAME(PN,FN2,FN) CHARACTER PN*40,FN2*4,FN*12 去掉 PN 中前面的空格 DO 10 I=1,40 IF(PN(I:I).EQ.' ') GOTO 10 IP=I10 20 30! 366. GOTO 20 367. 10 CONTINUE 368. 20 CONTINUE 369. FN(1:8)=PN(IP:IP+7) 370. ! 去掉 FN 中后面的空格 371. DO 30 I=8,1,-1 372. IF(FN(I:I).EQ.' ') GOTO 30 373. IL=I 374. GOTO 40 375. 30 CONTINUE 376. 40 CONTINUE 377. ! 生成文件名 FN=PN+FN2 378. FN(IL+1:IL+4)=FN2(1:4) 379. END 380. 381. SUBROUTINE PP(NP,NG,NT,T,XY,B) 382. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 383. DIMENSION MP(2,NP),ZP(NP),XY(2,NG),B(NT) 384. DO I=1,NP 385. READ(5,*) K,MP(1,I),MP(2,I),ZP(I) 386. ENDDO 387. WRITE(6,'(/A,I3)') &分布荷载数 NP=&,NP 388. WRITE(6,20) 389. DO I=1,NP 390. WRITE(6,30) I,MP(1,I),MP(2,I),ZP(I) 391. ENDDO 392. DO 10 I=1,NP 393. II=2*(MP(1,I)-1) 394. JJ=2*(MP(2,I)-1) 395. PX=0.5*ZP(I)*T*(XY(2,MP(1,I))-XY(2,MP(2,I))) 396. PY=0.5*ZP(I)*T*(XY(1,MP(2,I))-XY(1,MP(1,I))) 397. B(II+1)=B(II+1)+PX 398. B(II+2)=B(II+2)+PY 399. B(JJ+1)=B(JJ+1)+PX 400. B(JJ+2)=B(JJ+2)+PY 401. 10 CONTINUE 402. 20 FORMAT(1X, '序号', 4x, '首结点号', 4x, '末结点 号', 7x, '荷载值') 403. 30 FORMAT(1X, I3, 8x, I3, 9X, I3, 7X, F10.3) 404. END 405. 406. SUBROUTINE PG(NM,NE,NT,DENSITY,T,IJK,BCA,B) 407. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 408. DIMENSION IJK(3,NE),BCA(7,NE),MG(NM),B(NT) 409. 410. 411. 412. 413. 414. 415. 416. 417. 418. 419. 420. 421. 422. 423. 424. 425. 426. 427. 428. 429. 430. 431. 432. 433. 434. 435. 436. 437. 438. 439. 440. 441. 442. 443. 444. 445. 446. 447. 448. 449. 450. 451. 452.10 100IF (NM.GE.NE) THEN NM=NE DO I=1,NM MG(I)=I ENDDO ELSE READ(5,*) (MG(I),I=1,NM) ENDIF WRITE(6,'(/A,I3)') &计自重单元数 NM=&,NM WRITE(6,'(A)') &计自重的单元号为：& K1=0 K2=0 DO WHILE((NM-K2).GE.10) WRITE(6,100) (MG(K2+J),J=1,10) K1=K1+1 K2=10*K1 ENDDO WRITE(6,100) (MG(K2+J),J=1,NM-K2) DO 10 I=1,NM II=2*(IJK(1,MG(I))-1) JJ=2*(IJK(2,MG(I))-1) KK=2*(IJK(3,MG(I))-1) PE1=-DENSITY*T*BCA(7,MG(I))/3.0 B(II+2)=B(II+2)+PE1 B(JJ+2)=B(JJ+2)+PE1 B(KK+2)=B(KK+2)+PE1 CONTINUE FORMAT(10(1X,I3)) END SUBROUTINE INPUT(NG,NE,NB,IJK,XY,MB,ZB) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION IJK(3,NE),XY(2,NG),MB(2,NB),ZB(NB) DO I=1,NE READ(5,*) K,(IJK(J,I),J=1,3) ENDDO DO I=1,NG READ(5,*) K,XY(1,I),XY(2,I) ENDDO DO I=1,NB READ(5,*) K,MB(1,I),MB(2,I),ZB(I) ENDDO WRITE(6,15) DO K=1,NE 453. WRITE(6,20) K,(IJK(J,K),J=1,3) 454. ENDDO 455. WRITE(6,30) 456. DO K=1,NG 457. WRITE(6,35) K,XY(1,K),XY(2,K) 458. ENDDO 459. WRITE(6,40) 460. DO K=1,NB 461. WRITE(6,45) K,MB(1,K),MB(2,K),ZB(K) 462. ENDDO 463. 15 FORMAT(/1X, '单元号', 4X, 'I 结点号', 4x, 'J 结点号', 4x, 'K 结点号') 464. 20 FORMAT(2X, I3, 7X, I3, 8X, I3, 8X, I3) 465. 30 FORMAT(/1X, '结点号', 8X, 'X 坐标', 8X, 'Y 坐标 ') 466. 35 FORMAT(2X, I3, 6X, F8.3, 5X, F8.3) 467. 40 FORMAT(/1X, '序号', 4x, '支承结点号', 6x, '支 承方向', 6x, & 468. & '支承位移') 469. 45 FORMAT(2X, I2, 8X, I4, 11X, I2, 7X, F10.3) 470. END 471. 472. SUBROUTINE DBC(NB,ND,NT,NX,NX1,MB,ZB,SK,B) 473. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 474. DIMENSION MB(2,NB),ZB(NB),SK(NT,ND),B(NT) 475. DO 5 I=1,NB 476. N=2*MB(1,I)-MB(2,I) 477. Z=ZB(I) 478. IF (ABS(Z).LT.1.0E-10) THEN 479. IF (NX.NE.NX1) GOTO 30 480. SK(N,1)=1.0 481. DO J=2,ND 482. SK(N,J)=0.0 483. ENDDO 484. DO 10 K=2,ND 485. IF (N.LT.K) GOTO 20 486. M=N-K+1 487. SK(M,K)=0.0 488. 10 CONTINUE 489. 20 CONTINUE 490. 30 B(N)=0.0 491. ELSE 492. IF (NX.NE.NX1) GOTO 40 493. SK(N,1)=SK(N,1)*1.0E+15 494. 40 B(N)=SK(N,1)*Z 495. ENDIF 496. 5 CONTINUE 497. END 498. 499. SUBROUTINE CALND(NG,NE,IJK,ND) 500. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 501. DIMENSION IJK(3,NE) 502. ND=0 503. DO 10 I=1,NG 504. M1=0 505. DO 20 J=1,NE 506. DO 30 K=1,3 507. IF(IJK(K,J).EQ.I) THEN 508. M2=0 509. DO 40 L=1,3 510. M3=IJK(L,J)-I 511. IF(M3.GT.M2) THEN 512. M2=M3 513. ENDIF 514. 40 CONTINUE 515. GOTO 25 516. ENDIF 517. 30 CONTINUE 518. 25 IF(M2.GT.M1) M1=M2 519. 20 CONTINUE 520. IF(M1.GT.ND)ND=M1 521. 10 CONTINUE 522. ND=2*(ND+1) 523. WRITE(7,100) ND 524. 100 FORMAT(1X,'总刚矩阵的半带宽 ND 为：',I4) 525. END 526. 527. SUBROUTINE TRES(IP,NG,NE,NX,NT,B,STR,DELD,TOD,DELST,TOST) 528. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 529. DIMENSION B(NT),STR(3,NE),DELD(2,NG,NX),TOD(2,NG),DELST(3,NE,NX),TOST(3,NE) 530. DO I=1,NG 531. IO=2*(I-1) 532. DELD(1,I,IP)=B(IO+1) 533. DELD(2,I,IP)=B(IO+2) 534. TOD(1,I)=TOD(1,I)+DELD(1,I,IP) 535. TOD(2,I)=TOD(2,I)+DELD(2,I,IP) 536. ENDDO 537. DO I=1,NE 538. DO J=1,3 539. DELST(J,I,IP)=STR(J,I) 540. TOST(J,I)=TOST(J,I)+DELST(J,I,IP) 541. ENDDO 542. ENDDO 543. END 544. 545. 546. SUBROUTINE SUPT(IP,NG,NE,NB,NX,T,V,A1,A2,IJK,MB,BCA,DELD,DELSUP,TOTSUP) 547. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 548. DIMENSION IJK(3,NE),MB(2,NB),BCA(7,NE),DELD(2,NG,NX),EK(6,6) 549. DIMENSION D1(6),DELSUP(NB,NX),TOTSUP(NB) 550. DO 10 K=1,NB 551. I1=MB(1,K) 552. I2=MB(2,K) 553. RF=0.0 554. DO 20 I=1,NE 555. DO 30 J=1,3 556. I3=IJK(J,I) 557. IF (I1.NE.I3) THEN 558. GOTO 30 559. ELSE 560. I4=2*(J-1)+2-I2 561. CALL KE(I,NE,T,V,A1,A2,IJK,BCA,EK,0) 562. DO M=1,3 563. IO=IJK(M,I) 564. D1(2*(M-1)+1)=DELD(1,IO,IP) 565. D1(2*(M-1)+2)=DELD(2,IO,IP) 566. ENDDO 567. DO L=1,6 568. RF=RF+EK(I4,L)*D1(L) 569. ENDDO 570. GOTO 20 571. ENDIF 572. 30 CONTINUE 573. 20 CONTINUE 574. DELSUP(K,IP)=DELSUP(K,IP)+RF 575. TOTSUP(K)= TOTSUP(K)+DELSUP(K,IP) 576. 10 CONTINUE 577. END 578. 579. 580. 581. 582. 583. 584. 585. 复制代码SUBROUTINE CLEAR1(N,B) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION B(N) DO I=1,N B(I)=0.0 ENDDO ENDPFEM,EXE 平面问题有限元程序的数据文件输入格式： 基本输入数据文件名：*.DAT 计算结果输出数据文件名：*.OUT 中间参数输出数据文件名：*.OU1 下面按在文件中的前后顺序进行说明： 1 基本控制参数信息：NG,NE,MC,NX,NB,EO,VO,DENSITY ,T(共计 5 个整形数，4 个实型数) NG：结构的结点总数； NE：结构的单元总数； MC：平面问题的类型，MC=0，为平面应力，MC=1，为平面应变； NX：荷载工况数； NB：支承位移数； EO：材料弹性模量(Pa)； VO：材料泊松比； DENSITY ：容重(N/m3) T ：材料厚度(m)； 2 打印输出控制参数：NWA,NEW,NWK,NWP(4 个整形数) 等于 1 时，输出，否则不输出。 3 单元结点信息：(K,(IJK(I,K),I=1,3),K=1,NE) (每行 4 个整形数，共计 NE 行) K：单元号； IJK(1,K)：K 单元 I 结点编号； IJK(2,K)：K 单元 J 结点编号； IJK(3,K)：K 单元 K 结点编号； 4 结点坐标信息：((K,XY(1,K),XY(2,K)),K=1,NG)(每行 3 个整形数，共计 NG 行) K：结点号 XY(1,K)：K 结点 X 坐标； XY(2,K)：K 结点 Y 坐标； 5 支承信息：((K,MB(1,K),MB(2,K),ZB(K)),K=1,NB)(每行 3 个整形数，1 个实型数，共计 NB 行) K：支承位移序号； MB(1,K)：第 K 个支承位移所在的结点号； MB(2,K)：第 K 个支承位移的坐标方向； ZB(K)： 第 K 个支承位移的数值； 6 按 NX 荷载工况数输入荷载信息：每一荷载工况如下 (1) NF,NP,NM(3 个整型数) NF：集中荷载个数； NP：分布荷载个数； NM：计自重单元数； (2) 若 NF≠0，则输入下面数据 K,MF(1,K),MF(2,K),ZF(K)(每行 3 个整形数，1 个实型数，共计 NF 行) K：集中荷载序号； MF(1,K)：第 K 个集中荷载作用的结点号； MF(2,K)：第 K 个集中荷载的坐标方向； ZF(K)： 第 K 个集中荷载的数值； (3) 若 NP≠0，则输入下面数据 K,MP(1,K),MP(2,K),ZP(K)(每行 3 个整形数，1 个实型数，共计 NP 行) K：分布荷载序号； MP(1,K)：第 K 个分布荷载作用的结点号； MP(2,K)：第 K 个分布荷载的坐标方向； ZP(K)： 第 K 个分布荷载的数值； (4) 若 NM≠0，则输入下面数据 若 NM≥NE，则表示计所有单元的自重，不需输入计自重的单元号； 若 NM&NE，则需要输入计自重的单元号；
里面包含结构的计算简图以及源程序算例一计算简图及结果输出用平面刚架静力计算...Fortran语言 有限元程序... 18页 1下载券 Fortran 讲解-平面有限元... 17页...Fortran 讲解-平面有限元源... 17页 1财富值 Fortran语言 有限元程序分... ...有限元程序设计005 15页 免费如要投诉违规内容,请到百度文库投诉中心;如要提出功...Fortran语言编写的有限元结构程序_能源/化工_工程科技_专业资料。算例一计算简图...Fortran 讲解-平面有限元... 17页 1下载券 王勖成《有限单元法》源... 16页...结构力学刚架问题的有限元分析 1 有限元分析软件的开发 1.1 程序功能该程序为...Fortran_讲解-平面有限元... 17页 免费
平面弹性力学有限元源程... 暂无评价...读有限元Fortran程序笔记_计算机软件及应用_IT/计算机_专业资料。读有限元程序笔记 1. ALLOCATABLE::COORD(:,:),PROPS(:,:,:) PROPS(:,:,:)是三维数组。 2...有限元编程算例(fortran) - 本程序通过Fortran语言编写,程序在Intel Parallel Studio XE 2013 with VS2013中成功运行,程序为《计算力学》(龙...
All rights reserved Powered by
www.tceic.com
copyright ©right 。文档资料库内容来自网络，如有侵犯请联系客服。}
君，已阅读到文档的结尾了呢~~
2016年度精品--FORTRAN 90 程序编程规范
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
2016年度精品--FORTRAN 90 程序编程规范
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='http://www.docin.com/DocinViewer--144.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口fortran程序警告 nuused dummy argument NE，见图，怎么解决。_百度知道
fortran程序警告 nuused dummy argument NE，见图，怎么解决。
我有更好的答案
一般来说，Fortran的数组在定义时需要指定其大小。如果你要用real::h(i)那么i必须是确定的值，例如常数，虚参，COMMON中的变量。如果定义时尚不能确定大小，你需要使用可变的数组。这需要书写allocatable属性和allocate语句。
为您推荐：
其他类似问题
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。以下试题来自：
单项选择题阅读下列FORTRAN程序：
INTEGER X(7)
DATA X/4，5，7，8，10,12,11/
DO 10K=1，7
IF(MOD(X(K)，2).NE.0)THEN
END IF 10 CONTINUE
WRITE(*．20)S，T 20 FORMAT(2X，3F6.1)
程序运行后的输出结果是：A． 23.0 7.7B． 47.0 6.7C． 34.0
8.5D． 57.0 11.4
为您推荐的考试题库
你可能感兴趣的试题
1A． K4=28B． K4=27C． K4=30D． K4=292A． SIN(30°)B． SIN(π/6)C． SIN(0.523)D． SIN(π+30°)3A． Internet网B． 广域网、互连网和城域网C． 局域网、互连网和Internet网D． 广域网，局域网和城域网4A． (125)10B． (1101011)2C． (347)8D． (FF)165A． 硬盘中的一块区域B． 光盘中的一块区域C． 内存中的一块区域D． 优盘中的一块区域
热门相关试卷
最新相关试卷当前位置： >>
Fortran 讲解-平面有限元源程序
Fortran 讲解FORTRAN 是科学计算语言。很多年不碰它了（87 年用）。 。。你怎么不去 FORTRAN 论坛？ 试着说说： SUBROUTINE KE(IO,NE,NWE,T,A1,A2,V,EK,BCA) DIMENSION B(7),BCA(7,NE),EK(6,6) DO 10 I=1,7 B(I)=BCA(I,IO) 10 CONTINUE A=A1/B(7)*T DO 20 I=1,3 DO 20 J=I,3 I1=2*I J1=2*J EK(I1-1,J1-1)=A*(B(I)*B(J)+A2*B(I+3)*B(J+3)) *2 维数组通过表达式计算后赋值 EK(I1-1,J1)=A*(V*B(I)*B(J+3)+A2*B(I+3)*B(J)) EK(I1,J1-1)=A*(V*B(I+3)*B(J)+A2*B(I)*B(J+3)) EK(I1,J1)=A*(B(I+3)*B(J+3)+A2*B(I)*B(J)) 20 CONTINUE DO 30 I=3,6 DO 30 J=1,I EK(I,J)=EK(J,I) 30 CONTINUE IF(NWE.EQ.0) GOTO 60 WRITE(6,40) IO 制输出格式 40 FORMAT(/1X,'EK NE='I5) WRITE(6,50) EK 50 60 FORMAT(1X,6E11.4) RETURN END *格式说明，似乎 1X 前多个/ *输出 EK *同上 *返回操作系统 *程序结束 *循环未完成继续 *如果 NWE=0 转标号 60 处 *输出 IO，6=显示器或打印机，40 是表控格式，就是由 40 标号语句控 *循环未完成继续 *有是一个双重循环 *子例行程序 KE(可以 CALL 调用) *定义 1 个一维数组和 2 个两维数组*循环，10 是下面的标号 *给一维数组 B(I)赋值。但 BCA 函数我没见过，外部函数？ *未完成 7 次循环继续 *表达式计算结果赋值给 A *下面是双重循环――外循环 *内循环EK(I1-1,J1-1)=A*(B(I)*B(J)+A2*B(I+3)*B(J+3)) 假如 I=2,J=3 I1=2*I=4,J1=2*J=6，这样上面的语句就相当于： EK(3,5)=A*(B(2)*B(3)+A2*B(5)*B(6)) 呵呵，就这意思。一维、二维数组你总该知道吧？B(2)就是元素在数组中第二个，EK（3，5）第三行第五 列那个元素(行列式)。有限元离不开行列式，甚至需要用到张量（N 维数组） 。 FORTRAN 在计算上得天独厚。 平面弹性力学有限元源程序（FORTRAN） 1. $DEBUG 2. PROGRAM PLANE 3. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 4. ALLOCATABLE::IJK(:,:),XY(:,:),BCA(:,:),SK(:,:),STR( :,:),MB(:,:),ZB(:),B(:) 5. ALLOCATABLE::DELD(:,:,:),TOD(:,:),DELST(:,:,:),TOST (:,:),DELSUP(:,:),TOTSUP(:) 6. DIMENSION EK(6,6) 7. CHARACTER PN*40,FN*12 8. 9. WRITE(*,'(A)') ' 本程序为计算平面问题的有限元程序' 10. WRITE(*,'(A)') ' 特点：(1)采用三结点三角形单元；' 11. WRITE(*,'(A)') ' (2)采用等带宽存贮技术；' 12. WRITE(*,'(A)') ' (3)采用高斯消元法解线性方程组。' 13. WRITE(*,'(/A)') ' 输入计算问题名(PN):' 14. READ(*,'(A)') PN 15. CALL FNAME(PN,'.DAT',FN) 16. WRITE(*,'(2A)') ' 输入数据文件名为：',FN 17. OPEN(5,FILE=FN,STATUS='OLD') 18. CALL FNAME(PN,'.OUT',FN) 19. WRITE(*,'(/2A)') ' 结果输出数据文件名为： ',FN 20. OPEN(6,FILE=FN,STATUS='UNKNOWN') 21. CALL FNAME(PN,'.OU1',FN) 22. WRITE(*,'(/2A)') ' 参数输出数据文件名为： ',FN 23. OPEN(7,FILE=FN,STATUS='UNKNOWN') 24. 25. READ(5,*) NG,NE,MC,NX,NB,EO,VO,DENSITY,T 26. WRITE(6,120) NG,NE,MC,NX,NB 27. WRITE(6,130) EO,VO,DENSITY,T 28. READ(5,*) NWA,NWE,NWK,NWP 29. NT=2*NG 30. ALLOCATE (IJK(3,NE),XY(2,NG),BCA(7,NE),STR(3,NE),MB(2,NB),ZB(NB),B(NT)) 31. ALLOCATE (DELD(2,NG,NX),TOD(2,NG),DELST(3,NE,NX),TOST(3,NE),DELSUP(NB,NX),TO TSUP(NB)) 32. CALL CLEAR(2,NG,TOD) 33. CALL CLEAR(3,NE,TOST) 34. CALL CLEAR1(NB,TOTSUP) 35. IF (NG.EQ.0) THEN 36. STOP 000 37. ENDIF 38. 39. 40. 41. 42. 43. 44. 45. 46. 47. 48. 49. 50. 51. 52. 53. 54. 55. 56. 57. 100 58. 59. 110 60. 61. 62. 63. 64. 65. 66. 67. 68. 69. 70. 71. 72. 73. 74. 75. 76. 77. 78. 79. 80. 81.CALL INPUT(NG,NE,NB,IJK,XY,MB,ZB) CALL CALND(NG,NE,IJK,ND) ALLOCATE (SK(NT,ND)) IF(MC.EQ.0) THEN E=EO V=VO ELSE E=EO/(1-VO*VO) V=VO/(1-VO) ENDIF IP=0 NX1=NX A1=E/(1-V*V)/4.0 A2=0.5*(1-V) CALL ABC(NG,NE,IJK,XY,BCA,NWA) CALL CLEAR(NT,ND,SK) DO 100 K=1,NE CALL KE(K,NE,T,V,A1,A2,IJK,BCA,EK,NWE) CALL SUMK(K,NE,NT,ND,IJK,EK,SK) CONTINUE CALL CHECK(NT,ND,SK) CONTINUE IP=IP+1 WRITE(6,'(///1X,A,I4)') &荷载工况=&,IP READ(5,*) NF,NP,NM DO I=1,NT B(I)=0.0 ENDDO IF(NF.GT.0) THEN CALL PF(NF,NT,B) ENDIF IF(NP.GT.0) THEN CALL PP(NP,NG,NT,T,XY,B) ENDIF IF(NM.GT.0) THEN CALL PG(NM,NE,NT,DENSITY,T,IJK,BCA,B) ENDIF DO I=1,NB I1=2*(MB(1,I)-1)+2-MB(2,I) DELSUP(I,IP)=-B(I1) ENDDO IF(NF.NE.0.OR.NP.NE.0.OR.NM.NE.0) THEN CALL DBC(NB,ND,NT,NX,NX1,MB,ZB,SK,B) CALL GAUSS(NT,ND,NX1,NX,SK,B,NWK,NWP) 82. 83. 84.CALL STRESS(NE,NT,V,A1,A2,IJK,BCA,B,STR) CALL TRES(IP,NG,NE,NX,NT,B,STR,DELD,TOD,DELST,TOST) CALL SUPT(IP,NG,NE,NB,NX,T,V,A1,A2,IJK,MB,BCA,DELD,DELSUP,TOTSUP) 85. CALL OUTPUT(IP,NG,NE,NB,NX,MB,DELD,TOD,DELST,TOST,DELSUP,TOTSUP) 86. ELSE 87. WRITE(*,'(2x,A,I4,A)') &荷载工况=&,IP,&没有荷载!& 88. WRITE(6,'(2x,A,I4,A)') &荷载工况=&,IP,&没有荷载!& 89. ENDIF 90. NX1=NX1-1 91. IF (NX1.GT.0) GOTO 110 92. 120 FORMAT(1X, '结点总数=', I3, 1X, '单元总数=', I3, 1X, '问题类型=', & 93. & I1, 1X, '荷载工况数=', I2, 1X, '支承位移数=', I2) 94. 130 FORMAT(/1X, '弹性模量=', E10.3, 1X, '泊松比=', F5.3, 1X, '密度=', E10.3, & 95. & 1X, '厚度=', F6.3) 96. END 97. 98. 99. SUBROUTINE GAUSS(NT,ND,NX1,NX,SK,B,NWK,NWP) 100. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 101. DIMENSION SK(NT,ND),A(NT,NT),B(NT) 102. CALL CLEAR(NT,NT,A) 103. DO I=1,NT 104. DO J=1,ND 105. IF((I+J-1).LE.NT) THEN 106. A(I,I+J-1)=SK(I,J) 107. ENDIF 108. ENDDO 109. ENDDO 110. IF(NWK.EQ.1.AND.NX1.EQ.NX) THEN 111. WRITE(7,'(A)') &结构总刚(列出右上三角元素)& 112. DO I=1,NT 113. WRITE(7,100) (A(I,J),J=1,NT) 114. ENDDO 115. ENDIF 116. IF (NWP.EQ.1) THEN 117. WRITE(7,'(1X,A,I3,A)') &第&,NX-NX1+1,&荷载工况 的荷载列阵& 118. DO I=1,NT 119. WRITE(7,'(1x,E11.4)') B(I) 120. ENDDO 121. ENDIF 122. DO 50 K=1,NT-1 123. DO 60 I=K+1,NT 124. C1=A(K,I)/A(K,K) 125. DO 70 J=I,NT 126. A(I,J)=A(I,J)-C1*A(K,J) 127. 70 CONTINUE 128. B(I)=B(I)-C1*B(K) 129. 60 CONTINUE 130. 50 CONTINUE 131. B(NT)=B(NT)/A(NT,NT) 132. DO 80 I=NT-1,1,-1 133. DO 90 J=I+1,NT 134. B(I)=B(I)-A(I,J)*B(J) 135. 90 CONTINUE 136. B(I)=B(I)/A(I,I) 137. 80 CONTINUE 138. 100 FORMAT(1x, ) 139. END 140. 141. 142. SUBROUTINE CHECK(NT,ND,SK) 143. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 144. DIMENSION SK(NT,ND) 145. M=0 146. DO I=1,NT 147. IF(SK(I,1).LE.0.0) THEN 148. M=M+1 149. ENDIF 150. ENDDO 151. IF (M.GT.0) THEN 152. WRITE(*,*) &总刚矩阵的对角元素为负值!!& 153. STOP . ENDIF 155. END 156. 157. SUBROUTINE STRESS(NE,NT,V,A1,A2,IJK,BCA,B,STR) 158. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 159. DIMENSION D1(6),D2(3),B(NT),S(3,6),IJK(3,NE),BCA(7,NE),STR(3,NE) 160. CALL CLEAR(3,NE,STR) 161. DO 10 K=1,NE 162. DO I=1,3 163. II=IJK(I,K) 164. D1(2*(I-1)+1)=B(2*(II-1)+1) 165. D1(2*(I-1)+2)=B(2*(II-1)+2) 166. ENDDO 167. CALL CLEAR(3,6,S) 168. C1=2*A1/BCA(7,K) 169. DO 20 I=1,3 170. S(1,2*(I-1)+1)=C1*BCA(I,K) 171. S(1,2*(I-1)+2)=C1*V*BCA(I+3,K) 172. S(2,2*(I-1)+1)=C1*V*BCA(I,K) 173. S(2,2*(I-1)+2)=C1*BCA(I+3,K) 174. S(3,2*(I-1)+1)=C1*A2*BCA(I+3,K) 175. S(3,2*(I-1)+2)=C1*A2*BCA(I,K) 176. 20 CONTINUE 177. CALL MATMUL(3,6,1,S,D1,D2) 178. DO I=1,3 179. STR(I,K)=D2(I) 180. ENDDO 181. 10 CONTINUE 182. END 183. 184. SUBROUTINE OUTPUT(IP,NG,NE,NB,NX,MB,DELD,TOD,DELST,TOST,DELSUP,TOTSUP) 185. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 186. DIMENSION MB(2,NB),DELD(2,NG,NX),TOD(2,NG) 187. DIMENSION DELST(3,NE,NX),TOST(3,NE),DELSUP(NB,NX),TOTSUP(NB) 188. CHARACTER VE*6 189. WRITE(6,20) IP 190. WRITE(6,30) 191. DO I=1,NG 192. WRITE(6,40) I,DELD(1,I,IP),TOD(1,I),DELD(2,I,IP),TOD(2,I) 193. ENDDO 194. WRITE(6,50) 195. DO J=1,NE 196. WRITE(6,60) J,(DELST(L,J,IP),TOST(L,J),L=1,3) 197. ENDDO 198. WRITE(6,70) 199. DO J=1,NB 200. IF (MB(2,J).EQ.1) THEN 201. VE='x 方向' 202. ELSE 203. VE='Y 方向' 204. ENDIF 205. WRITE(6,80) MB(1,J),VE,DELSUP(J,IP),TOTSUP(J) 206. ENDDO 207. 20 FORMAT(/1X, '荷载工况=', I3, '的计算结果') 208. 30 FORMAT(/, 1X, '结点号', 5X, 'X 位移增量', 5X, 'X 累计位移', 5x, & 209. & 'Y 位移增量', 5X, 'Y 累计位移') 210. 40 FORMAT(1X, I3, 6X, F10.4, 4x, F10.4, 4X, F10.4, 4x, F10.4) 211. 50 FORMAT(/, 1X, '单元号', 6x, 'X 应力增量', 6x, 'X 累计应力', 6x, & 212. & 'Y 应力增量', 6x, 'Y 累计应力', 6x, '剪应力增量', 6x, & 213. & '累计剪应力') 214. 60 FORMAT(2X, I3, 6X, E10.3, 5X, E10.3, 5X, E10.3, 5X, E10.3, 6X, E10.3, & 215. & 6X, E10.3) 216. 70 FORMAT(/, 1X, '支座结点号', 6x, '支反力方向', 6x, '支反力增量', & 217. & 6x, '支反力累计量') 218. 80 FORMAT(5X, I3, 12X, A, 6x, E10.3, 8X, E10.3) 219. END 220. 221. 222. 223. 224. 225. 226. 227. 228. 229. 230. 231. 232. 233. 234. 235. 236. 237. 238. 239. 240. 241. 242.100SUBROUTINE MATMUL(M,N,L,A,B,C) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION A(M,N),B(N,L),C(M,L) DO 100 I=1,M DO 100 J=1,L C(I,J)=0.0 DO 100 K=1,N C(I,J)=C(I,J)+A(I,K)*B(K,J) END SUBROUTINE CLEAR(M,N,A) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION A(M,N) DO I=1,M DO J=1,N A(I,J)=0.0 ENDDO ENDDO END SUBROUTINE ABC(NG,NE,IJK,XY,BCA,NWA) 243. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 244. DIMENSION X(2,5),XY(2,NG),IJK(3,NE),B(7),BCA(7,NE) 245. DO 10 I=1,NE 246. DO 20 K=1,3 247. DO 20 J=1,2 248. X(J,K)=XY(J,IJK(K,I)) 249. 20 CONTINUE 250. DO 30 J=1,2 251. X(J,4)=X(J,1) 252. X(J,5)=X(J,2) 253. 30 CONTINUE 254. DO 40 K=1,3 255. B(K)=X(2,K+1)-X(2,K+2) 256. B(K+3)=X(1,K+2)-X(1,K+1) 257. 40 CONTINUE 258. B(7)=(B(1)*B(5)-B(4)*B(2))*0.5 259. IF (B(7).LE.0.0) THEN 260. WRITE(*,*) &单元号=&,I 261. WRITE(*,*) &单元的结点 I,J AND K 编号： &,(IJK(K,I),K=1,3) 262. STOP . ELSE 264. DO J=1,7 265. BCA(J,I)=B(J) 266. ENDDO 267. ENDIF 268. 10 CONTINUE 269. IF(NWA.EQ.1) THEN 270. WRITE(7,100) 271. DO I=1,NE 272. WRITE(7,110) I,(BCA(J,I),J=1,7) 273. ENDDO 274. ENDIF 275. 100 FORMAT(1X, '单元号', 3x, 'bi', 7x, 'bj', 7x, 'bk', 7x, 'ci', 7x, & 276. & 'cj', 7x, 'ck', 7x, '面积 A') 277. 110 FORMAT(2X, I3, 3X, F6.3, 3X, F6.3, 3X, F6.3, 3X, F6.3, 3X, F6.3, & 278. & 3X, F6.3, 3X, F6.3) 279. END 280. 281. SUBROUTINE SUMK(K,NE,NT,ND,IJK,EK,SK) 282. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 283. ,ND) 284. 285. 286. 287. 288. 289. 290. 291. 292. 293. 294. 295. 296. 297. 298. 299. 20 300. 10 301. 302. 303. 304. 305. 306. 307. 308. 309. 310. 311. 312. 313. 314. BCA(J+3,K) 315. ,K)*BCA(J,K) 316. )*BCA(J,K) 317. 318. )*BCA(J+3,K) 319. 20 320. 10 321.DIMENSIONIJ(3),EK(6,6),IJK(3,NE),SK(NTDO I=1,3 IJ(I)=IJK(I,K) ENDDO DO 10 I=1,3 DO 20 J=1,3 IF (IJ(I).GT.IJ(J)) GOTO 20 M=2*(IJ(I)-1)+1 N=2*(IJ(J)-1)+1-(2*(IJ(I)-1)+1)+1 MO=2*(I-1)+1 NO=2*(J-1)+1 SK(M,N)=SK(M,N)+EK(MO,NO) SK(M,N+1)=SK(M,N+1)+EK(MO,NO+1) SK(M+1,N)=SK(M+1,N)+EK(MO+1,NO+1) IF (IJ(I).EQ.IJ(J)) GOTO 20 SK(M+1,N-1)=SK(M+1,N-1)+EK(MO+1,NO) CONTINUE CONTINUE END SUBROUTINE KE(K,NE,T,V,A1,A2,IJK,BCA,EK,NWE) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION IJ(3),IJK(3,NE),BCA(7,NE),EK(6,6) CALL CLEAR(6,6,EK) DO I=1,3 IJ(I)=IJK(I,K) ENDDO DO 10 I=1,3 DO 20 J=I,3 IO=2*(I-1) JO=2*(J-1) EK(IO+1,JO+1)=BCA(I,K)*BCA(J,K)+A2*BCA(I+3,K)* EK(IO+1,JO+2)=V*BCA(I,K)*BCA(J+3,K)+A2*BCA(I+3 EK(IO+2,JO+2)=BCA(I+3,K)*BCA(J+3,K)+A2*BCA(I,K IF (I.EQ.J) GOTO 20 EK(IO+2,JO+1)=V*BCA(I+3,K)*BCA(J,K)+A2*BCA(I,K CONTINUE CONTINUE DO I=2,6 322. 323. 324. 325. 326. 327. 328. 329. 330. 331. 332. 333. 334. 335. 336. 337. 338. 339. 340. 341. 342. 343. 344. 345. 346. 347. 348. 349. 350. 351. 352. 353. 354. 355. 356. 357. 358. 359. 360. 361. 362. 363. 364. 365.40 50DO J=1,I-1 EK(I,J)=EK(J,I) ENDDO ENDDO DO I=1,6 DO J=1,6 EK(I,J)=EK(I,J)*A1*T/BCA(7,K) ENDDO ENDDO IF(NWE.EQ.1) THEN WRITE(7,40) K, &号单元的单刚& DO I=1,6 WRITE(7,50) (EK(I,J),J=1,6) ENDDO ENDIF FORMAT(1X, I3, A) FORMAT(1X, 6(E10.4, 4X)) END SUBROUTINE PF(NF,NT,B) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION MF(2,NF),ZF(NF),B(NT) DO I=1,NF READ(5,*) K,MF(1,I),MF(2,I),ZF(I) ENDDO WRITE(6,'(/A,I3)') &集中荷载个数 NF=&,NF WRITE(6,20) DO I=1,NF WRITE(6,30) MF(1,I),MF(2,I),ZF(I) ENDDO DO 10 I=1,NF II=2*(MF(1,I)-1)-MF(2,I)+2 B(II)=B(II)+ZF(I) CONTINUE FORMAT(1X, '结点号', 4x, '方向', 4x, '荷载值') FORMAT(2X, I3, 5X, I3, 4X, E10.4) END SUBROUTINE FNAME(PN,FN2,FN) CHARACTER PN*40,FN2*4,FN*12 去掉 PN 中前面的空格 DO 10 I=1,40 IF(PN(I:I).EQ.' ') GOTO 10 IP=I10 20 30! 366. GOTO 20 367. 10 CONTINUE 368. 20 CONTINUE 369. FN(1:8)=PN(IP:IP+7) 370. ! 去掉 FN 中后面的空格 371. DO 30 I=8,1,-1 372. IF(FN(I:I).EQ.' ') GOTO 30 373. IL=I 374. GOTO 40 375. 30 CONTINUE 376. 40 CONTINUE 377. ! 生成文件名 FN=PN+FN2 378. FN(IL+1:IL+4)=FN2(1:4) 379. END 380. 381. SUBROUTINE PP(NP,NG,NT,T,XY,B) 382. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 383. DIMENSION MP(2,NP),ZP(NP),XY(2,NG),B(NT) 384. DO I=1,NP 385. READ(5,*) K,MP(1,I),MP(2,I),ZP(I) 386. ENDDO 387. WRITE(6,'(/A,I3)') &分布荷载数 NP=&,NP 388. WRITE(6,20) 389. DO I=1,NP 390. WRITE(6,30) I,MP(1,I),MP(2,I),ZP(I) 391. ENDDO 392. DO 10 I=1,NP 393. II=2*(MP(1,I)-1) 394. JJ=2*(MP(2,I)-1) 395. PX=0.5*ZP(I)*T*(XY(2,MP(1,I))-XY(2,MP(2,I))) 396. PY=0.5*ZP(I)*T*(XY(1,MP(2,I))-XY(1,MP(1,I))) 397. B(II+1)=B(II+1)+PX 398. B(II+2)=B(II+2)+PY 399. B(JJ+1)=B(JJ+1)+PX 400. B(JJ+2)=B(JJ+2)+PY 401. 10 CONTINUE 402. 20 FORMAT(1X, '序号', 4x, '首结点号', 4x, '末结点 号', 7x, '荷载值') 403. 30 FORMAT(1X, I3, 8x, I3, 9X, I3, 7X, F10.3) 404. END 405. 406. SUBROUTINE PG(NM,NE,NT,DENSITY,T,IJK,BCA,B) 407. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 408. DIMENSION IJK(3,NE),BCA(7,NE),MG(NM),B(NT) 409. 410. 411. 412. 413. 414. 415. 416. 417. 418. 419. 420. 421. 422. 423. 424. 425. 426. 427. 428. 429. 430. 431. 432. 433. 434. 435. 436. 437. 438. 439. 440. 441. 442. 443. 444. 445. 446. 447. 448. 449. 450. 451. 452.10 100IF (NM.GE.NE) THEN NM=NE DO I=1,NM MG(I)=I ENDDO ELSE READ(5,*) (MG(I),I=1,NM) ENDIF WRITE(6,'(/A,I3)') &计自重单元数 NM=&,NM WRITE(6,'(A)') &计自重的单元号为：& K1=0 K2=0 DO WHILE((NM-K2).GE.10) WRITE(6,100) (MG(K2+J),J=1,10) K1=K1+1 K2=10*K1 ENDDO WRITE(6,100) (MG(K2+J),J=1,NM-K2) DO 10 I=1,NM II=2*(IJK(1,MG(I))-1) JJ=2*(IJK(2,MG(I))-1) KK=2*(IJK(3,MG(I))-1) PE1=-DENSITY*T*BCA(7,MG(I))/3.0 B(II+2)=B(II+2)+PE1 B(JJ+2)=B(JJ+2)+PE1 B(KK+2)=B(KK+2)+PE1 CONTINUE FORMAT(10(1X,I3)) END SUBROUTINE INPUT(NG,NE,NB,IJK,XY,MB,ZB) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION IJK(3,NE),XY(2,NG),MB(2,NB),ZB(NB) DO I=1,NE READ(5,*) K,(IJK(J,I),J=1,3) ENDDO DO I=1,NG READ(5,*) K,XY(1,I),XY(2,I) ENDDO DO I=1,NB READ(5,*) K,MB(1,I),MB(2,I),ZB(I) ENDDO WRITE(6,15) DO K=1,NE 453. WRITE(6,20) K,(IJK(J,K),J=1,3) 454. ENDDO 455. WRITE(6,30) 456. DO K=1,NG 457. WRITE(6,35) K,XY(1,K),XY(2,K) 458. ENDDO 459. WRITE(6,40) 460. DO K=1,NB 461. WRITE(6,45) K,MB(1,K),MB(2,K),ZB(K) 462. ENDDO 463. 15 FORMAT(/1X, '单元号', 4X, 'I 结点号', 4x, 'J 结点号', 4x, 'K 结点号') 464. 20 FORMAT(2X, I3, 7X, I3, 8X, I3, 8X, I3) 465. 30 FORMAT(/1X, '结点号', 8X, 'X 坐标', 8X, 'Y 坐标 ') 466. 35 FORMAT(2X, I3, 6X, F8.3, 5X, F8.3) 467. 40 FORMAT(/1X, '序号', 4x, '支承结点号', 6x, '支 承方向', 6x, & 468. & '支承位移') 469. 45 FORMAT(2X, I2, 8X, I4, 11X, I2, 7X, F10.3) 470. END 471. 472. SUBROUTINE DBC(NB,ND,NT,NX,NX1,MB,ZB,SK,B) 473. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 474. DIMENSION MB(2,NB),ZB(NB),SK(NT,ND),B(NT) 475. DO 5 I=1,NB 476. N=2*MB(1,I)-MB(2,I) 477. Z=ZB(I) 478. IF (ABS(Z).LT.1.0E-10) THEN 479. IF (NX.NE.NX1) GOTO 30 480. SK(N,1)=1.0 481. DO J=2,ND 482. SK(N,J)=0.0 483. ENDDO 484. DO 10 K=2,ND 485. IF (N.LT.K) GOTO 20 486. M=N-K+1 487. SK(M,K)=0.0 488. 10 CONTINUE 489. 20 CONTINUE 490. 30 B(N)=0.0 491. ELSE 492. IF (NX.NE.NX1) GOTO 40 493. SK(N,1)=SK(N,1)*1.0E+15 494. 40 B(N)=SK(N,1)*Z 495. ENDIF 496. 5 CONTINUE 497. END 498. 499. SUBROUTINE CALND(NG,NE,IJK,ND) 500. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 501. DIMENSION IJK(3,NE) 502. ND=0 503. DO 10 I=1,NG 504. M1=0 505. DO 20 J=1,NE 506. DO 30 K=1,3 507. IF(IJK(K,J).EQ.I) THEN 508. M2=0 509. DO 40 L=1,3 510. M3=IJK(L,J)-I 511. IF(M3.GT.M2) THEN 512. M2=M3 513. ENDIF 514. 40 CONTINUE 515. GOTO 25 516. ENDIF 517. 30 CONTINUE 518. 25 IF(M2.GT.M1) M1=M2 519. 20 CONTINUE 520. IF(M1.GT.ND)ND=M1 521. 10 CONTINUE 522. ND=2*(ND+1) 523. WRITE(7,100) ND 524. 100 FORMAT(1X,'总刚矩阵的半带宽 ND 为：',I4) 525. END 526. 527. SUBROUTINE TRES(IP,NG,NE,NX,NT,B,STR,DELD,TOD,DELST,TOST) 528. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 529. DIMENSION B(NT),STR(3,NE),DELD(2,NG,NX),TOD(2,NG),DELST(3,NE,NX),TOST(3,NE) 530. DO I=1,NG 531. IO=2*(I-1) 532. DELD(1,I,IP)=B(IO+1) 533. DELD(2,I,IP)=B(IO+2) 534. TOD(1,I)=TOD(1,I)+DELD(1,I,IP) 535. TOD(2,I)=TOD(2,I)+DELD(2,I,IP) 536. ENDDO 537. DO I=1,NE 538. DO J=1,3 539. DELST(J,I,IP)=STR(J,I) 540. TOST(J,I)=TOST(J,I)+DELST(J,I,IP) 541. ENDDO 542. ENDDO 543. END 544. 545. 546. SUBROUTINE SUPT(IP,NG,NE,NB,NX,T,V,A1,A2,IJK,MB,BCA,DELD,DELSUP,TOTSUP) 547. IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) 548. DIMENSION IJK(3,NE),MB(2,NB),BCA(7,NE),DELD(2,NG,NX),EK(6,6) 549. DIMENSION D1(6),DELSUP(NB,NX),TOTSUP(NB) 550. DO 10 K=1,NB 551. I1=MB(1,K) 552. I2=MB(2,K) 553. RF=0.0 554. DO 20 I=1,NE 555. DO 30 J=1,3 556. I3=IJK(J,I) 557. IF (I1.NE.I3) THEN 558. GOTO 30 559. ELSE 560. I4=2*(J-1)+2-I2 561. CALL KE(I,NE,T,V,A1,A2,IJK,BCA,EK,0) 562. DO M=1,3 563. IO=IJK(M,I) 564. D1(2*(M-1)+1)=DELD(1,IO,IP) 565. D1(2*(M-1)+2)=DELD(2,IO,IP) 566. ENDDO 567. DO L=1,6 568. RF=RF+EK(I4,L)*D1(L) 569. ENDDO 570. GOTO 20 571. ENDIF 572. 30 CONTINUE 573. 20 CONTINUE 574. DELSUP(K,IP)=DELSUP(K,IP)+RF 575. TOTSUP(K)= TOTSUP(K)+DELSUP(K,IP) 576. 10 CONTINUE 577. END 578. 579. 580. 581. 582. 583. 584. 585. 复制代码SUBROUTINE CLEAR1(N,B) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSION B(N) DO I=1,N B(I)=0.0 ENDDO ENDPFEM,EXE 平面问题有限元程序的数据文件输入格式： 基本输入数据文件名：*.DAT 计算结果输出数据文件名：*.OUT 中间参数输出数据文件名：*.OU1 下面按在文件中的前后顺序进行说明： 1 基本控制参数信息：NG,NE,MC,NX,NB,EO,VO,DENSITY ,T(共计 5 个整形数，4 个实型数) NG：结构的结点总数； NE：结构的单元总数； MC：平面问题的类型，MC=0，为平面应力，MC=1，为平面应变； NX：荷载工况数； NB：支承位移数； EO：材料弹性模量(Pa)； VO：材料泊松比； DENSITY ：容重(N/m3) T ：材料厚度(m)； 2 打印输出控制参数：NWA,NEW,NWK,NWP(4 个整形数) 等于 1 时，输出，否则不输出。 3 单元结点信息：(K,(IJK(I,K),I=1,3),K=1,NE) (每行 4 个整形数，共计 NE 行) K：单元号； IJK(1,K)：K 单元 I 结点编号； IJK(2,K)：K 单元 J 结点编号； IJK(3,K)：K 单元 K 结点编号； 4 结点坐标信息：((K,XY(1,K),XY(2,K)),K=1,NG)(每行 3 个整形数，共计 NG 行) K：结点号 XY(1,K)：K 结点 X 坐标； XY(2,K)：K 结点 Y 坐标； 5 支承信息：((K,MB(1,K),MB(2,K),ZB(K)),K=1,NB)(每行 3 个整形数，1 个实型数，共计 NB 行) K：支承位移序号； MB(1,K)：第 K 个支承位移所在的结点号； MB(2,K)：第 K 个支承位移的坐标方向； ZB(K)： 第 K 个支承位移的数值； 6 按 NX 荷载工况数输入荷载信息：每一荷载工况如下 (1) NF,NP,NM(3 个整型数) NF：集中荷载个数； NP：分布荷载个数； NM：计自重单元数； (2) 若 NF≠0，则输入下面数据 K,MF(1,K),MF(2,K),ZF(K)(每行 3 个整形数，1 个实型数，共计 NF 行) K：集中荷载序号； MF(1,K)：第 K 个集中荷载作用的结点号； MF(2,K)：第 K 个集中荷载的坐标方向； ZF(K)： 第 K 个集中荷载的数值； (3) 若 NP≠0，则输入下面数据 K,MP(1,K),MP(2,K),ZP(K)(每行 3 个整形数，1 个实型数，共计 NP 行) K：分布荷载序号； MP(1,K)：第 K 个分布荷载作用的结点号； MP(2,K)：第 K 个分布荷载的坐标方向； ZP(K)： 第 K 个分布荷载的数值； (4) 若 NM≠0，则输入下面数据 若 NM≥NE，则表示计所有单元的自重，不需输入计自重的单元号； 若 NM&NE，则需要输入计自重的单元号；
里面包含结构的计算简图以及源程序算例一计算简图及结果输出用平面刚架静力计算...Fortran语言 有限元程序... 18页 1下载券 Fortran 讲解-平面有限元... 17页...Fortran 讲解-平面有限元源... 17页 1财富值 Fortran语言 有限元程序分... ...有限元程序设计005 15页 免费如要投诉违规内容,请到百度文库投诉中心;如要提出功...Fortran语言编写的有限元结构程序_能源/化工_工程科技_专业资料。算例一计算简图...Fortran 讲解-平面有限元... 17页 1下载券 王勖成《有限单元法》源... 16页...结构力学刚架问题的有限元分析 1 有限元分析软件的开发 1.1 程序功能该程序为...Fortran_讲解-平面有限元... 17页 免费
平面弹性力学有限元源程... 暂无评价...读有限元Fortran程序笔记_计算机软件及应用_IT/计算机_专业资料。读有限元程序笔记 1. ALLOCATABLE::COORD(:,:),PROPS(:,:,:) PROPS(:,:,:)是三维数组。 2...有限元编程算例(fortran) - 本程序通过Fortran语言编写,程序在Intel Parallel Studio XE 2013 with VS2013中成功运行,程序为《计算力学》(龙...
All rights reserved Powered by
www.tceic.com
copyright ©right 。文档资料库内容来自网络，如有侵犯请联系客服。}
我要回帖
更多关于 不想做程序员了 的文章
更多推荐
- ·如何查看运费险从哪里查看赔没赔付已理赔多少
- ·政府一图看懂四套班子子是哪一图看懂四套班子子
- ·自助火锅烤肉一体店加盟培训地址？
- ·从事电商好做吗服务的中小企业如何打造自己的竞争优势？
- ·微信怎么召唤如何不充钱获得亚索皮肤？
- ·51电销机器人的弊端售后服务怎么样
- ·组长叫我他母亲的房子去打别人房子有罪吗
- ·周六平遥到五台山多少公里个人对公户打款
- ·关于支付宝的各种问题宝
- ·手机突然收到短信，非税代扣支付宝支出短信提醒人民币1200元，这是什么意思？为什么银行卡上无缘无故
- ·明星大侦探3网盘第一季百度云
- ·可以通过平台将蘑菇街开放平台网店转让出去吗
- ·艾薇资源百度云自取 要钱绕道 百度云最好 感谢
- ·空调遥控器温度显示不显示温度海尔空调遥控器温度显示不显示温度
- ·求助一下，switch怎么注册日服日服突然登不上了
- ·诺基亚7plus 香港报价，在香港买会不会便宜点
- ·更新到8.1后，手机6splus待机一夜耗电大概都是费多少电
- ·怎么在javajava中读取配置文件件中配置Hive的数据源
- ·陈建州超跑思域什么梗遭追尾是查查查钱网第一个报道吗？
- ·Fortran不想做程序员了如下，我想对NG,NE,MC等等赋值该怎么做
- ·谁能给我普及0911对战平台天梯地图的天梯分阿
- ·不能忘却的长征观后感联通4G 能用么
- ·霞姐是谁糖水怎么样，好不好的默认点评
- ·泊车是否有空位LED屏一目了然的了
- ·非自然死亡全集百度云
- ·电视剧《硬骨头之绝地归途电视剧》电脑_手机在线看_百度云超清资源缓存
- ·如何使用R语言中的MASS第三方包获取数据
- ·荣耀magicbook充电器的充电器能不能给荣耀8充电
- ·如何查看手机上查看微信小程序怎么做demo功能
- ·喜欢聊天好不好
- ·他43岁了，我29了，他说他爱我，发动画拥抱表情能随便发吗说你是我的小宝贝，我可以信他吗
- ·手机号就是微信号怎么改成手机号登陆说没有这个微信号
- ·红米5plus会迎来安卓8.0更新吗
- ·开始使用Maven组件插件的从麻烦别人开始问题，怎么解决
