C Rutina 103
C-------------------------------------------------------ITERA_INTERP---
C
C SUBRUTINA ITERA_INTERP
C
C ALGORITMO INCREMENTAL
C
C EN CASOS:
C
C 33 : Nodo en contacto que permanece en contacto
C 13 : Nodo libre que entra en contacto
C
C AJUSTA DELTA PARA QUE EFECTIVAMENTE EL NODO QUEDE EN CONTACTO
C
SUBROUTINE ITERA_INTERP(CNOD,CNTB,ELT,NOD,LAMBDA,BC,DIS,COD)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CNOD : Nodo
INTEGER CNOD
C CNTB : Contorno B de la zona de contacto
RECORD /S_contorno/ CNTB
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 LAMBDA : Coef. de pond. de la sol. del incr. de carga
REAL*8 LAMBDA
C BC : Condicion de contorno final
INTEGER BC
C DIS : Distancia para la que se produce el cruce
REAL *8 DIS
C COD = 1 : Nodo inicial
C = 2 : Nodo intermedio
C = 3 : Nodo final
INTEGER COD
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I : Contador
INTEGER I
C Contadores
INTEGER CELT,N1,N2
C xA,yA: Posicion del nodo
REAL*8 xA,yA
C x1B,x2B,y1B,y2B: Posicion de los extremos del elemento
REAL*8 x1B,x2B,y1B,y2B
C P,L :
REAL*8 P,L
C tx,ty,mod :
REAL*8 tx,ty
C
CCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCC
C
xA=NOD(cnod).XYZ(1)+NOD(cnod).D_SOLU(1)
yA=NOD(cnod).XYZ(2)+NOD(cnod).D_SOLU(2)
C
CCCCCCCCCCCCCCC
C RUTINA
CCCCCCCCCCCCCCC
C
C Bucle que recorre los elementos de B
C
i=0
DO WHILE (i.LT.CNTB.NELT)
i=i+1
celt=CNTB.ELT(i)
n1=ELT(celt).NI(1)
n2=ELT(celt).NI(2)
C
x1B=NOD(n1).XYZ(1)+NOD(n1).D_SOLU(1)
y1B=NOD(n1).XYZ(2)+NOD(n1).D_SOLU(2)
x2B=NOD(n2).XYZ(1)+NOD(n2).D_SOLU(1)
y2B=NOD(n2).XYZ(2)+NOD(n2).D_SOLU(2)
C
C Vector tangente: (tx,ty)
C
IF(x1B.EQ.x2B)THEN
IF(y1B.EQ.y2B)THEN
WRITE(F_W,11120) 1
WRITE (*,11120) 1
STOP
ELSE
tx=0.d0
ty=SIGN(1.d0,(y2B-y1B))
END IF
ELSE
IF(y1B.EQ.y2B)THEN
tx=SIGN(1.d0,(x2B-x1B))
ty=0.d0
ELSE
tx=(x2B-x1B)
ty=(y2B-y1B)
L=DSQRT(tx*tx+ty*ty)
tx=tx/L
ty=ty/L
END IF
END IF
C
C Distancia y posicion
C
p=((xA-X1B)*tx+(yA-y1B)*ty)/L
C
IF((p.GE.0.d0).AND.(p.LE.1.d0))THEN
C1111
dis=dis-(-(xA-X1B)*ty+(yA-y1B)*tx)
C1111
i=CNTB.NELT
END IF
C
END DO
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
11120 FORMAT(/,2x,'ERROR 11120 - ITERA_INTERP - Error Programa',
@ ' (#',i1,'#)')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
