| 查看: 1119 | 回复: 5 | |||
[交流]
Fortran平面桁架有限元程序
|
|||
|
program truss_2D use prep use solve implicit none integer :: i,j,m,n,nel,nne,nn,nodof,edof,gdof integer :: e2s(4) real :: L,FN real :: kl(4,4),kg(4,4),T(4,4) integer,allocatable :: elemNodes (:,:） ,nf(:,:） real ,allocatable :: Coords (:,:）,prop(:,:）, & KK(:,:）,loads(:,:） real ,allocatable :: nodedisp(:,:）,F(:） ,edg(:）, & fg(:）,fl(:）,delta(:） open ( 10,file = 'data.txt' ) open ( 11,file = 'out.txt' ) read (10,*) nel ! Number of elements read (10,*) nne ! Number of nodes per element allocate ( elemNodes(nel,nne),prop(nel,nne) ) read (10,*) nn ! Number of nodes read (10,*) nodof ! Number of degrees of freedom per node edof = nodof * nne allocate ( Coords(nn,nodof),nf(nn,nodof), & loads(nn,nodof),nodedisp(nn,nodof),edg(edof), & fg(edof),fl(edof) ) read (10,*) ( (elemNodes(i,j), j=1,nne), i=1,nel ) read (10,*) ( (prop(i,j), j=1,nne), i=1,nel ) read (10,*) ( (Coords(i,j), j=1,nodof), i=1,nn ) read (10,*) ( (nf(i,j), j=1,nodof), i=1,nn ) read (10,*) ( (loads(i,j), j=1,nodof), i=1,nn ) gdof = 0 do i = 1,nn do j = 1,nodof if ( nf(i,j)/=0 ) then gdof = gdof + 1 nf(i,j) = gdof end if end do end do allocate ( KK(gdof, gdof), F(gdof),delta(gdof) ) F = 0. call truss_F( m,n,nn,nodof,nf,gdof,loads,F) e2s = 0 KK = 0. do i = 1,nel call truss_T(i,nel,nne,nn,nodof,elemNodes,Coords,L,T) call truss_kl (i,nel,nne,L,prop,kl) call truss_kg (T,kl,kg) call truss_e2s(i,j,nn,nel,nne,nodof,nf,elemNodes,e2s) call form_KK (m,n,edof,gdof,kg,e2s,KK) end do call fem_Solver(KK,F,gdof,delta) ! 求解 nodedisp = 0. forall ( i = 1:nn,j = 1:nodof,nf(i,j)/=0 ) nodedisp(i,j) = delta( nf(i,j) ) end forall do i = 1,nel call truss_T(i,nel,nne,nn,nodof,elemNodes,Coords,L,T) call truss_kl (i,nel,nne,L,prop,kl) call truss_kg (T,kl,kg) call truss_e2s(i,j,nn,nel,nne,nodof,nf,elemNodes,e2s) edg = 0. do j = 1,edof if ( e2s(j) /= 0 ) then edg(j) = delta( e2s(j) ) end if end do fg = matmul( kg,edg ) fl = matmul( T,fg ) FN = fl(3) write (11,100) i write (11,200) FN end do 100 format (/,T10,'单元',I2 ) 200 format (T10,'轴力=',F18.4) end program truss_2D |
回复此楼» 猜你喜欢
材料考研调剂
已经有3人回复
-
材料调剂
已经有12人回复
-
英一数一408，总分284，二战真诚求调剂
已经有14人回复
-
085410 一志愿211 22408分数359求调剂
已经有4人回复
-
271求调剂
已经有19人回复
-
385分 生物学（071000）求调剂
已经有3人回复
-
一志愿安徽大学计算机科学与技术学硕，331分求调剂
已经有3人回复
-
318求调剂，计算材料方向
已经有8人回复
-
291求调剂
已经有25人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有6人回复
» 抢金币啦！回帖就可以得到:
-
汕头大学胡良胜课题组招收化学、化工、材料方向调剂研究生3名
+1/185
-
中国科学院过程工程研究所 能源环境课题组 博士后招聘启示
+1/184
-
陕西理工大学资源与环境、材料化工方向部分调剂指标-理工兼收
+1/181
-
中国科学院深圳先进院唐永炳研究员团队诚聘电化学/电池方向博士后
+1/43
-
河南工业大学化学化工学院 2026 年硕士调剂正在进行！
+1/40
-
新疆碳纳米材料制备与应用重点实验室2026年招聘启事
+1/34
-
【青岛大学】2026年生物与医药申请考核制博士生招生（含少数民族骨干人才）
+1/19
-
【哈工程】核学院某课题组招收全日制调剂研究生
+1/17
-
就业好-广东唯一石化院校——资源与环境专硕招生
+1/9
-
西北工业大学机电学院 复合材料加工制造方向招聘博士后
+1/8
-
课题组招收专业代码08开头的学硕！
+1/7
-
调剂招生
+1/7
-
苏州大学医学部纳米生物医学方向招收调剂生1名
+1/7
-
招聘电池方向博士
+1/7
-
招聘︱国家纳米科学中心丁宝全杰青团队诚聘博士后（特别研究助理）
+2/6
-
中国石油大学（北京）石油工程专业-开发方向课题组 2026年博士招生
+1/6
-
招收2026级研究生
+1/2
-
2026年北京石油化工学院环境学科-水污染控制工程方向课题组欢迎您
+1/2
-
Top-88悉尼科技大学数据科学/AI 招收2027年入学 校奖 博士生1到2名(国际和本地学生）
+1/2
-
英国剑桥好房出租
+1/1
简单回复
dsctg2楼
2017-02-24 11:48
回复
springer_(金币+1): 谢谢参与
2017-02-24 14:30
回复
springer_(金币+1): 谢谢参与
是 发自小木虫IOS客户端
2017-11-10 23:54
回复
2018-01-11 19:19
回复
springer_(金币+1): 谢谢参与
是 发自小木虫Android客户端
2018-01-27 19:42
回复
springer_(金币+1): 谢谢参与
顶 发自小木虫IOS客户端
