用Newton-Raphson法求解端部位移及单元应力，绘制力-位移曲线。测试F = 30kN时会发生什么，解释原因。
问题遇到的现象和发生背景
问题相关代码，请勿粘贴截图
运行结果及报错内容
我的解答思路和尝试过的方法
我想要达到的结果
牛顿-拉夫逊方法（Newton-Raphson method）
https://blog.csdn.net/xinm1001/article/details/52938113
牛顿-拉夫逊方法（Newton-Raphson method）
https://blog.csdn.net/xinm1001/article/details/52938113