【JY】各软件框架的弹塑性分析对比

文章详细比较了JYSAP、Perform3D、Etabs、Abaqus和SSG等软件在抗震性能化设计中的弹塑性分析,涉及周期、振型、位移分析和损伤评估,强调了弹塑性分析在设计中的关键作用以及不同软件之间的性能对比。
摘要由CSDN通过智能技术生成

e4ce5018a7665953de36801a33a550fe.png

写在前文

目前国家现行规范通常执行的是:进行多遇地震阶段的结构弹性设计,再进行罕遇地震下的构件承载力与变形校核。

罕遇地震下的抗震性能化设计一般采用弹塑性分析方法,是在多遇地震弹性设计基础上的变形复核:当不满足罕遇地震阶段的性能目标时,需要重新进行多遇地震下的结构布置与构件截面调整等设计内容(当轻微不满足时,可只提高弹塑性模型的配筋进行罕遇地震阶段的再评估),总的来说是一个迭代的过程。

当构件超过屈服承载力以后,虽然也可以根据等效弹性分析的内力结果进行有关极限承载力验算,甚至可能认为极限承载力验算可与变形验算的结果建立一定的关联性,但是当结构将出现发生非常显著的刚度退化和内力重分布,此时的结构内力将显著低于根据初始弹性刚度计算的结构内力,因此在采用接近于初始弹性刚度的所谓等效弹性刚度时,构件内力验算往往出现很大的失真并偏于保守。

因此弹塑性分析对于抗震性能化设计起至关重要的作用,并且需要一套准确稳定的计算内核才能保证整体设计的可靠性。

做一个小算例验证各个软件的梁柱单元的非线性性能对比。由于剪力墙的非线性行为较为复杂,本期先不讨论,后面有空再和大家分享~

做一个4×5的框架,采用600×600的柱子,和300×500的梁,采用YJK进行计算分析得到的计算配筋。分别导入 JYSAP、Etabs、Perfrom3d、Abaqus、SSG进行动力弹塑性分析计算。

(注:Abaqus采用自编子程序,非SSG导出子程序。)

带入一条人工波,以PGA=510Gal进行条幅,X:Y=1:0.85进行输入,并简要的对各软件的关键结果进行整理及总结。

73a2f72692ded3ad77f0663e8c7d2d75.png

JYSAPOS V1.0

注:JYSAP目前只是笔者闲时开发的一个小玩具,主要做做计算对比~ 目前为非商业工程软件。

9f45c9044397c698671cc7c3adeb32bd.png

周期与振型对比

周期表T1T2T3T4T5T6备注
JYSAP0.7110.6850.6200.2210.2140.194含配筋
Perform3D0.6920.6650.5940.2130.2070.185含配筋
SSG0.7080.6890.5870.2200.2150.184含配筋
Etabs0.7520.7300.6540.2310.2250.202不含配筋
Abaqus0.7210.6900.6280.2300.2250.196含配筋

可见各软件在线弹性下的周期,对比吻合度较高,可见各个软件线弹性行为基本一致~

从下各个振型图也可看出个软件的振型形态也基本一致。

JYSAPOS 振型

35d429a4892052c9d7e3f6bdb2017c56.png

Perfrom3D振型

4902d4f6f23299f931210ca36cc92ef6.png

Etabs振型

f3b245565877b6af27a9825d6c4e1cca.png

Abaqus振型

5f6e59e84a426bbd3381bd472db1d707.png

SSG振型

c1c9dd54111c0dde4a605c9f87643870.png

顶点位移时程对比

各软件对结构进行动力弹塑性分析,并取顶层位移进行双向的对比分析。

全部对比

将所有曲线放在一起会发现不论是X或者Y都显得十分凌乱,感觉各软件好像没得对比。为了方便对比,将分成两个类型的软件进行区分,好便于对比:

1、隐式计算+宏观单元类

2、显式计算+所谓的精细化单元类

4fb4d09f3b69ab88230bee18e09f23ec.png

图 X-方向

343addf862d41bf3ab2c45e8117a7e95.png

图 Y-方向

JY/Ets/P3d

关于隐式计算+宏观单元类,笔者的自编玩具JYSAP也是深受非线性大师 鲍威尔教授的影响,采用宏观单元的形式进行编写。可见这三个软件的计算结果趋势和相位都基本吻合。(梁单元均采用塑性铰单元、柱单元均采用纤维铰单元)

【JY】ETABS与Perform3D弹塑性分析功能对比示例

