求教用中心差分法求解非线性的问题

xsdlmpwbs

我用中心差分法编的程序，但是进行结构计算时出现了些问题：
1、在刚度不变时(应力应变曲线成直线），改变刚度的大小，隐式和显式积分对结构均无影响，皆收敛
2、在刚度改变时（应力应变曲线成折线），隐式积分影响很小，收敛；但是中心差分法找不到合适的刚度值使其收敛，我做的数据一直震荡发散。
问各位高手，非线性结构对显式积分是不是特别敏感啊，如何找到结构的最优刚度呢，求解释