[FEM][MATLAB][有限元][编程] 压杆稳定问题MATLAB有限元编程 (《有限单元法-编程与软件应用》章节节选)

实干、实践、积累、思考、创新。 接着博文《[力学][有限元][FEM]Basics of Buckling Analysis [曲屈分析基础]》继续介绍经典材料力学或结构力学课本上介绍的压杆稳定问题。该部分内容也是 书本 《有限单元法:编程与软件应用》屈曲分析章节的部分内容节选。 11.5 屈曲分析3:压杆稳定 作为屈曲分析的补充,本节讨论一下压杆稳定问题。 图 11‑10 压杆支座情况 算例结构为一根等截面轴心受压直杆,直杆材料为钢管,钢管外径100mm,管厚5mm,高5m。一共考虑了5种约束情况,分别为:1.两端铰接;2.一端铰接、一端嵌固;3.两端嵌固;4.一端嵌固、一端滑动;5.一端嵌固、一端自由。针对每种约束情况,分别将压杆划分为1个、2个、5个、10个、20个梁单元,进行屈曲分析。 11.5.1 MATLAB代码与注释 本节以底端嵌固、上端自由、划分20个单元的情况为例,给出进行屈曲分析所需的主要代码。 % Pressed Bar Buckling Analysis % Author : JiDong Cui(崔济东),Xuelong Shen(沈雪龙) % Website : …

[Midas2Abaqus] 钢柱屈曲分析 (Buckling Analysis of Steel Column)

实干、实践、积累、思考、创新。 底部固支的H型钢悬臂柱,顶部竖直向下作用轴压力,分析构件的屈曲模态。为了考虑翼缘或腹板的局部屈曲,可采用壳单元进行建模。本例分别采用Midas Gen及Abaqus软件进行模拟。 Midas Gen Model Midas Gen Results (a)一阶屈曲模态(屈曲因子:8.24E+005)(绕工字钢的弱轴) (b)二阶屈曲模态(屈曲因子:1.976E+006)(绕工字钢的强轴) (c)三阶屈曲模态(屈曲因子:1.976E+006)(绕工字钢轴向扭转) Midas2Abaqus Model 为了进行Abaqus分析,采用Midas2Abaqus软件进行模型转换。 Abaqus Model Abaqus Results (a)一阶屈曲模态(屈曲因子:8.26583E+005)(绕工字钢的弱轴) (b)二阶屈曲模态(屈曲因子:1.98061E+006)(绕工字钢的强轴) (c)三阶屈曲模态(屈曲因子:1.976E+006)(绕工字钢轴向扭转) 由以上两个软件分析结果可知,Midas Gen及Abaqus的分析结果基本一致,一阶屈曲模态为绕工字钢的弱轴,二阶屈曲模态为绕工字钢的强轴,三阶屈曲模态为绕工字钢轴向的扭转屈曲。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[力学][有限元][FEM]Basics of Buckling Analysis [曲屈分析基础] (《有限单元法-编程与软件应用》章节节选)

实干、实践、积累、思考、创新。 《有限单元法:编程与软件应用》 曲屈分析 章节的内容节选: 11.1 稳定问题分类 结构失稳(屈曲)是指在外力作用下结构的平衡状态开始丧失,稍有扰动变形便迅速增大,最后使结构发生破坏。稳定问题一般分为两类,第一类是理想化的情况,即达到某种荷载时,除结构原来的平衡状态存在外,还可能出现第二个平衡状态,所以又称平衡分岔失稳或分支点失稳 (Buckling, Bifuraction)(图1中OAB曲线),对应于数学中是求解特征值问题,故又称特征值屈曲,此类结构失稳时相应的荷载称为屈曲荷载。第二类是结构失稳时,变形将迅速增大,而不会出现新的变形形式,又称极值点失稳(Instability)(图1中OCD曲线),结构失稳时相应的荷载称为极限荷载。此外还有一种跳跃失稳(Snap-through)(图1中OEF曲线),当荷载达到某值时(E点),结构平衡状态发生一个明显的跳跃,突然过渡到非临近的另一个具有较大位移的平衡状态(F点),由于在跳跃时结构通常已经破坏,其后的状态不能被利用,所以可归入第二类失稳。 图 11‑1 稳定问题 11.2 最小势能原理 11.2.1 系统的平衡 若系统处于平衡状态,则一定存在一种位移状态,使得系统总势能为驻值,即               (11.2‑1) 公式中               …

