C
C-------------------------------------------------------------RESOL_T_TERM---
C
C SUBRUTINA RESOL_T_TERM
C
SUBROUTINE RESOL_T_TERM(PRBT,SOLT,MATT,CNTT,ELTT,NODT,
@ PNTT,ZCTOT,NPGT,SAijT,SFiT,SXjT,PCT,PTT)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES DEL PROBLEMA
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PRBT : Variables globales del problema
RECORD /S_problema/ PRBT
C SOLT : Solidos del problema
RECORD /S_solido/ SOLT(N_MAX_SOL)
C MATT : Materiales
RECORD /S_material/ MATT(N_MAX_SOL)
C CNTT : Contornos del problema
RECORD /S_contorno/ CNTT(N_MAX_CNT)
C ELTT : Elementos del problema
RECORD /S_elemento/ ELTT(N_MAX_ELT)
C ELTPT : Tipos de elemento
RECORD /S_tipelt/ ELTPT(N_MAX_TPEL)
C NODT : Nodos del problema
RECORD /S_nodo/ NODT(N_MAX_NOD)
C PNT : Puntos internos del problema
RECORD /S_pntint/ PNTT(N_MAX_PINT)
C PCT : Puntos de colocacion
RECORD /S_coloc/ PCT(N_MAX_PC)
C ZCTOT : Zonas de contacto del problema
RECORD /S_zncto/ ZCTOT(N_MAX_ZCTO)
C SAIJT(,) : Matriz del sistema.
REAL*8 SAijT(N_MAX_INC,N_MAX_INC)
C SXJ() : Vector de incógnitas.
REAL*8 SXjT(N_MAX_INC)
C SFIT() : Termino independiente.
REAL*8 SFiT(N_MAX_INC)
C NPGT : Numero de puntos de Gauss
INTEGER NPGT
C PTT : Paso de la iteración térmica
INTEGER PTT
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES INTERNAS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J : Contadores
INTEGER I,J
C
CCCCCCCCCCCCCCC
C
C RUTINA
C
CCCCCCCCCCCCCCC
C
C
C ANALIZAMOS SI ES EL PRIMER PASO
C
C
C CHEQUEO DE CONTINUIDAD Y CALCULO DE COORD. LOCALES
C
CALL CONECTIVIDAD(PRBT,SOLT,CNTT,PRBT.ELTP,ELTT,NODT)
C
C MONTAJE DE LAS CONDICIONES DE CONTORNO (incluye el contacto)
C
CALL MONTAJE_CDC_POT(PRBT,SOLT,ZCTOT,CNTT,PRBT.ELTP,
@ ELTT,NODT)
C
C La esctitura de los datos de entrada sólo la hacemos al principio.
C
SELECT CASE(PTT)
C
CASE(0)
C
C ESCRITURA DE DATOS DE ENTRADA EN FICHERO DE SALIDA
C
CALL ESC_ENT_POT(PRBT,SOLT,MATT,CNTT,ELTT,NODT,PNTT)
C
CASE(1)
C
CONTINUE
C
CASE DEFAULT
WRITE(F_W,01000)
WRITE(*,01000)
END SELECT
C
C
C
C ECUACIONES DE APOYO TERMICO (Problema de Neumann)
C
CALL APOYO_POT(PRBT,SOLT,NODT,SAijT,SFiT)
C
C MONTAJE DE LOS PUNTOS DE COLOCACION
C
CALL MONTAJE_PC_POT(PRBT,SOLT,MATT,CNTT,PRBT.ELTP,
@ ELTT,NODT,PCT)
C
C MONTAJE DE LAS ECUACIONES INTEGRALES DE CONTORNO
C
CALL MONTAJE_BIE_POT(PRBT,SOLT,MATT,CNTT,ELTT,NODT,PCT,
@ SAijT,SFiT,NPGT)
C
C
C
C
C
C
C ELECCION DEL ALGORITMO DE SOLUCION
C
SELECT CASE(PRBT.TF.TIP)
C
CASE(11,21)
C
C SOLUCION DIRECTA
C
CALL ALG_DIRECTO_POT(PRBT,SOLT,CNTT,ELTT,NODT,ZCTOT,
@ SAijT,SFiT,SXjT)
C
CASE(12,22)
C
C ALGORITMO INCREMENTAL CON ESCALADO DE LA SOLUCION
C
CALL ALG_INCREMENTAL_POT(PRBT,SOLT,CNTT,ELTT,NODT,ZCTOT,
@ SAijT,SFiT,SXjT)
CASE DEFAULT
WRITE(F_W,00002) PRBT.TF.TIP
WRITE (*,00002) PRBT.TF.TIP
STOP
END SELECT
C
C Borrado del sistema de ecuaciones
C
DO i=1,PRBT.NAP+PRBT.NBIE+PRBT.NCE
DO j=1,PRBT.NINC
SAijT(i,j)=0.d0
END DO
SFiT(i)=0.d0
END DO
C
C
CCCCCCCCCCCCCCCCCCCCCCCCC
C FIN DE LA RUTINA
CCCCCCCCCCCCCCCCCCCCCCCCC
C
RETURN
C
C FORMATOS DE ESCRITURA
C
01000 FORMAT(/,2x,'ERROR 01000 - RESOL_T_TERM - Error Programa ')
C
00002 FORMAT(/,2x,'ERROR 00002 - RESOL_T_TERM - Tipo de formulacion (',
@ i2,') erronea.')
00003 FORMAT(/,2x,'ERROR 00003 - RESOL_T_TERM - Tipo de problema (',
@ i2,') erronea.')
00004 FORMAT(/,2x,'ERROR 00004 - RESOL_T_TERM - Tipo de problema (',
@ i2,') erronea.')
10 FORMAT(i2)
C
RETURN
C FINAL DEL PROGRAMA
C
END
C
C GRUPO DE ELASTICIDAD Y RESISTENCIA DE MATERIALES.
C UNIVERSIDAD DE SEVILLA
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
