C Rutina 101
C-----------------------------------------------------------EQUIL_POT---
C
C SUBRUTINA EQUIL_POT
C
C MONTAJE DE LAS ECUACIONES DE EQUILIBRIO EN FORMA DEBIL
C CORRESPONDIENTES A UN FRAMENTO DE LA ZONA DE CONTACTO
C
C CONTACTO DEBIL
C
SUBROUTINE EQUIL_POT(ELTA,ELTB,NOD,x1A,x2A,x1B,x2B,
@ alfa,SAij,SFi,ceq)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C ELTA,ELTB : Elemento de los solidos A y B corresp. al fragmento
RECORD /S_elemento/ ELTA,ELTB
C NOD : Nodos del problema
RECORD /S_nodo/ NOD(N_MAX_NOD)
C x1A,x2A,x1B,x2B : Posicion local de comienzo y fin de la integracion
C en cada fragemento de la zona de contacto
REAL*8 x1A,x2A,x1B,x2B
C ALFA : Indice de la funcion de forma de la tension virtual
INTEGER ALFA
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 ecuacuones de contacto
INTEGER CEQ
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J,K,M : Contador
INTEGER I,J,K,M
C Contador
INTEGER CINC
C n1B,n2B : Nodos extremos de ELTB
INTEGER n1B,n2B
C LB: Longitud de ELTA
REAL*8 LB
C Ho*H(,) : Matriz de Integrales de las funciones de forma
REAL*8 Ho,H(2,5)
C FFA1,FFA2,FFB1,FFB2 : Valor de las funciones de forma en x1A,x2A,x1B,x2B
REAL*8 FFA1(2),FFA2(2),FFB1(2),FFB2(2)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
n1B=ELTB.NI(1)
n2B=ELTB.NI(2)
C
LB=(x2B-x1B)*DSQRT((NOD(n2B).XYZ(1)-NOD(n1B).XYZ(1))**2.d0+
@ (NOD(n2B).XYZ(2)-NOD(n1B).XYZ(2))**2.d0)/2.d0
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C ECUACIONES DE COMPATIBILIDAD
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
CALL FUNC_FORM(x1A,FFA1,1)
CALL FUNC_FORM(x2A,FFA2,1)
CALL FUNC_FORM(x1B,FFB1,1)
CALL FUNC_FORM(x2B,FFB2,1)
C
C **********
C SOLIDO A
C **********
C
C Para cada extremo del segmento
C
DO i=1,2
DO j=1,5
H(i,j)=0.d0
END DO
END DO
H(1,3) =1.d0
H(2,2) =1.d0
C
DO i=1,2
C
Ho=( FFB2(alfa)*( 2.d0*FFA1(i)+FFA2(i) ) +
@ FFB1(alfa)*( FFA1(i)+2.d0*FFA2(i) ) )*LB/6.d0
C
C Terminos que afectan a la matriz del sistema
C
DO j=1,NOD(ELTA.NI(i)).NINC
cinc=NOD(ELTA.NI(i)).INC(j)
DO m=2,3
SAij(ceq,cinc)=SAij(ceq,cinc)+
@ Ho*H(i,m)*NOD(ELTA.NI(i)).Gij(m,j)
END DO
END DO
C
C Terminos que afectan al vector de carga
C
DO j=1,NOD(ELTA.NI(i)).NBC
k=j+NOD(ELTA.NI(i)).NINC
DO m=2,3
SFi(ceq)=SFi(ceq)-
@ Ho*H(i,m)*NOD(ELTA.NI(i)).Gij(m,k)*
@ NOD(ELTA.NI(i)).VBC(j)
END DO
END DO
C
END DO
C
C **********
C SOLIDO B
C **********
C
C Para cada extremo del segmento
C
DO i=1,2
DO j=1,3
H(i,j)=0.d0
END DO
END DO
H(1,3) =1.d0
H(2,2) =1.d0
C
DO i=1,2
C
Ho=( FFB2(alfa)*( 2.d0*FFB2(i)+FFB1(i) ) +
@ FFB1(alfa)*( FFB2(i)+2.d0*FFB1(i) ) )*LB/6.d0
C
C Terminos que afectan a la matriz del sistema
C
DO j=1,NOD(ELTB.NI(i)).NINC
cinc=NOD(ELTB.NI(i)).INC(j)
DO m=2,3
SAij(ceq,cinc)=SAij(ceq,cinc)+
@ Ho*H(i,m)*NOD(ELTB.NI(i)).Gij(m,j)
END DO
END DO
C
C Terminos que afectan al vector de carga
C
DO j=1,NOD(ELTB.NI(i)).NBC
k=j+NOD(ELTB.NI(i)).NINC
DO m=2,3
SFi(ceq)=SFi(ceq)-
@ Ho*H(i,m)*NOD(ELTB.NI(i)).Gij(m,k)*
@ NOD(ELTB.NI(i)).VBC(j)
END DO
END DO
C
END DO
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
10400 FORMAT(/,2x,'ERROR 10400 - EQUIL_POT - Error Programa (#',
@ i1,'#)')
C
END
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
