matlab接触求解器开发

运用matlab做一个接触求解器（罚函数法，拉格朗日法均可但要有较为详细的程序说明），能够用平面单元计算一个圆筒与平板在接触过程中的接触刚度，过程中可能涉及到几何非线性与接触非线性，我自己目前是有一个求解器但是计算结果只能在接触力较小的时候准确，大接触力下得不到正确结果，因此希望所做求解器能够同时计算大接触力和小接触力的情况