[FEM][有限元][编程][Matlab][Code by myself] 8节点六面体单元(C3D8)(8-node linear brick)

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

程序作者 ( Author )

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

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的计算结果是完全一致的。

You 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! )