78 lines (77 with data), 2.2 kB

*MEMBRANE SECTION, ELSET=RECT_FEM-2D, MATERIAL=TENDON_RECT_FEM_MAT
1.0,
*MATERIAL, NAME=TENDON_RECT_FEM_MAT
** material definitions have STRESS (MPa) vs STRAIN (%)
*HYPERELASTIC, POISSON=0.498, TEST DATA INPUT, VAN DER WAALS, BETA=0
**
*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-8,
**
*DAMPING, ALPHA=100.
**
*MASS, ELSET=RECT_FEM_MASS
1E-7
**************RF DEFINITION***********************************
*PARAMETER
** RF_L=the nominal RFament length(mm)
RF_L=60.0
** RF_eref=the reference strain (+)=added slack, (-)=initial force
RF_eref=0.0
** RF_k=the nominal stiffness (N/unit strain); for (N/mm) divide RF_k by RF_L
RF_k = 50000.0
** RF_num_sp_par=number of springs in parallel representing whole RFament
RF_num_sp_par = 5.0
RF_num_sp_ser = 15.0
***************************
k_ratio=50000.0/RF_k
RF_sp_L = RF_L/RF_num_sp_ser
**
*PARAMETER
RF_Lref=RF_sp_L
RF_L0=RF_Lref/(RF_eref+1)
**
RF_d0=(RF_Lref-RF_L0)/RF_num_sp_ser
RF_d2=((RF_Lref+0.02*RF_Lref)-RF_L0)/RF_num_sp_ser
RF_d4=((RF_Lref+0.04*RF_Lref)-RF_L0)/RF_num_sp_ser
RF_d6=((RF_Lref+0.06*RF_Lref)-RF_L0)/RF_num_sp_ser
RF_d10=((RF_Lref+0.1*RF_Lref)-RF_L0)/RF_num_sp_ser
RF_d20=((RF_Lref+0.2*RF_Lref)-RF_L0)/RF_num_sp_ser
RF_d100=((RF_Lref+1.0*RF_Lref)-RF_L0)/RF_num_sp_ser
RF_d200=((RF_Lref+2.0*RF_Lref)-RF_L0)/RF_num_sp_ser
**
RF_f0=0.0*k_ratio
RF_f2=420.0/(RF_num_sp_par*RF_num_sp_ser*k_ratio)
RF_f4=830.0/(RF_num_sp_par*RF_num_sp_ser*k_ratio)
RF_f6=1500.0/(RF_num_sp_par*RF_num_sp_ser*k_ratio)
RF_f10=3600.0/(RF_num_sp_par*RF_num_sp_ser*k_ratio)
RF_f20=9000.0/(RF_num_sp_par*RF_num_sp_ser*k_ratio)
RF_f100=50000.0/(RF_num_sp_par*RF_num_sp_ser*k_ratio)
RF_f200=75000.0/(RF_num_sp_par*RF_num_sp_ser*k_ratio)
**
*SPRING, NONLINEAR, ELSET=RECT_FEM_SP

0.0 , -<RF_d100>
<RF_f0> , <RF_d0>
<RF_f2> , <RF_d2>
<RF_f4> , <RF_d4>
<RF_f6> , <RF_d6>
<RF_f10> , <RF_d10>
<RF_f20> , <RF_d20>
<RF_f100> , <RF_d100>
<RF_f200> , <RF_d200>
**
**
**
*SURFACE, NAME=SURF-RECT_FEM-POS, MAX RATIO=0.8
RECT_FEM-2D, SPOS
*SURFACE, NAME=SURF-RECT_FEM-NEG, MAX RATIO=0.8
RECT_FEM-2D, SNEG
*MPC
BEAM, RF_PAT_ATTACH, 3900000