【BiliBili】 【知乎】 【QQ空间】 【微博】


博学慎思,坚持创新,求知若飢,虛心若愚                       关键字     类型    
  • 博主
  • 全部
  • 日志
  • 研究
  • 论文
  • 自编程序
  • 工程
  • Dino笔记
  • DinoBox
  • OpenSEES
  • Dino公开课
  • 每周挑战
  • Dino爱编程
  • 留言板

DinoBox

ETE_Frame3D
结构分析程序
ETE_Frame3D
2022-01-19
ETE-section截面分析程序
ETE-section
2021-08-19
ETO_OpenSEES前后处理程序
程序通过邮件申请。
★注意,软件不能下载,是通过邮箱发送软件的。

申请ETO时,注意EMAIL 的写法:(格式一定要对,否则电脑收不到)
邮件题目: ETO (版本号) 下载
注册人: (真实姓名)
学校/单位: (大学或工作单位)
邮件回复,大约需要1~4天的时间。
ETO_程序
2014-03-26
上一篇 | 下一篇
自编程序 订阅所有自编程序的日志

[程序]ETE-Frame3D 程序增加壳元ShellElement

作者:dinochen1983 日期:2022-02-03

字体大小: 小 中 大

ETE-Frame3D 程序增加壳元ShellElement


微信公众号用户,点击下方【阅读原文】可以下载程序
2022-02-03,第2版发布
【程序下载】 程序: ETE-Frame3D有限元分析程序(版本2022-02-03)
【实例下载】 实例: Shell Element 实例1~实例5


ETE-FRAME3D上次增加了三维的杆系单元(Frame单元),这次增加壳单元(shell element),壳单元按节点数量可划分为三角形单元与四边形单元,根据受力可分为平面应力单元(Membrane Element)、板单元(Plate Element),壳单元(Shell Element),以下所有单元理论的综合就成了最终的FRAME-3D的壳单元。以下是FRAME-3D建立壳单元模型三维图(三角形单元建立圆筒模型)。以下是采用的单元介绍与实例验证记录。
FRAME-3D增加了壳单元以后,可以成为学习有限元理论的工具。
  
其实壳元=板元+膜元,壳元理论包括了T3,Q4,T6,Q8等参元的全部内容,王勖成老师的《有限单元法》里面有大量的介绍,书的后面附有源代码可供学习。
T3,Q4单元用于平面应力分析单元分析
T6,Q8单元理论用于板单元的分析
学会了平面单元与板单元理论,就可以组合成为壳单元了。
 
总结一下有限元程序编程的要点:
关于形函数与高斯积分可以专门做两期结构笔记来讲一下,现在的有限元教材讲得太复杂了。
(a) 了解不同单元的形函数及矩阵[B]的生成方法
(b) 了解单元的单刚自由度与总刚度矩阵之前的关系
(c) 三维单元都需要进行局部坐标与整体坐标的转换
(d) 了解高斯积分方法的计算
(e) 计算整体位移后,转换得到单元的变形
(f) 再通过[B][D]矩阵与单元变形计算单元的内力应力
平面单元,如壳元,板元及膜元的原因是类似的,单元刚度的形成有三个要素
(a) 能过形函数得到单元与等参元的几何关系,也就是求解[J]矩阵
(b) 对形函数求一次偏导二次偏导,得到[N]矩阵,用于形成[B]矩阵
(c) 对不同高斯积分点求解[B][D]矩阵后,积分得到单元矩阵,通过[J]矩阵的反变换把等参元的单刚变成原来坐标单元的单刚。
(1) 四边形平面应力单元Q4_Membrance测试
Q4_Membrance单元就是线性平面应力单元,线性代表插值函数为一次函数。
 
算例1描述:长度6m高度3m,厚度为100mm,材料Ec=25000MPa,柏松比为0.2,荷载为1.5kN*13。四边形平面应力单元模型采用线性4结点单元,4个高斯积分点。计算的S11应力竖向位移UZ如下图所示。Frame-3d与SAP2000的计算结果基本吻合。
 
FRAME3D应力S11云图绘制效果

sap2000应力S11云图绘制效果

sap2000与FRAME-3D的主应力对比图

sap2000与FRAME-3D的竖向位移对比图
 
(2) 四边形板单元Q4_Plate测试
Q4_Plate也就是四边形等参薄板单元,就是二次板单元,二次代表插值函数为二次函数,虽然只有4个结点,但是形函数是采用8结点的形函数。
算例2描述:长度6m高度3m,厚度为250mm,材料Ec=25000MPa,柏松比为0.2,荷载为1.5kN*13。四边形平面应力单元模型采用二次性8结点单元(外部为4个节点),4个高斯积分点。计算的M11应力,所有节点的竖向位移UZ如下图所示。Frame-3d与SAP2000的计算结果基本吻合。

SAP2000的M11内力云图

FRAME-3D的M11应力云图

FRAME-3D与SAP2000的竖向位移与应力的对比结果 
 
(3) 三角形平面应力单元T3_Membrance测试
 
T3也叫做常应变三角形单元,整个单元只有一个应变值。
算例3描述:长度6m高度3m,厚度为100mm,材料Ec=25000MPa,柏松比为0.2,荷载为15kN*7。三角形平面应力单元模型采用常应力3结点单元,1个高斯积分点。计算的Smax应力竖向位移UZ如下图所示。Frame-3d与SAP2000的计算结果基本吻合。

FRAME-3D的SMAX应力云图

FRAME-3D的SMAX应力云图

FRAME-3D与SAP2000的位移对比

