Продолжение приложения 2

i=0

ico=0

120 i=i+1

if(i.gt.limit)go to 700

if(nstor(i).le.0)go to 120

ico=ico+1

j=nstor(i)

c write node number x,y,z coordinates

write(7,12)j,x(j),y(j),z(j)

12 format(5x,i5,/,3e12.5)

go to 120

700 write(7,710)

710 format(9x,'0')

c write the element connection

itot=0

do 720 i=1,limit

if(icon(i,1).eq.0)go to 720

itot=itot+1

720 continue

c write total number of elements

write(7,511)itot

c ity is the element type

ity=11

izero=0

do 740 i=1,limit

if(icon(i,1).eq.0)go to 740

write(7,511)ity

write(7,741)icon(i,1),icon(i,2),icon(i,3),icon(i,4),izero

*,icon(i,5),icon(i,6),izero,icon(i,7)

741 format(5x,9i5)

740 continue

write(7,710)

return

end

Приложение 3 программа s1_3.F

parameter(nstm=18,np=6,ne=4,ncn=3,nszf=36,

* ndf=3,ndfg=6)

* вычисление и запись локальных матриц жесткости и масс треугольников

dimension rtb(nstm,nstm),cord(np,3),nop(ne,ncn),

* estifm(nstm,nstm),rtb(nstm,nstm),dmgl(nstm,nstm)

open(1,file='RES1_3')

open(8,FILE='ESTIFM',form='unformatted')

open(9,FILE='DMGL',form='unformatted')

open(10,file='NOP',form='unformatted')

open(13,file='DTE',form='unformatted')

read(13,*)em,ro,pr,thin

* ro=ro/9.81

write(1,*)' em=',em,' ro=',ro,' thin=',thin,' pr=',pr

do 777 i=1,np

do 777 j=1,3

777 cord(i,j)=0.

do 778 i=1,ne

do 778 j=1,ncn

778 nop(i,j)=0

do 445 i=1,np

445 read(13,*)j,(cord(j,j1),j1=1,3)

do 446 n=1,ne

446 read(13,*)k,(nop(n,i),i=1,3)

call NBND(nop,ne,nband,ndfg,ncn)

write(1,*)' nband=',nband

write(1,103)(j,(nop(j,i),i=1,ncn),j=1,ne)

103 format(2(i8,3i5))

write(1,889)(j,(cord(j,i),i=1,3),j=1,np)

889 format(2(i6,3f8.2))

do 2 nj=1,ne

print *' element',nj

Продолжение приложения 3

call sti24(nj,NOP,NE,NCN,THIN,Em,PR,np,ndf,cord,estifm,

* y12,y23,y31,x21,x32,x13,rtb,str)

call MinerS(str,ro,thin,y12,y23,y31,x21,x32,x13,dmgl,rtb)

WRITE(8) ((ESTIFM(I2,J2),I2=1,18),J2=1,18)

2 WRITE(9) ((dmgl(I2,J2),I2=1,18),J2=1,18)

WRITE(10) NOP,CORD

end file(10)

close(13)

end file(1)

close(1)

close(8)

close(9)

stop

end

Программа s2_3.F

parameter(NDFG=6,NCNM=3,NSTM=18,NP=6,NE=4

* ,NDF=6,NB=2,NBAND=24,IMAX=36,NSZF=36,NCN=3)

* пофазное формирование глобальной матрицы жесткости

dimension nop(ne,ncnm),estifm(nstm,nstm),

* cord(np,3),rq(nszf),nbc(nb),nfix(nb),

* s(nszf),a(imax,nband),

* npr(nszf),npr1(nszf),ndis(nszf)

open(16,file='RES2_3')

open(7,file='NPR',form='unformatted')

open(1,file='AGL',form='unformatted')

open(9,file='NOP')

read(9)nop,cord

close(9)

open(13,file='DTE2',form='unformatted')

do 445 n=1,nb

445 read(13,*) nbc(n),nfix(n)

close(13)

WRITE(16,9)(NBC(I),NFIX(I),I=1,NB)

Продолжение приложения 3

9 FORMAT(2i8)

call FRMEG(imax,nband,nszf,ne,ndfg,ncnm,

*nop,estifm,nstm,a,nb,nfix,nbc,ndis,knul)

open(33,file='NDIS',form='unformatted')

write(33)ndis

close(33)

call CORR(a,imax,nband,nszf,nszf1,ndis,knul,ik)

call PROF(a,imax,nband,npr,nszf1)

write(16,*)' nszf=',nszf,' nszf1=',nszf1

write(16,*)' nband=',nband,' imax=',imax,'ik=',ik,'knul=',knul

write(7)ik,(npr(i),i=1,nszf1)

end file(7)

close(7)

end file(16)

close(16)

close(1)

close(2)

stop

end