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

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

  • 程序作者 ( 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的计算结果是完全一致的。


8 0

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

WeChat_QRCode

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

Leave a Reply

Your email address will not be published. Required fields are marked *

:wink: :-| :-x :twisted: :) 8-O :( :roll: :-P :oops: :-o :mrgreen: :lol: :idea: :-D :evil: :cry: 8) :arrow: :-? :?: :!:

This site uses Akismet to reduce spam. Learn how your comment data is processed.