FRAME-3D与SAP2000的主应力SMAX的对比
 
 
(4) 三角形板单元T3_Plate测试
对于实际结构或构件,四边形的划分不一定满足边界需求,一般三角形与四边形是混合使用的。所以一定要增加三角形单元。
T3_Plate也就是三角等参薄板单元,就是二次板单元,二次代表插值函数为二次函数,虽然只有3个结点,但是形函数是采用6结点的形函数。
算例4描述:长度6m高度3m,厚度为100mm,材料Ec=25000MPa,柏松比为0.2,荷载为1.5kN*11。三角形板单元模型采用线性6结点单元(外部为3个节点),3个高斯积分点。计算的Smax应力竖向位移UZ如下图所示。Frame-3d与SAP2000的计算结果基本吻合。
 

FRAME-3D的SMAX应力云图

SAP2000的SMAX应力云图

SAP2000与FRAME-3D的应力与位移的对比结果 
(5) 圆孔平板拉伸问题测试

弹性力学有解析解的传统问题,圆孔平板拉伸问题。经典中的经典,通过对移原理,只建立1/4的模型
主要检测三角形与四边形单元协同使用。

 
算例5描述:长度6m高度3m,厚度为100mm,材料Ec=25000MPa,柏松比为0.2,荷载为1.5kN*11。计算的Smax应力对比如下图所示。Frame-3d与SAP2000的计算结果基本吻合。

FRAME-3D的应力SMAX计算结果

SAP2000的应力SMAX计算结果


SAP2000与FRAME-3D的SMAX应力对比
通过学习王勖成老师的《有限单元法》的相关理论,结合编程软件DELPHI与PYTHON编制出与SAP2000计算结果吻合的二维单元,以后可用于有限元理论的学习,工程优化的核心处理器,编程过程中与有限元理论还有很多细节,可以在以后的结构笔记上展开讨论。通过有限元编程,可以了解到SAP2000的壳元、板元膜元的计算核心,有一种豁然开朗的感觉,有限元工具再不是黑箱子了。附年可以下载以上5个实例的文件,往后增加一些不规则结构计算对比结果。
往后在DINOSTRU结构笔记,增加有限元的学习内容
====本节完====

请关注我的微信公众号:DINOSTRU


更多有趣的结构分析知识分享



[本日志由 dinochen1983 于 2022-02-03 11:37 PM 编辑]
文章来自: 本站原创
引用通告: 查看所有引用 | 我要引用此文章
Tags: 有限元有限元
评论: 0 | 引用: 0 | 查看次数: 1380
发表评论
昵 称:
密 码: 游客发言不需要密码.
验证码: 验证码
Dino:如果看不到验证码,再刷几下就可以了
内 容:
选 项:
虽然发表评论不用注册,但是为了保护您的发言权,建议您注册帐号.
字数限制 200 字 | UBB代码 关闭 | [img]标签 关闭

信息(Information)

DinoChen 视频频道

My GitHub Page



扫描以下二维码,加入公众号!


dinochen.com需要你的支持


日历(Calendar)

2022年7月
  • 日
  • 一
  • 二
  • 三
  • 四
  • 五
  • 六
  • 26
  • 27
  • 28
  • 29
  • 30
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6

查找(Search)

关键字
 
类 型

文章(Recent Article)

[OpenSEES]OpenSEES实例教程[汇总...[Dino笔记]粒子群算法在桁架结构中的应用[获奖]WSP亚洲2021年卓越奖获得一等奖与创新...[Dino笔记]基于遗传算法对剪力墙结构优化设计方...[Dino笔记]基于OpenSEES虚功原理的钢结...[OpenSEES] 新版ETO更新及相关下载方法...[新闻]祝韩老师荣获广东省科技进步一等奖[程序]ETE-section功能更新:可修改划分...[App开发]采用Swfit语言编制经典的推箱子游...[程序]ETE-section更新_考虑混凝土保护...

评论(Comments)

陈博士您好,这个小工具好像有点问题,有效翼板宽度b...打开的时候显示“你没有许可证”,怎么弄啊陈博士您好,这个是需要购买的软件还是可以试用的怎么查看文档计算过程呢,陈博士一直显示出现一个错误字符无法计算是什么原因呢?为什么一直出错计算不了?陈博士,能共享下文件嘛?第二版没有光盘陈博士您好,我换了好几台电脑ETE还是需要许可,陈博士您好,现在程序还是打不开,需要许可请问如何获得会员许可证?

留言(GuestBook)

Thanks regarding furnishi...陈博士您好,想请问一下零长度单元用在桥梁的支座上,...陈博,新版ETO下载以后界面打不开是什么原因呢,旧...陈博士您好,请问可以出一期画易损性曲线的教程吗陈博士您好,我在学习您4.33侧向多自由度简化模型...陈博士您好,計算wind2019的program ...陈博士,您好,我下载使用了您编写的压型钢板组合楼板...陈博士您好,我是opensee初学者,想向您请教两...陈博士你好,请问开洞混凝土截面,ETE-secti...陈博士您好,感谢您的回复,我的刚性隔板设置了转动惯...

TagBar

Loading...

用户面板(User Panel)

登录用户注册

统计(Statistics)

日志: 408 篇
评论: 893 个
留言: 413 个
访问: 14024138 次
在线: 7 人
建站时间: 2008-02-08

Powered By PJBlog3 V2.7 Build 05 CopyRight 2005 - 2011, 陈学伟的博客

Processed in 0.390625 second(s) , 6 queries 粤ICP备15071789号-1