一种高频底座力天平的动力校准方法与流程

文档序号:12064305阅读:来源:国知局

技术特征:

1.一种高频底座力天平的动力校准方法,其特征在于,包括以下步骤:

S1:对测量信号x(t)进行白化得到白化后信号z(t);

S2:寻求正交矩阵V,使得z(t)=Vq(t),进而得到分离信号q(t);

S3:在模态坐标下,对分离信号进行固有频率和模态阻尼比识别;

S4:根据识别得到的参数,对分离信号进行修正;

S5:由修正后的分离信号逆推,得到修正后的气动荷载。

2.根据权利要求1所述的高频底座力天平的动力校准方法,其特征在于,步骤S1中对测量信号x(t)进行白化的具体步骤是:

S1-1:计算x(t)的协方差矩阵Cx:

Cx=E[x(t)xT(t)]

S1-2:对Cx进行特征值分解,得到相应的特征向量矩阵E和特征值对角矩阵D;

S1-3:计算白化矩阵W:

W=D-1/2ET

假定Cx是正定的,其特征值均为正,从而保证白化矩阵W必然存在;

S1-4:计算白化后信号z(t):

z(t)=Wx(t)。

3.根据权利要求2所述的高频底座力天平的动力校准方法,其特征在于,步骤S2的具体步骤是:

S2-1:计算延时相关矩阵Rz(τ):

Rz(τ)=E[z(t)zT(t+τ)];

S2-2:根据Rz(τ)计算Given旋转角度,然后由Given旋转角度组成旋转矩阵;

S2-3:判断旋转矩阵的非对角元素平方是否小于一给定正数ε,结论为是,执行步骤S2-5;结论为否,执行步骤S2-4;

S2-4:利用得到的旋转矩阵更新Rz(τ)、V,V给定初值为单位矩阵,执行步骤S2-2;

S2-5:得到正交矩阵V,进而得到混合矩阵Φ=W-1V和分离矩阵B=Φ-1=VTW。

4.根据权利要求1所述的高频底座力天平的动力校准方法,其特征在于,步骤S3的具体步骤是:

S3-1:计算第i阶分离信号qi(t)的功率谱密度Sqi(f);

S3-2:计算无因次化的频响函数hi(f):

<mrow> <msub> <mi>h</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mi>f</mi> <msub> <mi>f</mi> <mi>i</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>2</mn> <msub> <mi>&zeta;</mi> <mi>i</mi> </msub> <mfrac> <mi>f</mi> <msub> <mi>f</mi> <mi>i</mi> </msub> </mfrac> <mi>i</mi> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow>

式中,fi和ζi分别为第i阶模态的固有频率和模态阻尼比,f为频率序列。

S3-3:计算第i阶模态气动力qr,i(t)的功率谱密度Sqr,i(f):

<mrow> <msub> <mi>S</mi> <msub> <mi>q</mi> <mrow> <mi>r</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>S</mi> <msub> <mi>q</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> </mrow> <msup> <mrow> <mo>|</mo> <msub> <mi>h</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>;</mo> </mrow>

S3-4:拟合曲线使得原始曲线与拟合曲线的残差最小;ki和αi为常数;Δf为谱线频率间隔,(n1Δf,n2Δf)为参加拟合的频率范围;

S3-5:得到四个参数值fi、ζi、ki和αi

5.根据权利要求3所述的高频底座力天平的动力校准方法,其特征在于,步骤S4中,根据得到的分离矩阵B及由各个解耦源信号识别得到的模态参数,可以得到模态坐标下修正后的模型基底气动力的功率谱密度矩阵为:

<mrow> <msub> <mi>S</mi> <msub> <mi>q</mi> <mi>r</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>h</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <msub> <mi>BS</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <msup> <mi>B</mi> <mi>T</mi> </msup> <msup> <mrow> <mo>(</mo> <msup> <mi>h</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>(</mo> <mi>f</mi> <mo>)</mo> <mo>)</mo> </mrow> <mi>H</mi> </msup> </mrow>

式中,h(f)为各主坐标的频率响应函数hi(f)构成的复对角矩阵,上标H表示矩阵的Hermite转置,Sx(f)表示最原始的测量信号的功率谱密度矩阵。

6.根据权利要求5所述的高频底座力天平的动力校准方法,其特征在于,步骤S5中,根据模态坐标与物理坐标之间的转换关系,得到:

<mrow> <msub> <mi>S</mi> <msub> <mi>x</mi> <mi>r</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&Phi;S</mi> <msub> <mi>q</mi> <mi>r</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <msup> <mi>&Phi;</mi> <mi>T</mi> </msup> </mrow>

是修正后的气动力的功率谱密度矩阵。

7.根据权利要求1所述的高频底座力天平的动力校准方法,其特征在于,针对非比例阻尼,在步骤S1之前要进行信号前处理:将测量信号x(t)相移90°得到x90(t),然后将二者组装成复信号;后续步骤都是对复信号进行处理。

8.根据权利要求7所述的高频底座力天平的动力校准方法,其特征在于,信号前处理的具体操作步骤如下:

将测得的长度为K的基底各分量信号xi(t)逆序排列,得到xir(t);建立临时序列di(t):

di(t)=-xir(t),1≤t≤K

di(t)=xi(t),K+1≤t≤2K;

di(t)=-xir(t),2K+1≤t≤3K

对di(t)进行Hilbert变换,得到di90(t),从di90(t)中分离出xi90(t),对每个基底分量的时间序列应用以上步骤,组装成矩阵,得到组装信号x(t)+ix90(t)。

9.根据权利要求6所述的高频底座力天平的动力校准方法,其特征在于,针对非比例阻尼,根据模态坐标与物理坐标之间的转换关系,得到:

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1