! Práctica: Cálculo de la matriz de rigidez de una estructura 
! ===========================================================
! Procedimiento: 
   ! Se restringen todos los gdl excepto uno
   ! Se aplica una fuerza (momento) Fu (Mu) en el gdl no restringido      
   ! Se calcula la estructura y se obtiene el desplazamiento
   ! Las rigideces en cada gdl son Fj/Di (Mj/Gi)
   !    siendo Fj,Mj las fuerzas aplicadas/reacciones en los gdl
   !    y Di, Dj los desplazamientos/giros en el gdl libre.  

*CONTROL DEL PROBLEMA
   IDPR    P1-CMRF
   TITU   "Cálculo de la matriz de rigidez de la estructura" 
   SISUNI  mkN
   TIPEST, EP_NR_XY

*PARAMETROS GENERALES
   iEstado = 4
   Fu = 1.0e+3
   Mu = 1.0e+3

*PARAMETROS DE ANALISIS
>COORDENADAS
! kN        X      Y      Z
   1      0.0    0.0    0.0
   2      0.0    3.0    0.0
   3      3.0    3.0    0.0
   4      8.0    0.0    0.0

>RESTRICCIONES
IF,iEstado,=,1,THEN
   1   DX DY DZ GX GY GZ
   2      DY DZ GX GY GZ
   3   DX DY DZ GX GY GZ
   4   DX DY DZ GX GY GZ
ELSEIF,iEstado,=,2,THEN
   1   DX DY DZ GX GY GZ
   2   DX    DZ GX GY GZ
   3   DX DY DZ GX GY GZ
   4   DX DY DZ GX GY GZ
ELSEIF,iEstado,=,3,THEN	
   1   DX DY DZ GX GY GZ
   2   DX DY DZ GX GY 
   3   DX DY DZ GX GY GZ
   4   DX DY DZ GX GY GZ 
ELSEIF,iEstado,=,4,THEN	
   1   DX DY DZ GX GY GZ
   2   DX DY DZ GX GY GZ
   3      DY DZ GX GY GZ
   4   DX DY DZ GX GY GZ 
ELSEIF,iEstado,=,5,THEN
   1   DX DY DZ GX GY GZ
   2   DX DY DZ GX GY GZ
   3   DX    DZ GX GY GZ
   4   DX DY DZ GX GY GZ
ELSEIF,iEstado,=,6,THEN	
   1   DX DY DZ GX GY GZ
   2   DX DY DZ GX GY GZ
   3   DX DY DZ GX GY 
   4   DX DY DZ GX GY GZ
ELSEIF,iEstado,=,7,THEN	
   1   DX DY DZ GX GY GZ
   2   DX DY DZ GX GY GZ
   3   DX DY DZ GX GY GZ
   4   DX DY DZ GX GY 
ENDIF

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

>PROPIEDADES
! kP        Ax    Ay      Az       Ix       Iy          Iz
   1   3000E-6   0.0     0.0      0.0      0.0     8000E-8
   2   5000E-6   0.0     0.0      0.0      0.0    16000E-8

>ELEMENTOS
   GRUPO 1 BNR
      1    1   2    1  1     0  0  CA    -1.0      0.0    0.0
      2    2   3    1  2     0  0  CA     0.0      5.0    0.0
      3    3   4    1  2     0  0  CA     4.0     10.0    0.0
    
>CARGAS
   Estado 1 "Fuerza unitaria gdl 1"
      CARGAS EN NUDOS        
      2  Fu  0.0  0.0  0.0  0.0  0.0  
   Estado 2 "Fuerza unitaria gdl 2"
      CARGAS EN NUDOS         
      2  0.0  Fu  0.0  0.0  0.0  0.0
   Estado 3 "Fuerza unitaria gdl 3"
      CARGAS EN NUDOS         
      2   0.0  0.0  0.0  0.0  0.0  Mu
   Estado 4 "Fuerza unitaria gdl 4"
      CARGAS EN NUDOS         
      3   Fu  0.0  0.0  0.0  0.0  0.0
   Estado 5 "Fuerza unitaria gdl 5"
      CARGAS EN NUDOS         
      3   0.0  Fu  0.0  0.0  0.0  0.0
   Estado 6 "Fuerza unitaria gdl 6"
      CARGAS EN NUDOS         
      3   0.0  0.0  0.0  0.0  0.0  Mu
   Estado 7 "Fuerza unitaria gdl 7"
      CARGAS EN NUDOS         
      4   0.0  0.0  0.0  0.0  0.0  Mu
*FIN
