[发明专利]一种结构损伤识别的不确定性定量分析方法及装置有效
| 申请号: | 202110634293.9 | 申请日: | 2021-06-08 |
| 公开(公告)号: | CN113297773B | 公开(公告)日: | 2023-01-24 |
| 发明(设计)人: | 吕中荣;尹智毅;汪利;刘济科 | 申请(专利权)人: | 中山大学 |
| 主分类号: | G06F30/23 | 分类号: | G06F30/23 |
| 代理公司: | 广州三环专利商标代理有限公司 44202 | 代理人: | 郭浩辉;颜希文 |
| 地址: | 510275 广东*** | 国省代码: | 广东;44 |
| 权利要求书: | 查看更多 | 说明书: | 查看更多 |
| 摘要: | |||
| 搜索关键词: | 一种 结构 损伤 识别 不确定性 定量分析 方法 装置 | ||
1.一种结构损伤识别的不确定性定量分析方法,其特征在于,包括:
建立结构的有限元模型,并根据所述有限元模型获得所述结构的计算特征数据,所述计算特征数据包括计算频率和计算模态;
采用Tikhonov正则化方法,结合所述结构的计算特征数据以及所述结构的实测特征数据,构造目标函数,其中,所述实测特征数据包括实测频率和实测模态,所述目标函数为:
其中,ΔR(α)=RM-RC(α)为所述结构的实测特征数据与所述结构的计算特征数据的差值项;RM为所述结构的实测特征数据,是一个特征数据矩阵,表示所述结构的实测频率,表示所述结构的实测模态,n为选用所述实测频率和所述实测模态的阶次;RC为所述结构的计算特征数据,是一个特征数据矩阵,表示所述结构的计算频率,表示所述结构的计算模态,m为选用所述计算频率和所述计算模态的阶次,m=n;W为权重矩阵,为引入的Tikhonov正则化约束项,λ为通过L曲线法选取的Tikhonov正则化系数,α为引入的损伤参数;
采用灵敏度方法求解所述目标函数,得到最优损伤参数和损伤参数的协方差矩阵,其中,所述最优损伤参数为:
式中,为上一次迭代的灵敏度矩阵,为损伤参数α的上一次迭代值,为的转置矩阵,I为单位矩阵,所述上一次迭代的结构的实测特征数据与所述结构的计算特征数据的差值项;
所述损伤参数的协方差矩阵为:
∑aR≈[ST(α*)WS(α*)+λI]-1
式中,S(α*)为最优损伤参数的灵敏度矩阵,ST(α*)为S(α*)的转置矩阵;
根据所述最优损伤参数和所述损伤参数的协方差矩阵进行不确定性定量分析,得到分析结果。
2.如权利要求1所述的结构损伤识别的不确定性定量分析方法,其特征在于,所述建立结构的有限元模型,具体为:
采用有限元方法,将所述结构离散为N个单元,以根据所述结构的动力学特征方程建立所述有限元模型。
3.如权利要求1所述的结构损伤识别的不确定性定量分析方法,其特征在于,所述采用灵敏度方法求解所述目标函数,得到最优损伤参数和损伤参数的协方差矩阵,具体为:
采用灵敏度方法迭代求解所述目标函数直至满足预设精度条件,得到所述最优损伤参数;
根据所述目标函数的海森矩阵,得到所述损伤参数的协方差矩阵。
4.如权利要求3所述的结构损伤识别的不确定性定量分析方法,其特征在于,所述采用灵敏度方法迭代求解所述目标函数直至满足预设精度条件,得到所述最优损伤参数,具体为:
采用灵敏度方法迭代求解所述目标函数,得到当前损伤参数,并计算所述当前损伤参数与上一损伤参数的相对误差;
在所述相对误差小于预设容许值时,将所述当前损伤参数作为所述最优损伤参数,结束迭代求解所述目标函数。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中山大学,未经中山大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202110634293.9/1.html,转载请声明来源钻瓜专利网。
- 上一篇:一种利用超声测量螺栓载荷的方法
- 下一篇:处理方法、移动终端及存储介质





