标题:版主求救!关于fortran有限元编程
取消只看楼主
lengend
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2008-10-21
 问题点数:0 回复次数:1 
版主求救!关于fortran有限元编程
最近在学习有限元,下面是书上fortran语言编写的一段程序,以前没学过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))
             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            FORMAT(1X,6E11.4)
60            RETURN
            END
搜索更多相关主题的帖子: 有限元 fortran 版主 
2008-10-21 17:57
lengend
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2008-10-21
得分:0 
回复 4# cnhanxiao 的帖子
非常有用,多谢!
2008-10-22 08:56



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-239667-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.959609 second(s), 8 queries.
Copyright©2004-2025, BCCN.NET, All Rights Reserved