[发明专利]基于刚体准则的平面框架结构静力非线性分析方法在审
申请号: | 201811119408.5 | 申请日: | 2018-09-25 |
公开(公告)号: | CN109299540A | 公开(公告)日: | 2019-02-01 |
发明(设计)人: | 杨永斌;陈朝晖;陶宇宸;李云飞;古洪;廖旻懋 | 申请(专利权)人: | 重庆大学 |
主分类号: | G06F17/50 | 分类号: | G06F17/50 |
代理公司: | 重庆鼎慧峰合知识产权代理事务所(普通合伙) 50236 | 代理人: | 周维锋 |
地址: | 400044 *** | 国省代码: | 重庆;50 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 刚体 平面框架结构 单元刚度矩阵 材料非线性 几何非线性 计算效率 耦合问题 静力 荷载位移曲线 非线性分析 梁单元模型 弹簧平面 刚度矩阵 工程应用 划分单元 极限荷载 几何刚度 矩阵形式 求解方程 线性矩阵 增量迭代 增量方程 塑性铰 求解 迭代 加载 更新 分析 | ||
1.一种基于刚体准则的平面框架结构静力非线性分析方法,其特征在于包括以下步骤:
步骤1、建立集中塑性铰弹簧平面梁单元模型;
步骤2、建立满足刚体准则的上述模型对应的单元刚度矩阵;
步骤3、基于更新拉格朗日格式与刚体准则,建立求解方程的增量迭代过程;
步骤4、迭代求解增量方程,得到结构加载过程的荷载位移曲线与极限荷载。
2.根据权利要求1所述的基于刚体准则的平面框架结构静力非线性分析方法,其特征在于,建立集中塑性铰平面梁单元模型,框架梁单元的屈服发生于零长度的塑性铰内,塑性铰之间的杆件为完全弹性。
3.根据权利要求1所述的基于刚体准则的平面框架结构静力非线性分析方法,其特征在于利用轴力与弯矩耦合的截面精细初始屈服与全屈服曲线,建立框架梁单元截面逐步屈服塑性铰弹簧刚度模型,塑性铰弹簧的刚度系数为
其中,EI/L为梁的弹性弯曲刚度,初始屈服函数为完全屈服函数为Fp为截面的轴向承载力,Mp为塑性极限弯矩,Fx和M分别为单元的轴力和弯矩。
4.根据权利要求1所述的基于刚体准则的平面框架结构静力非线性分析方法,其特征在于采用以下单元刚度矩阵:
[k]=[kg]+[kep]
其中,[kg]为单元几何刚度矩阵,[kep]为单元弹塑性刚度矩阵。
5.根据权利要求4所述的基于刚体准则的平面框架结构静力非线性分析方法,其特征在于利用刚体准则,即:对处于初始平衡状态的刚体单元,当发生刚体移动和转动时,由于单元没有发生变形,单元在原结点力作用下应保持平衡,原单元结点力的大小不变,力的方向随单元转动而转动,利用刚体准则和虚功原理建立的单元几何刚度矩阵如下:
当单元发生刚体位移{u}r时,单元虚应变能为零,几何刚度矩阵不会产生结点力增量,即:[kg]{u}r=0,其中{u}rT=<0 0 θr 0 Lθr θr>。
6.根据权利要求4所述的基于刚体准则的平面框架结构静力非线性分析方法,其特征在于:利用截面逐步屈服塑性铰框架梁单元结点弯矩的增量平衡方程,建立单元弹塑性刚度矩阵,形式如下
其中,
k36=2EIzSaSbL/(12E2Iz2+4EIzL(Sa+Sb)+SaSbL2);
k66=EIzSb(12EIz+4SaL)/(12E2Iz2+4EIzL(Sa+Sb)+SaSbL2);
其中,E为弹性模量,A为单元截面面积,L为单元长度,Iz为截面惯性矩,Sa为单元a端的弹簧刚度,Sb为单元b端的弹簧刚度。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于重庆大学,未经重庆大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201811119408.5/1.html,转载请声明来源钻瓜专利网。