*CONTROL DEL PROBLEMA
   IDPR,   Ej-3.11
   TITU,   "Ejercicio 3.11"
   SISUNI, mkN
   TIPEST, EP_NR_XY

*PARAMETROS GENERALES
   Rad   = 15
   NumP  = 49		! (mínimo 5)
   AngFv = 30		! (grados) 
   AngFg = 60		! (y 180º -60º) (grados)
   Cv    = 0.7
   Cg    = 1
   Pv    = 0.7*Rad*8
   Pg    = 1*2*Rad*8   
   pi    = 3.14159265
   Alfa  = 180/NumP
   NFv   = int(AngFv/Alfa+1)
   NFg   = int(AngFg/Alfa+1)

*PARAMETROS DE ANALISIS

>COORDENADAS
!   kN                X                                Y                    Z  
   DO,J,0,NumP-1
     J+1  , Rad*(1-cos(pi*(J/(NumP-1)))),  Rad*(sin(pi*(J/(NumP-1)))) ,    0.0
   ENDDO
   DO,J,0,NumP-1 
     J+1  , Rad*(1-cos(pi*(J/(NumP-1)))),  Rad*(sin(pi*(J/(NumP-1)))) ,    0.0
   ENDDO
 
>RESTRICCIONES
 !TODOS              DZ   GX    GY
   1       DX   DY   DZ   GX    GY       
   NumP    DX   DY   DZ   GX    GY 

>MATERIALES
! kM    E         G    nu      Rho     Alpha      SigmaE
   1  210.0e6   0.0   0.3    77.01   12.0e-6     275.0e3

>PROPIEDADES
!  kP  Familia  Perfil   Angulo
   1   HEB      400	    90.0

>ELEMENTOS
   GRUPO 1 BNR
      DO I,1,int(NumP/2)-1
      I I I+1 1 1  0   0 CA   Rad Rad/2 0
      ENDDO
      DO S,int(NumP/2)+2,NumP-1
      S S S+1 1 1  0   0 CA   Rad Rad/2 0
      ENDDO
      int(NumP/2)   , int(NumP/2)   ,int(NumP/2)+1,  1  1  0   0 CA   Rad Rad/2 0 LIBE  0 6
      int(NumP/2)+1 , int(NumP/2)+1 ,int(NumP/2)+2,  1  1  0   0 CA   Rad Rad/2 0 LIBE  0 0

>CARGAS
   ESTADO 1 "Cargas de viento"
      CARGAS EN NUDOS
      NFv     Pv   0.0   0.0   0.0   0.0   0     
   
   ESTADO 2 "Cargas gravitatorias"
      CARGAS EN NUDOS
      NFg           0.0  -Pg/2    0.0   0.0   0.0   0
      NumP-NFg+1    0.0  -Pg/2    0.0   0.0   0.0   0

>COMBINACIONES
   ESTADO 3 "1.0 * Cargas de viento + 1.0 * cargas gravitatorias"
      1  1.0   2  1.0

 *FIN DE DATOS
