一种含不确定性金属结构声振疲劳寿命的估算方法技术

本发明专利技术公开了一种含不确定性金属结构声振疲劳寿命的估算方法。该方法首先建立金属结构的有限元分析模型，考虑有限样本条件下材料与结构特性参数的不确定性效应，基于频率响分析理论得到系统的传递函数；然后根据随机噪声载荷的特点，针对典型宽带平稳各态历经的随机噪声载荷，用平方和开方(SRSS)方法得到各个关键结点位移、应力的均方根响应与应力功率谱密度函数(PSD)，再结合Dirlik模型与区间不确定传播二阶泰勒级数展开法与频域内振动疲劳的分析方法，建立雨流幅值概率密度函数和功率谱密度函数之间的关系，最后应用Miner线性累计损伤理论得到含不确定性金属结构的声振疲劳寿命区间范围。本发明专利技术在计算受典型随机噪声激励的结构寿命时充分考虑了结构与材料参数的分散性，因此得到的疲劳寿命更加合理。

【技术实现步骤摘要】
一种含不确定性金属结构声振疲劳寿命的估算方法
本专利技术涉及金属结构噪声振动疲劳
，特别涉及考虑不确定性作用下结构对于随机噪声载荷的动响应与声疲劳寿命估算方法。适用于在典型复杂噪声载荷下各种结构类型的声疲劳寿命高效预测，为工程抗声疲劳结构设计提供帮助。
技术介绍
噪声载荷对结构的作用本质是一种空间分布的，具有一定频率分布特性的动态随机压力载荷。当声压级量级超过140dB，它可以在结构上产生一定的分布应力响应，特别是当噪声的频率分布特性和它所作用的结构的动特性互相耦合时，结构便会产生显著的应力响应。声疲劳属于高周低应力问题，在这种动态应力的长时间作用下，结构的应力集中和缺陷部位便会形成疲劳裂纹，在工程上当裂纹成核并扩展到可检长度的阶段后即产生了疲劳裂纹形成寿命即所要研究的声振疲劳寿命。飞行器在服役过程中处于发动机喷气、气动力与壁面耦合作用、机体自身结构工作振动等产生的复杂噪声环境。随着飞行器飞行速度的提高、动力装置推力增大和作战要求性能的改善，飞行器运行中产生的各种噪声量值不断地增大。实践表明，不论是军用飞机还是民用飞机在使用中都会常常出现各种类型的声疲劳破坏现象。其中大多数表现为：各种翼面蒙皮及机身侧蒙皮裂纹、翼肋和机身环框裂纹、进气道内蒙皮裂纹、机尾在喷流作用下的各种破坏。传统的声疲劳寿命估算存在很多局限性，比如：在考虑动态应力效应时只考虑了低阶频率，有的只考虑了一阶固有频率，对频率分布宽能量高的载荷引起的效应有一定误差；在初步设计中设计者需清楚对象所用的材料和结构型式与已有曲线的是否相同，相同则可在得到均方根值应力后引用特定的结构型式及材料的声疲劳试验曲线来得到寿命，但不相同则需要重新进行结构参数设计，按照现有的曲线来等效转化与结果经验修正，此举工作量大且精度低；特定结构及材料的均方根值寿命曲线本身就很缺乏，目前已有的数量少，需要获得大量的试验数据支持。另外，制造加工工艺及材料非均质性所造成的初始缺陷和损伤不可避免，并在未来长期服役过程中于结构内部不断发展、蔓延、传播，严重影响着结构的力学行为及使用安全。建立以非概率理论框架为基础的不确定性表征技术、不确定性金属结构声振疲劳寿命估算技术具有显著的现实意义。
技术实现思路
本专利技术要解决的技术问题是：克服现有技术的不足，提供一种针对含不确定性金属结构声振疲劳寿命的估算方法。本专利技术充分考虑实际工程问题中普遍存在的不确定性因素，以非概率区间方法表征不确定量，引入不确定区间传播理论，所得到的设计结果更加符合真实情况，工程适用性更强。本专利技术采用的技术方案为：一种含不确定性金属结构声振疲劳寿命的估算方法，实现步骤如下：第一步：依据飞行任务剖面时间数据和状态声载荷测量带宽声压级，采用将其转换为随机噪声的频谱载荷，其中L＝Lb-10lg(Δf)，P0＝2×10-5Pa为参考声压，Lb为带宽声压级，Δf为频带宽度，G(f)为转化后的噪声功率谱密度值；第二步：引入区间向量x∈xI＝[E,a,b]合理定量贫信息、少数据条件下的结构的不确定参数，其中E为金属材料弹性模量，a与b分别表示结构不同部位的几何参数。则含结构参数的不确定性可以表示为：代表参量的取值上界，x＝[E,a,b]＝[Ec-Er,ac-ar,bc-br]，x代表参量的取值下界，其中上标c代表中心值，上标r代表半径；第三步：建立所关心结构的几何模型，分析构件的连接形式，合理划分网格施加边界条件形成有限元分析模型。在有限元分析软件中加载单位均匀面压力载荷，其中载荷的频率变化范围与第一步所得的噪声载荷的频率范围相同，对模型进行频率响应分析后提取结构各个结点的应力传递函数；第四步：应用平方和开方(SRSS)方法，在MSC.Patran随机振动模块中读取第三步得到的应力传递函数，并施加第一步所得的频谱噪声载荷得到结点的应力功率谱密度(PSD)曲线，其中随机振动分析输出文件还包括频率响应的功率谱密度、自相关函数、每单位时间的正斜率方向的零交点个数及应力响应的RMS值；第五步：随机噪声的典型特点是压力时间历程幅值是随机变化的，也就是无规律非衰减不能用解析函数表达的形式，这类噪声频率的变化范围宽，其谱值直到很高频率都是连续的。结合典型宽带随机过程的Dirlik模型与区间传播分析方法，以上一步得到的应力PSD曲线为输入，应用非概率区间过程理论中的泰勒级数展开法得到不确定性变量传播后各个关键结点的应力雨流幅值概率密度函数(PDF)变化范围；第六步：对于连续分布应力状态，将时间T内在应力范围(Si,Si+ΔSi)内的应力循环次数表示为ni＝vTp(Si)ΔSi。式中的v表示单位时间内的应力循环次数，由每秒的峰点数E[P]决定，即v＝E[P]，p(Si)表示应力级水平为Si时的幅值概率密度函数值，ΔSi为微小应力级变化范围；第七步：在MSC.Fatigue中选择对应金属材料的S-N曲线，其表达式为N(Si)＝K/Sm，应用Miner线性累计损伤理论D＝∑Di＝∑ni/Ni，当总体损伤度D＝1时得到构件发生疲劳破坏的时间疲劳寿命为：其中，K与m为材料常数，N为金属材料循环次数，D为损伤度。最终计算得到构件发生疲劳破坏时的时间寿命分布云图。根据非概率区间传播分析理论得到构件的寿命时间范围。进一步的，所述第二步中区间不确定性参数向量x可以表示为：其中，xc＝(Ec,ac,bc)，xr＝(Er,ar,br)，e∈Ξ3，Ξ3定义为所有元素包含在[-1,1]内的3维向量集合，符号“×”定义为两个向量各对应元素，乘积仍为3维向量。进一步的，所述第五步中应用的Dirlik模型适用于典型宽带随机过程的模拟，对于窄带或者带有尖峰(即脉冲幅值)特性的随机过程其结果准确性供参考。进一步的，所述第五步中非概率区间传播分析方法使用的是二阶泰勒级数展开法，其响应区间的上下界可表示为：其中，与Ψi表示i个不确定量分别取上下界与xi的响应结果的上下界，Ψi(xc)表示不确定变量取中值xc时结构所关心量响应值，表示Ψi(xc)在第i个不确定参数中点的一阶导数展开值，表示不确定变量区间中点处的二阶展开，Δxi与Δxj分别表示第i个与第j个不确定变量的区间半径。本专利技术与现有技术相比的优点在于：本专利技术提供了声疲劳寿命估算的新思路，弥补和完善了传统理论的局限性。传统估算声疲劳的方法只能针对特定的结构形式，此方法能计算任意结构形式；传统方法需要一条均方根值应力寿命曲线，此方法只需要传统疲劳寿命曲线。在考虑材料与结构分散性时该方法不需要知道不确定参数的概率分布形式，只需要知道金属材料以及结构几何参数的上下界限就可以方便的求解含不确定性金属材料任意结构形式的声疲劳寿命区间范围，在工程实用方面更加方便可靠。附图说明图1是本专利技术针对含不确定性金属结构声振疲劳寿命的估算方法流程图；图2是本专利技术使用的有限元软件分析的流程示意图；图3是本专利技术针对飞行器空腔结构进行简化并施加边界条件后的有限元模型；图4是本专利技术计算得到模型的关键结点的功率谱密度函数示意图；图5是本专利技术针对7075-HV-T6金属材料引用的S-N曲线；图6是本专利技术计算得到的金属材料声振疲劳寿命分布云图；图7是本专利技术对含不确定参数的结构得到的结点最小寿命区间范围。具体实施方式下面结合附图以及具体实施例进一步说明本专利技术。本文档来自技高网...

