好久没用过MATLAB了，最近发现MATLAB的界面都和以前的不一样了。于是装个2014版本来学习一下，顺便复习一下有限元。
用MATLAB做个最简单的2D桁架分析并与SAP2000对比。
桁架信息：
弹性模量 E ：200000
截面面积 A ：4535
MATLAB计算结果：
SAP2000计算结果：
结果对比：
| DOF | MATLAB | SAP2000 |
| 1 | 0 | 0 |
| 2 | 0 | 0 |
| 3 | 0.2481 | 0.248071 |
| 4 | -1.586 | -1.586026 |
| 5 | 0.6615 | 0.661521 |
| 6 | -1.586 | -1.586026 |
| 7 | 0.9096 | 0.909592 |
| 8 | 0 | 0 |
| 9 | 0.7856 | 0.785557 |
| 10 | -1.0862 | -1.086157 |
| 11 | 0.4548 | 0.454796 |
| 12 | -1.9205 | -1.920515 |
| 13 | 0.124 | 0.124035 |
| 14 | -1.0862 | -1.086157 |
从上面的节点位移对比可知，MATLAB和SAP2000计算结果是一致的。
时间冲忙，就整理到这里吧，发现就这么简单的整理一下，备忘一下都会花很多时间。
PS. 希望后面可以继续有时间和激情整理MATLAB用于有限元的一些小东西，希望可以整理成资料放上来和大家分享，对于初学者就很有用了。
- 注释 ( Comments )
源代码已经收录到新书《有限单元法：编程与软件应用》，详见链接：http://www.jdcui.com/?page_id=9731
- 微信公众号 ( Wechat Subscription)
欢迎关注 “结构之旅” 微信公众号
崔老师，您好，2D Truss我按照书上的代码在MATLAB里面写了，可以运行，结果也一样,但好像有些问题，就是总的刚度矩阵不满足奇异性，也就是行列式不为0，
det(stiffness)
ans =
12.1527e+027
我感觉这好像不满足书里的奇异性，望崔老师解答疑问，我在3D4面体单元TET4运行时虽然满足总刚行列式为0，求解位移结果出现了
Displacement
警告: 矩阵为奇异工作精度。
> In D3_cantilever_TET4 (line 111)
我不清楚什么原因
约束条件加上去没。
第一个代码没调试成功
代码不开源？可以每篇文章附带代码吗，这样便于学习哈
都收录在即将出本的书本上了：http://www.jdcui.com/?page_id=9731，书本上都开源了。到时候有兴趣可以看看书。
崔博士你好，能不能麻烦你把matlab 3d桁架有限元的程序发我一下，我最近正在学习matlab的有限元编程 😉