C
C----------------------------------------------------APOYO_POT---
C
C SUBRUTINA APOYOTERMICO
C
C MONTAJE DE LAS ECUACIONES PARA APOYAR PROBLEMA TÉRMICO
C
SUBROUTINE APOYO_POT(PRB,SOL,NOD,SAij,SFi)
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 SOL : Solidos del problema
RECORD /S_solido/ SOL(N_MAX_SOL)
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
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,K,Q : Contadores
INTEGER i,k,q
C CINC,CNOD,CSOL : Incognita,nodo y solido actual
INTEGER CINC,CNOD,CSOL
C NSOL : Numero de solidos del problema
INTEGER NSOL
C NAPO : Numero de apoyos del problema
INTEGER NAPO
C CEC : Contador del numero de ecuaciones
INTEGER CEC
C NECS : Numero de ecuaciones de apoyo térmico que hay que montar
INTEGER NECS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
NSOL=PRB.NSOL
napo=PRB.NAP
cec=PRB.NAP
C
CCCCCCCCCCCCCCCC
C RUTINA
CCCCCCCCCCCCCCCC
C
C
C BUCLE QUE RECORRE LOS SOLIDOS
C
DO i = 1 , NSOL
csol=SOL(i).NUM
IF((SOL(csol).APOPOT.APO).EQ.4)THEN
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C COMPROBACION DE QUE NODO DE APOYO ES UN NODO DEL SOLIDO
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
IF ((SOL(csol).NOD(1)).LE.(SOL(csol).APOPOT.NOD).AND.
@ (SOL(csol).APOPOT.NOD).LE.
@ (SOL(csol).NOD(SOL(csol).NNOD))) THEN
CONTINUE
ELSE
WRITE(F_W,01204) SOL(csol).APOPOT.NOD,SOL(csol).NUM
WRITE(*, 01204) SOL(csol).APOPOT.NOD,SOL(csol).NUM
STOP
END IF
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C MONTAJE ECUACION ADICIONAL Tapo=DATO
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
PRB.NINC=PRB.NINC+SOL(csol).APOPOT.NAP
cec=cec+SOL(csol).APOPOT.NAP
napo=napo+SOL(csol).APOPOT.NAP
SOL(csol).APOPOT.NINC=PRB.NINC
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C MONTAJE FILA EN EL SISTEMA DE ECUACIONES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C Bucle que recorre los nodos/incógnitas del solido (FILA de SAIJ)
C
DO k=1,SOL(csol).NNOD
cnod=SOL(csol).NOD(k)
IF(cnod.EQ.SOL(csol).APOPOT.NOD)THEN
cinc=NOD(cnod).INC(1)
SAIJ(cec,cinc)=1
ELSE
cinc=NOD(cnod).INC(1)
SAIJ(cec,cinc)=0
END IF
C
C Fin bucle que recorre los nodos/incognitas del solido
C
END DO
C
C FILA de SFi
C
SFi(napo)=SOL(csol).APOPOT.TEMP
C
END IF
C
C FIN DEL BUCLE QUE RECORRE LOS SOLIDOS
C
END DO
C
NECS=napo-PRB.NAP
q=0
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C MONTAJE COLUMNAS EN EL SISTEMA DE ECUACIONES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
C BUCLE QUE RECORRE LOS SOLIDOS
C
DO i = 1 , NSOL
csol=SOL(i).NUM
IF((SOL(csol).APOPOT.APO).EQ.4)THEN
q=q+1
C
C Bucle que recorre los nodos/incógnitas del solido (FILA de SAIJ)
C
DO k=1,SOL(csol).NNOD
cinc=NECS+k
SAIJ(cinc,PRB.NINC-napo+q)=1
C
C Fin bucle que recorre los nodos/incognitas del solido
C
END DO
C
NECS=NECS+SOL(csol).NNOD
C
END IF
C
C DO i =NECS, NSOL
C
C FIN DEL BUCLE QUE RECORRE LOS SOLIDOS
C
END DO
C
C Actualizo PRB.NAPO para introducir debajo de las ecuaciones
C de apoyo las BIEs y las NCEs
C
PRB.NAP=napo
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
01204 FORMAT(/,2x,'ERROR 01202 - APOYO_POT',/,
@ ' Nodo de apoyo',i4,' del solido ',i2,' fuera de rango.')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
