C Rutina 101
C--------------------------------------------------------FRICCION_DEB---
C
C SUBRUTINA FRICC_DEB
C
C MONTAJE DE LAS ECUACIONES DE FRICCION EN FORMA DEBIL
C CORRESPONDIENTES A UN NODO
C
SUBROUTINE FRICC_DEB(PRB,NOD,MU,COD,DIS,SAij,SFi,ceq)
C
IMPLICIT NONE
C
C CONTORNO FICTICIO
C Y
C CONTACTO DEBIL DISCONTINUO
C
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PRB : Tipo de problema
RECORD /S_problema/ PRB
C NOD : Nodo
RECORD /S_nodo/ NOD
C MU : Coeficiente de friccion
REAL*8 MU
C COD = 1 : Nodo inicial
C = 2 : Nodo intermedio (==> tensiones continuas)
C = 3 : Nodo final
INTEGER COD
C DIS : Distancia del nodo al contorno opuesto
REAL*8 DIS
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 CEQ : Contador de ecuaciones de contacto
INTEGER CEQ
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C ECUACION DE FRICCION
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
SELECT CASE(NOD.CBC(1))
C
CASE(31,21)
C
CALL FRICCION_DEB(NOD,MU,COD,DIS,SAij,SFi,ceq)
C
CASE(13,12)
C
CALL FRICCION_DEB(NOD,MU,COD,DIS,SAij,SFi,ceq)
C
CASE(11)
C
CALL FRICCION_DEB(NOD,MU,COD,DIS,SAij,SFi,ceq)
C
CASE DEFAULT
WRITE(F_W,10200) 1
WRITE( *,10200) 1
STOP
END SELECT
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
10200 FORMAT(/,2x,'ERROR 10200 - FRICC_DEB - Error Programa (#',
@ i1,'#)')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
