*CONTROL DEL PROBLEMA
   IDPR,   Ej13.1
   TITU,   "Ejercicio 13.1: Estructura simétrica con carga puntual asimétrica"
   SISUNI, mkN
   TIPEST, EP_NR_XY

*PARAMETROS GENERALES
   Modelo = 2   ! 0 <> Estructura completa
                ! 1 <> Media estructura con carga simétrica
                ! 2 <> Media estructura con carga antisimétrica
   H   = 4.0
   L1  = 4.0
   L2  = 4.0
   F   = 40.0

   Ar  = 20.0     ! para obtener errores en deformaciones axiales pequeños
	
*PARAMETROS DE ANALISIS

>COORDENADAS
! nudo      X      Y        Z   (metros)
   1      0.0     0.0     0.0
   2       L1     0.0     0.0
   4      0.0      H      0.0
IF,Modelo,=,0,THEN
   5       L1       H     0.0
   3     L1+L2     0.0    0.0
   6     L1+L2      H     0.0
ELSE
   3       L1      H      0.0
ENDIF
    
>RESTRICCIONES
   TODOS         DZ  GX GY
   1      DX DY            GZ
   2      DX DY            GZ  ! puedo quitar esta restricción en el caso simétrico, pero quizá también interesa ver que la barra no trabaja.
IF,Modelo,=,0,THEN
   3      DX DY            GZ
ELSEIF,Modelo,=,1,THEN
   3      DX  DY           GZ  ! porque en el libro no permite alargamiento en barras no??   
ELSEIF,Modelo,=,2,THEN
ENDIF

>MATERIALES
! kM        E     G     Nu    Rho    Alpha        Se
   1  210.0e6   0.0    0.3   77.01     0.0   260.0e3


>PROPIEDADES
! kP    A    Ay    Az    Ix    Iy        Iz     
   1   Ar   0.0   0.0   1.0   1.0   5000.e-8     
   2   Ar   0.0   0.0   1.0   1.0  10000.e-8    

>ELEMENTOS
   GRUPO 1 BNR
   ! kE   kI  kJ    kM  kP	
      1    1   4     1   1    0   0
 
IF,Modelo,=,0,THEN
      2    4   5     1   2    0   0
      3    2   5     1   2    0   0 
      4    5   6     1   2    0   0
      5    3   6     1   1    0   0
ELSE
      2    4   3     1   2    0   0
      3    2   3     1   1    0   0
ENDIF  
     
>CARGAS
   ESTADO 1 "Cargas en nudos (kN)"
      CARGAS EN ELEMENTOS
      IF,Modelo,=,0,THEN
         PUNTUAL H/2
         1   0.0   -F   0.0   0.0   0.0   0.0
      ELSE
         PUNTUAL H/2
         1   0.0   -F/2   0.0   0.0   0.0   0.0

      ENDIF

*FIN



