- 1、本文档共26页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
線弹性小变形平面板壳静力有限元计算程序
元计算有限元自动生成系统所开发源代码系列
线弹性小变形平面板壳静力有限元计算程序
简介
元计算(www.ectec.asia)公司所开发的并行有限元程序自动生成系统(pFEPG)可根据用户需要开发出各种有限元计算程序源代码。该源代码系列即为pFEPG所开发出来的求解各学科典型问题的有限元计算程序。
该组程序为线弹性小变形平面板壳静力有限元计算程序。
starta.for,对位移场的数据进行初始化;
implicit real*8 (a-h,o-z)
character*12 fname,filename(20)
common /aa/ ia(250000000)
common /bb/ ib(125000000)
c.... open disp0 file to get the numbers of nodes and degree of freedom
c.... knode .... number of nodes, kdgof .... number of d.o.f.
open(1,file= ,form=unformatted)
read(1) knode,kdgof
close(1)
kvar=knode*kdgof
write(*,*) knode,kdgof,kvar =
write(*,(1x,4i7)) knode,kdgof,kvar
kvar1=kvar+1
kcoor=3
kelem knb1=kdgof*knode*1
if (knb1/2*2 .lt. knb1) knb1=knb1+1
kna4=kcoor*knode*2
kna1=kdgof*knode*2
kna2=kdgof*knode*2
kna3=kdgof*knode*2
kna5=knode*1
if (kna5/2*2 .lt. kna5) kna5=kna5+1
knb4=kelem*1
if (knb4/2*2 .lt. knb4) knb4=knb4+1
knb2=kvar1*1
if (knb2/2*2 .lt. knb2) knb2=knb2+1
knb3=kvar1*1
if (knb3/2*2 .lt. knb3) knb3=knb3+1
kna0=1
kna1=kna1+kna0
kna2=kna2+kna1
kna3=kna3+kna2
kna4=kna4+kna3
kna5=kna5+kna4
if (kna5-1.gt.250000000) then
write(*,*) exceed memory of array ia
write(*,*) memory of ia = 250000000
write(*,*) memory needed = ,kna5, in prgram start
stop 55555
endif
knb0=1
knb1=knb1+knb0
knb2=knb2+knb1
knb3=knb3+knb2
knb4=knb4+knb3
if (knb4-1.gt.125000000) then
write(*,*) exceed memory of array ib
write(*,*) memory of ib = 125000000
write(*,*) memory needed = ,knb4, in prgram start
stop 55555
endif
call start(knode,kdgof,kcoor,kvar,
*kelem,maxt,kvar1,ia(kna0),ia(kna1),ia(kna2),
*ia(kna3),ia(kna4),ib(knb0),ib(knb1),ib(knb2),
*ib(knb3),
*filename)
end
subroutine start(knode,kdgof,kcoor,kvar,
*kelem,maxt,kvar1,u0
您可能关注的文档
最近下载
- Blender软件二次开发:Blender二次开发综合案例分析.pdf VIP
- 国标图集05J909《工程做法》—标准图集.pdf
- Blender软件二次开发:Blender游戏引擎与交互式内容开发.pdf VIP
- 孔子诞辰日致辞稿(共7篇) .doc
- 第一单元第四课《上学路上》课件 湘美版(2024)一年级上册.pptx
- 浙教版信息科技五上 第4课 算法中的数据 课件.pptx
- 经济全球化对国家经济安全的影响研究.pptx
- 聊城大学音乐学院中国音乐史课件 西周春秋战国音乐.ppt
- Blender软件二次开发:BlenderPython脚本优化与调试.pdf VIP
- 第二课 美的本质和特征课件 --中职技校《美育》第四版.pptx
文档评论(0)