C Rutina 100
C-----------------------------------------------MONTAJE_CONTACTO_POT---
C
C SUBRUTINA MONTAJE_CONTACTO_POT
C
C MONTAJE DE LAS ECUACIONES LINEALES DE CONTACTO (PARA TODOS
C LOS CONTACTOS) Y LAS COLUMNAS ADICIONALES QUE INTRODUCEN
C LAS ECUACIONES DE APOYO
C
SUBROUTINE MONTAJE_CONTACTO_POT(PRB,ZCTO,CNT,ELT,NOD,SAij,SFi,SOL)
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 del problema
RECORD /S_zncto/ ZCTO(N_MAX_ZCTO)
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 SAIJ(,) : Matriz del sistema.
REAL*8 SAij(N_MAX_INC,N_MAX_INC)
C SFi() : Termino independiente.
REAL*8 SFi(N_MAX_INC)
C SOL: SOLIDOS
RECORD /S_solido/ SOL(N_MAX_SOL)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I : Contador
INTEGER I
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
PRB.NCE=0
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C MONTAJE ECUACIONES SOLIDO
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C Bucle que recorre las zonas de contacto
C
DO i = 1 , PRB.NZCTO
C
C Calculo de las posiciones de los nodos de A
C en el contorno de B
C
CALL POS_REL_A_POT(CNT,PRB.ELTP,ELT,NOD,ZCTO(i))
C
C Calculo de las posiciones de los nodos de B
C en el contorno de A
C
CALL POS_REL_B(CNT,ELT,NOD,ZCTO(i))
C
C Escritura de las condiciones de contacto
C
CALL ESCRIBE_CONTACTO_POT(CNT,PRB.ELTP,ELT,NOD,ZCTO(i))
C
C Imposicion de las condiciones de contacto
C
CALL CONTACTO_DEBIL_POT(PRB,CNT,PRB.ELTP,ELT,NOD,
@ ZCTO(i),SAij,SFi,SOL)
C
C Fin bucle que recorre las zonas de contacto
C
END DO
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
10001 FORMAT(/,2x,'ERROR 10001 - MONTAJE_CONTACTO_POT',/,
@ ' Zona de contacto: ',i4,'. Tipo de zona (',i2,
@ ') no implementada.')
10004 FORMAT(/,2x,'ERROR 10004 - MONTAJE_CONTACTO_POT',/,
@ ' Nodo (',i4,').',
@ /,' Pertenece a dos zonas de contacto y esta libre',
@ /,' Caso no implementado.')
10005 FORMAT(/,2x,'ERROR 10005 - MONTAJE_CONTACTO_POT -',
@ ' Numero de incognitas (',i4,')'
@ ,/,' distinto del numero de ecuaciones (',i4,').')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
