! Procedimiento: 
   ! Para comprobar que lo obtenido en el libro está correcto tomamos
   ! 	comprob = 0
   ! Para determinar con el programa la fuerza aplicada tomamos
   ! 	comprob = 1. Las componentes de la reacción en 4 serán 
   !    los valores de Fx y Fy respectivamente.
   !Se comprueba que esas componentes hacen nulo el desplazamiento 
   !	de 4 y se obtienen el módulo y el ángulo de la fuerza.  
! ===========================================================

*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    
   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 del libro   
   F   = 87.94       !Fuerza aplicada en 4
   alfa  =  0.4597   !Dirección (rad) de F   
   
 
   Fx  =   172.245
   Fy  =  -136.22
   ang = atan(-Fy/Fx)      ! angulo alfa obtenido
   mod = sqrt(Fx^2+Fy^2)   ! F obtenido

   comprob = 2   ! 0 = comprobamos los parámetros del libro
                 ! 1 = Obtenemos Fx y Fy
                 ! 2 = comprobamos que está bien

!________________________________
!Parámetros necesarios para el cálculo
   car1 =  (2*nm+3*(int(1/3*nm)+1)) !inicio de fuerza de empotramiento
   car2 =  (2*nm+3*(int(2/3*nm)+1)) !fin de fuerza de empotramiento
   gen = int((car2-car1)/3)                  !barras que genera

*PARAMETROS DE ANALISIS

>COORDENADAS
! nudo      X      Y      Z
   1       0.0    0.0     0.0	L	nm+1	1	lm	0.0	0.0	    
   nm+2    nm*lm   h      0.0	L	nm+1	1     -lm	0.0	0.0                                      
  
>MATERIALES
! kM        E     G      nu    Pesp   Alpha     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   2*nm+1   1   1   1
   2*nm+2  1  2*nm+2  1 1   0  0  
   2*nm+3  1  2*nm+1  1 1   0  0   L   nm-1    3   1  -1
   2*nm+4  2  2*nm+1  1 1   0  0   L   nm-1    3   1  -1
   2*nm+5  3  2*nm+1  1 1   0  0   L   nm-1    3   1  -1
      

            
>RESTRICCIONES
  TODOS                  DZ  GX  GY  GZ
   1              DX  DY
   int(2/3*nm)+1      DY
 IF,comprob,=,1,THEN
   nm+1   DX  DY
 ELSE
 ENDIF 

>CARGAS
ESTADO 1 "Fuerza F con dirección alpha en nudo 4"
       CARGAS EN NUDOS
    IF,comprob,=,0,THEN
          nm+1   F*cos(alfa)    -F*sin(alfa)   0.0    0.0   0.0   0.0 
      
    ELSEIF,comprob,=,2,THEN  
          nm+1   Fx    Fy   0.0    0.0   0.0   0.0 
    ELSE
    ENDIF

       CARGAS EN ELEMENTOS
        Fuerzas de empotramiento
        car1 E*Ax/(2*sqrt(2))*al 0 0 0 0 0 -E*Ax/(2*sqrt(2))*al 0 0 0 0 0 L gen 3 
        
*FIN