[发明专利]一种基于伴随变量法的含有多不确定参数的结构区间响应传播分析方法有效
申请号: | 201711349714.3 | 申请日: | 2017-12-15 |
公开(公告)号: | CN108121865B | 公开(公告)日: | 2018-12-21 |
发明(设计)人: | 王晓军;罗振先;王磊;任强;倪博文;丁旭云 | 申请(专利权)人: | 北京航空航天大学 |
主分类号: | G06F17/50 | 分类号: | G06F17/50 |
代理公司: | 北京科迪生专利代理有限责任公司 11251 | 代理人: | 杨学明;卢纪 |
地址: | 100191*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种基于伴随向量法的含有多不确定参数的结构区间响应传播分析方法,本发明对结构区间响应进行Taylor级数展开,在不确定参数区间内得到系统响应的一阶近似模型,基于此进而求得区间响应的上下界。针对含有多不确定参数的系统,在Taylor级数展开时,基于伴随变量法求解响应关于不确定参数的偏导数,可以减少系统重分析的次数得到响应关于所有不确定参数的偏导数,进而确定系统响应不确定区间的上下界。该方法在保证计算精度的前提下,极大地提高了计算效率。因此可高效地预测结构的响应边界。 | ||
搜索关键词: | 响应 级数展开 系统响应 偏导数 上下界 参数区间 计算效率 减少系统 近似模型 预测结构 分析 向量法 求解 一阶 传播 保证 | ||
【主权项】:
1.一种基于伴随变量法的含有多不确定参数的结构区间响应传播分析方法,其特征在于,该方法实现步骤如下:第一步:建立结构的有限元模型,包括建立几何模型,赋予材料属性,施加荷载和边界条件;第二步:确定结构中的不确定性参数及其不确定区间,可记作
其中α和
分别是不确定区间的下、上界;第三步:建立结构的静力学平衡方程:K(α)U=F(α)其中,
为结构不确定性参数向量,含有m个分量αi(i=1,2,...,m),刚度矩阵K和载荷F都是关于不确定参数的函数,故结构位移U也是不确定性参数的函数,可表示为:
第四步:基于Taylor级数展开法对结构位移进行一阶近似,可以得到位移关于不确定性参数函数的近似表达:
这里,![]()
其中,
为结构位移关于不确定性参数的函数,α0是不确定性参数的中心值向量,
是结构位移函数
在不确定性参数取中心值α0时的值,δ是不确定性参数的半径向量;g是结构位移函数关于各不确定性参数在中心值处的偏导数向量;第五步:计算步骤四中位移函数
对不确定量α的偏导,因此可对步骤三中的静力平衡方程两侧求偏导得到:
进一步整理可得:
第六步:求解响应对不确定性参数的偏导,在静力学分析中,位移、应力以及载荷的第j个响应φj可以表示成位移的函数
则响应φj对不确定性参数的偏导数为:
将步骤五中的方程代入上述方程中,可得到,
其中,
第七步:计算结构响应的区间半径ΔU,根据一阶Taylor级数展开公式可以得到结构响应的一阶增量:
第八步:计算结构响应区间的上下界
由步骤三中的静力平衡方程求得结构响应的中心值Uc,结合步骤七中求得的结构响应的区间半径ΔU求得结构响应区间的上下界:U=Uc‑ΔU,
第九步:结束。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京航空航天大学,未经北京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201711349714.3/,转载请声明来源钻瓜专利网。