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

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

  • 程序作者 ( Author )

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

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


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

  • 微信公众号 ( Wechat Subscription)

WeChat_QRCode

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

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

  1. Zeyu Reply

    崔博士您好,请问这个程序是否可以获得学习呢?

  2. Ade Reply

    崔博士 ,您好!读了您的文章,受益匪浅。正巧,我们最近有个作业需要编制程序。我想询问下,如何编制MALTAB程序用高斯积分实现自重节点等效荷载?

    • CJD Post authorReply

      参考相关书籍把,自重,其实就是一种体力,很多书都有给出体力的等效节点荷载计算过程的。

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.