C
C LA SUBRUTINA INICOMA INICIALIZA LAS VARIABLES COMPLEJAS QUE C
C INTERVENDRAN EN LAS INTEGRACIONES ANALITICAS C
C C
C-------------------------------------------------------------------------C
SUBROUTINE INICOMA(B11,B12,B16,B22,B26,B66,MU1,MU2,P,Q,CUK,
@ QP,PP,QQP,QQQ)
C
IMPLICIT NONE
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS (ENTRADA)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C NTIP ................... tipo de problema:
C 0 ==> deformacion plana
C 1 ==> tension plana
C INTEGER NTIP
C E(3) ........... modulos de elasticidad E11, E22, E33
C POIS(3) ........ coeficiente de Poison POIS12, POIS13, POIS23
C G12 ........... modulo de cortadura G12
C REAL*8 E(3),POIS(3),G12
C
C B11 ........... coeficiente B11 de la ley de comportamiento
C B12 ........... coeficiente B12 de la ley de comportamiento
C B22 ........... coeficiente B22 de la ley de comportamiento
C B66 ........... coeficiente B66 de la ley de comportamiento
C
REAL*8 B11,B12,B22,B16,B26,B66
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS (SALIDA)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
complex*16 P(2,2),Q(2,2),CUK(2)
complex*16 MU1,MU2
complex*16 QP(2,2,2),PP(2,2,2)
complex*16 QQQ(3,2,2),QQP(3,2,2)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
INTEGER I,J,L,N
REAL*8 PI
COMPLEX*16 IC
C PI .... sin comentario
C IC .... unidad imaginaria
IC=(0.D00,1.D00)
PI=4.D00*ATAN(1.D00)
C
CCCCCCCCCCCCCCC
C RUTINA
CCCCCCCCCCCCCCC
C
C Cálculo de las raices del polinomio característico de cuarto grado.
C
CALL dxzrhqr(B11,B22,B12,B16,B26,B66,MU1,MU2)
C
C Si las raices del polinomio caracteristico son diferentes
C
c MU1=-MU1
c MU2=-MU2
DO I=1,2
DO J=1,2
P(I,J)=0.D00
Q(I,J)=0.D00
END DO
END DO
P(1,1)=B11*MU1**2.D00+B12-B16*MU1
P(1,2)=B11*MU2**2.D00+B12-B16*MU2
P(2,1)=B12*MU1+B22/MU1-B26
P(2,2)=B12*MU2+B22/MU2-B26
Q(1,1)=-MU1
Q(1,2)=-MU2
Q(2,1)=(1.d00,0.d00)
Q(2,2)=(1.d00,0.d00)
DO I=1,2
CUK(I)=0.D00
END DO
CUK(1)=1.D00/(2.D00*(-B11*MU1**3.D00+B16*MU1**2.D00+
@B22/MU1-B26))
CUK(2)=1.D00/(2.D00*(-B11*MU2**3.D00+B16*MU2**2.D00+
@B22/MU2-B26))
DO I=1,2
DO J=1,2
DO L=1,2
QP(I,J,L)=0.D00
PP(I,J,L)=0.D00
END DO
END DO
END DO
DO I=1,2
DO J=1,2
DO L=1,2
QP(J,I,L)=Q(J,L)*P(I,L)*CUK(L)/(PI*IC)
PP(I,J,L)=P(I,L)*P(J,L)*CUK(L)/(PI*IC)
END DO
END DO
END DO
DO J=1,3
DO N=1,2
DO L=1,2
QQP(J,N,L)=0.D00
QQQ(J,N,L)=0.D00
END DO
END DO
END DO
C
C Nucleos para el calculo del tensor de tensiones en puntos internos
C *** Caso de raices diferentes ***
C
DO N=1,2
C
C N indica la coordenada segun los ejes principales del material
C
DO L=1,2
C
C L indica la raiz del polinomio caracteristico
C
QQP(1,N,L)=Q(1,L)*Q(1,L)*P(N,L)*CUK(L)/(PI*IC)
QQP(2,N,L)=Q(1,L)*Q(2,L)*P(N,L)*CUK(L)/(PI*IC)
QQP(3,N,L)=Q(2,L)*Q(2,L)*P(N,L)*CUK(L)/(PI*IC)
QQQ(1,N,L)=Q(1,L)*Q(1,L)*Q(N,L)*CUK(L)/(PI*IC)
QQQ(2,N,L)=Q(1,L)*Q(2,L)*Q(N,L)*CUK(L)/(PI*IC)
QQQ(3,N,L)=Q(2,L)*Q(2,L)*Q(N,L)*CUK(L)/(PI*IC)
END DO
END DO
RETURN
END
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
