C Rutina 064
C-------------------------------------------------------------EIJ_NUM2---
C
C SUBRUTINA EIJ_NUM
C
C CALCULO DE LAS CONSTANTES DE INTEGRACION DESDE UN P.C.
C A UN ELEMENTO INTEGRACION NUMERICA
C
C EL PC NO PERTENECE AL ELEMENTO
C
SUBROUTINE EIJ_NUM2(PC,PRB,MAT,ELT,NOD)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS (ENTRADA)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PC : Punto de colocacion
RECORD /S_coloc/ PC
C PRB : Variables globales del solido
RECORD /S_problema/ PRB
C MAT : Material
RECORD /S_material/ MAT
C ELT : Elemento
RECORD /S_elemento/ ELT
C NOD() : Nodos del problema
RECORD /S_nodo/ NOD(N_MAX_NOD)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C Aij,Bij : Matrices clasicas de constantes de integracion
REAL*8 Aij(2,4),Bij(2,4)
C nd1,nd2 : Extremos del elemento
INTEGER nd1,nd2
C Lk,Tk,N1,N2,Nmod : Longitud, angulo y normal del elemento
REAL*8 Lk,Tk,N1,N2,Nmod
C X() : Coordenadas globales de los puntos de Gauss
REAL*8 X(2)
C SF_U(,),SF_T(,) : Nucleos de la integral de contorno
C SF_U(i,j) : Uij
C SF_T(i,j) : Tij
REAL*8 SF_U(2,2),SF_T(2,2)
C FF() : FUNCIONES DE FORMA
REAL*8 FF(2)
C TP : Tipo de problema
INTEGER TP
C ELTP : Tipo de elemento
INTEGER ELTP
C I,j : Contadores
INTEGER I,J
C
CCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCC
C
DO i=1,2
DO j=1,4
Aij(i,j)=0.d0
Bij(i,j)=0.d0
END DO
END DO
C
TP=PRB.TP.TIP
C
CCCCCCCCCCCCCCC
C RUTINA
CCCCCCCCCCCCCCC
C
ELTP=ELT.TIP
C
C CALCULO DE Aij y Bij
C
SELECT CASE(ELTP)
C
C ELEMENTOS LINEAL CONTINUO
C
CASE(1)
C Nodos extremos
nd1=ELT.NG(1)
nd2=ELT.NG(2)
C Longitud del elemento
Lk= (NOD(nd2).XYZ(1)-NOD(nd1).XYZ(1))**2.d0
Lk=Lk+(NOD(nd2).XYZ(2)-NOD(nd1).XYZ(2))**2.d0
Lk=DSQRT(Lk)
C
C Normal al elemento
C
N1=(NOD(nd2).XYZ(2)-NOD(nd1).XYZ(2))
N2=(NOD(nd1).XYZ(1)-NOD(nd2).XYZ(1))
Nmod=DSQRT(N1*N1+N2*N2)
N1=N1/Nmod
N2=N2/Nmod
C
C Angulo entre los ejes locales y globales del elemento
C
CALL ANGULO(Tk,N1,N2)
C
DO i = 1 , PRB.TI.NPG
C Coordenadas correspondientes a los puntos de Gauss
X(1)=NOD(nd1).XYZ(1)*(1.d0-PRB.TI.S(i))/2.d0
@ +NOD(nd2).XYZ(1)*(1.d0+PRB.TI.S(i))/2.d0
X(2)=NOD(nd1).XYZ(2)*(1.d0-PRB.TI.S(i))/2.d0
@ +NOD(nd2).XYZ(2)*(1.d0+PRB.TI.S(i))/2.d0
C Funciones de forma en los puntos de Gauss
CALL FUNC_FORM(PRB.TI.S(i),FF,ELTP)
C Solucion fundamental en los puntos de Gauss
C Nota: Se llama con '0' porque el PC esta fuera del elemento
CALL SOL_FUND2(PC.XYZ,X,MAT.PRP,SF_U,SF_T,N1,N2,TP)
C Incrementa Aij,Bij
CALL Aij_Bij(Aij,Bij,PRB.TI.W(i),SF_U,SF_T,FF,Lk)
END DO
CASE DEFAULT
WRITE(F_W,08511) ELTP
WRITE (*,08511) ELTP
STOP
END SELECT
C
C Matriz Eij
C
CALL MAT_Eij(PC.ELT,PC.LOC,PC.Eij(ELT.NUM),Aij,Bij,PC.Cij,Tk,
@ ELT.NUM)
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
08510 FORMAT(/,2x,'ERROR 08510 - EIJ_NUM2 - Error Programa (#',i1,'#)')
08511 FORMAT(/,2x,'ERROR 08511 - EIJ_NUM2 - Tipo de elemento (',i2,
@ ') no implementado')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
