下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
1、Fotran90版一平面剛架有限元分源程序代碼program mainreal , allocatable :ks(:,:)allocatable lnd(:,:)allocatable crd(:,:)allocatable ea(:)allocatable ei(:)allocatable jcs(:,:)allocatable pj(:,:)allocatable bl(:)allocatable p(:)open (5,file= inputdates.in )read (5,*) ne,nj,ns,npjnj3=3*njallocate (lnd(ne,2)allocate (crd
2、(nj,2)allocate (ea(ne)allocate (ei(ne)allocate (jcs(ns,4)allocate (pj(npj,4)allocate (bl(ne)allocate (ks(nj3,nj3)allocate (p(nj3)write (*, (1x,plane fram structural analysis/) )write (*, (1x,structural parameters/) )write (*, (/1x,total number of)write (*, (1x,element=,i5/1x,joints=,i5/1x,constructe
3、d joints=,i5/1x,loads=,i5/)&ne,nj,ns,npjcall readin(ne,nj,ns,npj,lnd,crd,ea,ei,jcs,pj,bl)call formf(nj3,npj,pj,p)call cks(nj3,ne,ea,ei,bl,lnd,crd,nj,ks)call dealbc(ns,jcs,nj3,ks,p)call solve(nj,nj3,ks,p)close (5)stopend !subroutine readin(ne,nj,ns,npj,lnd,crd,ea,ei,jcs,pj,bl) dimensionlnd(ne,2),crd(
4、nj,2),ea(ne),ei(ne),jcs(ns,4),pj(npj,4),bl(ne)read (5,*) (lnd(i,j),j=1,2),i=1,ne)write (*, (/1x,element dates,/1x,element,4x,conection,8x) )write (*, (1x,i5,2x,i5,3x,to,i5,3x) (i,(lnd(i,j),j=1,2),i=1,ne)read (5,*) (crd(i,j),j=1,2),i=1,nj)write (*, (1x,nodel coordinates/3x,node,6x,x-coordinates,7x,y-
5、coordinates) )write (*, (1x,i5,5x,f10.4,10x,f10.4) ) (i,(crd(i,j),j=1,2),i=1,nj)read (5,*) (ea(i),i=1,ne)read (5,*) (ei(i),i=1,ne)write (*, (/1x,material parameters,/1x,element,11x,ea,13x,ei) )write (*, (1x,i5,5x,2e15.6) ) (i,ea(i),ei(i),i=1,ne)read (5,*) (jcs(i,j),j=1,4),i=1,ns)write (*, (/1x,const
6、rained nodes,/3x,nodes,1x,X,4x,Y,4x,R) )write (*, (4i5) ) (jcs(i,j),j=1,4),i=1,ns)read (5,*) (pj(i,j),j=1,4),i=1,npj)write (*, (/1x,joints of londs,/1x,joint,5x,PX,8x,PY,8x,Mxy) )write (*, (1x,f5.0,3f10.4) ) (pj(i,j),j=1,4),i=1,npj)do ie=1,nei=lnd(ie,1)j=lnd(ie,2)dx=crd(j,1)-crd(i,1)dy=crd(j,2)-crd(
7、i,2)bl(ie)= sqrt (dx*2+dy*2)end dowrite (*, (/1x,the length of the elements,/1x,element numers,5x,numembers length)write (*, (1x,i5,10x,f10.4) ) (i,bl(i),i=1,ne)returnend !subroutine turn(ne,lnd,crd,bl,t,ie,nj)dimension lnd(ne,2),crd(nj,2),bl(ne),t(6,6)i1=lnd(ie,1)j1=lnd(ie,2) dx=crd(j1,1)-crd(i1,1)
8、dy=crd(j1,2)-crd(i1,2)si=dy/bl(ie)co=dx/bl(ie)do i=1,6t(i,1:6)=0.0end dot(1,1)=cot(1,2)=sit(2,1)=-sit(2,2)=co t(3,3)=1.0 do i=1,3 do j=1,3 t(i+3,j+3)=t(i,j) end do end do !write (*,(/1x,the dates of t) !write (*,(/1x,6f10.4) (t(i,j),j=1,6),i=1,6) return end! subroutine stif(ne,ea,ei,bl,kd,ie) dimens
9、ion ea(ne),ei(ne),bl(ne) real kd(6,6) a1=ea(ie) e1=ei(ie) s=bl(ie) do i=1,6kd(i,1:6)=0.0 end do kd(1,1)=a1/s kd(2,2)=12.0*e1/s*3 kd(3,2)=-6.0*e1/s*2 kd(3,3)=4.0*e1/s kd(4,1)=-kd(1,1) kd(4,4)=kd(1,1) kd(5,2)=-kd(2,2) kd(5,3)=-kd(3,2) kd(6,3)=2.0*e1/s kd(6,5)=-kd(3,2) kd(6,6)=kd(3,3) do i=1,6 do j=1,i
10、kd(j,i)=kd(i,j)end do end do return end ! subroutine stie(ne,ea,ei,bl,kd,lnd,crd,t,ke,ie,nj) dimension ea(ne),ei(ne),bl(ne),lnd(ne,2),crd(nj,2),t(6,6),ek(6,6) real kd(6,6),ke(6,6) call stif(ne,ea,ei,bl,kd,ie)call turn(ne,lnd,crd,bl,t,ie,nj)do i=1,6do j=1,6ek(i,j)=0.0do k=1,6 ek(i,j)=ek(i,j)+kd(i,k)*
11、t(k,j) end do end do end do do i=1,6 do j=1,6ke(i,j)=0.0do k=1,6 ke(i,j)=ke(i,j)+t(k,i)*ek(k,j) end do end do end do return end!subroutine formf(nj3,npj,pj,p) dimension p(nj3),pj(npj,4) p=0.0do i=1,npjj=pj(i,1)p(3*j-2)=p(3*j-2)+pj(i,2)p(3*j-1)=p(3*j-1)+pj(i,3)p(3*j)=p(3*j)+pj(i,4)end do!write(*,(/1x
12、,f10.4) (p(i),i=1,nj3)returnend !subroutine cks(nj3,ne,ea,ei,bl,lnd,crd,nj,ks)dimension ea(ne),ei(ne),bl(ne),lnd(ne,2),crd(nj,2),& t(6,6),ek(6,6)real ks(nj3,nj3),kd(6,6),ke(6,6) ks(1:nj3,1:nj3)=0.0 do ie=1,necall stie(ne,ea,ei,bl,kd,lnd,crd,t,ke,ie,nj)do i=1,2do ii=1,3 ir=3*(i-1)+ii iw=3*(lnd(ie,i)-
13、1)+iido j=1,2 do jj=1,3jr=3*(j-1)+jj jw=3*(lnd(ie,j)-1)+jj ks(iw,jw)=ks(iw,jw)+ke(ir,jr) end do end do end do end do end do!write (*,(/1x,dates of ks)!write (*,(/1x,9e10.3) (ks(i,j),j=1,nj3),i=1,nj3) return end! subroutine dealbc(ns,jcs,nj3,ks,p) dimension jcs(ns,4),p(nj3) real ks(nj3,nj3) do i=1,ns
14、i1=jcs(i,1) do j=1,3j1=jcs(i,j+1) if (j1=0) then elseiw=3*(i1-1)+j do k=1,nj3 if (k=iw) thenks(iw,iw)=1.0 elseks(iw,k)=0.0 ks(k,iw)=0.0 end if end dop(iw)=0.0end if end do end do !write (*,(/1x,the dates of bc) !write (*,(/1x,9e15.6,5x,e15.6) (ks(i,j),j=1,nj3),p(i),i=1,nj3) return end! subroutine so
15、lve(nj,nj3,ks,p) dimension p(nj3)real ks(nj3,nj3)do k=1,nj3-1do i=k+1,nj3c=ks(i,k)/ks (k,k)do j=k,nj3ks(i,j)=ks(i,j)-c*ks(k,j)end dop(i)=p(i)-c*p(k)end doend dop(nj3)=p(nj3)/ks(nj3,nj3)do i=nj3-1,1,-1do j=i+1,nj3p(i)=p(i)-ks(i,j)*p(j)end dop(i)=p(i)/ks(i,i)end dowrite (*, (/1x,displacement of the jo
16、ints,/1x,jointsnumbers,3x,X-dis,10x,Y-dis,10x,Z-rot)write (*, (1x,i5,5x,3e15.6) (i,p(3*i-2),p(3*i-1),p(3*i),i=1,nj)returnend!結(jié)果輸出如下:z| r pg。起 XI _ziI h | |l l| I I I I I 0 I I I I I I I I I M |l l| I I I I I H I I I I I I I I I h I |l I I I Isubroutine readin (ne njs n石 npj* Ind crd, ea, eibl Jdi me
17、ns i on Ind (nep 2) crd (nj# 2)r ea (ne.)ei (ne), j cs (ns, 4.)pj (npjf 4 J P bl (ne)read(5, (Clndd, jX j-L 2) i=L writeelement dates:, rite (f, (li,2k, 15,i :read (黑 (電機i. j)sEirite (也(1歸ncidel coardinate,write (*, llz, i5f 5k, flO. 4j 10k, fl(read(catC, i=l3ne)read(5嚴)(ei (i)* i=l,ne)wrile (*, (/I
18、sjf materi al pararet:write t*.f i5j 5k, 2el5T 6)read (5,+Jj)4),i=lwrile 住,constrained nodewrite (也氣4i5)#) (jcs(if j)F jpljread (5f *) (pj j) j=lj 4)A i=lj Jsrrite (氣氣/hj joints of Londswrite (叫 v(lzp f5. 0, 3HO. t) v) (p.eV *1: TortrilbTurkplac.TrjrcdaTryeoa-b-.gTrycd. no-e l eiHieintconstrain B-dnodes
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 遺傳病遠程診斷技術(shù)的創(chuàng)新與應(yīng)用-洞察分析
- 2025年人教版(2024)八年級地理上冊階段測試試卷
- 2025年新科版九年級生物上冊階段測試試卷
- 2025年度美團商家入駐合同及運營支持細則4篇
- 維修服務(wù)模式創(chuàng)新-洞察分析
- 2025年人教新起點九年級物理下冊月考試卷含答案
- 2025年上教版七年級生物上冊月考試卷
- 2025年外研版三年級起點九年級歷史上冊階段測試試卷含答案
- 銀幣收藏投資策略-洞察分析
- 鞋履新材料應(yīng)用趨勢分析-洞察分析
- 物業(yè)民法典知識培訓(xùn)課件
- 2023年初中畢業(yè)生信息技術(shù)中考知識點詳解
- 2024-2025學(xué)年山東省德州市高中五校高二上學(xué)期期中考試地理試題(解析版)
- 《萬方數(shù)據(jù)資源介紹》課件
- 麻風(fēng)病病情分析
- 《急診科建設(shè)與設(shè)備配置標(biāo)準》
- 第一章-地震工程學(xué)概論
- JJF(陜) 063-2021 漆膜沖擊器校準規(guī)范
- 《中國糖尿病防治指南(2024版)》更新要點解讀
- TSGD7002-2023-壓力管道元件型式試驗規(guī)則
- 2024年度家庭醫(yī)生簽約服務(wù)培訓(xùn)課件
評論
0/150
提交評論