塑性材料的有限元模拟你是不是也搞过？矩阵位移法的 MATLAB 例程，我最近发现了一个蛮实用的，叫plasticFEM，主要用来做塑性，特点是结构清晰、注释详尽，适合用来当做二次开发的起点。
runFEM.m和runReturnMap.m两个脚本负责整体流程，其他五个函数各司其职，像bmatHex8.m计算的是八节点六面体的 B 矩阵，returnMap.m负责弹塑性应力更新。嗯，虽然只支持位移控制（不施加力），但对于理解材料非线性挺友好的。
代码风格也不错，不是那种让人看不懂的“科研风”，每个函数开头都有清楚的输入输出，比如globTanStiff.m里同时算了刚度矩阵、内力、应力应变，一步到位，效率还挺高。要是你在用 MATLAB 搞有限元，不妨试试看。
参考文献也给得实在，像 Simo 和 Hughes 那本《Computational Inelasticity》，还有金的新手向书籍，都对理解背景模型蛮有的。
如果你想了解其他类似资源，可以看看SlopeFEA这个项目，或者简易有限元程序实现这种入门向代码，结合着学，效率会更高哦～