102 lines (100 with data), 2.7 kB

*MEMBRANE SECTION, ELSET=PFL-LAT-2D, MATERIAL=LIG_PFL-LAT_MAT
1.0,
*MATERIAL, NAME=LIG_PFL-LAT_MAT
** material definitions have STRESS (MPa) vs STRAIN (%)
*HYPERELASTIC, POISSON=0.45, TEST DATA INPUT, VAN DER WAALS, BETA=1
**
*UNIAXIAL TEST DATA
0.0062,  0.07
0.0144,  0.15
0.0288, 0.25
0.0504, 0.36
0.0728, 0.48
0.1003, 0.67
0.1185, 0.79
*DENSITY			
5E-7,
**
*DAMPING, ALPHA=200.
**
*MASS, ELSET=PFL-LAT_MASS
1E-5
**************PFL DEFINITION***********************************
*PARAMETER
** LPRL_L=the nominal LPRLament length(mm)
LPRL_L=50.0
** LPRL_eref=the reference strain (+)=added slack, (-)=initial force
LPRL_eref= -0.05
** LPRL_k=the nominal stiffness (N/unit strain); for (N/mm) divide LPRL_k by LPRL_L
LPRL_k = 1000.0
** LPRL_num_sp_par=number of springs in parallel representing whole LPRLament
LPRL_num_sp_par = 5.0
LPRL_num_sp_ser = 18.0
***************************
**
*PARAMETER
**
k_ratio=1000.0/LPRL_k
LPRL_sp_L = LPRL_L/LPRL_num_sp_ser
LPRL_Lref=LPRL_sp_L
LPRL_L0=LPRL_Lref/(LPRL_eref+1)
**
LPRL_d0=(LPRL_Lref-LPRL_L0)/LPRL_num_sp_ser
LPRL_d2=((LPRL_Lref+0.02*LPRL_Lref)-LPRL_L0)/LPRL_num_sp_ser
LPRL_d4=((LPRL_Lref+0.04*LPRL_Lref)-LPRL_L0)/LPRL_num_sp_ser
LPRL_d6=((LPRL_Lref+0.06*LPRL_Lref)-LPRL_L0)/LPRL_num_sp_ser
LPRL_d10=((LPRL_Lref+0.1*LPRL_Lref)-LPRL_L0)/LPRL_num_sp_ser
LPRL_d20=((LPRL_Lref+0.2*LPRL_Lref)-LPRL_L0)/LPRL_num_sp_ser
LPRL_d100=((LPRL_Lref+1.0*LPRL_Lref)-LPRL_L0)/LPRL_num_sp_ser
**
LPRL_f0=0.0*k_ratio
LPRL_f2=3.0/(LPRL_num_sp_par*LPRL_num_sp_ser*k_ratio)
LPRL_f4=12.0/(LPRL_num_sp_par*LPRL_num_sp_ser*k_ratio)
LPRL_f6=30.0/(LPRL_num_sp_par*LPRL_num_sp_ser*k_ratio)
LPRL_f10=70.0/(LPRL_num_sp_par*LPRL_num_sp_ser*k_ratio)
LPRL_f20=180.0/(LPRL_num_sp_par*LPRL_num_sp_ser*k_ratio)
LPRL_f100=1000.0/(LPRL_num_sp_par*LPRL_num_sp_ser*k_ratio)
**
*SPRING, NONLINEAR, ELSET=PFL-LAT-SP

0.0 , -<LPRL_d100>
<LPRL_f0> , <LPRL_d0>
<LPRL_f2> , <LPRL_d2>
<LPRL_f4> , <LPRL_d4>
<LPRL_f6> , <LPRL_d6>
<LPRL_f10> , <LPRL_d10>
<LPRL_f20> , <LPRL_d20>
<LPRL_f100> , <LPRL_d100>
**
**
*SURFACE, NAME=SURF-PFL-LAT, MAX RATIO=0.8
PFL-LAT-2D
**
*CONNECTOR SECTION, ELSET=MUSC-PFL-LAT
TRANSLATOR
LPFL_CONN,
*ORIENTATION, NAME=LPFL_CONN, DEFINITION=NODES
90200, 1900000, 90300
**
*NODE
11100001, 0.0, 0.0, 0.0
**
**
*MPC
BEAM,   1090018,     90300
BEAM,   1090086,     90300
BEAM,   1090016,     90300
BEAM,   1090087,     90300
BEAM,   1090015,     90300
**
**
*MPC
BEAM,     90200,    1900000
**
** BEAM ANTERIOR NODES TO PAT RBRN
*MPC
BEAM, 1090000, 3900000
BEAM, 1090001, 3900000
BEAM, 1090002, 3900000
BEAM, 1090003, 3900000
BEAM, 1090004, 3900000
**