[FEM][有限元][编程][Matlab][Code by myself] FEM Analysis: 2D Euler Beam Element [有限元分析: 2D欧拉梁单元]

(  有空和小伙伴一起写写有限元程序

  • 程序作者 ( Author)

JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2

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

  • 问题描述( Problem Description)

节点1、5处为固定支座,节点4处受到+x方向P=200kN的集中力作用。结构中各杆件采用相同的材料,弹性模量E=30000MPa,梁、柱截面面积分别为0.08m2和0.16m2,梁、柱截面惯性矩分别为0.0128/12 m4和0.0256/12 m4

基于MATLAB编程实现该框架结构的弹性静力分析,所有构件采用欧拉梁单元进行模拟,并将MATLAB的计算结果与SAP2000、midas Gen分析结果进行对比。

  • MATLAB 编程结果

  • SAP2000 计算结果

MATLAB vs. SAP2000 Displacement

节点 1节点 2节点 3节点 4节点 5节点 6节点 7节点8
MATLAB011.665333.039553.7172011.666633.037253.5102
SAP2000011.665333.039553.7172011.666633.037253.5102
相对偏差( % )00.000200.000100.0002-0.00010
MATLAB00.15850.26550.31140-0.1585-0.2655-0.3114
SAP200000.15850.26550.31140-0.1585-0.2655-0.3114
相对偏差( % )00-0.0002000-0.00020.0001
MATLAB0-0.0054-0.0065-0.00490-0.0054-0.0065-0.0049
SAP20000-0.0054-0.0065-0.00490-0.0054-0.0065-0.0049
相对偏差( % )0000.0001000.00010
  • Midas Gen计算结果

MATLAB vs. MidasGen Displacement

响应类型节点 1节点 2节点 3节点 4节点 5节点 6节点 7节点8
DxMATLAB011.665333.039553.7172011.666633.037253.5102
Midas Gen011.665333.039553.7172011.666633.037253.5102
相对偏差( % )00000000
DzMATLAB00.15850.26550.31140-0.1585-0.2655-0.3114
Midas Gen00.15850.26550.31140-0.1585-0.2655-0.3114
相对偏差( % )00000000
RyMATLAB0-0.0054-0.0065-0.00490-0.0054-0.0065-0.0049
Midas Gen0-0.0054-0.0065-0.00490-0.0054-0.0065-0.0049
相对偏差( % )00000000

You already voted!

  • 注释 ( Comments )

源代码已经收录到新书《有限单元法:编程与软件应用》,详见链接:http://www.jdcui.com/?page_id=9731

  ( 如果您发现有错误,欢迎批评指正。邮箱: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

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

3 thoughts on “[FEM][有限元][编程][Matlab][Code by myself] FEM Analysis: 2D Euler Beam Element [有限元分析: 2D欧拉梁单元]

  1. Yang Reply

    Hi Dr. Cui,
    I wanna buy your new book, but I am now in USA, so I am wondering if there is a way that I can get it. Thank you. My email is attached to make communication easier.

    Best,
    ZY

Leave a Reply

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

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