c17d599563c495bec2249d9861a132aa.png

图 X-方向

c7f2c2994eb19f59122f8fc9997b3c32.png

图 Y-方向

Abaqus/SSG

关于显式计算+所谓的精细化单元类,采用国内结构工程师保姆级软件SSG,或采用各个工程师自编各具特色的Abaqus的子程序均可进行分析,但是由于显式计算存在蝴蝶效应,因此对结果的近似性或者唯一性很难把握,计算分析结果如下:

d6ead73016de1898e836cc04a6d21784.png

图 X-方向

918ae788b13881e457012ca028e43791.png

图 Y-方向

位移角对比

将所有位移角曲线进行对比整理,如下图所示。

2ea53ceac1abff1f48f5ff8f4869930b.png

损伤评估

Etabs默认用的美标加固规范的性能水准评价,IO、LS、CP,详情可看:

【JY】浅析基于性能的抗震分析方法——性能设计

Abaqus的结果没有用任何标准,仅仅展示的是混凝土的受压损伤。

JYSAPOS、Perfrom3D的评价标准调整为T/CECA 20024-2022 建筑结构抗震性能化设计标准中混凝土构件位移角的性能标准。

e72bbb960d239bc18da435f66d65c605.png ed224f0355b3a9b1bdce33e4afd16379.png

SSG也调整为T/CECS 906-2021 建筑结构非线性分析技术规程中采用材料层面进行评价。(若采用混凝土构件位移角的性能标准,SSG给出的答案是毫无损伤!!)

可见Perform3d、JYSAP可以非常吻合的对比上,证明了JYSAP初步领悟的鲍威尔教授的精髓。

虽然不同的性能评价标准结果趋势可能存在差异,但是Etabs的损伤趋势也可与上述Perform3d和JYSAP吻合得上。

Abaqus和SSG的损伤趋势从表明上估计,仍然勉强可以对的上,但基于材料层面的性能评价显示上显得过于严苛,对于评价是否正确还有待考究。

4f6f250fef3cba06683fcb6de2ce4779.png befe596cff262119e0571bffe526c443.png

小结与JY观点

结构的抗震性能能力在很大程度上取决于模型假设。每种建模的方法可能导致结构产生存在差异的倒塌行为模式。因此能力评估是一件极其困难的事情,并不是简单的钢筋+混凝土=钢筋混凝土,从目前的技术手段来看,精细化分析只是个噱头,并不能做到仿真层面的精细化,即使结构很小也很难建立精确的分析模型,尤其是考虑循环变形时。

01cc14d46c4236c43d948f67036cc457.png

从基于构件变形的宏观校核,包括针对梁、柱、剪力墙等受弯、压弯或拉弯构件的弹塑性位移角指标,是通过大量构件抗震试验统计得到,涵盖了钢砼构件内部各部分的相互作用(热膨胀、徐变、滑移)、力学上有机结合等等。因此,基于构件变形的模型经常采用构件的力—变形模型,包括弯矩一曲率模型、弯矩一转角模型、剪力一转角模型、剪力一位移模型、轴力-位移模型等,且考虑轴力与弯矩的相关性(PMM)等这类宏观力学模型对结构进行计算分析。

从基于混凝土应变或损伤因子、钢筋和钢材应变的材料层次细观校核。由于该评价方式并不能通过直观统计的方式进行验证,通常建议作为备选的评价方式予以采用;对于复杂受力/新型的构件,未通过大量构件抗震试验的,材料应变评价将是唯一的手段。

总的来说,位移角结果在构件试验研究中容易获得从而很自然地成为一项重要的评价指标;甚至说塑性铰或截面曲率的评价聚焦于构件的塑性变形段,更便于真实反映构件的塑性变形情况。因此基于构件变形的宏观校核应常作为优先选用的评价指标。因此,从评价效果看,基于塑性铰转角特别是截面曲率的评价指标是较为理想的评价指标但从与试验结果的对比或验证关系看,位移角更简便直观的。

然而基于材料应变的评价具有一定的唯一性,但它无法反映构件全截面的综合受力情况,特别是将混凝土和钢筋的评价予以分离而难以进行构件变形的综合评判,因此在构件变形模式相对简单时往往不作为优先选用的评价指标。 实际上,基于材料应变与基于构件变形二者在本质上应属表征相同的损伤特征,目前有学者在研究且建立二者相关关系。

爱德华·威尔逊 说过:

