Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Hi, i have this open source written in Frotran language
- this program is a 3d tree-dimensional nonlinear finite element of reinforced concrete
- this program has:
- different elements types (2/3/6/8/20/36 noded elements link and brick elements)
- nonlinear propertis of concrete/steel
- capability of cracking/crushing of concrete
- many post-cracking models
- .... and other capabilities
- please if you want help contact me at
- be.abou.raid@gmail.com
- the full list of the fortran file can be downloaded here:
- https://www.up-4ever.org/8pj5pusp5b4o
- some information and theoris
- https://tamtolearn.blogspot.com/
- ------------- here small portion of the program----
- PROGRAM P3DNFEA
- IMPLICIT REAL*8 (A-H,O-Z)
- C DIM1=====TOTDOF,DIM2====DOFEL,DIM3====NQP
- C DIM4=====TOTNOD,DIM5====TOTGAS,DIM6====TOTGASC
- C DIM7=====TOTELS,DIM8====NBAR,DIM9====HBAND
- C NOTE
- C 1. THE MAXIMUM DIMENSIONS OF ARRAYS RELATED TO ELEMENT
- C FORMULATION ARE SET ACCORDING TO THE REQUIREMENTS OF
- C THE 20 NODED BRICK ELEMENT. WHEN THE 32 NODED BRICK
- C IS USED THE MAXIMUM DIMENSIONS OF THESE ARRAYS SHOULD
- C BE MODIFIED ACCORDING TO NUMBER OF NODES PER ELEMENT.
- C 2. A PARABOLIC UNIAXIAL STRESS-STRAIN RELATION IS USED IN
- C THIS PROGRAM WHICH IS SUITABLE FOR THE ANALYSIS OF R.C
- C MEMBERS. FOR THE ELASTO-PLASTIC ANALYSIS OF HOMOGENEOUS
- C MEMBERS, SUITABLE RELATIONS ARE GIVEN IN A PROGRAM NAMED
- C **STEEL**
- CHARACTER *72 UNITS,TITLE
- INTEGER DIM1,DIM2,DIM3,DIM4,DIM5,DIM6,DIM7,DIM8,DIM9
- cc changes: dim1=3*totnod ccccccccccccccccccccccccc
- PARAMETER(DIM1=153,DIM2=60,DIM3=27,DIM4=51
- *,dim5=204,DIM6=108,DIM7=4,dim8=48,DIM9=120)
- INTEGER DIF, DIMEN, DOFNOD, ELNUM, ELTOP, ELTYP,
- * HBAND,RESNOD,RESTR,TOTGAS, BARNUM,ELTOPB,
- * TOTDOF, TOTELS, TOTNOD,TOTGASC
- REAL *8 LOADS,NU,INLOAD,COPRK,frac
- LOGICAL FIRST
- DIMENSION ERR(DIM1),
- * WGHTR(2),DE(3,3),
- * NWORK(20),RESTR(DIM4,4),
- * WGHT(DIM3),WORK(3),FACT(3),
- * ORLOAD(DIM1),PTDISP(DIM1),DLOADS(DIM1)
- COMMON sysk(dim1,dim9),TDISP(DIM1), INLOAD(DIM1),
- * COORD(DIM4,3),D(6,6),CURY(DIM6),
- * LOADS(DIM1),FTFT(6),FCFC(6),
- * DDC(DIM6,6,6),DCOSN(DIM6,3,3),DAMM(DIM6,2),
- * RLNTH(DIM8,3) ,STRNG(6,DIM5),HH(6),DDN(DIM6,3,3),
- * FYFY(6),EE(6), AA(6), POS(6),H1H1(DIM8,3),
- * STRSG(6,DIM5),PLOAD(DIM1),EFFST(DIM5),ECRK(3,DIM6),
- * DDET(DIM7,DIM3),
- * EPSTN(DIM5),
- * DCOSS(DIM6,6,6),
- * TRSTRN(6,DIM6),TPSTRN(3,DIM6),
- * ABSLNC(3,DIM6),SCR(3,DIM6),
- * TRSTRS(6,DIM6),TPSTRS(3,DIM6),
- * INDCT(3,DIM5),ELTOP(DIM7,22),
- * MPROP(DIM7+DIM8),ELTOPB(DIM8,2),NOUTP(2),NF(DIM4,3)
- DATA ID /6/,JD/6/,JELTOP/22/,JNF/3/,
- * IRESTR/DIM4/,IELTOP /DIM7/,JRESTR/4/,
- * ILOADS /DIM1/, INF /DIM4/,IGAS/DIM5/,
- * ISYSK/DIM1/,JSYSK/DIM1/,IGASC/DIM6/
- Cc for pure shear (nin=18) (shear+tension (nin=19) open (tensinsh) ccccccc
- DATA NIN /18/, NOUT /11/, NTEMP /13/,NOUT1/12/,NTEMP1/14/
- *,NSYSK/10/,lout/19/
- C
- C SET ITEST FOR FULL CHECKING
- C
- ccc ITEST = 0
- open(lout,file='table',status='replace')
- rewind(lout)
- OPEN(NOUT,FILE='SHEARDN1.OUT',status='replace')
- REWIND(NOUT)
- OPEN(NOUT1,FILE='TAPE11')
- REWIND(NOUT1)
- c ********************************************************************
- OPEN(NTEMP,FILE='WARM1',FORM='UNFORMATTED',status='replace')
- c *****************************************************************
- c OPEN(NTEMP,FILE='WARM1',form='unformated',status='UNKNOWN')
- REWIND(NTEMP)
- C *************************************
- C OPEN(NTEMP1,FILE='WARD1',FORM='UNFORMATTED',status='replace')
- C **************************************************
- OPEN(NTEMP1,FILE='WARD1',form='unformatted',status='replace')
- c OPEN(NTEMP1,FILE='WARD1')
- C ,STATUS='NEW',ACCESS='DIRECT',FORM='UNFORMATTED')
- REWIND(NTEMP1)
- open(nsysk,file='sysk1',form='unformatted')
- c OPEN(NSYSK,FILE='SYSK1.MAT',FORM='UNFORMATTED',
- c * ACCESS='DIRECT',RECL=8)
- REWIND(NSYSK)
- C* **********************
- C* * *
- C* * INPUT DATA SECTION *
- C* * *
- C* **********************
- OPEN(NIN,FILE='ph1.dat',STATUS='OLD')
- REWIND(NIN)
- READ(NIN,*)COPRK
- CCC CALL GETTIM(IH,IM,IS,IT)
- ccc WRITE(NOUT,60000) 'TIME = ',IH,IM,IS,IT
- 60000 FORMAT (1X,A7,I2,':',I2,':',I2,':',I2)
- READ(NIN,*)TITLE
- WRITE(NOUT,*)TITLE
- WRITE(LOUT,*)TITLE
- WRITE(*,*)TITLE
- READ(NIN,*)UNITS
- WRITE(NOUT,*)UNITS
- WRITE(*,*)UNITS
- C
- C INPUT OF NODAL GEOMETRY
- C
- WRITE (NOUT,9010)
- READ (NIN,8010) TOTNOD, DIMEN
- WRITE (NOUT,9020) TOTNOD, DIMEN
- DO 1010 I=1,TOTNOD
- READ (NIN,8020) NODNUM, (WORK(J),J=1,DIMEN)
- WRITE (NOUT,9030) NODNUM, (WORK(J),J=1,DIMEN)
- DO 1020 J=1,DIMEN
- 1020 COORD(NODNUM,J) = WORK(J)
- 1010 CONTINUE
- C
- C INPUT OF ELEMENT TOPOLOGY
- C
- WRITE (NOUT,9040)
- C
- C REINFORCEMENT ELEMENTS ARE TYPES 2,3&4
- C
- READ (NIN,8010) ELTYP, TOTELS, NODEL
- WRITE (NOUT,9020) ELTYP,TOTELS, NODEL
- DO 1040 I=1,TOTELS
- READ (NIN,8040) ELNUM, (NWORK(J),J=1,NODEL), MPROP(ELNUM)
- WRITE (NOUT,9090) ELNUM, (NWORK(J),J=1,NODEL), MPROP(ELNUM)
- ELTOP(ELNUM,1) = ELTYP
- ELTOP(ELNUM,2) = NODEL
- DO 1050 J=1,NODEL
- 1050 ELTOP(ELNUM,J+2) = NWORK(J)
- 1040 CONTINUE
- READ(NIN,8010)NBAR
- WRITE(NOUT,9020)NBAR
- DO 1030 I=1,NBAR
- READ(NIN,8025)BARNUM,(NWORK(J),J=1,2),MPROP(BARNUM+TOTELS),
- * (WORK(K),K=1,DIMEN)
- WRITE(NOUT,9035)BARNUM,(NWORK(J),J=1,2),MPROP(BARNUM+TOTELS),
- * (WORK(K),K=1,DIMEN)
- DO 1900 J=1,DIMEN
- IF(J.LT.3)ELTOPB(BARNUM,J)=NWORK(J)
- 1900 RLNTH(BARNUM,J)=WORK(J)
- 1030 CONTINUE
- C
- C INPUT OF MATERIAL PROPERTIES AND
- C CONSTRUCTION OF ELASTIC STRESS-STRAIN MATRIX D
- C
- WRITE (NOUT,9050)
- READ(NIN,8060) ALPH,CHI,FRAC
- WRITE(NOUT,9170) ALPH,CHI,FRAC
- WRITE(lout,9170) ALPH,CHI,FRAC
- READ(NIN,8060) GAMA,SHI,SHIM
- WRITE(NOUT,9170) GAMA,SHI,SHIM
- WRITE(lout,9170) GAMA,SHI,SHIM
- READ(NIN,8010)NMATS
- WRITE(NOUT,9020)NMATS
- DO 1060 IMAT = 1,NMATS
- READ (NIN,8060) NU, E, FY , A , H ,FT, FC
- POS(IMAT) = NU
- EE(IMAT) = E
- FYFY(IMAT) = FY
- AA(IMAT) = A
- HH(IMAT)=H
- FTFT(IMAT)=FT
- FCFC(IMAT)=FC
- 1060 WRITE (NOUT,9170) NU, E, FY , A , H , FT ,FC
- WRITE (lout,9170) NU, E, FY , A , H , FT ,FC
- NU=POS(1)
- E=EE(1)
- itest=0
- CALL DISO(D, ID, JD, E, NU, NUMSS, ITEST)
- CALL MATNUL(DE,3,3,3,3,ITEST)
- FACTOR=E/(1.0-NU*NU)
- DE(1,1)=FACTOR
- DE(2,2)=DE(1,1)
- DE(1,2)=NU*FACTOR
- DE(2,1)=DE(1,2)
- DE(3,3)=E/(2.0+2.0*NU)
- C
- C INPUT OF NUMBER OF DEGREES OF FREEDOM
- C PER NODE, INPUT OF RESTRAINED NODE
- C DATA AND CONSTRUCTION OF NODAL FREEDOM
- C ARRAY NF
- C
- write(lout,*)
- write(lout,*)'incrms thetaP (deg) AvGama xy(%) Sh Strs(GPa)
- *ExSh Strs(GPa) ExAx Strs(GPa)'
- 1979 format(i5,5f13.6)
- WRITE (NOUT,9070)
- READ (NIN,8010) DOFNOD,RESNOD
- WRITE (NOUT,9020) DOFNOD,RESNOD
- K = DOFNOD + 1
- DO 1070 I=1,RESNOD
- READ (NIN,8010) (RESTR(I,J),J=1,K)
- 1070 WRITE (NOUT,9020) (RESTR(I,J),J=1,K)
- CALL FORMNF(RESTR, IRESTR, JRESTR, RESNOD, TOTNOD, dofnod,
- * NF, INF, JNF, TOTDOF)
- C
- C LOADING DATA INPUT
- C
- WRITE (NOUT,9080)
- READ (NIN, 8010)LODCAS
- WRITE(NOUT,9020)LODCAS
- READ (NIN,8030)(FACT(I), I=1,LODCAS)
- WRITE (NOUT,9060)(FACT(I),I=1,LODCAS)
- CALL VECNUL(ORLOAD,ILOADS,TOTDOF,ITEST)
- READ (NIN, 8010)LODNOD
- WRITE(NOUT,9020)LODNOD
- DO 1080 I=1,LODNOD
- READ (NIN,8020) NODNUM, (WORK(J),J=1,DOFNOD)
- WRITE (NOUT,9030) NODNUM, (WORK(J),J=1,DOFNOD)
- DO 1090 J=1,DOFNOD
- K = NF(NODNUM,J)
- 1090 IF (K.NE.0) ORLOAD(K)=WORK(J)
- 1080 CONTINUE
- READ (NIN, 8010)NGAUS,NCRIT,TOTGAS,TOTGASC
- WRITE(NOUT,9020)NGAUS,NCRIT,TOTGAS,TOTGASC
- READ(NIN,8010)NALGO,KSERCH,NINCS,NITER
- WRITE(NOUT,9020)NALGO,KSERCH,NINCS,NITER
- READ(NIN,8050)TOLER,ICONVER,NOUTP(1),NOUTP(2)
- 8050 FORMAT(F12.0,3I5)
- WRITE(NOUT,9140)TOLER,ICONVER,NOUTP(1),NOUTP(2)
- 9140 FORMAT(F7.2,3I5)
- C9140 FORMAT(' CONVERGENCE TOLERANCE = ',F7.2,
- C * ' CONVERG. CRITERION= ',I2,
- C * ' INITIAL OUTPUT PARAMETER= ',I2,
- C * ' FINAL OUTPUT PARAMETER= ',I2)
- C
- C CALCULATION OF SEMI-BANDWIDTH
- C
- DIF=0
- DO 1100 NELE=1,TOTELS
- CALL FREDIF(NELE, ELTOP, IELTOP, JELTOP, NF, INF, JNF,
- * DOFNOD, nodel, DIF)
- 1100 CONTINUE
- HBAND = DIF + 1
- WRITE (*,*) 'HBAND =',HBAND
- WRITE(NOUT,9020)HBAND,TOTDOF
- CALL VECNUL(LOADS,ILOADS,TOTDOF,ITEST)
- CALL VECNUL(INLOAD,ILOADS,TOTDOF,ITEST)
- CALL VECNUL(TDISP,ILOADS,TOTDOF,ITEST)
- CALL VECNUL(EPSTN,IGAS,TOTGAS,ITEST)
- CALL MATNUL(STRSG,NUMSS,IGAS,NUMSS,TOTGAS,ITEST)
- CALL VECNUL(EFFST,IGAS,TOTGAS,ITEST)
- CALL VECNUL(CURY,IGASC,TOTGASC,ITEST)
- DO 111 J=1,TOTGAS
- DO 111 I=1,3
- 111 INDCT(I,J)=0
- CALL MATNUL(STRNG,NUMSS,IGAS,NUMSS,TOTGAS,ITEST)
- CALL MATNUL(ABSLNC,3,IGASC,3,TOTGASC,ITEST)
- CALL MATNUL(SCR,3,IGASC,3,TOTGASC,ITEST)
- CALL MATNUL(ECRK,3,IGASC,3,TOTGASC,ITEST)
- CALL MATNUL(TRSTRS,NUMSS,IGASC,NUMSS,TOTGASC,ITEST)
- CALL MATNUL(TRSTRN,NUMSS,IGASC,NUMSS,TOTGASC,ITEST)
- CALL MATNUL(TPSTRN,3,IGASC,3,TOTGASC,ITEST)
- CALL MATNUL(TPSTRS,3,IGASC,3,TOTGASC,ITEST)
- DO 1760 K=1,TOTGASC
- DAMM(K,1)=1.0
- DAMM(K,2)=1.0
- DO 1760 I=1,NUMSS
- DO 1760 J=1,NUMSS
- DDC(K,I,J)=0.0
- IF(I.LE.3.AND.J.LE.3)THEN
- DCOSN(K,I,J)=0.0
- DDN(K,I,J)=0.0
- ENDIF
- 1760 DCOSS(K,I,J)=0.0
- CALL VECNUL(ERR,ILOADS,TOTDOF,ITEST)
- C
- C LOOP OVER EACH INCREMENT
- C
- KRESL=2
- MOUT=0
- open (mout,file='panel.OUT',status='REPLACE')
- write (mout,*)title
- write (mout,*)alph,chi
- write (mout,*)'E=',e,'Fy=',fy, 'ft=',ft,'fc=',fc
- write(mout,*)
- DO 1110 IINCS=1,NINCS
- WRITE(*,70000)'IINCS = ',IINCS
- 70000 FORMAT (1X,A9,1X,I3)
- write(mout,*)'increm.=',iincs
- WRITE(NOUT1,*)'****IINCS = ',IINCS
- WRITE(NOUT,9100)IINCS
- 9100 FORMAT(1H /17H NO OF INCREMENT=,I3)
- DO 1120 K=1,TOTDOF
- IF(IINCS.LT.11)THEN
- LOADS(K)=FACT(1)*ORLOAD(K)
- INLOAD(K)=INLOAD(K)+LOADS(K)
- ELSEIF(IINCS.LT.50)THEN
- LOADS(K)= FACT(2)*ORLOAD(K)
- INLOAD(K)=INLOAD(K)+LOADS(K)
- ENDIF
- 1120 CONTINUE
- C WRITE (*,70001) 'INLOAD = ',INLOAD(1)
- C70001 FORMAT (1X,A9,1X,E10.3)
- C
- C LOOP OVER EACH ITERATION
- C
- C MSERCH IS COUNTER OF LINE SEARCH WITNIN LOAD INCREMENT
- MSERCH=0
- DO 1130 IITER = 1,NITER
- cc write(*,*)'iiter=== niter===',iiter,niter
- IF(KRESL.NE.1)THEN
- KRESL=2
- IF(NALGO.EQ.1.AND.IINCS.EQ.1.AND.IITER.EQ.1)KRESL=1
- IF(NALGO.EQ.2)KRESL=1
- IF(NALGO.EQ.3.AND.IITER.EQ.1)KRESL=1
- IF(NALGO.EQ.4.AND.IINCS.EQ.1.AND.IITER.EQ.1)KRESL=1
- IF(NALGO.EQ.4.AND.IITER.EQ.2)KRESL=1
- CCC IITER=== 1 , 6 , 11
- IF(NALGO.EQ.5)THEN
- I=(IITER-1)/5
- K=IITER-1
- IF(K.EQ.5*I)KRESL=1
- ENDIF
- CCC IITER=== 2 , 7 , 12
- IF(NALGO.EQ.6.AND.IINCS.EQ.1.AND.IITER.EQ.1)KRESL=1
- IF(NALGO.EQ.6.AND.IITER.GE.2)THEN
- I=(IITER-2)/10
- K=IITER-2
- IF(K.EQ.10*I)KRESL=1
- ENDIF
- ENDIF
- C CHECK WETHER THE EVALUATION OF STIFNESS MATRIX
- C IS REQUIRED
- C
- IF(KRESL.EQ.1)THEN
- C* *********************
- C* * *
- C* * STIFFNESS MATRIX *
- C* * *
- C* *********************
- ccc write(*,*)'cal stiffk 302 totels=',totels
- ccccccccccccc add sysk cccccccccccccccccc
- CALL STIFFK(DOFNOD,NGAUS,TOTELS,
- * DIMEN,IINCS,NUMSS,IITER,
- * NCRIT,NTEMP,NSYSK,WGHT,
- * HBAND,TOTDOF,NQP,NBAR,WGHTR,NQPR)
- cc write(*,*)'dofnod=',dofnod,'ngaus=',ngaus,'totels=',totels,'dim
- cc *en=',dimen,'iincs=',iincs,'numss=',numss,'iiter=',iiter,'ncrit=',
- cc *ncrit,'ntemp=',ntemp,'nsysk=',nsysk,'wght=',wght,'hband=',hband,
- cc *'totdof=',totdof,'nqp=',nqp,'nbar=',nbar,'wightr=',wightr,
- cc *'nqpr=',nqpr
- C* *********************
- C* * *
- C* * EQUATION SOLUTION *
- C* * *
- C* *********************
- C
- C SOLUTION OF SYSTEM MATRIX FOR THE
- C NODAL DISPLACEMENTS
- C
- IRES=0
- ccc write(*,*)'isysk=',isysk,'jsysk=',jsysk,'totdof=',totdof,'hband'
- cc *,hband,'band=',band,'loads=',loads,'ires=',ires,'nsysk=',nsysk
- CALL SOLVE (ISYSK, JSYSK, TOTDOF, HBAND,sysk,LOADS,IRES)
- KRESL=2
- ELSE
- KRESL=2
- IRES=1
- CALL SOLVE (ISYSK, JSYSK, TOTDOF, HBAND,sysk,lOADS,IRES)
- ccc write(*,*)'solve=====solve=reuls========solve=========='
- ENDIF
- C NSERCH ISCOUNTER OF LINE SEARCH WITHIN EACH ITERATION
- NSERCH=0
- ETRO=1.000
- DO 1780 K=1,TOTDOF
- PTDISP(K)=TDISP(K)
- 1780 DLOADS(K)=LOADS(K)
- C BEFORE USING LINE SERCH ACTIVATE THE FOLLOWING STATEMENT
- IF (KSERCH.EQ.0)GO TO 7390
- c if (kserch.eq.1)go to 7390
- C PREVIOUS STATE OF STRUCTURE
- ccc write(*,*)'===========ntemp1==========',ntemp1
- ccc write(*,*)'ntemp1=',ntemp1
- REWIND(NTEMP1)
- DO 1790 K=1,TOTGAS
- WRITE(NTEMP1)EPSTN(K),EFFST(K)
- WRITE(NTEMP1)(INDCT(I,K),I=1,3)
- WRITE(NTEMP1)(STRSG(I,K),I=1,6)
- WRITE(NTEMP1)(STRNG(I,K),I=1,6)
- ccc write(*,*)'epstn(',k,')=',epstn(k)
- 1790 CONTINUE
- DO 1972 K=1,TOTGASC
- WRITE(NTEMP1)CURY(K),(DAMM(K,I),I=1,2)
- WRITE(NTEMP1)(ECRK(I,K),I=1,3)
- WRITE(NTEMP1)(SCR(I,K),I=1,3)
- WRITE(NTEMP1)(ABSLNC(I,K),I=1,3)
- WRITE(NTEMP1)(TPSTRN(I,K),I=1,3)
- WRITE(NTEMP1)(TPSTRS(I,K),I=1,3)
- WRITE(NTEMP1)(TRSTRN(I,K),I=1,6)
- WRITE(NTEMP1)(TRSTRS(I,K),I=1,6)
- DO 1972 J=1,6
- IF(J.LE.3)THEN
- WRITE(NTEMP1)(DDN(K,J,I),I=1,3)
- WRITE(NTEMP1)(DCOSN(K,J,I),I=1,3)
- ENDIF
- WRITE(NTEMP1)(DDC(K,J,I),I=1,6)
- WRITE(NTEMP1)(DCOSS(K,J,I),I=1,6)
- 1972 CONTINUE
- REWIND(NTEMP1)
- SO=0.0
- DO 1800 K=1,TOTDOF
- 1800 SO=SO+ERR(K)*DLOADS(K)
- 7380 CONTINUE
- IF(ETRO.NE.1.0)THEN
- NSERCH=NSERCH+1
- MSERCH=MSERCH+1
- WRITE(NOUT,*)'MSERCH=',MSERCH,' NSERCH=',NSERCH,' IITER=',IITER
- ENDIF
- DO 1810 K=1,TOTDOF
- 1810 LOADS(K)=ETRO*DLOADS(K)
- 7390 CONTINUE
- C
- CALL VECNUL(ERR,ILOADS,TOTDOF,ITEST)
- C* *********************
- C* * *
- C* * RESIDUAL FORCES *
- C* * *
- C* *********************
- ccc write(*,*)'cal residu 387 totels=',totels
- REWIND(NTEMP)
- CALL RESIDU(NCRIT,NQP,
- * NUMSS,DIMEN,TOTELS,DOFNOD,
- * WGHT,IITER,KRESL,IINCS,
- * TOTDOF,NOUT,NTEMP,DE,
- * ALPH,GAMA,CHI,SHI,FRAC,SHIM,
- * NBAR,WGHTR,NQPR,coprk)
- DO 1820 K=1,TOTDOF
- 1820 ERR(K)=INLOAD(K)-PLOAD(K)
- ccc write(*,*)'iincs=',iincs,'iiter=',iiter,'nserch=',nserch
- IF(KSERCH.EQ.0)GOTO 7370
- c if(kserch.eq.1)goto 7370
- PETRO=ETRO
- IF(IINCS.EQ.1)GOTO 7370
- IF(IITER.EQ.1)GOTO 7370
- IF(NSERCH.EQ.1) GOTO 7370
- C IF(NSERCH.EQ.3) GOTO 7370
- C* *********************
- C* * *
- C* * LINE SEARCHES *
- C* * *
- C* *********************
- SN=0.0
- DO 1830 K=1,TOTDOF
- 1830 SN=SN+ERR(K)*DLOADS(K)
- ccc write(*,*)'nout====',nout
- ccc write(*,*)'so=',so,'sn=',sn,'iiter=',iiter
- ccc WRITE(*,70002)'SO =',SO,'SN =',SN,'IITER =',IITER
- ccc70002 FORMAT (1X,A4,1X,E10.3,5X,A4,1X,E10.3,5X,A7,1X,I3)
- WRITE(NOUT,*)'SO=',SO,' SN=',SN,' IITER=',IITER
- cccc write(nout,*)so,sn,iiter
- ABSSO=ABS(SO)
- IF(ABSSO.EQ.0.0)ABSSO=0.0000001
- SRTO=ABS(SN)/ABSSO
- IF(SRTO.LT.0.6)GOTO 7370
- SOSN=SO-SN
- IF(SOSN.EQ.0.0)SOSN=0.0000001
- ETRO=SO*PETRO/SOSN
- ccc WRITE(*,*)'SRTO=',SRTO,' ETRO=',ETRO
- WRITE(NOUT,*)'SRTO=',SRTO,' ETRO=',ETRO
- IF(ETRO.GT.10.0)ETRO=10.0
- IF(ETRO.LT.0.1)ETRO=0.1
- C RECALL THE PREVIOUS STATE OF STRUCTURE
- REWIND(NTEMP1)
- DO 1840 K=1,TOTGAS
- READ(NTEMP1)EPSTN(K),EFFST(K)
- READ(NTEMP1)(INDCT(I,K),I=1,3)
- READ(NTEMP1)(STRSG(I,K),I=1,6)
- READ(NTEMP1)(STRNG(I,K),I=1,6)
- 1840 CONTINUE
- DO 1842 K=1,TOTGASC
- READ(NTEMP1)CURY(K),(DAMM(K,I),I=1,2)
- READ(NTEMP1)(ECRK(I,K),I=1,3)
- READ(NTEMP1)(SCR(I,K),I=1,3)
- READ(NTEMP1)(ABSLNC(I,K),I=1,3)
- READ(NTEMP1)(TPSTRN(I,K),I=1,3)
- READ(NTEMP1)(TPSTRS(I,K),I=1,3)
- READ(NTEMP1)(TRSTRN(I,K),I=1,6)
- READ(NTEMP1)(TRSTRS(I,K),I=1,6)
- DO 1842 J=1,6
- IF(J.LE.3)THEN
- READ(NTEMP1)(DDN(K,J,I),I=1,3)
- READ(NTEMP1)(DCOSN(K,J,I),I=1,3)
- ENDIF
- READ(NTEMP1)(DDC(K,J,I),I=1,6)
- READ(NTEMP1)(DCOSS(K,J,I),I=1,6)
- 1842 CONTINUE
- REWIND(NTEMP1)
- GOTO 7380
- 7370 CONTINUE
- DO 1450 K=1,TOTDOF
- 1450 TDISP(K)=TDISP(K)+LOADS(K)
- ------------------ please find the full listing in the link above----
- ****************
- fiber concrete finite element method
- concrete structures finite element method
- structural concrete finite-element analysis for limit-state design
- finite element for reinforced concrete
- concrete finite element analysis
- reinforced concrete finite element analysis
- concrete centre finite element analysis
- prestressed concrete finite element analysis
- finite element design of concrete structures
- finite element design of concrete structures rombach pdf
- finite element analysis of concrete structures
- finite element design of concrete structures pdf
- finite element analysis of concrete gravity dam
- finite element design of concrete structures rombach
- finite element modeling of concrete expansion and confinement
- finite element analysis of concrete box culverts
- finite element analysis of concrete fracture
- concrete slabs design with finite element methods
- simulation of concrete batch plant production
- simulation of concrete fracture using abaqus
- simulation of concrete structures
- simulation of concrete beam
- simulation of concrete gravity dams
- simulation on concrete
- about finite element analysis
- about finite element method
- about finite element
- about finite element analysis pdf
- finite element methods 1970’s and beyond
- finite element by mit
- finite element in analysis and design
- finite element in matlab
- finite element in python
- finite element in structural analysis
- finite element for fluid mechanics
- finite element for dummies
- finite element for structural analysis
- finite element for analysis
- finite element for engineer
- finite element in plasticity
- finite element in fluid mechanics
- finite element in geotechnical engineering
- finite element in geotechnical engineering pdf
- finite element in civil engineering
- finite element in engineering pdf
- finite element like method
- finite element of beam
- finite of element
- finite element analysis was carried out
- finite element post processing
- finite element post processing software
- postdoc in finite element
- finite element analysis post processing
- finite element analysis post-tensioned
- finite element analysis post
- finite element post
- finite element endodontic post
- finite element preprocessor
- finite element analysis pro and cons
- finite element analysis staad pro
- finite element in staad pro
- finite element analysis using staad pro
- finite element method in staad pro
- finite element analysis materials science
- finite element method sciencedirect
- finite element method material science
- finite element method science journal
- computational materials science finite element analysis
- computational materials science finite element method
- through finite element method
- finite element method under pressure
- finite element startup
- finite element versus finite difference
- finite element difference
- finite element difference method
- finite element versus finite volume
- finite difference element volume methods
- finite element analysis concepts via solidworks
- finite element vs finite difference
- finite element vs finite volume
- finite element vs finite difference vs finite volume
- finite element vs finite difference methods
- finite element vs discrete element
- finite element triangle vs quad
- finite element with matlab
- finite element with python
- finite element with matlab pdf
- finite element with mathematica
- finite element with c++
- finite element with julia
- finite element with
- concrete analysis of tcp and udp
- concrete analysis of measurement
- concrete analysis of
- concrete analysis in abaqus
- concrete buildings analysis for safe construction
- analysis in concrete beam
- concrete petrographic analysis near me
- concrete analysis of concrete conditions
- concrete analysis and testing marathon fl
- concrete analysis and design software
- concrete analysis of the concrete situation
- concrete beams with openings analysis and design
- concrete modeling in abaqus
- concrete modeling in ansys
- concrete model in abaqus
- concrete for modelling
- reinforced concrete modeling in abaqus
- fiber concrete modeling in abaqus
- modeling concrete slab in revit
- reinforced concrete beam modeling in abaqus
- concrete modeling and simulation
- modeling of concrete performance
- modeling of concrete
- modeling of concrete flow
- modeling of concrete shear wall
- modeling of concrete beam
- modeling of concrete frames
- modeling on concrete
- modeling concrete strength
Add Comment
Please, Sign In to add comment