【技术保护点】
一种含不确定性金属结构声振疲劳寿命的估算方法，其特征在于实现步骤如下：第一步：依据飞行任务剖面时间数据和状态声载荷测量带宽声压级，采用将其转换为频谱载荷，其中L＝Lb‑10lg(△f)，P0＝2×10‑5Pa为参考声压，Lb为带宽声压级，△f为频带宽度，G(f)为转化后的噪声功率谱密度值；第二步：引入区间向量x∈xI＝[E,a,b]合理定量贫信息、少数据条件下结构的不确定参数，其中E为金属材料弹性模量，a与b分别表示结构的几何参数，含结构参数的不确定性可以表示为：x‾=[E‾,a‾,b‾]=[Ec+Er,ac+ar,bc+br],]]>代表参量的取值上界，x＝[E,a,b]＝[Ec‑Er,ac‑ar,bc‑br]，x代表参量的取值下界，其中上标c代表中心值，上标r代表半径；第三步：建立所关心结构的几何模型，分析构件的连接形式，合理划分网格施加边界条件形成有限元分析模型，在有限元分析软件中加载单位均匀面压力载荷，其中载荷的频率变化范围与第一步所得的噪声载荷的频率范围相同，对模型进行频率响应分析后提取结构各个结点的应力传递函数；第四步：应用平方和开方(SRSS)方法，在MSC.Patran随机振动模块中读取第四步得到的系统传递函数，并施加第一步所得的频谱噪声载荷得到结点的应力功率谱密度(PSD)曲线，其中随机振动分析输出文件还包括频率响应的自相关函数、每单位时间的正斜率方向的零交点个数及应力响应的均方值(RMS值)；第五步：随机噪声的典型特点是压力时间历程幅值是随机变化的，也就是无规律非衰减不能用解析函数表达的形式，这类噪声频率的变化范围宽，其谱值直到很高频率都是连续的，结合模拟此类典型宽带随机过程的Dirlik模型与区间传播分析方法，以上一步得到的应力PSD曲线为输入，引入非概率区间过程理论中泰勒展开分析法得到不确定变量扩张后各个关键结点的应力雨流幅值概率密度函数(PDF)变化范围；第六步：对于连续分布应力状态，将时间T内在应力范围(Si,Si+△Si)内的应力循环次数表示为ni＝vTp(Si)△Si，式中的v表示单位时间内的应力循环次数，由每秒的峰点数E[P]决定，即v＝E[P]，p(Si)表示应力级水平为Si时的幅值概率密度函数值，△Si为微小应力级变化范围；第七步：在MSC.Fatigue中选择对应金属材料的S‑N曲线，其曲线的特征表达式为N(Si)＝K/Sm，应用Miner线性累计损伤理论D＝∑Di＝∑ni/Ni，当总体损伤度D＝1时得到构件发声疲劳破坏的时间疲劳寿命为：T=KE[P]∫SmP(S)dS]]>其中，K与m为材料常数，N为金属材料循环次数，D为损伤度，最终计算得到构件疲劳破坏时的时间寿命分布云图与关键结点最小寿命时间范围。...

