[发明专利]一种考虑结构对称性的张拉整体结构非线性找形方法在审
申请号: | 201810081302.4 | 申请日: | 2018-01-26 |
公开(公告)号: | CN108345728A | 公开(公告)日: | 2018-07-31 |
发明(设计)人: | 陈耀;孙求知;冯健 | 申请(专利权)人: | 东南大学 |
主分类号: | G06F17/50 | 分类号: | G06F17/50;G06F17/16 |
代理公司: | 南京苏高专利商标事务所(普通合伙) 32204 | 代理人: | 柏尚春 |
地址: | 211100 江*** | 国省代码: | 江苏;32 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种考虑结构对称性的张拉整体结构非线性找形方法,针对任意张拉整体结构,通过分析结构对称坐标系下平衡矩阵第一分块子矩阵,综合运用矩阵舒尔分解、奇异值分解理论求得一组节点坐标和整体可行自应力模态。主要步骤为:首先输入结构拓扑信息,明确结构所属对称群,求出对称坐标系下与结构外荷载向量和杆件内力向量相关联的全对称子空间,定义一组初始力密度,形成初始力密度矩阵,计算对称坐标系下平衡矩阵第一分块子矩阵,求解该分块子矩阵零空间,并求得整体自应力模态,更新力密度矩阵,通过迭代最终得到一组整体可行自应力模态,并求出满足力密度矩阵秩亏条件和结构稳定条件的节点坐标和杆件力密度。 | ||
搜索关键词: | 矩阵 张拉整体结构 分块子矩阵 自应力模态 密度矩阵 结构对称性 杆件 向量 对称 奇异值分解 节点坐标 结构对称 结构稳定 输入结构 拓扑信息 荷载 对称群 零空间 全对称 子空间 组节点 内力 求解 迭代 平衡 秩亏 关联 分解 更新 分析 | ||
【主权项】:
1.一种考虑结构对称性的张拉整体结构非线性找形方法,其特征在于,该方法包括以下步骤:步骤1基于待求解结构的拓扑信息,形成拓扑矩阵C,明确结构所属对称群;步骤2求出与结构外荷载向量相关联的全对称子空间
和与杆件内力向量相关联的全对称子空间
步骤3根据结构压杆受压为正、拉索受拉为负定义一组初始力密度q0,通过公式D0=CT*diag(q0)*C求出初始力密度q0对应的力密度矩阵D0,令迭代计数变量i=0,其中,diag()表示矩阵的对角线元素为括号中的元素,其余元素为0,CT为拓扑矩阵转置;步骤4对第i次的力密度矩阵Di进行舒尔分解,并利用矩阵舒尔分解剔除使结构所有节点某一方向坐标取值均为同一数值的解,得出一组坐标,并求出该组坐标对应的平衡矩阵Ai;步骤5通过下式计算对称坐标系下平衡矩阵Hi的第一分块子矩阵![]()
其中
表示
的转置;步骤6对所述对称坐标系下平衡矩阵Hi的第一分块子矩阵
进行奇异值分解,求出其零空间Ni,通过下式求出结构整体自应力模态Si:
步骤7通过步骤6得出的结构整体自应力模态求出新的杆件力密度qi+1,通过公式Di+1=CT*diag(qi+1)*C更新力密度矩阵,令迭代计步变量i=i+1,其中Di+1为更新后的力密度矩阵;步骤8判断所述步骤7中更新的力密度矩阵Di+1是否满足秩亏要求以及待求解结构的切线刚度矩阵是否正定,若上述两个条件均满足,则本方法流程结束;若不满足,则返回步骤4。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于东南大学,未经东南大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201810081302.4/,转载请声明来源钻瓜专利网。