发明名称 一种可评估结构非概率可靠度的安全系数设计方法
摘要 本发明公开了一种可评估结构非概率可靠度的安全系数设计方法,首先利用结构参数的有限样本数据,利用非统计度量方法合理表征不确定参数的非概率特性;近似求解、Monte Carlo模拟法等确定结构应力、强度的不确定分布特性;构建结构的应力‑强度非概率集合干涉模型,提出相应的非概率集合理论可靠度评估公式;结合非概率集合理论可靠度评估公式,建立与非概率集合理论可靠度、强度变异系数以及应力变异系数等相关的安全系数的解析表达式;针对工程结构,根据相应可靠度要求计算对应的安全系数,进行安全系数结构设计,获得满足一定可靠度的结构最优设计方案。本发明兼顾了设计本身安全性和经济性,保留了传统安全系数设计方法概念简单,容易实行的特点。
申请公布号 CN105844060A 申请公布日期 2016.08.10
申请号 CN201610255291.8 申请日期 2016.04.22
申请人 北京航空航天大学 发明人 邱志平;陈潇;王晓军;王磊;孙佳丽
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 成金玉;卢纪
主权项 一种可评估结构非概率可靠度的安全系数设计方法,其特征在于实现步骤如下:第一步:将结构参数,包括载荷、材料弹性常数、强度指标以及结构尺寸的有限样本数据写成原始数据矩阵<img file="FDA0000971560480000011.GIF" wi="588" he="295" />其中x<sub>1</sub>(1),x<sub>1</sub>(2),…x<sub>m</sub>(n)是来源于试验或等精度测量的原始数据,m为结构参数的总个数,n为每个参数样本数据的个数;利用非统计度量方法对有限样本数据进行筛选与评估,得到结构参数的合理不确定区间表征向量<img file="FDA0000971560480000012.GIF" wi="443" he="63" />其中<img file="FDA0000971560480000013.GIF" wi="52" he="70" />为第p个结构参数的区间表征,I表示区间,T表示矩阵转置;第二步:利用第一步得到的区间表征向量<img file="FDA0000971560480000014.GIF" wi="447" he="63" />进行区间运算,得到结构参数的不确定信息,包括载荷、材料弹性常数、强度指标以及结构尺寸的上、下界与中心值、区间半径的关系式,有:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>X</mi><mi>U</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mn>1</mn><mi>U</mi></msubsup><mo>,</mo><msubsup><mi>x</mi><mn>2</mn><mi>U</mi></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>x</mi><mi>m</mi><mi>U</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mn>1</mn><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>x</mi><mn>1</mn><mi>r</mi></msubsup><mo>,</mo><msubsup><mi>x</mi><mn>2</mn><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>x</mi><mn>2</mn><mi>r</mi></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>x</mi><mi>m</mi><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>x</mi><mi>m</mi><mi>r</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup></mrow>]]></math><img file="FDA0000971560480000015.GIF" wi="1061" he="71" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>X</mi><mi>L</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mn>1</mn><mi>L</mi></msubsup><mo>,</mo><msubsup><mi>x</mi><mn>2</mn><mi>L</mi></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>x</mi><mi>m</mi><mi>L</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mn>1</mn><mi>c</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mn>1</mn><mi>r</mi></msubsup><mo>,</mo><msubsup><mi>x</mi><mn>2</mn><mi>c</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mn>2</mn><mi>r</mi></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>x</mi><mi>m</mi><mi>c</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>m</mi><mi>r</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup></mrow>]]></math><img file="FDA0000971560480000016.GIF" wi="1045" he="63" /></maths>其中X<sup>U</sup>为结构参数的上界表达式,X<sup>L</sup>为结构参数的下界表达式,上标U代表变量的取值上界,上标L代表变量的取值下界;上标c代表区间中心值,上标r代表区间半径;第三步:将第二步得到的结构参数的不确定信息引入到工程结构应力S和强度R的计算表达式S(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>m</sub>)和R(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>m</sub>)中,引入非概率不确定传播理论和方法,确定结构应力S、强度R的不确定区间<img file="FDA0000971560480000017.GIF" wi="117" he="62" />和<img file="FDA0000971560480000018.GIF" wi="125" he="62" />及不确定分布特性S<sup>c</sup>,S<sup>r</sup>和R<sup>c</sup>,R<sup>r</sup>,有:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>S</mi><mi>c</mi></msup><mo>=</mo><mfrac><mrow><mover><mi>S</mi><mo>&OverBar;</mo></mover><mo>+</mo><munder><mi>S</mi><mo>&OverBar;</mo></munder></mrow><mn>2</mn></mfrac><mo>,</mo><msup><mi>S</mi><mi>r</mi></msup><mo>=</mo><mfrac><mrow><mover><mi>S</mi><mo>&OverBar;</mo></mover><mo>-</mo><munder><mi>S</mi><mo>&OverBar;</mo></munder></mrow><mn>2</mn></mfrac><mo>,</mo><msup><mi>R</mi><mi>c</mi></msup><mo>=</mo><mfrac><mrow><mover><mi>R</mi><mo>&OverBar;</mo></mover><mo>+</mo><munder><mi>R</mi><mo>&OverBar;</mo></munder></mrow><mn>2</mn></mfrac><mo>,</mo><msup><mi>R</mi><mi>r</mi></msup><mo>=</mo><mfrac><mrow><mover><mi>R</mi><mo>&OverBar;</mo></mover><mo>-</mo><munder><mi>R</mi><mo>&OverBar;</mo></munder></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000971560480000019.GIF" wi="901" he="126" /></maths>其中<u>S</u>,<u>R</u>分别为结构应力S和强度R的下界;<img file="FDA00009715604800000110.GIF" wi="94" he="63" />分别为结构应力S和强度R的上界,S<sup>c</sup>,R<sup>c</sup>分别为结构应力S和强度R的中心值,S<sup>r</sup>,R<sup>r</sup>分别为结构应力、强度的区间半径;第四步:根据第三步得到的结构应力S、强度R的不确定分布特性,包括中心值S<sup>c</sup>,R<sup>c</sup>和区间半径S<sup>r</sup>,R<sup>r</sup>,建立非概率可靠度的功能函数方程:M(R,S)=R‑S对结构应力S、强度R进行标准区间变换,得到标准化变量空间下的极限状态方程为:M(δ<sub>R</sub>,δ<sub>S</sub>)=R<sup>c</sup>‑S<sup>c</sup>+R<sup>r</sup>δ<sub>R</sub>‑S<sup>r</sup>δ<sub>S</sub>=0其中δ<sub>R</sub>∈[‑1,1],δ<sub>S</sub>∈[‑1,1]为标准化的强度R、应力S区间变量;利用应力‑强度非概率集合干涉模型,得到非概率集合理论可靠度R<sub>集合</sub>:<img file="FDA0000971560480000021.GIF" wi="789" he="317" />第五步:根据第三步得到的结构应力S、强度R的不确定分布特性,包括中心值S<sup>c</sup>,R<sup>c</sup>和区间半径S<sup>r</sup>,R<sup>r</sup>,定义区间理论下的中心安全系数n<sub>m</sub>为结构强度R、应力S的中心值之比,即n<sub>m</sub>=R<sup>c</sup>/S<sup>c</sup>;定义强度变异系数C<sub>R</sub>和应力变异系数C<sub>S</sub>分别为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>C</mi><mi>R</mi></msub><mo>=</mo><mfrac><msup><mi>R</mi><mi>r</mi></msup><msup><mi>R</mi><mi>c</mi></msup></mfrac><mo>,</mo><msub><mi>C</mi><mi>S</mi></msub><mo>=</mo><mfrac><msup><mi>S</mi><mi>r</mi></msup><msup><mi>S</mi><mi>c</mi></msup></mfrac></mrow>]]></math><img file="FDA0000971560480000022.GIF" wi="356" he="119" /></maths>得到区间理论下的中心安全系数n<sub>m</sub>、强度变异系数C<sub>R</sub>、应力变异系数C<sub>S</sub>及非概率集合理论可靠度R<sub>集合</sub>的关系式:<img file="FDA0000971560480000023.GIF" wi="990" he="135" />第六步:针对工程结构选择合适的可靠度R<sub>集合</sub>;利用得到的强度变异系数C<sub>R</sub>、应力变异系数C<sub>S</sub>,求解与非概率集合理论可靠度R<sub>集合</sub>对应的中心安全系数n<sub>m</sub>,并进行安全系数结构设计,获得满足一定可靠度的最优设计方案。
地址 100191 北京市海淀区学院路37号
您可能感兴趣的专利