CN112326188B - 一种基于doe和响应面模型的六分力应变天平公式的系数标定方法及计算方法 - Google Patents
一种基于doe和响应面模型的六分力应变天平公式的系数标定方法及计算方法 Download PDFInfo
- Publication number
- CN112326188B CN112326188B CN202011243088.1A CN202011243088A CN112326188B CN 112326188 B CN112326188 B CN 112326188B CN 202011243088 A CN202011243088 A CN 202011243088A CN 112326188 B CN112326188 B CN 112326188B
- Authority
- CN
- China
- Prior art keywords
- surface model
- coefficient
- response surface
- matrix
- doe
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
- G01M9/062—Wind tunnel balances; Holding devices combined with measuring arrangements
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明公开一种基于DOE和响应面模型的六分力应变天平公式的系数标定,包括如下步骤:步骤一、采用实验设计方法生成样本点,根据样本点在六个分力方向进行加载,并同时采集六个分力方向上对应于加载量的输出电压信号值增量;步骤二、采用近似模型标定一次干扰系数和二次干扰系数:步骤三、通过天平校准公式导出系数矩阵。本发明还公开了一种基于DOE和响应面模型的六分力应变天平公式的计算方法,采用所述的基于DOE和响应面模型的六分力应变天平公式的系数标定方法获得系数矩阵,并得到六分力应变天平的六分力:F=A‑1Δn;式中,A‑1为系数矩阵的逆矩阵,Δn为输出电压信号增量矩阵,F为六分力矩阵。本发明求解速度快、精度高。
Description
技术领域
本发明涉及汽车风洞试验技术领域,更具体的是,本发明涉及一种基于 DOE和响应面模型的六分力应变天平公式的系数标定方法及计算方法。
背景技术
随着国内汽车行业的不断发展,汽车空气动力学性能备受重视,而风洞试验作为汽车空气动力学开发过程中的重要一环,其不仅为理论分析和数值计算的基础,并且还可以用来检验理论结果的正确性和可靠性。所以世界各大车企以及高校都纷纷投资建设汽车风洞,用以更好的开展汽车空气动力学性能研究。而国内目前已建成汽车风洞相对较少,针对全尺寸汽车风洞六分力应变天平公式的标定及求解经验不足。
《风洞天平》一书中所述,对于一般六分力天平,不考虑一次非对称干扰项与三次立方干扰项的天平校准通式的隐式为
与该天平校准通式对应的天平工作公式为
书中公式不够直观,应用起来比较复杂。且从公式可知六分力天平公式的系数矩阵最少为6×6阶矩阵,天平公式系数标定需求解36个未知系数且方程组只有六个方程,系数标定十分复杂。若要提高精度考虑二次项干扰,系数矩阵即为6×27阶矩阵,含有162个未知数,复杂程度更高。
发明内容
本发明设计开发了一种基于DOE和响应面模型的六分力应变天平公式的系数标定方法,可采用一阶响应面模型只标定一次干扰系数,也可以采用二阶响应面模型同时标定一次干扰系数和二次干扰系数。
本发明设计开发了一种基于DOE和响应面模型的六分力应变天平公式的计算方法,求解速度快,使用简便,求解精度高。
本发明提供的技术方案为:
一种基于DOE和响应面模型的六分力应变天平公式的系数标定方法,包括如下步骤:
步骤一、采用实验设计方法生成样本点,根据样本点在六个分力方向进行加载,并同时采集六个分力方向上对应于加载量的输出电压信号值增量;
步骤二、采用近似模型标定一次干扰系数和二次干扰系数:
步骤三、通过天平校准公式导出系数矩阵。
优选的是,所述输出电压信号值增量矩阵的响应实际值与近似值之差满足:
ε0=ε+β0;
式中,ε为输出电压信号值增量矩阵的响应实际值与近似值之间的随机误差,β0为拟合响应面模型的常数项。
优选的是,所述近似模型选择响应面模型,满足:
优选的是,一阶响应面模型中的响应近似值满足:
式中,β0,β1…βM为多项式系数,x1,…xM为输入变量。
优选的是,所述一阶响应面模型所需最少样本点数满足:
S1=M+1;
式中,S1为一阶响应面模型所需最少样本点数,M为输入变量个数。
优选的是,二阶响应面模型中的二阶响应近似值满足:
式中,β0,β1…βM,βM+1,…β2M,∑i≠jβijxixj为多项式系数。
优选的是,所述二阶响应面模型所需最少样本点数满足:
S2=(M+1)(M+2)/2;
式中,S2为二阶响应面模型所需最少样本点数。
一种基于DOE和响应面模型的六分力应变天平公式的计算方法,采用所述的基于DOE和响应面模型的六分力应变天平公式的系数标定方法获得系数矩阵,并得到六分力应变天平的六分力:
F=A-1Δn;
式中,A-1为系数矩阵的逆矩阵,Δn为输出电压信号增量矩阵,F为六分力矩阵。
本发明所述的有益效果:
(1)本发明提供的基于DOE和响应面模型的六分力应变天平公式的系数标定方法,可采用一阶响应面模型只标定一次干扰系数,也可以采用二阶响应面模型同时标定一次干扰系数和二次干扰系数。
(2)本发明提供的基于DOE和响应面模型的六分力应变天平公式的系数标定方法,采用DOE方法生成样本点,样本点数量少,标定简便且精度高。
(3)本发明提供的基于DOE和响应面模型的六分力应变天平公式的计算方法,求解速度快,使用简便,求解精度高。
附图说明
图1为本发明所述六分力天平测量系统示意图。
图2为本发明所述基于DOE和响应面模型的六分力应变天平公式的系数标定方法流程图。
具体实施方式
下面结合附图对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
本发明提供了一种基于DOE和响应面模型的六分力应变天平公式的系数标定方法,如图2所示,具体包括:
拟定天平校准公式如下:
Δn=AF;
式中,Δn为输出电压信号增量矩阵,A为系数矩阵,F为六分力矩阵。
其中,具体的:
1)只考虑一次干扰系数:
A=(aij)6×6;
F=[X Y Z MX MY MZ];
2)同时考虑一次干扰系数和二次干扰系数:
A=(aij)6×27;
可知,天平校准公式为六个方程构成的方程组,式中:
天平公式标定即通过标准加载六分力矩阵并采集输出电压信号增量矩阵来求取系数矩阵。
其中,近似模型方法是通过数学模型的方法逼近一组输入变量与输出变量的方法,通过建立经验公式,获得输入输出之间的量化关系,常用的近似模型主要包括:响应面模型、神经网络模型、正交多项式模型、克里格模型。
在本实施例中,采用响应面模型的方法,所述响应面模型利用多项式函数拟合设计空间,其优势在于:通过较少的试验在局部范围内比较精确的逼近函数关系,并用简单的代数表达式展现出来,计算简单;可以拟合复杂的响应关系,鲁棒性良好;数学理论基础扎实,系统性、实用性强,适用范围广。
近似模型将输入变量和输出响应之间的关系描述为:
对于响应面模型,响应近似值可以是一阶、二阶、三阶和四阶,一阶响应面模型中的响应近似值满足:
式中,β0,β1…βM为多项式系数,x1,…xM为输入变量;
二阶响应面模型中的二阶响应近似值满足:
式中,β0,β1…βM,βM+1,…β2M,∑i≠jβijxixj为多项式系数;
其中,试验设计方法(DOE)是以概率论与数理统计为理论基础、合理安排试验的一种方法,包括三个步骤:试验计划、执行试验和结果分析。DOE 方法中常用的采样方法有:参数试验、全因子设计、正交数组、拉丁超立方设计等,可根据实际情况具体选择,本实施例中主要应用DOE方法中提供的参数试验采样方法生成样本点用以建立近似模型,样本点数量少,标定简便,建立一阶响应面模型所需最少样本点数满足:
S1=M+1;
式中,S1为一阶响应面模型所需最少样本点数,M为输入变量个数。
所述二阶响应面模型所需最少样本点数满足:
S2=(M+1)(M+2)/2;
式中,S2为二阶响应面模型所需最少样本点数。
对比一阶响应面模型中的响应近似值、二阶响应面模型中的二阶响应近似值与天平校准公式,以输出电压信号增量矩阵为响应,以六分力矩阵为输入变量,即可采用一阶响应面模型拟合一次干扰系数,也可采用二阶响应面模型同时拟合一次干扰系数和二次干扰系数。
输入变量和输出响应的关系可以表示为:
ε0=ε+β0;
式中,ε为输出电压信号值增量矩阵的响应实际值与近似值之间的随机误差,β0为拟合响应面模型的常数项,天平标定时β0为极小量,可忽略不计,通过校验ε0大小来验证近似模型可信度,最后导出系数矩阵即可。
本发明提供的基于DOE和响应面模型的六分力应变天平公式的系数标定方法,可采用一阶响应面模型只标定一次干扰系数,也可以采用二阶响应面模型同时标定一次干扰系数和二次干扰系数。
本发明还提供一种基于DOE和响应面模型的六分力应变天平公式的计算方法,采用上述的基于DOE和响应面模型的六分力应变天平公式的系数标定方法获得系数矩阵,运用标定好的天平公式在风洞测力实验中根据采集到的六分量的输出电压信号值增量求解六分力,本发明提供采用逆矩阵法求解:
对天平校准公式两侧分别左乘系数矩阵A的逆矩阵,即可得出:
F=A-1Δn;
式中,A-1为系数矩阵的逆矩阵,Δn为输出电压信号增量矩阵,F为六分力矩阵。
代入六分量的输出电压信号增量即可求出六分力,求解便捷。
本发明提供的基于DOE和响应面模型的六分力应变天平公式的计算方法,求解速度快,使用简便,求解精度高。
实施例
本实施例采用的六分力天平测量系统如图1所示,以一阶响应面模型标定的一次干扰系数进行说明:
考虑天平六分力之间的交互作用不显著,故此采用DOE方法中提供的参数试验法生成样本点,若想进一步提高拟合精度,可选用其他能够提供交互作用信息的样本生成方法,如正交数组法、拉丁超立方法等,由于只考虑一次干扰系数,近似模型选用一阶响应面模型。
天平校准中心取为天平校准加载系统的加载平面中心,使用天平校准标准砝码通过天平校准加载系统对天平加载,每个分力方向从0到70%设计量程选取合适间距进行加载并重复三次,采集各分量输出电压信号,并对电压信号进行去坏值和均值处理,最后得到六个分力方向上对应于每个加载量的输出电压信号值增量矩阵。
通过一阶响应面模型拟合得到的系数矩阵为:
另外进行三组组合加载以校验近似模型可信度,数据如表一和表二所述:
表一 三组组合标准加载值
标准加载 | X/N | Y/N | Z/N | M<sub>X</sub>/Nm | M<sub>Y</sub>/Nm | M<sub>Z</sub>/Nm |
第一组 | 1000 | 0 | 400 | 0 | 239.12 | 0 |
第二组 | 500 | 0 | 200 | 147 | 119.56 | 166.6 |
第三组 | 1000 | 200 | 400 | -294 | 0 | 0 |
表二 三组组合求解载荷值
求解载荷 | X | Y | Z | M<sub>X</sub> | M<sub>Y</sub> | M<sub>Z</sub> |
第一组 | 999.2432 | -2.6991 | 399.8281 | 2.5391 | 242.7378 | -2.4466 |
第二组 | 499.1344 | 0.5868 | 200.5781 | 152.8381 | 125.1637 | 172.3349 |
第三组 | 998.6436 | 202.2807 | 400.8745 | -287.1814 | 0.7044 | 0.4014 |
由系数矩阵求解得到的六分力结果可知,其与各方向实际加载载荷基本一致,近似模型具有一定的可信度,本发明提供的系数标定方法具有合理性。同时也可采用其他能够提供交互作用信息的样本生成方法以及同时拟合一次干扰系数和二次干扰系数的二阶响应面模型进一步提高近似模型的拟合精度。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
Claims (6)
1.一种基于DOE和响应面模型的六分力应变天平公式的系数标定方法,其特征在于,包括如下步骤:
步骤一、采用实验设计方法生成样本点,根据样本点在六个分力方向进行加载,并同时采集六个分力方向上对应于加载量的输出电压信号值增量;
步骤二、采用近似模型标定一次干扰系数和二次干扰系数:
其中,所述输出电压信号值增量矩阵的响应实际值与近似值之差满足:
ε0=ε+β0;
式中,ε为输出电压信号值增量矩阵的响应实际值与近似值之间的随机误差,β0为拟合响应面模型的常数项;
所述近似模型选择响应面模型,满足:
步骤三、通过天平校准公式导出系数矩阵。
3.如权利要求2所述的基于DOE和响应面模型的六分力应变天平公式的系数标定方法,其特征在于,所述响应面模型中的样本点通过试验设计方法获得,所述一阶响应面模型所需最少样本点数满足:
S1=M+1;
式中,S1为一阶响应面模型所需最少样本点数,M为输入变量个数。
5.如权利要求4所述的基于DOE和响应面模型的六分力应变天平公式的系数标定方法,其特征在于,所述二阶响应面模型所需最少样本点数满足:
S2=(M+1)(M+2)/2;
式中,S2为二阶响应面模型所需最少样本点数。
6.一种基于DOE和响应面模型的六分力应变天平公式的计算方法,其特征在于,采用如权利要求1-5中任一项所述的基于DOE和响应面模型的六分力应变天平公式的系数标定方法获得系数矩阵,并得到六分力应变天平的六分力:
F=A-1Δn;
式中,A-1为系数矩阵的逆矩阵,Δn为输出电压信号增量矩阵,F为六分力矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011243088.1A CN112326188B (zh) | 2020-11-10 | 2020-11-10 | 一种基于doe和响应面模型的六分力应变天平公式的系数标定方法及计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011243088.1A CN112326188B (zh) | 2020-11-10 | 2020-11-10 | 一种基于doe和响应面模型的六分力应变天平公式的系数标定方法及计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112326188A CN112326188A (zh) | 2021-02-05 |
CN112326188B true CN112326188B (zh) | 2022-02-25 |
Family
ID=74315516
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011243088.1A Active CN112326188B (zh) | 2020-11-10 | 2020-11-10 | 一种基于doe和响应面模型的六分力应变天平公式的系数标定方法及计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112326188B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112800633B (zh) * | 2021-04-06 | 2021-07-13 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种多元校准风洞天平数据的处理方法 |
CN117109666B (zh) * | 2023-10-24 | 2024-01-05 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种风洞天平应变计粘贴工艺质量的检验装置和检验方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462395A (zh) * | 2017-08-16 | 2017-12-12 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种天平分量间干扰系数校准方法 |
CN109029908A (zh) * | 2018-09-03 | 2018-12-18 | 吉林大学 | 一种全尺寸汽车风洞六分力应变天平公式的系数标定方法及计算方法 |
CN111046505A (zh) * | 2019-11-28 | 2020-04-21 | 河海大学 | 基于响应面模型的轴流泵辐条参数优化设计方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB612749A (en) * | 1946-05-31 | 1948-11-17 | Siemens Brothers & Co Ltd | Improvements relating to torsionmeters |
US2708368A (en) * | 1950-12-21 | 1955-05-17 | Continental Silver Co Inc | Equipment for rating by volume, weight, and zone |
US3084253A (en) * | 1960-09-06 | 1963-04-02 | Barnes Eng Co | Continuously self-calibrating differential detection system |
US20180289577A1 (en) * | 2015-01-18 | 2018-10-11 | Calabrian High Tech Srl | Six Degrees of Freedom Self Balanced Hybrid Serial-Parallel Robotic Systemfor Rehabilitation of Upper and Lower Limbs, Left and Right |
CN112800633B (zh) * | 2021-04-06 | 2021-07-13 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种多元校准风洞天平数据的处理方法 |
-
2020
- 2020-11-10 CN CN202011243088.1A patent/CN112326188B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462395A (zh) * | 2017-08-16 | 2017-12-12 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种天平分量间干扰系数校准方法 |
CN109029908A (zh) * | 2018-09-03 | 2018-12-18 | 吉林大学 | 一种全尺寸汽车风洞六分力应变天平公式的系数标定方法及计算方法 |
CN111046505A (zh) * | 2019-11-28 | 2020-04-21 | 河海大学 | 基于响应面模型的轴流泵辐条参数优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112326188A (zh) | 2021-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112326188B (zh) | 一种基于doe和响应面模型的六分力应变天平公式的系数标定方法及计算方法 | |
CN101398311B (zh) | 基于灰色系统理论的重复动态测量数据处理方法 | |
CN112800633B (zh) | 一种多元校准风洞天平数据的处理方法 | |
CN109933925B (zh) | 一种金属板材的冲压成形性能预测方法 | |
CN113537319B (zh) | 一种针对类型化模型的风洞天平终身学习校准方法及系统 | |
CN112414668B (zh) | 一种风洞试验数据静气弹修正方法、装置、设备及介质 | |
CN116090111B (zh) | 一种基于深度学习模型的汽车钢板弹簧疲劳寿命预测方法 | |
CN115809504A (zh) | 一种多源气动数据不确定度量化评估方法和系统 | |
Ma et al. | Parameter reduction of composite load model using active subspace method | |
CN114583767B (zh) | 一种数据驱动的风电场调频响应特性建模方法及系统 | |
CN116412995A (zh) | 一种分区间确定风洞应变天平测量不确定度的方法 | |
Wang et al. | An efficient hybrid sequential approximate optimization method for problems with computationally expensive objective and constraints | |
Wang et al. | Modified Tikhonov regularization in model updating for damage identification | |
CN112489734A (zh) | 一种内燃机替代燃料二甲醚燃烧反应机理模型的简化方法 | |
Sindia et al. | Parametric fault testing of non-linear analog circuits based on polynomial and V-transform coefficients | |
CN114993604B (zh) | 一种基于深度学习的风洞天平静态校测方法 | |
CN114254533B (zh) | 考核疲劳振动对产品组部件固定角度影响和预测的方法 | |
CN107587955A (zh) | 基于深信度网络的火箭发动机推力偏移量的标定方法 | |
CN115186385A (zh) | 车身刚度学科的后处理方法、装置、设备及存储介质 | |
CN107203677B (zh) | 一种电子系统多退化进程研究方法 | |
Lee | Stochastic optimization for design under uncertainty with dependent random variables | |
Ulbrich | Iterative strain-gage balance calibration data analysis for extended independent variable sets | |
CN110765560A (zh) | 一种基于时变阻尼机械机构振动预测方法 | |
Akbar et al. | Calibration of Piezoelectric Dynamometer based on Neural Networks | |
Zhang et al. | Sensitivity Analysis of Weather Factors Affecting PV Module Output Power Based on Artificial Neural Network and Sobol Algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |