CN107066663-CN201611252501-一种基于满应力约束准则的桁架结构非概率可靠性拓扑优化方法.pdfVIP

(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 CN 107066663 A (43)申请公布日 2017.08.18 (21)申请号 201611252501.4 (22)申请日 2016.12.30 (71)申请人 北京航空航天大学 地址 100191 北京市海淀区学院路37号 (72)发明人 王磊　刘易斯　王晓军　蔡逸如　 刘东亮　夏海军　耿新宇　 (74)专利代理机构 北京科迪生专利代理有限责 任公司 11251 代理人 杨学明　顾炜 (51)Int.Cl. G06F 17/50(2006.01) 权利要求书3页 说明书7页 附图6页 (54)发明名称 一种基于满应力约束准则的桁架结构非概 率可靠性拓扑优化方法 (57)摘要 本发明公开了一种基于满应力约束准则的 桁架结构非概率可靠性拓扑优化方法。首先建立 一个桁架结构的基结构，即在设计域内所有的节 点上都连上杆件。利用顶点组合法得到应力的上 下界，再基于非概率集合可靠性指标，得到应力 的非概率集合可靠度，将其作为约束条件带入更 新设计变量的迭代方程中，使得桁架结构各杆的 横截面积不断变化，直至达到收敛条件为止。各 杆的横截面积在变化的过程中需要满足如下的 条件：如果某单元没有达到其尺寸容许的上下 限，那么该单元至少在某一工况下达到容许应力 A 的可靠性上限，也即满可靠性应力状态。收敛条 3 件为：如果前后两次加载各杆的横截面积变化小 6 6 6 于容差，即可认为收敛。本发明基于该方法拓扑 6 0 7 优化得到的结构更为安全。 0 1 N C CN 107066663 A 权　利　要　求　书 1/3 页 1.一种基于满应力约束准则的桁架结构非概率可靠性拓扑优化方法，其特征在于，该 方法实现步骤如下： 第一步：推导确定性的桁架结构横截面积的迭代公式，首先根据桁架结构的一般特点， 建立含优化模型的数学表达式，基于满应力准则的桁架结构拓扑优化的数学表达列式如 下： s.t.: 0＜A≤A ,i＝1,...,N i Bf＝F 其中i为杆的编号，N为基结构中杆件的总数目，ρ为材料密度，A为杆的截面积，A为杆横 截面积的下界，f为各杆所受的轴力，σ为各杆所受的应力， 为各杆所能承受的极限应力，B 为角度矩阵，F为在直角坐标系下各节点受的力，以上方程组的含义为：在满足应力不超过 极限应力、横截面积不超过容许上下限以及满足力平衡的条件下，桁架结构的总质量最小， 基于以上的优化模型方程组建立拉格朗日方程，方程如下所示： 其中α ,μ ,λ是拉格朗日乘子，将拉格朗日方程对其每一个变量求导可以得到KKT条件， i i 其具体表述如下方程组所示： μ(A-A)＝0,μ ≥0,i＝1,...,N i i i 0＜A≤A ,i＝1,...,N i Bf＝F 再对各杆的横截面积以及应力状况进行分析。假设第i根杆的横截面积大于它的下界， 即A＜A ，根据以上的KKT条件我们可以得到μ ＝0，KKT条件中的第二个方程可以做如下简 i i 化： 进一步假设第i根杆的应力小于它的极限应力，即 可以得到α ＝0，KKT条件