! ================================================================
! Procedimiento: 
!    Para comprobar resultados del libro, utilizar:
! 	comprob = 0
!    Para determinar las componentes de la fuerza a aplicar, 
!    imponiendo desplazamientos cero al nudo 4. 
!    Las componentes de la fuerza tienen igual módulo,
!    dirección, y sentido contrario a las reacciones obtenidas.
!    Utilizar:
! 	comprob = 1
!    Para comprobar las fuerzas obtenidas en comprob = 1, utilizar: 
!      comprob = 2  
! ================================================================

*CONTROL DEL PROBLEMA
   IDPR,   Ej-8.1
   TITU,   "Ejercicio 8.1. Estructura articulada hiperestática"
   SISUNI, mkN
   TIPEST, EP_NA_XY

*PARAMETROS GENERALES
   h   =  2.0
   lm  =  2.0
   nm  =  3           !número de módulos (proyecto)
   
   al   =  -0.01      !alargamiento inicial de la barra 12 (m)
   Ax  =  1.00e-3     !área de las barras (m^2)
   E   =  200.0e6     !módulo de elasticidad (kPa)

! Parámetros para comprobación de resultados del libro   
   F   = 87.94       !Fuerza aplicada en 4
   alfa  =  0.4597   !Dirección (rad) de F   
   
! Parámetros para obtener el módulo y dirección de la fuerza para desplazamiento nulo en 4 
   Fx  =  78.798 
   Fy  =  -39.046
   ang = atan(-Fy/Fx)      ! angulo alfa obtenido
   mod = sqrt(Fx^2+Fy^2)   ! F obtenido

   comprob = 2   ! 0 = Se comprueban los resultados del libro
                 ! 1 = Se obtiene Fx y Fy
                 ! 2 = Se comprueba que las fuerzas Fx Fy son correctas

*PARAMETROS DE ANALISIS

>COORDENADAS
! nudo      X      Y      Z
   1       0.0    0.0     0.0   L   4	  1   lm   0.0	  0.0	    
   5       nm*lm   h      0.0	  L   4	  1  -lm   0.0	  0.0                                      
  
>MATERIALES
! kM        E     G     æ      Pesp   Cdil     SigmaE
  1         E    0.0    0.3   77.01     0.0    260.0e3 

>PROPIEDADES GEOMETRICAS
! kP      A     Ay    Az    Ix    Iy     Iz   
  1      Ax    0.0   0.0   0.0   0.0.   0.0

>ELEMENTOS
  GRUPO 1 BNA
      1  1  2  1 1   0  0   L   3   1   1   1
      4  4  5  1 1   0  0
      5  5  6  1 1   0  0   L   3   1   1   1
      8  1  8  1 1   0  0  
      9  1  7  1 1   0  0
     10  2  7  1 1   0  0 
     11  3  7  1 1   0  0
     12  2  6  1 1   0  0
     13  3  6  1 1   0  0
     14  4  6  1 1   0  0  
           
>RESTRICCIONES
  TODOS        DZ  GX  GY  GZ ! Estructura plana de nudos articulados
   1   DX  DY
   3       DY
 IF,comprob,=,1,THEN
   4   DX  DY
 ENDIF 

>CARGAS
ESTADO 1 "Fuerza F con dirección alpha en nudo 4"
   CARGAS EN NUDOS
   IF,comprob,=,0,THEN
      4   F*cos(alfa)    -F*sin(alfa)   0.0 
   ELSEIF,comprob,=,2,THEN  
      4   Fx              Fy            0.0 
   ENDIF
   
   CARGAS EN ELEMENTOS
      FUERZAS DE EMPOTRAMIENTO
         12   E*Ax/(2*sqrt(2))*al 0 0 0 0 0 -E*Ax/(2*sqrt(2))*al 0 0 0 0 0

*FIN