等几何分析，壳结构非线性问题，matlab代码，现阶段已有线性部分，需要加入非线性内容及压电层合内容