C Rutina 013
C-------------------------------------------------------LEE_CONTACTO---
C
C SUBRUTINA LEE_CONTACTO
C
C LECTURA DE LAS ZONAS DE CONTACTO DE UN PROBLEMA
C
SUBROUTINE LEE_CONTACTO(PRB,ZCTO,NZCTO,CNT)
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 ZCTO : Zonas de contacto
RECORD /S_zncto/ ZCTO(N_MAX_ZCTO)
C NZCTO : Numero de zonas de contacto
INTEGER NZCTO
C CNT : Contornos del problema
RECORD /S_contorno/ CNT(N_MAX_CNT)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES LOCALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I : Contadores
INTEGER I,J
C CNTA : Contorno que controla los desplazamientos
INTEGER CNTA
C CNTB : Contorno que controla las tensiones
INTEGER CNTB
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C LECTURA ZONAS DE CONTACTO
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
READ(F_E,*) NZCTO
C
IF(NZCTO.GT.N_MAX_ZCTO)THEN
WRITE(F_W,01701) NZCTO,N_MAX_ZCTO
WRITE (*,01701) NZCTO,N_MAX_ZCTO
STOP
END IF
C
C ZONAS DE CONTACTO
C
DO i=1,NZCTO
READ(F_E,*) ZCTO(i).NUM,ZCTO(i).TIP,
@ ZCTO(i).CNT(1),ZCTO(i).CNT(2),ZCTO(i).MU
C
cntA=0
cntB=0
j=0
DO WHILE((cntA.EQ.0).OR.(cntB.EQ.0))
j=j+1
IF(j.GT.N_MAX_CNT)THEN
IF(cntA.EQ.0)THEN
WRITE(F_W,01703) NZCTO,ZCTO(i).CNT(1)
WRITE (*,01703) NZCTO,ZCTO(i).CNT(1)
STOP
END IF
IF(cntB.EQ.0)THEN
WRITE(F_W,01703) NZCTO,ZCTO(i).CNT(2)
WRITE (*,01703) NZCTO,ZCTO(i).CNT(2)
STOP
END IF
ELSE IF((CNT(j).NUM).EQ.(ZCTO(i).CNT(1)))THEN
cntA=j
ELSE IF((CNT(j).NUM).EQ.(ZCTO(i).CNT(2)))THEN
cntB=j
END IF
END DO
ZCTO(i).CNT(1)=cntA
ZCTO(i).CNT(2)=cntB
C
SELECT CASE(ZCTO(i).TIP)
C
CASE(11)
C
C CONTORNO FICTICIO
C
ZCTO(i).MU=1.0d30
C
CASE(21)
C
C CONTORNO DE CONTACTO
C
SELECT CASE(PRB.TF.TIP)
CASE(12,13,14,22,23,24)
CONTINUE
CASE(11,21)
WRITE(*,01707) PRB.TF.TIP
borrar CASE(11)
borrar WRITE(*,01705)
borrar CALL TIP_FORM(PRB.TF,21)
borrar CASE(21)
borrar WRITE(*,01705)
borrar CALL TIP_FORM(PRB.TF,22)
CASE DEFAULT
WRITE(F_W,01706) PRB.TF.TIP
WRITE( *,01706) PRB.TF.TIP
STOP
END SELECT
C
CASE DEFAULT
WRITE(F_W,01704) i,ZCTO(i).TIP
WRITE( *,01704) i,ZCTO(i).TIP
STOP
END SELECT
C
IF(CNT(cntA).NELT.GT.CNT(cntB).NELT)THEN
WRITE(*,01702) i,cntA,cntB
END IF
END DO
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA. SUBRUTINA ENTRADA
C
01700 FORMAT(/,2x,'ERROR 01700 - LEE_CONTACTO - Error Programa (#',
@ i1,'#)')
01701 FORMAT(/,2x,'ERROR 01701 - LEE_CONTACTO - ',
@ 'Numero de zonas de contacto (',i2,
@ ') fuera del rango permitido (',i2,').')
01702 FORMAT(/,2x,'WARNING 01702 - LEE_CONTACTO - Zona de contacto',i2
@ ,/,' El contorno A (',i2,') tiene una discretizacion mas',
@ ' fina que el contorno B (',i2,').',/,
@ ' Para obtener una solucion mas precisa se sugiere',
@ ' cambiar ambos contornos.')
01703 FORMAT(/,2x,'ERROR 01703 - LEE_CONTACTO - Zona de contacto',i2
@ ,/,' El contorno especificado (',i2,') no existe.')
01704 FORMAT(/,2x,'ERROR 01704 - LEE_CONTACTO',/,
@ ' Zona de contacto: ',i4,'. Tipo de zona (',i2,
@ ') no implementada.')
borrar01705 FORMAT(/,2x,'WARNING 01705 - LEE_CONTACTO',/,
borrar @ ' El problema es de contacto debil y por tanto se',
borrar @ ' usara',/,' el algoritmo incremental con escalado de',
borrar @ ' la solucion')
01706 FORMAT(/,2x,'ERROR 01706 - LEE_CONTACTO',/,
@ ' Tipo de formulacion (',i2,') no implementada.')
01707 FORMAT(/,2x,'WARNING 01707 - LEE_CONTACTO',/
@ ,2x,'El problema tiene zonas de contacto,',/
@ ,2x,'pero se ha definido un algoritmo'
@ ,' de solucion directa (TF=',i2,').',/
@ ,2x,'Por tanto no se comprobara la validez de las'
@ ,' condiciones de contacto definidas.')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
