[FEM][MATLAB][有限元] FEM Modal Analysis Programming with MATLAB (Frame Elements) (框架单元模态分析编程)

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 接着博文《[动力学] 模态分析基础 (Basics of modal analysis)》继续介绍框架单元用于模态分析的方法。该部分内容也是 书本 《有限单元法:编程与软件应用》模态分析章节的部分内容节选。 10.3.1 梁单元质量矩阵 集中质量矩阵(局部坐标) 设杆材料密度为 ,单元长度为 ,截面积为 ,每个节点分担单元1/2的平动质量,无转动惯量,则单元质量矩阵            (10.3‑1) 10.3.2 算例:2D框架模态分析 算例采用与2D框架结构静力分析中相同的结构,采用欧拉梁单元,材料密度为2.5493e-9t/mm3,采用集中质量矩阵,将单元质量集中于两端节点,且不考虑集中质量后的节点转动质量。 由于算例模态分析的Matlab代码与前面章节中静力分析的Matlab代码大体相同,因此这里仅给进行模态分析所需的新增代码。 % Truss 2D modal analysis …

[FEM][MATLAB][有限元] FEM Modal Analysis Programming with MATLAB (Truss Element) (桁架单元模态分析编程)

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 接着博文《[动力学] 模态分析基础 (Basics of modal analysis)》继续介绍桁架单元用于模态分析的方法。该部分内容也是 书本 《有限单元法:编程与软件应用》模态分析章节的部分内容节选。 10.2 2D桁架结构模态分析 模态分析需要的参数主要有结构的刚度矩阵和质量矩阵,其中的刚度矩阵与静力分析中的刚度矩阵相同,这里仅给出进行模态分析所需的质量矩阵等其他条件。 10.2.1 桁架单元质量矩阵 集中质量矩阵(局部坐标) 设杆材料密度为 ,单元长度为 ,截面积为 ,每个节点分担单元1/2质量,则单元质量矩阵                 (10.2‑1) 整体坐标下的质量矩阵转换 整体坐标下单元一致质量矩阵为    …

[程序] 2001版及2010版抗震设计规范CQC振型耦联系数对比(Cross-Modal coefficients of CQC method)

实干、实践、积累、思考、创新。 程序图标 ( Program LOGO ) 程序介绍 ( Program Introduction) CQC扭转耦联系数计算工具。Cross-Modal coefficients of CQC method in Chinese Seismic Code GB 50011-2001 and GB 50011-2010。 (1)振型扭转耦联系数与振型的周期比及阻尼比有关,在阻尼比一定的情况下,随着周期比的减小,耦联系数减小,并趋于0。即振型的周期越接近,扭转耦联系数越大,当振型周期(频率)相等时,扭转耦联系数为1。 (2)2001版及2010版抗震设计规范CQC振型耦联系数对比。2001版CQC扭转耦联系数是2010版扭转耦联系数的在阻尼比相等情况下的简化。当振型频率比较接近,振型之间阻尼比相差比较大时,两者相差比较大,即2001版规范的简化公式误差比较大。 程序图例 ( Program Gallery ) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] Abaqus Disp Post: Abaqus 后处理工具(Disp篇)

20180628写下的题目,忙完项目现在终于有时间整理了。 实干、实践、积累、思考、创新。 程序图标 ( Program LOGO ) 程序介绍 ( Program Introduction) Abaqus 后处理小工具之一, Abaqus Disp Post,结合python及.NET 开发的工程设计工具,提取结构的位移相关的结果,并自动化形成便于报告的图表。Abaqus分析利器。 程序图例 ( Program Gallery ) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] 隔震设计简化计算(抗规附录L)(Seismic Isolation)

实干、实践、积累、思考、创新。 程序介绍 ( Program Introduction) 抗规附录L计算减隔震,实际没啥用。和小伙伴做的小软件,那时候小伙伴刚开始学习编程,主要用于编程训练,学习规范。 程序图例 ( Program Gallery ) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[结构设计][规范] 等效侧向力法计算非结构构件地震力 (抗规13.2节)

实干、实践、积累、思考、创新。 程序图标 ( Program LOGO ) 程序介绍 ( Program Introduction) 抗规13.2节计算非结构构件地震力,实际没啥用。和小伙伴做的小软件,2015年,那时候小伙伴刚开学编程,主要用于编程训练,学习规范。 程序图例 ( Program Gallery ) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] Abaqus Fre Post: Abaqus频率数据处理

