吧友们我用python下载openseespy，根据PEER的数据库建了一个很简单的墩柱模型 粘结滑移效应什么都没考虑直接fix底端 但是拟合数据库的拟静力试验一直不收敛，尝试了很多方法改善收敛性但都没办法跑完全程，改善方法包括但不限于：
1.将nonlinearBeamColumn单元换成dispBeamColumn
2.将算法从Newton改为KrylovNewton
3.减小步长（该方法在0位移处荷载反向时起作用，但在后续加载中不知为何再次报错）
4.降低容差到10^-5 m4.将积分点从2个依次增加至5个
5.划分墩柱单元至5个
6.增加截面纤维数量
上述方法都会出现 -3 error报错
我在youtube上观看了一位youtuber的pushover分析在重力和轴向荷载分析完成后，在循环荷载的pattern下多加了sp指令指定墩顶强制位移，但是出现了 -2 error报错
我自己写的代码如下，包括单独写了一份根据试验测量位移加载位移的函数：




循环加载的函数如下：
loadlist是试验测量的位移列表

基础模型跑不通的话后续模型更正和材料本构修改都没法进行
希望吧友大神能够给予一些指导 谢谢各位