功率谱密度曲线的解析表达式

功率谱密度曲线的解析表达式 摘要 本文主要论述了在用有限单元法对结构进行随机振动分析时传统功率谱密度值计算 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 的缺点，用分段函数给出了功率谱密度曲线的 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 解析表达式和近似解析表达式，使得功率谱密度值的计算简单易行。 关键词 功率谱 随机振动 解析表达式 The Analytic Expression of the Power Spectrum Density(PSD)Curve. Chen Gengchao (The North West China Research Institute of Electronic Equipment. Xi’an,710065,China) Abstract This paper mainly discusses the defects of the traditional PSD value calculating method when analyzing the structure's random vibration by FEA method. Put forward the standard analytic expression of the PSD curve and the approximate analytic expression by segmental function. This makes the PSD value calculation simply and easily. Keywords power spectrum random vibration analytic expression 1 引言 在对结构做随机振动试验或用有限元软件对结构做随机振动分析时，施加在结构上的随机激励载荷条件(即功率谱密度曲线)通常是用不同频率范围内单位倍频程所对应的分贝数来表达的。典型的随机激励载荷的表达形式如表1所示(g为重力加速度)，其在对数坐标系下功率谱密度值随频率变化的曲线如图1所示。 表1:典型随机激励载荷条件 谱值 10,200 Hz +6dB/oct 随 2/Hz 200,1500 Hz 0.12g机 激 1500,2000 Hz -12dB/oct 励 总均方根加速度 14.4g 频率(Hz) 用有限元软件对结构做随机振动分析 时，需要以表格的形式按离散点给出不同 图1:典型随机激励载荷的功率谱密度曲线 频率所对应的功率谱密度值来拟合出相应 的功谱密度曲线。用传统的方法来计算不同频率所对应的功率谱密度值，比较繁琐，效率很低。为此，我们用分段函数给出了功率谱密度曲线的标准解析表达式和方便计算的近似解析表达式，使得功率谱密度值的计算简单易行。 2 传统功率谱密度值的计算方法 假设频率对应的功率谱密度值为，频率对应的功率谱密度值为d，用来表示对的fffdfmnnn111 倍频程数，则有关系式: fnf1 式(1) m,log2 用表示单位倍频程使功率谱密度值增加的分贝数，相当于对数功率谱密度曲线中直线段的斜率，x 对应表1中的“+6”或“-12”，则有关系式: dnd1 式(2) xm,10log10 如果知道某频率对应的功率谱密度值，则可由式(1)和式(2)可得到频率对应的功率dffnn1谱密度值。如果知道某频率对应的功率谱密度值，则可由式(1)和式(2)可得到频率对fdfdn111 应的功率谱密度值。 dn 以表1所给的载荷条件为例。 2已知200Hz时的功率谱密度值为0.12g，求10Hz时的功率谱密度值。由式(1)可得200Hz对d120.12g200d101m,log,4.3226，4.322,10log10Hz的倍频程数为:，由式(2)则有，可得到10Hz210时的功率谱密度值为0.0294。 d12已知1500Hz时的功率谱密度值为0.12g，求2000Hz时的功率谱密度值。由式(1)可得2000Hzdn d2000n15001500m,log,0.415对1500Hz的倍频程数为:，由式(2)则有，可0.415，(,12),10log210得到2000Hz时的功率谱密度值为3.661。 dn 表2:通用随机激励载荷条件 用此方法可求出各所需频率点的功率谱密度值。 3 功率谱密度曲线的解析表达式, Hz dB/oct faf随 12 2 , Hz g/Hz bf机 f32 假设通用的随机激励载荷条件以表2的形式给出。激 ,Hz dB/oct ffc34 f通过分析频率和功率谱密度值d的关系，对表2给出的励 总均方根加速度 g s载荷条件，我们用分段函数给出的功率谱密度曲线的表达式如 下: flog10a2 Hz 式(3) d,A(10)f,f,f122 为常值 Hz 式(4) f,f,fd,bg23 flog10c2 Hz 式(5) f,f,fd,B(10)34 2把频率、和它们所对应的功率谱密度值分别代人式(3)和式(5)，可确定常数A和B。fbgf32 2bg 式(6) A,f2log10a2(10) 2bg 式(7) B,f3log10c2(10) 对表1所给载荷条件，=6、b=0.12、=-12、=200、=1500，由式(6)和式(7)可得到ffac32,413A,2.9876，10B,5.2787，10,。则用分段函数给出的功率谱密度曲线的表达式如下: flog4610,210,f,200 Hz 式(8)d,2.9876，10，(10) 2200,f,1500 为常值 Hz 式(9)d,0.12g flog131210,21500,f,2000 Hz 式(10)d,5.2787，10，(10) )可得到10 Hz和2000Hz时的功率谱密度值分别为0.0294和3.661，与传统由式(8)和式(10 计算方法得到的值相同。这也验证了解析表达式的正确性。 4 功率谱密度曲线的近似解析表达式 在式(3)中，假设常数为3的整数倍，即(为正或负的整数)。a,3mam 3m100.3mm，所以式(3)可近似表达为: 由于10,(10),2 ffloglogmmma322。 d,A(2),A(2),Af,Af 同理，在式(5)中，如果常数为3的整数倍，即(为正或负的整数)。c,3ncn 则式(5)可近似表达为: ffloglognnnc322。 d,A(2),A(2),Af,Af 故对表2给出的载荷条件，当常数和为3的整数倍时，我们用分段函数给出的功率谱密度曲ac 线的近似表达式如下: a3 Hz，常数为3的整数倍 式(11)d,Aff,f,fa12 2 为常值 Hz 式(12)f,f,fd,bg23 c3 Hz，常数为3的整数倍 式(13)f,f,fd,Afc34 2把频率、和它们所对应的功率谱密度值分别代人式(11)和(13)，可确定常数A和B。fbgf32 2a3 式(14) A,bgf2 2c3B,bgf 式(15) 3 对表1所给载荷条件，=6、=0.12、=-12、=200、=1500，由式(14)和式(15)可得bffac32,413A,2.8812，10B,5.8344，10到,。则用分段函数给出的功率谱密度曲线的表达式如下: ,4210,f,200 Hz 式(16) d,2.8812，10，f 2200,f,1500 为常值， Hz 式(17) d,0.12g 13,41500,f,2000 Hz 式(18) d,5.8344，10，f 由式(16)可得到10 Hz时的功率谱密度值为0.0288，与传统计算方法得到的值0.0294相比，误 差为2.04%。由式(18)可得到2000Hz时的功率谱密度值为3.647，与传统计算方法得到的值3.661 相比，误差为0.382%。所以，近似表达式的计算精度是足够的。 由于功率谱密度曲线的近似表达式是幂函数，可很方便地通过积分得到曲线下的面积，从而得到随机激励的总均方根值。对表1所给的随机激励载荷条件，按近似表达式(16)、(17)、(18)进行积分，得到总均方根加速度为14.1g，与14.4g的准确值的误差为2.08%，可满足计算精度要求。 7 总结 在用有限单元法对结构进行随机振动分析时，用传统的方法来计算不同频率所对应的功率谱密度值，比较繁琐，效率很低。本文结合实例，用分段函数给出了功率谱密度曲线的标准解析表达式，并进一步给出了特殊情况下的近似解析表达式，使得功率谱密度值的计算简单易行。