C Rutina 020
C--------------------------------------------------------CONECTIVIDAD---
C
C SUBRUTINA CONECTIVIDAD
C
C CHEQUEO DE LA CONTINUIDAD ENTRE ELEMENTOS
C Y CALCULO DE LAS COORDENADAS LOCALES DE LOS NODOS
C
SUBROUTINE CONECTIVIDAD(PRB,SOL,CNT,ELTP,ELT,NOD)
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 CNT : Contornos del problema
RECORD /S_contorno/ CNT(N_MAX_CNT)
C ELTP : Tipos de elemento del problema
RECORD /S_tipelt/ ELTP(N_MAX_TPEL)
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
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C TIP : Tipo de problema
INTEGER TIP
C NSOL : Numero de solidos
INTEGER NSOL
C I,J,K : Contadores
INTEGER I,J,K
C Contadores
INTEGER CSOL,CCNT,CELTP,CELT,UELT,PNOD,ANOD,CNOD,UNOD
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
TIP=PRB.TP.TIP
NSOL=PRB.NSOL
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CONTINUIDAD ENTRE CONTORNOS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
pnod=0
anod=0
C
C Bucle que recorre los solidos
C
DO csol = 1 , NSOL
C
C Bucle que recorre los contornos del sólido
C
i=1
DO WHILE (i.LE.SOL(csol).NCNT)
ccnt=SOL(csol).CNT(i)
C
C Primer nodo del contorno
C
celt=CNT(ccnt).ELT(1)
cnod=ELT(celt).NG(1)
C
C Ultimo nodo del contorno
C
uelt=CNT(ccnt).ELT(CNT(ccnt).NELT)
unod=ELT(uelt).NG(PRB.ELTP(ELT(uelt).TIP).NNG)
C
IF(pnod.EQ.0)THEN
C Contorno "inicial"
IF(cnod.EQ.unod)THEN
C Contorno cerrado
i=i+1
ELSE
C Contorno no cerrado
pnod=cnod
anod=unod
i=i+1
END IF
ELSE
C Contorno "intermedio"
IF(cnod.EQ.ANOD)THEN
C Contorno continuo con el anterior
IF(unod.EQ.pnod)THEN
C Ultimo contorno de un conjunto cerrado
pnod=0
i=i+1
ELSE
C Contorno "intermedio"
anod=unod
i=i+1
END IF
ELSE
C El contorno anterior era abierto
WRITE (*,03103) pnod,anod
pnod=0
END IF
END IF
C
C Fin bucle que recorre los contornos del sólido
C
END DO
C
C Fin bucle que recorre los solidos
C
END DO
C
IF(pnod.NE.0.AND.pnod.NE.anod)THEN
C El ultimo contorno era abierto
WRITE (*,03103) pnod,anod
END IF
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CONTINUIDAD DENTRO DE CADA CONTORNO
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C Bucle que recorre los solidos
C
DO csol = 1 , NSOL
C
C Bucle que recorre los contornos del solido
C
DO i=1,SOL(csol).NCNT
ccnt=SOL(csol).CNT(i)
anod=0
C
C Bucle que recorre los elementos de un contorno
C
DO j = 1 , CNT(ccnt).NELT
C
C COMPROBACION DE CONTINUIDAD
C
celt=CNT(ccnt).ELT(j)
cnod=ELT(celt).NG(1)
unod=ELT(celt).NG(PRB.ELTP(ELT(celt).TIP).NNG)
IF(anod.NE.0.AND.anod.NE.cnod)THEN
WRITE(F_W,03101) csol,celt,cnod,anod
WRITE (*,03101) csol,celt,cnod,anod
STOP
ELSE
anod=unod
END IF
C
C CALCULO DE LAS COORDENADAS LOCALES DE LOS NODOS
C Y ASIGNACION DE LOS TIPOS DE NODO
C
celtp=ELT(celt).TIP
SELECT CASE(celtp)
CASE (1)
C ELEMENTO LINEAL CONTINUO
ELT(celt).NGLOC(1)= -1.d0
ELT(celt).NGLOC(2)= 1.d0
ELT(celt).NILOC(1)= -1.d0
ELT(celt).NILOC(2)= 1.d0
x DO k=1,ELTP(celtp).NNI
x IF(NOD(ELT(celt).NI(k)).TIP.EQ.0)THEN
x NOD(ELT(celt).NI(k)).TIP=ELTP(celtp).TNI(k)
x ELSE IF(NOD(ELT(celt).NI(k)).TIP.NE.
x @ ELTP(celtp).TNI(k))THEN
x WRITE(F_W,03104) NOD(ELT(celt).NI(k)).NUM
x WRITE (*,03104) NOD(ELT(celt).NI(k)).NUM
x STOP
x END IF
x END DO
DO k=1,ELTP(celtp).NNI
NOD(ELT(celt).NI(k)).TIP=ELTP(celtp).TNI(k)
END DO
CASE DEFAULT
WRITE(F_W,03102) csol,celt,ELT(celt).TIP
WRITE (*,03102) csol,celt,ELT(celt).TIP
STOP
END SELECT
C
C Fin bucle que recorre los elementos de un contorno
C
END DO
C
C Fin bucle que recorre los contornos 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
03101 FORMAT(/,2x,'ERROR 03101 - CONECTIVIDAD - Solido: ',i2,
@ '. Elemento:',i4,
@ /,' El nodo ',i4,' no coincide con el ultimo nodo del elemento',
@ /,' anterior (',i4').')
03102 FORMAT(/,2x,'ERROR 03102 - CONECTIVIDAD - Solido: ',i2,
@ '. Elemento:',i4,
@ /,' Tipo de elemento (',i2,') no implementado.')
03103 FORMAT(/,2x,'WARNING 03103 - CONECTIVIDAD',/,
@ ' Contorno abierto entre los nodos',i4,' y ',i4, '.',/)
03104 FORMAT(/,2x,'ERROR 03104 - CONECTIVIDAD - Nodo: ',i4,
@ /,' Incompatibilidad en el tipo de nodo. Chequear los',
@ ' elementos concurrentes al nodo')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
