Customized RMR solver with glenohumeral stability formulations

[imihasan]

In this customized version, we proposed different formulations to model glenohumeral stability using:

• Inequality constraints, where the direction of the glenohumeral contact force was constrained within different stability border approximations:
• A point
• An ellipse
• A polynomial fit of empirical values obtained from concavity compression tests

Penalty terms in the objective function as:

• Conditional penalty: that triggers if the direction of the contact force exceeds a specified circular stability border
• Planar penalty: increases as the direction of the glenohumeral contact force points further away from the glenoid cavity border. Here, the glenoid cavity is considered as a flat plane. This plane is the plane that contains the glenoid cavity border.
• Curve penalty: increases as the direction of the glenohumeral contact force points further away from the glenoid cavity border. Here, the curvature of the glenoid cavity is approximated as spherical.