C Rutina 062
C-----------------------------------------------------CALC_EIJ_POT---
C
C SUBRUTINA CALC_EIJ_POT
C
C CALCULO DE LAS CONSTANTES DE INTEGRACION Y LA MATRIZ EIJ
C DESDE UN P.C. A LOS ELEMENTOS DE UN SOLIDO (INCLUYENDO
C EL TERMINO LIBRE)
C
SUBROUTINE CALC_EIJ_POT(PRB,PC,MAT,SOL,CNT,ELT,NOD,NPG)
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 PC : Punto de colocacion
RECORD /S_coloc/ PC
C MAT : Material
RECORD /S_material/ MAT
C SOL : Solido
RECORD /S_solido/ SOL
C CNT : Contornos
RECORD /S_contorno/ CNT(N_MAX_CNT)
C ELT : Elementos
RECORD /S_elemento/ ELT(N_MAX_ELT)
C NOD() : Nodos del problema
RECORD /S_nodo/ NOD(N_MAX_NOD)
C NPG : Numero de puntos de Gauss
INTEGER NPG
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J : Contadores
INTEGER I,J
C CCNT,CELT,CELTP : Contadores
INTEGER CCNT,CELT,CELTP
C TI : Tipo de integracion
INTEGER TI
C TP : Tipo de problema
INTEGER TP
C COD = 1 : P.C. no pertenece al elemento
C 2 : P.C. en el origen del elemento
C 3 : P.C. en el interior del elemento
C 4 : P.C. en el extremo del elemento
INTEGER COD
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
TI=PRB.TI.TIP
TP=PRB.TP.TIP
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALC. CONST. INTEGRACION
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C Bucle que recorre los contornos del solido
C
DO i = 1 , SOL.NCNT
ccnt=SOL.CNT(i)
C
C Bucle que recorre los elementos
C
DO j = 1 , CNT(ccnt).NELT
celt=CNT(ccnt).ELT(j)
celtp=ELT(celt).TIP
C
SELECT CASE(celtp)
CASE(1)
C
C ELEMENTO LINEAL CONTINUO
C
CALL PC_ELT(PC,ELT(celt),NOD,COD)
C
CASE DEFAULT
WRITE(F_W,08312) celtp
WRITE (*,08312) celtp
STOP
END SELECT
C
C CALCULA CTES DE INTEGRACION
C
SELECT CASE (TI)
C
CASE(2)
C
C INTEGRACION ANALITICA
C
SELECT CASE (TP)
C
CASE(41)
C
C Potencial 2D
C
CALL Eij_ANA_POT(PC,MAT,ELT(celt),NOD,COD)
CASE DEFAULT
WRITE(F_W,08314) TI , TP
WRITE (*,08314) TI , TP
STOP
END SELECT
C
CASE(3)
C
C INTEGRACION MIXTA
C
SELECT CASE (TP)
CASE(41)
C
C Potencial 2D plano
C
CALL Eij_MIX_POT(PRB,MAT,PC,ELT(celt),NOD,COD,NPG)
C
CASE DEFAULT
WRITE(F_W,08314) TI , TP
WRITE (*,08314) TI , TP
STOP
END SELECT
C
CASE DEFAULT
WRITE(F_W,08313) TI
WRITE (*,08313) TI
STOP
END SELECT
C
C
C Fin bucle que recorre los elementos
C
END DO
C
C Fin bucle que recorre los contornos del solido
C
END DO
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
08310 FORMAT(/,2x,'ERROR 08310 - CALC_EIJ_POT - Error Programa (#',
@ i1,'#)')
08311 FORMAT(/,2x,'ERROR 08311 - CALC_EIJ_POT - Tipo de punto de',
@ ' colocacion ',i2,' no implementado.')
08312 FORMAT(/,2x,'ERROR 08312 - CALC_EIJ_POT - Tipo de elemento ',
@ i2,' no implementado.')
08313 FORMAT(/,2x,'ERROR 08313 - CALC_EIJ_POT',/,' Tipo de integracion '
@ ,i2,' no implementada para ningun problema.')
08314 FORMAT(/,2x,'ERROR 08314 - CALC_EIJ_POT',/,' Tipo de integracion '
@ ,i2,' no implementada para el problema tipo: ',i2,'.')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