实干、实践、积累、思考、创新。 好吧,一年前的东西半年前说整理,现在终于整理了。 程序图标 ( Program LOGO ) 程序介绍 ( Program Introduction) Abaqus 后处理小工具之一, Abaqus Fre Post。 程序图例 ( Program Gallery ) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] FileTimeManager: A Tool for Managing Files by Time [按时间管理文件的工具]

程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 按时间管理文件。搜索指定时间范围的文件。并将其复制到某个位置。 用于一些特殊情况。当时间对你很重要的时候。 程序界面( Program Interface ) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[编程][算法][游戏] 编程训练: 3 5 7 游戏开发

程序图标 ( Program Icon )      程序介绍 ( Program Introduction) 3 5 7扑克牌游戏,基本规则是,每次从一行中抽取任意张牌,最后一张牌,谁抽到谁输。 前一排年会上,容总发起了 3 5 7游戏挑战,我连续下了两盘,都输了,所以回去研究了一下,游戏十分有意思,也是很有规律的,并编写了这个小程序。现在看这个游戏,还是很简单的。感兴趣的小伙伴可以试试。 程序界面 ( Program Interface ) 下载 ( Download ) GAME 357: Game357_byCJD_20180929.7z 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] Pushover Performance Point Solution Program [Pushover 分析性能点求解程序][Based on ATC 40]

说好的做完项目后更新……半年前写的,现在终于有时间更新一下了。 程序图标 ( Program Icon )      程序介绍 ( Program Introduction) 用于pushover分析性能点的求解,提供了基于ATC-40及FEMA 440的求解方法。也是之前用Abaqus进行pushover分析编写的程序。对这方面感兴趣的小伙伴可以看一看。 程序界面 ( Program Interface ) 注释 ( Comments ) 如果您需要这个软件,请在这给网站捐助:http://www.jdcui.com/?page_id=4813,捐助金额不低于200RMB,捐助后在评论区留下您的评论,软件会通过管理员QQ:3014479529给您发送,对软件有任何疑问也欢迎联系,敬请尊重劳动成果 !!! 如果原软件功能无法满足您的需求,您可以和我联系,可以针对您的需求,做局部修改和调整。 相关话题 ( Related Topics) [01]  [Tool] Pushover Performance Point Solution Program [Pushover …

[日记][FEM] 新书《有限单元法: 编程与软件应用》马上就要出版啦!!

新书《有限单元法: 编程与软件应用》马上就要出版啦。 对这方面内容感兴趣的小伙伴 欢迎多多关注 多多支持!!!先来看看Demo吧。 PS. 我们精心设计的封面很不错哦 ~ ~ ~ 🙂   🙂   🙂 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[软件案例] 倾斜外框引起的扭转效应 (ENGT超高层应用案例1)

实干、实践、积累、思考、创新。 在竖向荷载作用下,超高层倾斜外框会在楼面产生水平分力,若外框架统绕核心筒一向倾斜,则这水平分力将产生扭矩。进而在剪力墙中产生反向扭矩。而这一扭矩要靠楼板传递。基本原理如下图,(PS: 看不到屏幕的手写板实在是一个废物。)需要避免过大的扭转效应。 以下是一组方案核心筒的扭矩对比图。 方案2顶部区域扭矩较方案1大。从抵抗扭转的负担来看,方案1优于方案2。 采用 ENGT: Engineering Tookit [建筑结构辅助设计工具集成系统] 可以快速提取指定工况下剪力墙的扭矩,如下图所示。 关于 ENGT ( About ENGT) Link: ENGT: Engineering Tookit [建筑结构辅助设计工具集成系统] ENGT 应用案例 ( Application Examples of ENGT) [01] [结构][设计][高层建筑][笔记] 倾斜外框引起的扭转效应 (ENGT 超高层应用案例1) [02] …

[VBA][编程] VBA使用“TextFrame2”设置坐标轴属性出错

用VBA开发EXCEL,为修改坐标轴刻度线字体的大写,通过录制宏看到如下代码: ActiveChart.Axes(xlValue).Select With Selection.Format.TextFrame2.TextRange.Font   .Size = 12   .Bold = True End With 可以从 宏 中看见,代码是通过 TextFrame2 属性对刻度线姿态进行控制的。 但将上述代码在程序中编译,则发现报错,提示  TextFrame2 无法识别,不存在。 不知道是Bug还是什么,最后通过 TickLabels 属性解决。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] DataToExcel: A Tool for Converting Data to Excel Charts [将数据转换为Excel图表]

