| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 2263 人关注过本帖
标题:版主求救!关于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
快速回复:版主求救!关于fortran有限元编程
数据加载中...
 
   



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

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