C
C-------------------------------------------------------------LECTURA---
C
C SUBRUTINA LECTURA
C
C LECTURA DE DATOS DE ENTRADA
C
SUBROUTINE LECTURA(TP,TF,NSOL,TI,PRB,ETAP,SOL,MAT,CNT,ELT,NOD,
@ PINT,ZCTO,NPG)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES EXTERNAS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C tp,tf,nsol,ti : Variables auxiliares
INTEGER TP,TF,NSOL,TI
C PRB : Variables globales del problema
RECORD /S_problema/ PRB
C ETAP : Etapa de carga
INTEGER ETAP
C SOL : Solidos del problema
RECORD /S_solido/ SOL(N_MAX_SOL)
C MAT : Materiales
RECORD /S_material/ MAT(N_MAX_SOL)
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 PINT : Puntos internos del problema
RECORD /S_pntint/ PINT(N_MAX_PINT)
C ZCTO : Zonas de contacto del problema
RECORD /S_zncto/ ZCTO(N_MAX_ZCTO)
C NPG : Numero de puntos de Gauss
INTEGER NPG
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES INTERNAS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J : Contadores
INTEGER I,J
C CSOL,CCNT,CELT,CNOD,CCDC :Contadores
INTEGER CSOL,CCNT,CELT,CNOD,CPINT,CCDC
C Nombre de los ficheros
CHARACTER*50 NAME
C
CCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VALORES INICIALES
C
CCCCCCCCCCCCCCCCCCCCCCCCCC
C
i=0
j=0
csol=0
ccnt=0
celt=1
cnod=0
cpint=0
ccdc=0
C
CCCCCCCCCCCCCCC
C
C RUTINA
C
CCCCCCCCCCCCCCC
C
C APERTURA DEL FICHERO DE SALIDA
C
WRITE(NAME,*) TRIM(PRB.NAME),"_SAL.TXT"
WRITE(*,*) "Archivo de salida: ",NAME
OPEN(F_S,file=NAME,status='replace')
C
C APERTURA DEL FICHERO DE ERRORES
C
WRITE(NAME,*) TRIM(PRB.NAME),"_ERR.TXT"
WRITE(*,*) "Archivo de errores: ",NAME
OPEN(F_W,file=NAME,status='replace')
C
C APERTURA DEL FICHERO DE CONTACTO
C
IF(W_C.GE.1)THEN
WRITE(NAME,*) TRIM(PRB.NAME),"_CTO.TXT"
WRITE(*,*) "Archivo de contacto: ",NAME
OPEN(F_C,file=NAME,status='replace')
END IF
C
C APERTURA DEL FICHERO DE BASE DE DATOS
C
IF(W_D.GE.1)THEN
WRITE(NAME,*) TRIM(PRB.NAME),"_BDD.TXT"
WRITE(*,*) "Archivo de base de datos: ",NAME
OPEN(F_D,file=NAME,status='replace')
END IF
C
C DATOS DEL TIPO DE FORMULACION
C
CALL TIP_FORM(PRB.TF,tf)
C
C DATOS DEL TIPO DE PROBLEMA
C
CALL TIP_PROBL(PRB,tp)
C
C NUMERO DE SOLIDOS
C
IF(nsol.LT.1.OR.nsol.GT.N_MAX_SOL) THEN
WRITE(F_W,01001)
WRITE(* ,01001)
STOP
ELSE
PRB.NSOL=nsol
END IF
C
C TIPO DE INTEGRACION
C
CALL TIP_INT(PRB.TI,ti)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C LECTURA VARIABLES SOLIDO
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C BUCLE QUE RECORRE LOS SOLIDOS
C
DO csol = 1 , nsol
C
C LECTURA DE DATOS DEL SOLIDO Y DEL MATERIAL
C
CALL LEE_SOLIDO(csol,SOL(csol),MAT(csol),PRB.MATP,PRB.TF.TIP)
C
C CÁLCULO DE LAS CONSTANTES DEL MATERIAL ELÁSTICO COMPLEJAS !MT000225-1
C !MT000225-1
CALL MATERIAL(TP,MAT(csol)) !MT000225-1
!MT000225-1
C LECTURA DE LOS NODOS
C
CALL LEE_NODOS(ETAP,SOL(csol),NOD,PRB.TP.DIM)
C
C BUCLE QUE RECORRE LOS CONTORNOS DEL SOLIDO
C
DO i = 1 , SOL(csol).NCNT
C
ccnt=ccnt+1
SOL(csol).CNT(i)=ccnt
C
C LECTURA DE DATOS DE UN CONTORNO
C
CALL LEE_CONTORNO(ccnt,CNT(ccnt))
CNT(ccnt).SOL=SOL(csol).NUM
C
C LECTURA DE LOS ELEMENTOS DE UN CONTORNO
C
CALL LEE_ELEMENTOS(CNT(ccnt),PRB.ELTP,ELT)
C
C FIN DEL BUCLE QUE RECORRE LOS CONTORNOS DEL SOLIDO
C
END DO
C
C LECTURA DE LAS CONDICIONES DE CONTORNO
C
CALL LEE_CDC(ETAP,SOL(csol),NOD)
C
C LECTURA DE LOS PUNTOS INTERNOS
C
IF(SOL(csol).NPINT.GT.0) THEN
CALL LEE_PINTS(SOL(csol),PINT,PRB.PINTP)
END IF
C
C FIN DEL BUCLE QUE RECORRE LOS SOLIDOS
C
END DO
C !xx_contacto
C LECTURA DE LAS ZONAS DE CONTACTO !xx_contacto
C !xx_contacto
CALL LEE_CONTACTO(PRB,ZCTO,PRB.NZCTO,CNT) !xx_contacto
C
C FIN DE LA RUTINA LECTURA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
01000 FORMAT(/,2x,'ERROR 01000 - LECTURA - Error Programa (#',i1,'#)')
01001 FORMAT(/,2x,'ERROR 01001 - LECTURA - Num. solidos mayor de lo ',
@ 'permitido')
01002 FORMAT(/,2x,'ERROR 01002 - LECTURA - Tipo de elemento no valido')
C
10 FORMAT(a200)
20 FORMAT(a15)
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
