[FEM][有限元][编程][Matlab][Code by myself] 平面8节点二次“完全积分”单元(CPS8)

  (业余时间 和小伙伴一起写写有限元程序 ）

• 程序作者 ( Author )

JiDong Cui (崔济东) ¹, XueLong Shen (沈雪龙)²

1.广州容柏生建筑结构设计事务所； 2.华南理工大学建筑设计研究院

• 基本概念 ( Concept )

平面8节点二次“完全积分”单元，用Abaqus中的命名规则，即CPS8。

• 问题描述( Problem Description)

XZ平面内的悬臂梁，悬臂长度2.0m，梁高0.5m，梁宽0.2m。梁左端嵌固，右端受到-z方向的集中力1000kN。材料弹性模量E=200000MPa，材料泊松比为0.3。

MATLAB编程实现该悬臂梁的弹性静力分析，采用Q4平面单元进行模拟，并将MATLAB的计算结果与Abaqus分析结果进行对比。

• MATLAB 编程

disp

stress

• Abaqus

Model

stress

单元应力对比 (MATLAB vs. Abaqus)

单元编号	节点编号	ABAQUS结果	MATLAB结果	相对偏差（%）
49	52	105.4220	105.4217	0.000
	53	125.7960	125.7956	0.000
	69	273.3910	273.3905	0.000
	70	226.6090	226.6091	0.000
	169	109.9160	109.9162	0.000
	201	169.1640	169.1643	0.000
	202	244.3070	244.3074	0.000
	203	182.3680	182.368	0.000

可以看出，MATLAB编程计算结果是和 Abaqus计算结果完全一致的。

---

Vote Up +19 Vote Down -3You already voted!

• 注释 ( 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)

欢迎关注 “结构之旅” 微信公众号