[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) 整体坐标下的质量矩阵转换 整体坐标下单元一致质量矩阵为    …

[程序资料] 多方案曲线结果分析对比 (YJK Multi-Model Compare 功能介绍)

给 YJK Multi-Model Compare 增加多模型多种整体曲线对比的功能,用于方案对比及工程备份。   About Multi-Model Compare 关于软件,可在这个链接查看,链接:[程序] Multi-Model Compare: 基于盈建科的多模型参数优化比选工具 [超限结构设计助手] 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[结构设计] 超高层伸臂位置的敏感性分析 (Analysis of optimal outriggers location in high-rise building structures)

实干、实践、积累、思考、创新。 通过定性的考虑可的伸臂位置,进行分组模型对比分析,可以定性评估伸臂位置对结构位移及周期的影响,定性把握伸臂的合理位置。 以下某个工程的参数分析 伸臂位置对位移角的影响: 由上图可见,伸臂布置于2/3/8/9/10位置对于提高结构刚度,减小位移角效率相对较大。 伸臂位置对周期的影响: 由上图可见,伸臂布置于3/5/7位置对于控制结构的周期,提高结构整体刚度较为有效。 做个简单笔记备忘,后续再利用计算机优化算法系统研究…… 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] midas2Abaqus: midas Gen Model to Abaqus Model [midasGen模型转Abaqus]

The page you current visit has been moved to: http://www.jdcui.com/?page_id=10470 最新内容已移步到这个页面: http://www.jdcui.com/?page_id=10470 程序图标 ( Program LOGO ) 程序介绍 ( Program Introduction) Midas Gen 模型转 Abaqus模型的工具 (Midas Gen Model to Abaqus Model) 程序图例 ( Program Gallery ) …

[软件笔记] 记录盈建科一个错误:无法生成楼板

调整模型,无法进入前处理。最后检查发现总是卡在 生成楼板一步。仔细检查发现个别楼层,无法生成楼板,点击生成楼板就会死机。 检查无果,最后核心筒构件删除,重新再建核心筒就可以生成楼板了。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序案例] BeamEvelopes应用案例:移动荷载下车道侧壁内力包络分析

车道侧壁与地下室外墙通过车道的楼板相连,将楼板等效成支座,则力会施加到车道侧壁。 车道侧壁的弯矩受力分析类似一个竖向的连续梁。考虑到车道板位置的不确定性,大部分车道板是倾斜的,为后期留点安全富余,因此需要考虑车道板位置的变化对车道侧壁受力的影响,即进行移动荷载的包络分析。一下分析均采用www.jdcui.com 开发的软件 BeamEvelopes 进行,具体链接:([软件][编程][开发][结构设计] BeamEvelopes: 连续梁包络分析工具(Continuous Beam Envelope Analysis Program)) (1)受力简图: (2)仅考虑荷载施加在跨中的弯矩图及剪力图如下: 由上图可见,跨中加载,以第一跨为例,最大弯矩为:3.304e+9,最大剪力为1.623e+6; (3)弱考虑集中荷载施在跨内位置可变动,求解弯矩包络图及剪力包络图如所示: 由上图可见,包络分析的结果,以第一跨为例,最大弯矩为:4.196e+09,最大剪力为3.114e+006;可见均比跨中加载大了很多,因此设计师需要结合实际情况谨慎考虑。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] BeamEvelopes: Continuous Beam Envelope Analysis Program [连续梁包络分析工具]

程序研发组给设计师打造的一款小工具,可以用于地下室车道侧臂的受力分析设计。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 连续梁静力分析,连续梁移动荷载包络分析。可以用于地下室侧壁相邻的车道侧壁的受力分析。Continuous beam static analysis, continuous beam moving load envelope analysis. It can be used for the force analysis of the side wall …

[程序案例] 伸臂及腰桁架对巨柱轴力的影响 (ENGT 超高层应用案例3)

