C Rutina 064
C-------------------------------------------------------------ANA3_TERM---
C
C SUBRUTINA ANA3_TERM
C
C CALCULO DE LAS CONSTANTES DE INTEGRACION DESDE UN P.C.
C A UN ELEMENTO INTEGRACION ANALITICA
C
C EL PC PERTENECE AL ELEMENTO(EN EL EXTREMO)
C
SUBROUTINE ANA3_TERM(MATS,PCS,ELTS,NODS,ALFA,MATT,PRBS)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS (ENTRADA)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C PRBS : Variables globales del problema
RECORD /S_problema/ PRBS
C
C PCS : Punto de colocacion
RECORD /S_coloc/ PCS
C MATS,MATT : Material
RECORD /S_material/ MATS
RECORD /S_material/ MATT
C ELTS : Elemento
RECORD /S_elemento/ ELTS
C NODS() : Nodos del problema
RECORD /S_nodo/ NODS(N_MAX_NOD)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C nd1,nd2 : Extremos del elemento
INTEGER nd1,nd2
C Lk,Tk,N1,N2,Nmod,PI : Longitud, angulo y normal del elemento,PI
REAL*8 Lk,Tk,N1,N2,Nmod,PI
C E,NU,G,ALFA,GAM,K,CTE,ALF : Constantes del material para el sólido
REAL*8 E,NU,G,ALFA,GAM,K,CTE,ALF
C C1,C2,C3,D1,D2
REAL*8 C1,C2,C3,D1,D2
C EIJT : Nodos del problema
C RECORD /S_matriz_EijTERM/ EijT
C ELTP : Tipo de elemento
INTEGER ELTP
C Cij(I,J),Dij(I,J)
REAL*8 Cij(2,2)
REAL*8 Dij(2,2)
C
C
CCCCCCCCCCCCCCC
C RUTINA
CCCCCCCCCCCCCCC
C
ELTP=ELTS.TIP
C
SELECT CASE(ELTP)
C
CASE(1)
C CALCULO DE LAS CONSTANTES NECESARIAS
E=MATS.PRP(1)
NU=MATS.PRP(2)
G=E/(2.d0*(1+NU))
SELECT CASE(PRBS.TP.TIP)
CASE(22,32)
ALF=ALFA
CASE(21,31)
ALF=ALFA*(1.d0+NU)/(1+2.d0*NU)
NU=NU/(1.d0+NU)
C Modificado por: Diógenes Morales Sierra Fecha:13/03/05
C Objeto de la modificación:
C Falta el valor de E para que realice bien TP
C
E=E*(1-NU**2)
CASE DEFAULT
WRITE(F_W,08412) PRBS.TP.TIP
WRITE (*,08412) PRBS.TP.TIP
STOP
END SELECT
GAM=(E*ALF)/(1.d0-(2.d0*NU))
CTE=(GAM*(1.d0-2.d0*NU))/(2.d0*(1-NU))
Pi=DACOS(-1.d0)
K=MATT.PRP(1)
C
C Nodos extremos
nd1=ELTS.NG(1)
nd2=ELTS.NG(2)
C Longitud del elemento
Lk= (NODS(nd2).XYZ(1)-NODS(nd1).XYZ(1))**2.d0
Lk=Lk+(NODS(nd2).XYZ(2)-NODS(nd1).XYZ(2))**2.d0
Lk=DSQRT(Lk)
C
C Normal al elemento
C
N1=(NODS(nd2).XYZ(2)-NODS(nd1).XYZ(2))
N2=(NODS(nd1).XYZ(1)-NODS(nd2).XYZ(1))
Nmod=DSQRT(N1*N1+N2*N2)
N1=N1/Nmod
N2=N2/Nmod
C
C CALCULO DE EIJT PARA ELTS(desde pcs)
C
C
C Calculo de cij
C
c1=(cte*Lk)/(8.d0*PI*G)
c2=c1*DLOG(Lk)
c3=c1-c2
C Almacenamos en EIJT.EIJT(,),EIJT.CIJT(,)
C
Cij(1,1)=-c2*n1
Cij(1,2)=c3*n1
Cij(2,1)=-c2*n2
Cij(2,2)=c3*n2
PCS.EIJT(ELTS.NUM).EijT(1,1)=Cij(1,1)
PCS.EIJT(ELTS.NUM).EijT(1,2)=Cij(1,2)
PCS.EIJT(ELTS.NUM).EijT(2,1)=Cij(2,1)
PCS.EIJT(ELTS.NUM).EijT(2,2)=Cij(2,2)
C
C Calculo de dij
C
d1=(cte*Lk*Lk)/(8.d0*PI*G*K)
d2=d1*DLOG(Lk)
C Almacenamos en EIJT.EIJT(,),EIJT.CIJT(,)
C
Dij(1,1)=((-1.d0/9.d0)*d1-(d2/3.d0)*2.d0)*n2
Dij(1,2)=((1.d0/9.d0)*d1-(d2/3.d0))*n2
Dij(2,1)=((1.d0/9.d0)*d1+(d2/3.d0)*2.d0)*n1
Dij(2,2)=((-1.d0/9.d0)*d1+(d2/3.d0))*n1
PCS.EIJT(ELTS.NUM).EijT(1,3)=Dij(1,1)
PCS.EIJT(ELTS.NUM).EijT(1,4)=Dij(1,2)
PCS.EIJT(ELTS.NUM).EijT(2,3)=Dij(2,1)
PCS.EIJT(ELTS.NUM).EijT(2,4)=Dij(2,2)
C
CASE DEFAULT
WRITE(F_W,08411) ELTP
WRITE (*,08411) ELTP
STOP
END SELECT
C Actualizamos PCS.EIJT
C
C PCS.EIJT(ELTS.NUM)=EIJT
C
C
C
C
C FIN DE LA RUTINA
C
RETURN
C
08411 FORMAT(/,2x,'ERROR 08411 - ANA3_TERM - Tipo de elemento (',i2,
@ ') no implementado')
08412 FORMAT(/,2x,'ERROR 08412 - ANA3_TERM - Tipo de problema (',i2,
@ ') no implementado')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