程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 读取数据绘图,并保存为 EXCEL 图表。Excel都懒得用。 (更强大的 Data To Excel Graph V2019版已经更新了,有需要的小伙伴请移步这个网址:http://www.jdcui.com/?p=12685) 程序界面 ( Program Interface ) (更强大的 Data To Excel Graph V2019版已经更新了,有需要的小伙伴请移步这个网址:http://www.jdcui.com/?p=12685) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] Midas Gen分组节点提取工具 [Midas Gen Grouping Node Extraction Tool]

小工具,不一定有用,但是也可能很有用,也许做项目就用到了。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) Midas Gen分组节点提取工具。不知道为什么会有这个东西。Any way,也许你会用到呢。。。 主要将Midas Gen分组中的节点提取出来,保存到Excel。 程序界面 ( Program Interface ) 程序下载( Program Download) 如果您需要这个软件,请在这给网站捐助:http://www.jdcui.com/?page_id=4813,捐助金额不低于100RMB,捐助后在评论区留下您的评论,软件将通过邮件给您发送,敬请尊重劳动成果 !!!。) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] LoopModifier 2018: Hysteretic Loop Data Correction [滞回曲线修正工具]

新版的滞回曲线修正软件已经更新了,移步这个页面:[试验][软件][科研] LoopModifier 2024: Hysteretic Loop Modification Program[滞回曲线修正工具 2024版] 博士论文做试验的时候写的 试验数据修正 小工具。分享给有需要的小伙伴。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 修正试验滞回曲线上的一些明显错误。 总之,一言难尽。用过比较差的国产设备做过试验的小伙伴可能深有体会。 程序界面 ( Program Interface ) 应用案例 ( Example) 在校的时候,做过一个试验,当时试验仪器及设备较差,结果临时除了问题,好不容易昨晚试验,一画滞回曲线,结果傻了眼,根本没办法用。后面仔细分析了数据,数据还是可以用的,但是需要想写办法来修正。于是当晚写了这程序,挽救了这个试验结果,也不用重新再做,都是心血啊,节省了时间、精力及金钱。以下具体的效果,堪比整容,但实际上修正没那么复杂,也不是随意改的: 程序下载( Program Download) 如果您需要这个软件,请在这给网站捐助:http://www.jdcui.com/?page_id=4813,捐助金额不低于 …

[程序] MTS_PLOT: MTS作动筒滞回曲线绘图工具

博士论文做试验的时候,写的小工具。分享给有需要的小伙伴。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 以前做试验的时候,一开始采用 MTS作动筒时,系统还没完全配置好。无法实时输出 构件的滞回曲线,进行实时查看。于是写了个小程序,实时查看 MTS 作动筒的力-位移滞回曲线。可设置软件自动刷新的时间间隔。并且将数据输出来。 程序界面 ( Program Interface ) 程序下载( Program Download) (  若需要程序,如果您需要这个软件,请在这给网站捐助:http://www.jdcui.com/?page_id=4813,捐助不少于50RMB,请在评论区留下您的评论,软件将通过邮件给您发送,敬请尊重劳动成果 !!!) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[FEM][有限元][编程][Matlab][Code by myself] 4节点四面体单元 (C3D4)

  (业余时间 和小伙伴一起写写有限元程序 ) 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所;2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 4节点四面体线性“完全积分”单元,用Abaqus中的命名规则,即 C3D4。 问题描述( Problem Description) 悬臂梁,悬臂长度2.0m,梁高0.5m,梁宽0.2m。梁左端嵌固,受重力作用。材料弹性模量E=200000MPa,材料泊松比为0.3。 基于MATLAB编程实现该悬臂梁的弹性静力分析,采用C3D4单元进行模拟,并将MATLAB的计算结果与Abaqus、midas Gen分析结果进行对比。 MATLAB 编程 model disp stress …

[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 …

[FEM][有限元][编程][Matlab][Code by myself] 平面4节点线性“完全积分”单元(CPS4)

  (业余时间 和小伙伴一起写写有限元程序 ) 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所; 2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 平面4节点线性“完全积分”单元,用Abaqus中的命名规则,即CPS4。 问题描述( Problem Description) XZ平面内的悬臂梁,悬臂长度2.0m,梁高0.5m,梁宽0.2m。梁左端嵌固,右端受到-z方向的集中力1000kN。材料弹性模量E=200000MPa,材料泊松比为0.3。 MATLAB编程实现该悬臂梁的弹性静力分析,采用Q4平面单元进行模拟,并将MATLAB的计算结果与SAP2000、midas Gen分析结果进行对比。 MATLAB 编程 disp stress SAP2000 …