(业余时间 和小伙伴一起写写有限元程序 ）
- 程序作者 ( Author )
JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2
1.广州容柏生建筑结构设计事务所；2.华南理工大学建筑设计研究院
- 基本概念 ( Concept )
8节点六面体线性“完全积分”单元，用Abaqus中的命名规则，即 C3D8。
- 问题描述( Problem Description)
悬臂梁，悬臂长度2.0m，梁高0.5m，梁宽0.2m。梁左端嵌固，受重力作用。材料弹性模量E=200000MPa，材料泊松比为0.3。
基于MATLAB编程实现该悬臂梁的弹性静力分析，采用C3D8单元进行模拟，并将MATLAB的计算结果与SAP2000、midas Gen分析结果进行对比。
- MATLAB 编程
disp
stress
- SAP2000
model
disp
stress
单元应力对比
| 单元编号 | 节点号 | SAP2000结果 | MATLAB结果 | 相对偏差（%） |
| 1 | 1 | -1.7462 | -1.7462 | 0 |
| 5 | -1.3304 | -1.3304 | 0 | |
| 6 | 0.128 | 0.128 | 0 | |
| 7 | 0.00E+00 | 0.00E+00 | 0 | |
| 28 | -1.7462 | -1.7462 | 0 | |
| 29 | -1.3304 | -1.3304 | 0 | |
| 31 | 0.128 | 0.128 | 0 | |
| 30 | 0.00E+00 | 0.00E+00 | 0 | |
| 7 | 12 | -0.5314 | -0.5314 | 0 |
| 15 | -0.596 | -0.596 | 0 | |
| 16 | 0.0348 | 0.0348 | 0 | |
| 13 | 0.0547 | 0.0547 | 0 | |
| 37 | -0.5314 | -0.5314 | 0 | |
| 40 | -0.596 | -0.596 | 0 | |
| 41 | 0.0348 | 0.0348 | 0 | |
| 38 | 0.0547 | 0.0547 | 0 | |
| 15 | 24 | -0.0105 | -0.0105 | 0 |
| 2 | -0.0173 | -0.0173 | 0 | |
| 27 | 5.55E-04 | 5.55E-04 | 0 | |
| 25 | 0.0027 | 0.0027 | 0 | |
| 49 | -0.0105 | -0.0105 | 0 | |
| 52 | -0.0173 | -0.0173 | 0 | |
| 53 | 5.55E-04 | 5.55E-04 | 0 | |
| 50 | 0.0027 | 0.0027 | 0 |
可以看出，MATLAB编程计算结果是和 SAP2000 及 Midas Gen的计算结果是完全一致的。
- 注释 ( Comments )
( 如果您发现有错误，欢迎批评指正。邮箱：jidong_cui@163.com . 如果您喜欢这篇博文，请在上面给我 点个赞 吧！ 🙂 🙂
( If you found any mistakes in the post, please let me know. Email : jidong_cui@163.com. If you like this posts, please give me a “thumbs up“ rating on the above button! )
- 微信公众号 ( Wechat Subscription)
欢迎关注 “结构之旅” 微信公众号