C
C-------------------------------------------------------------RESOL_S_TERM---
C
C SUBRUTINA RESOL_S_TERM
C
SUBROUTINE RESOL_S_TERM(PRB,PRBS,SOLS,MATS,CNTS,ELTS,NODS,PNTS,
@ ZCTOS,NPGS,SAijS,SFiS,SXjS,PCS,NODT,PRBT,ALFA,MATT,PS)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES DEL PROBLEMA
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PRB,PRBS,PRBT : Variables globales del problema
RECORD /S_problema/ PRB
RECORD /S_problema/ PRBS
RECORD /S_problema/ PRBT
C SOLS : Solidos del problema
RECORD /S_solido/ SOLS(N_MAX_SOL)
C MATS,MATT : Materiales
RECORD /S_material/ MATS(N_MAX_SOL)
RECORD /S_material/ MATT(N_MAX_SOL)
C CNTS : Contornos del problema
RECORD /S_contorno/ CNTS(N_MAX_CNT)
C ELTS : Elementos del problema
RECORD /S_elemento/ ELTS(N_MAX_ELT)
C NODS,NODT : Nodos del problema
RECORD /S_nodo/ NODS(N_MAX_NOD)
RECORD /S_nodo/ NODT(N_MAX_NOD)
C PNTS : Puntos internos del problema
RECORD /S_pntint/ PNTS(N_MAX_PINT)
C PCS : Puntos de colocacion
RECORD /S_coloc/ PCS(N_MAX_PC)
C ZCTOS : Zonas de contacto del problema
RECORD /S_zncto/ ZCTOS(N_MAX_ZCTO)
C SAIJS(,) : Matriz del sistema.
REAL*8 SAijS(N_MAX_INC,N_MAX_INC)
C SXJS() : Vector de incógnitas.
REAL*8 SXjS(N_MAX_INC)
C SFIS() : Termino independiente.
REAL*8 SFiS(N_MAX_INC)
C NPGS : Numero de puntos de Gauss
INTEGER NPGS
C ALFA : Coeficientes de dilatacion termica
RECORD /S_alfaTERM/ ALFA
C PS : Paso de la iteracion
INTEGER PS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES INTERNAS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J : Contadores
INTEGER I,J
C SVT : Ecuación correspondiente a los apoyos puntuales
REAL*8 SVT(3*N_MAX_SOL,N_MAX_INC)
C SW : Proyeccion de los MSR sobre las incognitas
REAL*8 SW(N_MAX_INC,3*N_MAX_SOL)
C
CCCCCCCCCCCCCCC
C
C RUTINA
C
CCCCCCCCCCCCCCC
C
C
C
C ANALIZAMOS SI ES EL PRIMER PASO
C
c SELECT CASE(PS)
C
c CASE(0)
C
C
PRBS.NETAP=1
C
C
C Chequeo de continuidad y calculo de coord. locales
C
CALL CONECTIVIDAD(PRBS,SOLS,CNTS,PRBS.ELTP,ELTS,NODS)
C
C Montaje de las condiciones de contorno
C
CALL MONTAJE_CDC(PRBS,SOLS,MATS,CNTS,PRBS.ELTP,ELTS,NODS)
C !xx_contacto
C Montaje de las condiciones de contacto !xx_contacto
C !xx_contacto
CALL CONTACTO_CDC(PRBS,SOLS,MATS,ZCTOS,CNTS,ELTS,NODS) !xx_contacto
C
C Escritura de datos de entrada en fichero de salida
SELECT CASE(PS)
C
CASE(0)
C
CALL ESC_ENT(PRBS,SOLS,MATS,CNTS,ELTS,NODS,PNTS)
CASE(1)
C
CONTINUE
C
CASE DEFAULT
WRITE(F_W,01000)
WRITE(*,01000)
END SELECT
C !xx_apoyo
C Ecuaciones para evitar movimientos de sólido rígido !xx_apoyo
C !xx_apoyo
CALL CALC_APOYO(PRBS,SOLS,NODS,SVT,SW) !xx_apoyo
C
C Montaje de los puntos de colocacion
C
CALL MONTAJE_PC(PRBS,SOLS,MATS,CNTS,PRBS.ELTP,ELTS,NODS,PCS)
C
C Montaje de las ecuaciones integrales de contorno
C
CALL MONTAJE_BIE(PRBS,SOLS,MATS,CNTS,ELTS,NODS,PCS,
@ SAijS,SFiS,SVT,SW)
C
C
C
C Chequeo del algoritmo de solucion
C
IF((PRBS.NETAP.GT.1).AND.(PRBS.TF.TIP.EQ.11))THEN
CALL TIP_FORM(PRBS.TF,12)
END IF
IF((PRBS.NETAP.GT.1).AND.(PRBS.TF.TIP.EQ.21))THEN
CALL TIP_FORM(PRBS.TF,22)
END IF
C
C
C Calculo de las tensiones termicas
C
CALL TENS_TERM(PRB,PRBS,SOLS,MATS,CNTS,ELTS,NODS,PCS,SAijS,SFiS,
@ NODT,PRBT,ALFA,MATT)
C
C
C CALL ESC_SAIJ(PRBS,SAIJS,SFIS)
C
C
C ELECCION DEL ALGORITMO DE SOLUCION !xx_contacto
C !xx_contacto
SELECT CASE(PRBS.TF.TIP) !xx_contacto
C !xx_contacto
CASE(11,21) !xx_contacto
C
C SOLUCION DIRECTA
C
CALL ALG_DIRECTO(PRBS,SOLS,MATS,CNTS,ELTS,NODS,ZCTOS,
@ SAijS,SFiS,SXjS,SVT,SW)
C !xx_contacto
CASE(12,22) !xx_contacto
C !xx_contacto
C ALGORITMO INCREMENTAL CON ESCALADO DE LA SOLUCION !xx_contacto
C !xx_contacto
CALL ALG_INCREMENTAL(PRBS,SOLS,MATS,CNTS,ELTS,NODS,ZCTOS, !xx_contacto
@ SAijS,SFiS,SXjS,SVT,SW) !xx_contacto
C !xx_contacto
CASE(13,23) !xx_contacto
C !xx_contacto
C ALGORITMO INCREMENTAL CON ESCALADO DE LA SOLUCION !xx_contacto
C Y MONTAJE DE TODAS LAS ECUACIONES (BIE+CTO) EN CADA PASO !xx_contacto
C !xx_contacto
CALL ALG_ACTUAL(PRBS,SOLS,MATS,CNTS,ELTS,NODS,PCS,ZCTOS, !xx_contacto
@ SAijS,SFiS,SXjS,SVT,SW) !xx_contacto
C !xx_contacto
CASE(14,24) !xx_contacto
C !xx_contacto
C ALGORITMO ITERATIVO !xx_contacto
C !xx_contacto
CALL ALG_ITERATIVO(PRBS,SOLS,MATS,CNTS,ELTS,NODS,ZCTOS, !xx_contacto
@ SAijS,SFiS,SXjS,SVT,SW) !xx_contacto
C !xx_contacto
CASE DEFAULT !xx_contacto
WRITE(F_W,00002) PRBS.TF.TIP !xx_contacto
WRITE (*,00002) PRBS.TF.TIP !xx_contacto
STOP !xx_contacto
END SELECT !xx_contacto
C
C Borrado del sistema de ecuaciones
C
DO i=1,PRBS.NAP+PRBS.NBIE+PRBS.NCE
DO j=1,PRBS.NINC
SAijS(i,j)=0.d0
END DO
SFiS(i)=0.d0
END DO
C
C
CCCCCCCCCCCCCCCCCCCCCCCCC
C FIN DE LA RUTINA
CCCCCCCCCCCCCCCCCCCCCCCCC
C
C
C FORMATOS DE ESCRITURA
C
01000 FORMAT(/,2x,'ERROR 01000 - RESOL_S_TERM - Error Programa ')
C
C
00002 FORMAT(/,2x,'ERROR 00002 - Resol_s_term - Tipo de formulacion (',
@ i2,') erronea.')
C
RETURN
C FINAL DEL PROGRAMA
C
END
C
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