实干、实践、积累、思考、创新。 模型对比分析,研究伸臂及腰桁架对巨柱力的影响。 采用 ENGT: Engineering Tookit [建筑结构辅助设计工具集成系统] 可以快速提取指定巨柱的内力,图表形式展示相应的数据,如下图所示。 可见伸臂及腰桁架的刚度对巨柱的轴力均有影响,中下部腰桁架的影响更大。 关于 ENGT ( About ENGT) Link: ENGT: Engineering Tookit [建筑结构辅助设计工具集成系统] ENGT 应用案例 ( Application Examples of ENGT) [01] [结构][设计][高层建筑][笔记] 倾斜外框引起的扭转效应 (ENGT 超高层应用案例1) [02] [结构设计][超高层][ENGT] 巨柱倾角对外框剪力分担比的影响 …

[程序案例] 巨柱倾角对外框剪力分担比的影响 (ENGT 超高层应用案例2)

实干、实践、积累、思考、创新。 盈建科默认的外框剪力分担比的计算柱、墙剪力采用内的是CQC的剪力,总剪力采用的是规定水平力,因此,是不满足静力平衡的,即 外框剪力 + 墙剪力 不等于 总剪力。 可以将规定水平力按风荷载的方式输入,进而可以按静力平衡的方式求解外框剪力分担比。 采用 ENGT: Engineering Tookit [建筑结构辅助设计工具集成系统] 可以快速提取外框剪力分担比,并可以通过简易算法,评估巨柱倾斜角度对外框剪力分担比的影响,如下图所示。 关于 ENGT ( About ENGT) Link: ENGT: Engineering Tookit [建筑结构辅助设计工具集成系统] ENGT 应用案例 ( Application Examples of ENGT) [01] [结构][设计][高层建筑][笔记] 倾斜外框引起的扭转效应 …

[程序] midas Gen节点处理工具(如何避免盈建科斜杆丢失)

程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 基于midas Gen对模型节点处理,避免YJK模型杆件丢失。PS:不知道你有没有遇到YJK斜杆丢失的问题。 程序界面( Program Interface ) 例子( Example ) 处理前,构件丢失。 处理完后,构件不再丢失。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[结构分析] Pushover Analysis by Abaqus [Abaqus Pushover 分析]

大半年前做的分析(20171208),时间过得真快,有点忘记了。当时太忙忘记整理。不过既然题目写了下来,那就还是大概整理一下吧。当时写的小软件还在,所以,有些小程序这个好,相当于对一些东西的备忘,笔记。以下是采用ABAQUS进行pushover分析的简单步骤: (1)打开程序([01]  [工具][结构分析][编程][Abaqus] Pushover Performance Point Solution Program [Pushover 分析性能点求解程序][Based on ATC 40]),导入振型参数,输出ABAQUS pushover分析用的 .inp 参数文件 (2)进行Abaqus Push Over分析,获得X、Y两个方向塔楼的基底剪力-顶点位移曲线 (3)导入相关参数,用程序迭代求解性能点,并输出相关参数,整理成报告 (a)X向性能点求解过程 (b)Y向性能点求解过程 表1.1 X、Y向性能点处结构指标 X向 Y向 剪力V(KN) 154835 148871 剪重比 …

[程序] 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] …

[软件笔记] YJK简化模型楼板导荷注意问题

上图的简化方式建的楼板,板2 和 板1的楼板荷载会重复,多算一个板1的荷载。至少需要设置一个梁将板1边线和板2边线连接,如下图所示。这个在建简化模型的时候需要注意。 PS. 幸亏 常博士 指出问题。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[程序] YJKWv02qPost: 盈建科0.2q0数据处理 [超限结构设计助手]

程序图标 ( Program Icon )      程序介绍 ( Program Introduction) [超限结构设计助手] 系列程序之一。 程序界面 ( Program Interface ) 相关软件 ( Related Program ) [01] ENGT: Engineering Toolkit [建筑结构辅助设计工具集成系统] [02] [风洞试验][结构设计][软件] RWDI风洞试验荷载数据处理工具 [03] [结构设计][软件][Program] YJK风洞荷载试验数据处理工具[ A …

[程序] 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) 欢迎关注 “结构之旅” 微信公众号

[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] 平面6节点二次“完全积分”单元(CPS6)

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