Closed dashdotrobot closed 5 years ago
Instead of hunting for singular matrix configurations near T_c_theor, the buckling tension should be obtained from the generalized eigenvalue problem:
K_rim_matl + K_spk_matl + T_bar*(K_spk_matl - K_rim_matl) = 0
This should be considerably more robust than iteratively tuning T_bar until the condition number is maximized.
Instead of hunting for singular matrix configurations near T_c_theor, the buckling tension should be obtained from the generalized eigenvalue problem:
K_rim_matl + K_spk_matl + T_bar*(K_spk_matl - K_rim_matl) = 0
This should be considerably more robust than iteratively tuning T_bar until the condition number is maximized.