C Rutina 090
C-----------------------------------------------------MONTAJE_PNTS_INT---
C
C SUBRUTINA MONTAJE_PNTS_INT
C
C CALCULO DE LA SOLUCION EN LOS PUNTOS INTERNOS
C
SUBROUTINE MONTAJE_PNTS_INT(PRB,SOL,MAT,CNT,ELT,NOD,PNT)
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 SOL : Solidos del problema
RECORD /S_solido/ SOL(N_MAX_SOL)
C MAT : Materiales
RECORD /S_material/ MAT(N_MAX_SOL)
C CNT : Contornos del problema
RECORD /S_contorno/ CNT(N_MAX_CNT)
C ELT : Elementos del problema
RECORD /S_elemento/ ELT(N_MAX_ELT)
C NOD : Nodos del problema
RECORD /S_nodo/ NOD(N_MAX_NOD)
C PINT : Puntos internos del problema
RECORD /S_pntint/ PNT(N_MAX_PINT)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J,M,N,P,Q,R,S : Contadores
INTEGER i,j,m,n,p,q,r,s
C CSOL,CCNT,CELTP,CELT,CNOD,CMAT,CPNT,CNDTP,CPNTP :Contadores
INTEGER CSOL,CCNT,CELTP,CELT,CNOD,CMAT,CPNT,CNDTP,CPNTP
C Hi(,): Matriz de las ctes de integracion en el calculo de tensiones
C Hi(1,*)=Hrr*
C Hi(2,*)=Hrz*
C Hi(3,*)=Hzr*
C Hi(4,*)=Hzz*
C Hi(5,*)=Hoo*
REAL*8 Hi(5,8)
C Eij(,): Matriz de las ctes de integracion en el calculo de desplazamientos
RECORD /S_matriz_Eij/ Eij
C x(): Valor de las variables del elemento
REAL*8 x(N_MAX_VARELT)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALCULO DE LOS DESPLAZAMIENTOS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C BUCLE QUE RECORRE LOS SOLIDOS
C
DO i = 1 , PRB.NSOL
csol=SOL(i).NUM
cmat=MAT(csol).TIP
C
C BUCLE QUE RECORRE LOS PUNTOS INTERNOS DEL SOLIDO
C
DO j = 1 , SOL(csol).NPINT
cpnt=SOL(csol).PINT(j)
cpntp=PNT(cpnt).TIP
DO q=1,PRB.PINTP(cpntp).NVAR
PNT(cpnt).SOLU(q)=0.d0
END DO
C
C Bucle que recorre los contornos del solido
C
DO m = 1 , SOL(csol).NCNT
ccnt=SOL(csol).CNT(m)
C
C Bucle que recorre los elementos
C
DO n = 1 , CNT(ccnt).NELT
celt=CNT(ccnt).ELT(n)
celtp=ELT(celt).TIP
C
C Punto de integracion fuera del elemento
C
DO p=1,2
DO q=1,8
Eij.Eij(p,q)=0.d0
END DO
END DO
C
CALL CALC_CTES_INT_DESP(PRB,PNT(cpnt),MAT,ELT(celt),NOD,
@ Eij)
C
DO p=1,5
DO q=1,8
Hi(p,q)=0.d0
END DO
END DO
C
CALL CALC_CTES_INT_TENS(PRB,PNT(cpnt),MAT,ELT(celt),NOD,
@ Hi)
C
C BUCLE QUE RECORRE LOS NODOS DEL ELEMENTO
C
DO r=1 , PRB.ELTP(celtp).NVE
X(r)=0.d0
END DO
C
DO p = 1 , PRB.ELTP(celtp).NNI
cnod=ELT(celt).NI(p)
cndtp=PRB.ELTP(celtp).TNI(p)
DO r = 1 , PRB.ELTP(celtp).NVE
DO s = 1 , PRB.NDTP(cndtp).NINC
X(r)=X(r)+PRB.ELTP(celtp).Rij(p,r,s)*
@ NOD(cnod).SOLU(s)
END DO
END DO
C
C FIN BUCLE QUE RECORRE LOS NODOS DEL ELEMENTO
C
END DO
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C C
C CALCULO DE LOS DESPLAZAMIENTOS C
C C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
DO q=1,PRB.TP.DIM
DO p=1,8
PNT(cpnt).SOLU(q)=PNT(cpnt).SOLU(q)-Eij.Eij(q,p)*X(p)
END DO
END DO
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C C
C CALCULO DE LAS TENSIONES C
C C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
DO q=PRB.TP.DIM+1,PRB.PINTP(cpntp).NVAR
DO p=1,8
s=q-PRB.TP.DIM
PNT(cpnt).SOLU(q)=PNT(cpnt).SOLU(q)+Hi(s,p)*X(p)
END DO
END DO
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 BUCLE QUE RECORRE LOS PUNTOS INTERNOS DEL SOLIDO
C
END DO
C
C FIN BUCLE QUE RECORRE LOS SOLIDOS
C
END DO
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
07201 FORMAT(/,2x,'ERROR 07201 - MONTAJE_PNTS_INT',/,
@ ' Tipo de punto interno ',i2,'. No implementado.')
07202 FORMAT(/,2x,'ERROR 07202 - MONTAJE_PNTS_INT',/,
@ ' Tipo de material ',i2,'. No implementado.')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
