C C---------------------------------------------------------ALMAC_S_TERM--- C C SUBRUTINA ALMAC_S_TERM C C ALMACENAR LA SOLUCION ELASTICA C SUBROUTINE ALMAC_S_TERM(NSOL,SOLS,NODS,NODSR,NODTU,PRBT,ITER) C IMPLICIT NONE C INCLUDE '_param.for' INCLUDE '_struct.for' C CCCCCCCCCCCCCCCCCCCCCCCCCCC C C VARIABLES EXTERNAS C CCCCCCCCCCCCCCCCCCCCCCCCCCC C C SOLS : Solidos del problema RECORD /S_solido/ SOLS(N_MAX_SOL) RECORD /S_problema/ PRBT C NODS,NODSR : Nodos del problema RECORD /S_nodo/ NODS(N_MAX_NOD) RECORD /S_nodo/ NODSR(N_MAX_NOD) RECORD /S_nodo/ NODTU(N_MAX_NOD) C CCCCCCCCCCCCCCCCCCCCCCCCCCC C C VARIABLES EXTERNAS C CCCCCCCCCCCCCCCCCCCCCCCCCCC C C NSOL : Indicadores INTEGER NSOL INTEGER ITER C C CCCCCCCCCCCCCCCCCCCCCCCCCCC C C VARIABLES INTERNAS C CCCCCCCCCCCCCCCCCCCCCCCCCCC C C CSOL,CNOD,NNOD,I,J : Indicadores INTEGER CSOL,CNOD,NNOD,I,J,N C C 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 C WRITE(F_L,*) '----------------' C WRITE(F_L,*) 'ITERACION ',ITER C WRITE(F_L,*) '----------------' C DO csol = 1, nsol C i=1 C nnod=SOLS(csol).NNOD DO WHILE(i.LE.nnod) cnod=SOLS(csol).NOD(i) NODSR(cnod).NUM=cnod DO j = 1, 6 NODSR(cnod).SOLU(j)=NODS(cnod).SOLU(j) END DO c Inicializamos NOS(cnod).SOLU(j) c DO j = 1, 6 c NODS(cnod).SOLU(j)=0.d0 c END DO i=i+1 END DO C END DO C C DO i = 1 , PRBT.NSOL csol=SOLS(i).NUM C WRITE(F_L,*) C WRITE(F_L,10) ' SOLIDO ',csol C WRITE(F_L,*) ' **************************' C WRITE(F_L,*) C WRITE(F_L,16) 'Nodo ','PosX','PosY','T','Qa','Qb' C @ ,'DesplazamientoX' C @ ,'DesplazamientoY','t1a','t2a','t1b','t2b' C @ ,'D1','D2','D3' C WRITE(F_L,*) DO j = 1 , SOLS(csol).NNOD cnod=SOLS(CSOL).NOD(j) C WRITE(F_L,15) cnod,(NODS(cnod).XYZO(n),n=1,2) C @ ,(NODTU(cnod).SOLU(n),n=1,3),(NODSr(cnod).SOLU(n),n=1,6) C @ ,(NODS(cnod).SOLU(n),n=8,10) 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 C e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.


ANÁLISIS DEL PROBLEMA TERMOELÁSTICO DE CONTACTO MEDIANTE EL MÉTODO DE LOS ELEMENTOS DE CONTORNO CON MALLAS CONFORMES Y NO CONFORMES

: MORALES SIERRA, DIÓGENES
: Ingeniería Industrial
Contenido del proyecto: