C Rutina 080
C-------------------------------------------------------------ESC_DEF---
C
C SUBRUTINA ESC_ITR
C
C ESCRITURA DE LA DEFORMADA
C
SUBROUTINE ESC_ITR_POT(PRB,SOL,CNT,ELT,NOD,X)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PRB : Variables globales del problema
RECORD /S_problema/ PRB
C SOL : Solidos
RECORD /S_solido/ SOL(N_MAX_SOL)
C CNT : Contornos del problema
RECORD /S_contorno/ CNT(N_MAX_CNT)
C ELT : Elementos del problema
RECORD /S_elemento/ ELT(N_MAX_ELT)
C NOD : Nodos del problema
RECORD /S_nodo/ NOD(N_MAX_NOD)
C X : Factor de escala para el que se viola la primera cond. de contacto
REAL*8 X
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J,K,M,N : Contadores
INTEGER I,J,K,M,N
C CCNT,CELT,ANOD,CNOD : Contadores
INTEGER CCNT,CELT,ANOD,CNOD,PELT,PNOD
C CELTP : Contador
INTEGER CELTP, CSOL, NNOD
C r0,rx,r1 : Posicion inicial, intermedia y final de los nodos
REAL*8 r0(2),rx(2),r1(2)
C u0,ux,u1 : Desplazamiento inicial, intermedia y final de los nodos
REAL*8 u0(2),ux(2),u1(2)
C t0,tx,t1 : Tensiones inicial, intermedia y final de los nodos
REAL*8 t0(2),tx(2),t1(2)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
anod=0
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCC
C ESCRITURA DEFORMADA
CCCCCCCCCCCCCCCCCCCCCCCCCCCC
CCCCCCCCCCCCCCC
C
C RUTINA
C
CCCCCCCCCCCCCCC
C
i=1
C
C CALCULAMOS LA SUMA PARA TODOS LOS NODOS
C
C Bucle que recorra todos los solidos
C
WRITE(F_H,*) '----------------'
WRITE(F_H,*) 'NUEVA ITERACION '
WRITE(F_H,*) '----------------'
C
C
C
DO i = 1 , PRB.NSOL
csol=SOL(i).NUM
WRITE(F_H,*)
WRITE(F_H,10) ' SOLIDO ',csol
WRITE(F_H,*) ' **************************'
WRITE(F_H,*)
WRITE(F_H,16) 'Nodo ','PosX','PosY'
@ ,'DesplazamientoX'
@ ,'DesplazamientoY','t1a','t2a','t1b','t2b'
@ ,'D1','D2','D3'
WRITE(F_H,*)
DO j = 1 , SOL(csol).NNOD
cnod=SOL(CSOL).NOD(j)
WRITE(F_H,15) cnod,(NOD(cnod).XYZO(n),n=1,2)
@ ,(NOD(cnod).D_SOLU(n),n=1,7)
@ ,x
END DO
C FIN DEL BUCLE QUE RECORRE LOS SOLIDOS
END DO
10 FORMAT(2x,a,i2,a,a)
12 FORMAT(4x,i4,10(' - ',a3,': ',e22.16))
14 FORMAT(4x,i4,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16)
15 FORMAT(4x,i4,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16
@ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16
@ ,4x,e22.16,4x,e22.16,4x,e22.16)
16 FORMAT(5x,a,12x,a,21x,a,21x,a,21x,a,21x,a,21x,a,21x,a,21x,a,21x,a
@ ,21x,a,21x,a,21x,a,21x,a,21x,a)
17 FORMAT(19x,i4,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16)
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
