Downloads: 1
Diff of
/Knee Extension Testing/props_PatellarTendon.inp
[000000]
..
[c02310]
Switch to unified view
| a | b/Knee Extension Testing/props_PatellarTendon.inp | ||
|---|---|---|---|
| 1 | *MEMBRANE SECTION, ELSET=PAT_LIG-2D, MATERIAL=TENDON_PAT_MAT |
||
| 2 | 1.0, |
||
| 3 | *MATERIAL, NAME=TENDON_PAT_MAT |
||
| 4 | ** material definitions have STRESS (MPa) vs STRAIN (%) |
||
| 5 | *HYPERELASTIC, POISSON=0.498, TEST DATA INPUT, VAN DER WAALS, BETA=0 |
||
| 6 | ** |
||
| 7 | *UNIAXIAL TEST DATA |
||
| 8 | 0.0062, 0.07 |
||
| 9 | 0.0144, 0.15 |
||
| 10 | 0.0288, 0.25 |
||
| 11 | 0.0504, 0.36 |
||
| 12 | 0.0728, 0.48 |
||
| 13 | 0.1003, 0.67 |
||
| 14 | 0.1185, 0.79 |
||
| 15 | *DENSITY |
||
| 16 | 5E-8, |
||
| 17 | ** |
||
| 18 | *DAMPING, ALPHA=100. |
||
| 19 | ** |
||
| 20 | *MASS, ELSET=PAT_LIG_MASS |
||
| 21 | 1E-5 |
||
| 22 | **************PL DEFINITION*********************************** |
||
| 23 | *PARAMETER |
||
| 24 | ** PL_L=the nominal FIBERament length(mm) |
||
| 25 | PL_L=61.6 |
||
| 26 | ** PL_eref=the reference strain (+)=added slack, (-)=initial force |
||
| 27 | PL_eref=0.0 |
||
| 28 | ** PL_k=the nominal stiffness (N/unit strain); for (N/mm) divide FIBER_k by FIBER_L |
||
| 29 | PL_k = 50000.0 |
||
| 30 | ** PL_num_sp_par=number of springs in parallel representing whole FIBERament |
||
| 31 | PL_num_sp_par = 6.0 |
||
| 32 | PL_num_sp_ser = 16.0 |
||
| 33 | ** |
||
| 34 | k_ratio=50000.0/PL_k |
||
| 35 | PL_sp_L = PL_L/PL_num_sp_ser |
||
| 36 | ** |
||
| 37 | *PARAMETER |
||
| 38 | PL_Lref=PL_sp_L |
||
| 39 | PL_L0=PL_Lref/(PL_eref+1) |
||
| 40 | ** |
||
| 41 | PL_d0=(PL_Lref-PL_L0)/PL_num_sp_ser |
||
| 42 | PL_d2=((PL_Lref+0.02*PL_Lref)-PL_L0)/PL_num_sp_ser |
||
| 43 | PL_d4=((PL_Lref+0.04*PL_Lref)-PL_L0)/PL_num_sp_ser |
||
| 44 | PL_d6=((PL_Lref+0.06*PL_Lref)-PL_L0)/PL_num_sp_ser |
||
| 45 | PL_d10=((PL_Lref+0.1*PL_Lref)-PL_L0)/PL_num_sp_ser |
||
| 46 | PL_d20=((PL_Lref+0.2*PL_Lref)-PL_L0)/PL_num_sp_ser |
||
| 47 | PL_d100=((PL_Lref+1.0*PL_Lref)-PL_L0)/PL_num_sp_ser |
||
| 48 | PL_d200=((PL_Lref+2.0*PL_Lref)-PL_L0)/PL_num_sp_ser |
||
| 49 | ** |
||
| 50 | PL_f0=0.0*k_ratio |
||
| 51 | PL_f2=420.0/(PL_num_sp_par*PL_num_sp_ser*k_ratio) |
||
| 52 | PL_f4=830.0/(PL_num_sp_par*PL_num_sp_ser*k_ratio) |
||
| 53 | PL_f6=1500.0/(PL_num_sp_par*PL_num_sp_ser*k_ratio) |
||
| 54 | PL_f10=3600.0/(PL_num_sp_par*PL_num_sp_ser*k_ratio) |
||
| 55 | PL_f20=9000.0/(PL_num_sp_par*PL_num_sp_ser*k_ratio) |
||
| 56 | PL_f100=50000.0/(PL_num_sp_par*PL_num_sp_ser*k_ratio) |
||
| 57 | PL_f200=100000.0/(PL_num_sp_par*PL_num_sp_ser*k_ratio) |
||
| 58 | ** |
||
| 59 | *SPRING, NONLINEAR, ELSET=PAT_LIG-SP |
||
| 60 | |||
| 61 | 0.0 , -<PL_d100> |
||
| 62 | <PL_f0> , <PL_d0> |
||
| 63 | <PL_f2> , <PL_d2> |
||
| 64 | <PL_f4> , <PL_d4> |
||
| 65 | <PL_f6> , <PL_d6> |
||
| 66 | <PL_f10> , <PL_d10> |
||
| 67 | <PL_f20> , <PL_d20> |
||
| 68 | <PL_f100> , <PL_d100> |
||
| 69 | <PL_f200> , <PL_d200> |
||
| 70 | ** |
||
| 71 | *SURFACE, NAME=SURF-PAT_LIG, MAX RATIO=0.8 |
||
| 72 | PAT_LIG-2D, SPOS |
||
| 73 | ** |
||
| 74 | *NODE |
||
| 75 | 290, 25, 15, -50 |
||
| 76 | 291, 25, 15, -53 |
||
| 77 | ** |
||
| 78 | *ORIENTATION, DEFINITION=NODES, NAME=PAT_TEND_ORI |
||
| 79 | 290, 2900000, 291 |
||
| 80 | *ELEMENT, TYPE=CONN3D2, ELSET=PAT_TEND_PULL |
||
| 81 | 290, 291, 290 |
||
| 82 | ** |
||
| 83 | *CONNECTOR SECTION, ELSET=PAT_TEND_PULL |
||
| 84 | TRANSLATOR, |
||
| 85 | PAT_TEND_ORI, |
||
| 86 | ** |
||
| 87 | *MPC |
||
| 88 | BEAM, TEND_PAT_ATTACH, 3900000 |
||
| 89 | BEAM, TEND_TIB_ATTACH, 290 |
||
| 90 | BEAM, 291, 2900000 |