非线性壳编程，弧长法编程
可以实现和复现文章内容最好，目前的程序只做了单层壳单元的内容，缺少多层壳和弧长法后屈曲路径