我大学一年级的物理教师曾教条地警告学生:“绝不使用你推导不出来的方程。”他还这样说,“如果一个人用 5分钟时间来解决一个与人性命做关的问题,他应该花3分钟的时间去研究这个问题,并清楚地理解它。”在过去的 40 年里,这些简单而实际的教导一直在指导着我的工作,而且我希望这种哲学思想已经传给了我的学生,就现代结构工程而言,可以把这些话重新叙述为 “除非你完全理解了程序采用的理论和近似方法,否则绝不使用任何一个结构分析程序”和“在没有清楚地定义荷载、材料属性以及边界条件之前,不要创建计算机模型”!

400de81030e4ce82d6b7de8eb18bc802.gif

更多精彩,关注建源学堂!

【往期精彩】

# 性能分析

【JY】基于性能的抗震设计浅析(一)

【JY】基于性能的抗震设计浅析(二)

【JY】浅析消能附加阻尼比

【JY】近断层结构设计策略分析与讨论

【JY】浅析各动力求解算法及其算法数值阻尼(人工阻尼)理念

【JY|体系】结构概念设计之(结构体系概念)

【JY|理念】结构概念设计之(设计理念进展)

【JY|减震】结构概念之(消能减震黏滞阻尼器)

【JY|隔震】结构概念设计之(隔震概念设计)

【JY】有限单元分析的常见问题及单元选择

【JY】结构动力学之显隐式

【JY】浅谈结构设计

【JY】浅谈混凝土损伤模型及Abaqus中CDP的应用

【JY】浅谈混凝土结构/构件性能试验指标概念(一)

【JY】浅谈混凝土结构/构件性能试验指标概念(二)

【JY】橡胶系支座/摩擦系支座全面解析

【JY】《公路桥梁抗震性能评价细则》宣贯PPT

【JY】减隔震设计思考:隔震篇

【JY】浅析基于性能的抗震分析方法——性能设计

【JY】这个房子应该做抗震or减隔震?

【JY】减隔震设计思考Ⅱ

# 概念机理

【JY】基于Ramberg-Osgood本构模型的双线性计算分析

【JY】结构动力学初步-单质点结构的瞬态动力学分析

【JY】从一根悬臂梁说起

【JY】反应谱的详解与介绍

【JY】结构瑞利阻尼与经济订货模型

【JY】主成分分析与振型分解

【JY】浅谈结构多点激励之概念机理(上)

【JY】浅谈结构多点激励之分析方法(下)

【JY】板壳单元的分析详解

【JY】橡胶支座的简述和其力学性能计算

【JY】振型求解之子空间迭代

【JY】橡胶支座精细化模拟与有限元分析注意要点

【JY】推开土木工程振型求解之兰索斯法(Lanczos法)的大门

【JY】基于OpenSees和Sap2000静力动力计算案例分析

【JY】建筑结构施加地震波的方法与理论机理

【JY】力荐佳作《结构地震分析编程与应用》

【JY】ABAQUS子程序UEL的有限元原理与应用

【JY】钢筋混凝土正截面极限承载力设计的基本原理和快速计算方法

【JY】消能减震黏滞阻尼器的力学原理与应用

【JY】力荐 | 区域建筑地震安全性有限元分析示例

【JY】ABAQUS正交各向异性弹性本构模型

# 软件讨论

【JY】浅析时程分析中的阻尼设置

【JY】减隔震元件计算表格分享

【JY】复合材料分析利器—内聚力单元

【JY】SDOF计算教学软件开发应用分享

【JY】Abaqus案例—天然橡胶隔震支座竖(轴)向力学性能

【JY】Abaqus6.14-4如何关联fortran?

【JY】如何利用python来编写GUI?

【JY】如何解决MATLAB GUI编程软件移植运行问题?

【JY】浅谈结构分析与设计软件

【JY|STR】求解器之三维结构振型分析

【JY】SignalData软件开发应用分享

【JY】基于Matlab的双线性滞回代码编写教程

【JY】动力学利器 —— JYdyn函数包分享与体验

【JY】混凝土分析工具箱:CDP模型插件与滞回曲线数据

【JY】结构工程分析软件讨论(上)

【JY】结构工程分析软件讨论(下)

【JY】各类有限元软件计算功能赏析与探讨

【JY】ETABS与Perform3D弹塑性分析功能对比示例

【JY】自建9000条高质量人工波可供下载

【JY】代码|极简反应谱分析

#其他

【JY】位移角还是有害位移角?

【JY】如何利用python来编写GUI?

【JY】今日科普之BIM

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值