【技术特征摘要】
1.一种含不确定性金属结构声振疲劳寿命的估算方法，其特征在于实现步骤如下：第一步：依据飞行任务剖面时间数据和状态声载荷测量带宽声压级，采用将其转换为频谱载荷，其中L＝Lb-10lg(Δf)，P0＝2×10-5Pa为参考声压，Lb为带宽声压级，Δf为频带宽度，G(f)为转化后的噪声功率谱密度值；第二步：引入区间向量x∈xI＝[E,a,b]合理定量贫信息、少数据条件下结构的不确定参数，其中E为金属材料弹性模量，a与b分别表示结构的几何参数，含结构参数的不确定性可以表示为：代表参量的取值上界，x＝[E,a,b]＝[Ec-Er,ac-ar,bc-br]，x代表参量的取值下界，其中上标c代表中心值，上标r代表半径；第三步：建立所关心结构的几何模型，分析构件的连接形式，合理划分网格施加边界条件形成有限元分析模型，在有限元分析软件中加载单位均匀面压力载荷，其中载荷的频率变化范围与第一步所得的噪声载荷的频率范围相同，对模型进行频率响应分析后提取结构各个结点的应力传递函数；第四步：应用平方和开方方法，在MSC.Patran随机振动模块中读取第三步得到的应力传递函数，并施加第一步所得的频谱噪声载荷得到结点的应力功率谱密度曲线，其中随机振动分析输出文件还包括频率响应的自相关函数、每单位时间的正斜率方向的零交点个数及应力响应的均方值；第五步：随机噪声的典型特点是压力时间历程幅值是随机变化的，也就是无规律非衰减不能用解析函数表达的形式，这类噪声频率的变化范围宽，其谱值直到很高频率都是连续的，结合模拟此类典型宽带随机过程的Dirlik模型与区间传播分析方法，以上一步得到的应力PSD曲线为输入，引入非概率区间过程理论中泰勒展开分析法得到不确定变量扩张后各个关键结点的应力雨流幅值概率密度函数变化范围；第六步：对于连续分布应力状态，将时间T内在应力范...

【专利技术属性】
技术研发人员：邱志平，苏欢，王磊，王晓军，田靖军，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11