删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Kriging序贯设计方法在滑动轴承优化中的应用

本站小编 哈尔滨工业大学/2019-10-24

Kriging序贯设计方法在滑动轴承优化中的应用

张泽斌,张鹏飞,郭红,李永

(郑州大学 机械工程学院, 郑州 450001)



摘要:

利用 Kriging代理模型提供目标函数无偏预测值和理论置信区间的优势,比较传统试验设计方法和基于Kriging模型的序贯加点方法对模型的影响, 结合设计空间的全局搜索和最优解临近区间的局部搜索,引入并行加点准则及相应的收敛条件,得到精度、效率高的代理模型,用2个经典优化测试函数进行验证和评估. 结果表明,与传统试验设计方法相比,基于Kriging的序贯加点方法得到的模型全局精度更高且能更快地收敛到优化问题的真实最优解. 最后,以动静压滑动轴承为优化设计对象,单位承载力下摩擦功耗为目标函数,考虑几何结构及工况等约束条件,采用传统试验设计方案和Kriging序贯加点方案分别建立目标函数的模型,分别进行优化设计,并同传统的复合形优化结果进行对比. 3种方案对比结果显示,Kriging加点方案在有限迭代步数下,对于降低单位承载力下的摩擦功耗效果最为显著,验证了该方法快速收敛的特性.

关键词:  Kriging  代理模型  序贯设计  动静压轴承  优化设计

DOI:10.11918/j.issn.0367-6234.201810147

分类号:TH122TH117.2

文献标识码:A

基金项目:国家自然科学基金资助项目(51575498)



Implementation of Kriging model based sequential design on the optimization of sliding bearing

ZHANG Zebin,ZHANG Pengfei,GUO Hong,LI Yong

(School of Mechanical Engineering,Zhengzhou University,Zhengzhou 450001,China)

Abstract:

A built-in prime benefit of Kriging surrogate model resides in its unbiased prediction and the associated confidence intervals, comparisons have been done between the classical DoE methods and Kriging model based sequential design. This latter considers both the design space global exploration and optimum-neighborhood exploitation. A parallel point-adding training strategy and its corresponding convergence criterion were introduced to improve the model precision. The approach had been illustrated with two classical optimization test functions. Results show not only a more accurate model but also a possibility to reduce the number of sampling points. Finally, the training strategy was experimented for the optimal design of a sliding bearing. Friction power loss per unit load capacity was taken as the objective function to be modeled. Two surrogate model based optimizations had been per-formed based on classical DoE such as Orthogonal Latin Hypercube and Kriging sequential strategy respectively. These two optimization results are compared with a previously optimization using Complex-optimization method. Within a limited number of iterations, the Kriging model based training strategy showed the best convergence to the global optimum among those 3 methods.

Key words:  Kriging  Surrogate model  Sequential Design  Hybrid Bearing  Optimization


张泽斌, 张鹏飞, 郭红, 李永. Kriging序贯设计方法在滑动轴承优化中的应用[J]. 哈尔滨工业大学学报, 2019, 51(7): 178-183. DOI: 10.11918/j.issn.0367-6234.201810147.
ZHANG Zebin, ZHANG Pengfei, GUO Hong, LI Yong. Implementation of Kriging model based sequential design on the optimization of sliding bearing[J]. Journal of Harbin Institute of Technology, 2019, 51(7): 178-183. DOI: 10.11918/j.issn.0367-6234.201810147.
基金项目 国家自然科学基金资助项目(51575498) 作者简介 张泽斌(1982—),男,讲师 通信作者 张泽斌,zebin.zhang@zzu.edu.cn 文章历史 收稿日期: 2018-10-25



Contents            -->Abstract            Full text            Figures/Tables            PDF


Kriging序贯设计方法在滑动轴承优化中的应用
张泽斌, 张鹏飞, 郭红, 李永    
郑州大学 机械工程学院, 郑州 450001

收稿日期: 2018-10-25
基金项目: 国家自然科学基金资助项目(51575498)
作者简介: 张泽斌(1982—),男,讲师
通信作者: 张泽斌,zebin.zhang@zzu.edu.cn


摘要: 利用Kriging代理模型提供目标函数无偏预测值和理论置信区间的优势,比较传统试验设计方法和基于Kriging模型的序贯加点方法对模型的影响, 结合设计空间的全局搜索和最优解临近区间的局部搜索,引入并行加点准则及相应的收敛条件,得到精度、效率高的代理模型,用2个经典优化测试函数进行验证和评估.结果表明,与传统试验设计方法相比,基于Kriging的序贯加点方法得到的模型全局精度更高且能更快地收敛到优化问题的真实最优解.最后,以动静压滑动轴承为优化设计对象,单位承载力下摩擦功耗为目标函数,考虑几何结构及工况等约束条件,采用传统试验设计方案和Kriging序贯加点方案分别建立目标函数的模型,分别进行优化设计,并同传统的复合形优化结果进行对比. 3种方案对比结果显示,Kriging加点方案在有限迭代步数下,对于降低单位承载力下的摩擦功耗效果最为显著,验证了该方法快速收敛的特性.
关键词: Kriging    代理模型    序贯设计    动静压轴承    优化设计    
Implementation of Kriging model based sequential design on the optimization of sliding bearing
ZHANG Zebin, ZHANG Pengfei, GUO Hong, LI Yong    
School of Mechanical Engineering, Zhengzhou University, Zhengzhou 450001, China


Abstract: A built-in prime benefit of Kriging surrogate model resides in its unbiased prediction and the associated confidence intervals, comparisons have been done between the classical DoE methods and Kriging model based sequential design. This latter considers both the design space global exploration and optimum-neighborhood exploitation. A parallel point-adding training strategy and its corresponding convergence criterion were introduced to improve the model precision. The approach had been illustrated with two classical optimization test functions. Results show not only a more accurate model but also a possibility to reduce the number of sampling points. Finally, the training strategy was experimented for the optimal design of a sliding bearing. Friction power loss per unit load capacity was taken as the objective function to be modeled. Two surrogate model based optimizations had been per-formed based on classical DoE such as Orthogonal Latin Hypercube and Kriging sequential strategy respectively. These two optimization results are compared with a previously optimization using Complex-optimization method. Within a limited number of iterations, the Kriging model based training strategy showed the best convergence to the global optimum among those 3 methods.
Keywords: Kriging    Surrogate model    Sequential Design    Hybrid Bearing    Optimization    
近20年来,仿真计算技术迅速发展,很多工程实际问题的计算量可通过试验设计的方法进行控制.传统试验设计包括正交设计、均匀试验设计(UD)[1]、拉丁超立方设计(LHS)[2]、最优拉丁超立方设计(OLH)[3]等.这些方法都一次性给定试验点数目,使试验点对设计空间有较好的填充特性.但这些方法没有考虑目标函数特性,不同设计问题,对试验点的分布要求不同,不能一概而论.此外,这些方法都一次性固定试验点数量,加点或者取点不足将不能满足其取点规则.

为解决上述问题,利用Kriging模型[4]不仅能给出对目标函数的预估值,还能给出预估值的理论误差及其分布情况的独特优势. Jones等人提出的基于Kriging模型和期望值改进准则的有效全局优化(EGO)算法[5-6]已在航空航天[7-8]、结构优化[9-10]等领域广泛应用.将基于“加点准则优化”的Kriging序贯设计[11]引入到动静压滑动轴承[12]优化设计问题中,利用已有函数信息进一步优化试验点的布局,细化研究目标函数的优选区域并兼顾非优选区域已建立函数的精确度,用以提高轴承的综合性能、降低其摩擦功耗,从而提升旋转机械的传动效率.

旋转机械高效、稳定运转的首要问题之一是处理好承载能力和摩擦功耗之间的关系.单位承载力下摩擦功耗能有效反映这项设计指标,常被作为滑动轴承设计的优化目标.然而,传统优化方法[13]需要大样本的仿真模拟导致优化效果不甚理想.构建考虑几何结构影响的单位承载力下摩擦功耗的单目标优化模型.利用序贯设计的优势,更加有效地发挥试验点对于优化设计的指导作用[14].通过优化加点获取有代表性的样本点并选择新训练样本点进行迭代求解,以提高模型代理精度.最终快速准确获取全局最优解的Kriging模型.

1 Kriging代理模型介绍Kriging方法是一种基于统计学理论的插值方法,其模型包含均值项和随机项.前者为回归函数,后者为符合一定预设统计学规律的随机变量,即

$\hat{y}(x)=\boldsymbol{F}^{\mathrm{T}}(x) \beta+z(x).$

式中:$\hat y$(x)为未知点的函数估计;β为回归系数;FT(x)为关于x的多项式,用于模拟随机过程的期望值函数,通常使用常数、一次或二次多项式; z(x)为随机项,表达模型预测的不确定性,其随机分布规律决定了模型理论误差的分布规律.

随机项的函数模型可以通过事先预设的数据间的相关函数模型来确定.采样点的协方差如下所示:

$E[z(\boldsymbol{w}) z(\boldsymbol{x})]=\sigma^{2} \boldsymbol{R}\left(\theta, \boldsymbol{w}_{j}, \boldsymbol{x}_{j}\right).$

式中R(θ, wj, xj)为样本点wjxj之间的相关函数,工程中通常选用高斯相关函数

$\boldsymbol{R}\left(\theta, \boldsymbol{w}_{j}-\boldsymbol{x}_{j}\right)=\prod\limits_{j=1}^{n} \exp \left(-\theta\left(\boldsymbol{w}_{j}-\boldsymbol{x}_{j}\right)^{2}\right).$

式中θ是高斯相关函数的关键参数,通过优化θ,调节设计点的相关性.

给定m个设计点X=[x1, $\ldots $, xm]T, XRm×n和响应Y=[y1, $\ldots $, ym]T, YRm×1,则任意待测点x的线性无偏估计值${\hat y}$(x)为

$\hat{y}(\boldsymbol{x})=\boldsymbol{c}^{\mathrm{T}} \boldsymbol{Y}=\boldsymbol{F}^{\mathrm{T}}(\boldsymbol{x}) \hat{\boldsymbol{\beta}}+\boldsymbol{r}^{\mathrm{T}}(\boldsymbol{x}) \boldsymbol{R}^{-1}(\boldsymbol{Y}-\boldsymbol{F} \hat{\boldsymbol{\beta}}).$

式中:${\hat \beta }$的最小二乘估计值以及方差估计σ2可分别为

$\hat{\beta}=\left(\boldsymbol{F}^{\mathrm{T}} \boldsymbol{R}^{-1} \boldsymbol{F}\right)^{-1} \boldsymbol{F}^{\mathrm{T}} \boldsymbol{R}^{-1} \boldsymbol{Y},$

$\sigma^{2}=\frac{1}{m}(\boldsymbol{Y}-\boldsymbol{F} \hat{\boldsymbol{\beta}})^{\mathrm{T}} \boldsymbol{R}^{-1}(\boldsymbol{Y}-\boldsymbol{F} \hat{\boldsymbol{\beta}}).$

此外,R是相关矩阵,其表达式为

$\boldsymbol{R}=\left[\begin{array}{ccc}{\boldsymbol{R}\left(\boldsymbol{w}_{1}, \boldsymbol{x}_{1}\right)} & {\cdots} & {\boldsymbol{R}\left(\boldsymbol{w}_{1}, \boldsymbol{x}_{n}\right)} \\ {\vdots} & {\ddots} & {\vdots} \\ {\boldsymbol{R}\left(\boldsymbol{w}_{n}, \boldsymbol{x}_{1}\right)} & {\cdots} & {\boldsymbol{R}\left(\boldsymbol{w}_{n}, \boldsymbol{x}_{n}\right)}\end{array}\right].$

rT(x)是随机模型相关系数矩阵,表示样本点与预测点之间的相关性,可以表示为

$\boldsymbol{r}^{\mathrm{T}}(\boldsymbol{x})=\left(\boldsymbol{R}\left(\boldsymbol{w}_{1}, \boldsymbol{x}_{1}\right), \boldsymbol{R}\left(\boldsymbol{w}_{2}, \boldsymbol{x}_{1}\right), \ldots, \boldsymbol{R}\left(\boldsymbol{w}_{n}, \boldsymbol{x}_{1}\right),\right)^{\mathrm{T}}.$

Kriging模型预估计值的均方差为

${\rm{MSE}}\left\{ {\hat y\left( \mathit{\boldsymbol{x}} \right)} \right\} = {\sigma ^2}\left\{ {1 - {\mathit{\boldsymbol{r}}^{\rm{T}}}{\mathit{\boldsymbol{R}}^{ - 1}}\mathit{\boldsymbol{r}} + {{\left( {1 - {\mathit{\boldsymbol{F}}^{\rm{T}}}{\mathit{\boldsymbol{R}}^{ - 1}}\mathit{\boldsymbol{r}}} \right)}^2}/{\mathit{\boldsymbol{F}}^{\rm{T}}}{\mathit{\boldsymbol{R}}^{ - 1}}\mathit{\boldsymbol{F}}} \right\}.$ (1)

相关参数θ通过最大似然估计计算得到:

$\psi \left( \theta \right) = - \frac{{m\ln \left( {{{\hat \sigma }^2}} \right)}}{2} - \frac{{\ln \left| \mathit{\boldsymbol{R}} \right|}}{2}.$

通过求解ψ(θ)函数的极大值可得到最优θ值,使得Kriging模型达到最佳精度.

2 基于序贯设计的模型训练本文采用试验设计以提高Kriging全局代理精度和寻找全局最优点,其核心在于新训练样本点的选择以及Kriging模型收敛判据的确定.序贯设计具体步骤如下:

1) 根据工程优化设计要求,确定研究对象并建立相应数学模型、约束条件,并确认设计变量及取值范围.

2) 初始样本选取.通过OLH获取能较好反映设计空间的空间特性的初始样本,并数值求解得到对应的响应,从而得到Kriging模型的启动样本.

3) 建立代理模型.选择二次回归函数和高斯相关函数,并根据步骤2)获取的初始样本集,运用DACE工具箱[15]构建初始Kriging代理模型.

4) 加点准则.通过选择合理的加点法则循环选择新增训练样本点进行迭代求解,用以提高Kriging模型精度.

5) 收敛判定.根据代理模型精度(最优点附近的相对误差)作为收敛准则.如果模型满足收敛条件,则代理模型满足精度要求,建模过程即终止,可以基于模型进入优化阶段.否则进入下一次循环返回步骤4).

2.1 加点准则建立代理模型后,通过一定加点法则选择新的样本点.基于Kriging模型的算法[16]中常用的加点准则包括改善期望准则[17](EI)、均方差准则[18](MSE)和最小化代理模型预测准则(MSP)等.

MSE准则是直接运用Kriging代理模型提供的均方差估计(式(1))最大处作为新样本点的位置.该方法对提高模型全局精度,避免优化结果收敛于局部最优有较好的效果,但收敛的速度偏慢.

EI准则是综合考虑预测值以及方差加权的一种高效全局优化(EGO)方法.设当前最优点响应值为Ymin,且Kriging代理模型预测值满足均值${\hat y}$(x)和均方差σ2的正态分布,其概率密度为

$P\left( \mathit{\boldsymbol{x}} \right) = \frac{1}{{\sqrt {2\pi } \sigma \left( \mathit{\boldsymbol{x}} \right)}}\exp \left[ { - \frac{{{{\left( {{Y_{\min }} - I - \hat y\left( \mathit{\boldsymbol{x}} \right)} \right)}^2}}}{{2{\sigma ^2}\left( \mathit{\boldsymbol{x}} \right)}}} \right].$

针对最小化优化问题,目标函数响应值增量为

$I\left( \mathit{\boldsymbol{x}} \right) = {Y_{\min }} - y\left( \mathit{\boldsymbol{x}} \right).$

则响应值目标提高的期望值为

$\left[ {I\left( \mathit{\boldsymbol{x}} \right)} \right] = \sigma \left( \mathit{\boldsymbol{x}} \right)\left[ {u\mathit{\Phi }\left( u \right) + \varphi \left( u \right)} \right],$

$u = \frac{{{Y_{\min }} - \hat y\left( \mathit{\boldsymbol{x}} \right)}}{{\sigma \left( \mathit{\boldsymbol{x}} \right)}}.$

式中:Φ$φ$分别表示标准正态分布函数和标准正态概率密度函数.

通过计算目标函数改善的期望,把期望最大值点作为新训练样本点,从而寻找全局最优解的优化方法.该方法有助于提高优化的局部精度,在全局最优点所在区域确定的前提下,这种方法可提高优化结果的精度.

综合考虑模型精度和优化结果的精度,以及单独使用一种加点准则收敛速度较慢、效率不高的缺点, 使用“MSE+EI”结合的“并行加点准则[8]”,即一次优化迭代增加2个样本点.不仅可以保证代理模型全局精度,还可以更有效逼近全局最优解.

2.2 收敛判定根据代理模型在最优点附近的相对误差作为优化终止的收敛条件:

$\left| {\hat f\left( {{x_k}} \right) - f\left( {{x_k}} \right)} \right| \le {\varepsilon ^ * }f\left( {{x_k}} \right).$ (2)

式中:k为优化迭代次数; ${\hat f}$(xk)为当前代理模型预测值; f(xk)为当前精确数值模拟响应值; ε*是人为给定的收敛精度,阀值根据采样点仿真计算的数值误差来确定,当模型在最优点处的相对误差小于设定ε*=0.001时,即认为代理模型满足精度要求.

3 序贯设计方案的理论测试使用2个经典的优化测试函数(六峰值驼背函数和Rosenbrock函数),对LHS、OLH以及均匀设计和优化加点后的Kriging模型的全局精度以及寻优结果进行比较.

3.1 测试函数六峰值驼背函数(Six-Hump Camel Back function):

$f\left( {{x_1},{x_2}} \right) = x_1^2\left( {4 - 2.1x_1^2 + \frac{{x_1^4}}{3}} \right) + {x_1}{x_2} + x_2^2\left( { - 4 + 4x_2^2} \right).$

图 1(a)为六峰值驼背函数等值线图,此函数由一个四阶和一个六阶的一维多项式函数联立而成,并含有二阶混合项. 2个全局最优点:(0.0898,-0.7127)和(-0.0898 0.7127),4个局部最优点,最小响应值fmin=-1.029 8.

Fig. 1
图 1 测试函数等值线图 Fig. 1 Test function contour map


Rosenbrock函数:

$f\left( {{x_1},{x_2}} \right) = {\left( {1 - {x_1}} \right)^2} + 100{\left( {{x_2} - x_1^2} \right)^2}.$

与六峰值驼背函数等值线图不同,Rosenbrock函数图像如图 1(b)关于x1=0并不对称.此函数仅有一个全局最优点,位于(1, 1),对应响应值fmin=0,全局最优点位于一个近似抛物线形的梯度值极小的“深谷”中.此函数对于检验优化算法的寻优效率有很好的效果.

3.2 方案评价标准采用均方根误差对不同试验设计方案下建立的Kriging代理模型进行全局误差分析.将二维设计空间等分为n×n的网格,选取网格节点为待测点,并建立如下形式的均方根误差模型评价代理模型误差:

${\rm{RMSE}} = \sqrt {\sum\limits_{i = 1}^{{n_t}} {{{\left( {{e^{\left( i \right)}}} \right)}^2}/{n_i}} } ,{e^{\left( i \right)}} = \left\| {\left( {\hat y_{\rm{t}}^{\left( i \right)} - y_{\rm{t}}^{\left( i \right)}} \right)} \right\|.$ (3)

式中:nt为检测样本点的个数,${\hat y}$t(i)为Kriging在测试点i的预测值,yt(i)为该点真实输出响应.比较函数在不同试验方法下均方根误差(RMSE),来衡量Kriging代理模型的精度.

3.3 不同方案结果分析在测试函数建模过程中,只改变试验设计方法,而样本点数、加点准则以及误差分析方法均保持一致.首先利用DACE工具箱进行初始Kriging建模,再使用“并行加点准则”不断更新训练样本点优化模型,当代理模型满足收敛条件(式(2))时,试验设计加点后样本点分布情况如图 2、图 3所示,其中,“实心圆”点表示初始样本点,“空心圆”点表示通过“并行加点准则”新添加的样本点,“五角星”表示实际全局最优点.

Fig. 2
图 2 六峰值驼背函数序贯设计样本点分布情况 Fig. 2 Sample distribution of sequential design for Six-Hump Camel Back function


Fig. 3
图 3 Rosenbrock函数序贯设计样本点分布情况 Fig. 3 Sample distribution of sequential design for Rosenbrock function


2个测试函数均采用OLH设计构建初始点集,采用并行加点准则增加新样本点,直到代理模型收敛.利用式(3)评价各测试函数下模型精度,结果如表 1、2所示.

表 1
表 1 各试验设计方案Kriging模型均方根误差汇总表(六峰值驼背函数) Tab. 1 RMSE of Kriging model for different DoE methods(Six-Hump Camel Back function) 设计方法方案样本点数目

17 21 25 29

LHS 1.089 5 1.046 5 0.939 5 0.818 9

OLH 0.749 1 0.735 0 0.718 7 0.664 0

均匀设计 0.707 3 0.773 3 0.715 6 0.633 5

序贯设计 0.749 1 0.762 2 0.652 1 0.547 3



表 1 各试验设计方案Kriging模型均方根误差汇总表(六峰值驼背函数) Tab. 1 RMSE of Kriging model for different DoE methods(Six-Hump Camel Back function)


表 2
表 2 各试验设计方案Kriging模型均方根误差汇总表(Rosenbrock函数) Tab. 2 RMSE of Kriging model for different DoE methods(Rosenbrock function) 设计方法方案样本点数目

17 23 29 37

LHS 243.306 8 150.500 4 120.451 4 61.703 7

OLH 133.430 7 51.754 0 42.401 6 27.469 5

均匀设计 195.896 2 110.470 4 78.455 7 115.052 3

序贯设计 133.430 7 47.620 4 18.293 4 1.993 6



表 2 各试验设计方案Kriging模型均方根误差汇总表(Rosenbrock函数) Tab. 2 RMSE of Kriging model for different DoE methods(Rosenbrock function)


分析表 1可知,初始样本采用OLH方案,因此,二者有相同的误差.由表 1中误差数据可知,虽然在增加4个点后,相对于OLH设计,序贯设计的误差有所升高,这种现象可能的原因是测试函数波动较大,在梯度较大的区域需要多个设计点来捕捉函数值的变化,导致某些增加的试验点距离较近,进而改变了相关矩阵R的条件数,影响矩阵求逆运算的精度,最终导致模型精度有一定程度降低.但随着继续加点,模型精度则明显提高,而增至29点时,相对于其他设计的优势可明显体现出来.而OLH和均匀设计的模型精度变化规律类似,随着点数增加,精度略微改善.相较之下,并行加点准则由于考虑了函数对于自变量的响应特性,加点位置更加有针对性,对于模型精度的提高更有效. LHS、OLH和均匀设计的模型精度变化规律类似,随着点数增加精度有所改善,但改善程度有限.

而Rosenbrock函数同六峰值驼背函数的显著不同是函数值变化平缓,梯度值小,特别是最优点所在区域位于一个梯度极小的谷中,因此序贯设计所增加的点的分布相对较均匀,有利于相关矩阵的运算,从空间相关性上更有利于捕捉函数的全局特性.由表 2可知,样本点数由17增加到37个时,模型精度呈线性升高的趋势,序贯设计的精度明显优于其它设计方案.

当Kriging代理模型满足收敛条件时,比较4种设计方法下代理模型预测全局最优解能力,结果如表 3、4所示.对于六峰值驼背函数,基于LHS设计的模型全局最优解距离真实最优解较远,而基于OLH、均匀设计和序贯设计的模型预测的全局最优点位置接近,后3种设计中,序贯设计预测最优点对应的响应值更接近真实最优解.对于Rosenbrock函数可得类似的结论.

表 3
表 3 各试验设计方案Kriging模型全局最优解汇总表(六峰值驼背函数) Tab. 3 Global optimal solution of Kriging model for different DoE methods(Six-Hump Camel Back function) 设计方法 全局最优点 响应值

LHS (-0.153 8,0.743 6) -1.334 3

OLH (0.051 3,-0.692 3) -1.003 1

均匀设计 (0.051 3,-0.692 3) -0.996 5

序贯设计 (0.051 3,-0.692 3) -1.023 3

理论最优点 (-0.089 8,0.712 7) -1.029 8



表 3 各试验设计方案Kriging模型全局最优解汇总表(六峰值驼背函数) Tab. 3 Global optimal solution of Kriging model for different DoE methods(Six-Hump Camel Back function)


表 4
表 4 各试验设计方案Kriging模型全局最优解汇总表(Rosenbrock函数) Tab. 4 Global optimal solution of Kriging model for different DoE methods(Rosenbrock function) 设计方法 全局最优点 响应值

LHS (0.564 1,0.354 1) -1.458 3

OLH (1.282 1,1.589 7) -0.763 4

均匀设计 (1.487 2,2.000 0) -30.481 3

序贯设计 (0.871 8,0.769 2) 0.024 90

理论最优点 (1, 1) 0



表 4 各试验设计方案Kriging模型全局最优解汇总表(Rosenbrock函数) Tab. 4 Global optimal solution of Kriging model for different DoE methods(Rosenbrock function)


综合分析得,序贯设计考虑设计空间的空间填充特性以及目标函数的函数特性,加点位置更加有针对性,对模型精度的提高更有效.所以序贯设计的模型误差降低速度明显高于传统试验设计.由各试验设计方案的模型全局最优解情况可知,序贯设计有更好地收敛到优化问题的真实最优解的能力.这是由于使用“并行加点准则”克服了寻优算法容易陷入局部最优的缺点,改善了对设计空间搜索的全局性.

4 动静压滑动轴承优化设计随着滑动轴承的应用范围不断向着高速、重载方向发展,需要处理好轴承承载能力和摩擦功耗之间的关系.轴承的摩擦功耗不仅带来能量损失,其导致的润滑油温度增加直接影响润滑油粘度,大大降低轴承的综合性能.当前滑动轴承领域中常采用随机方向搜索法、复合形法等常规优化方法对轴承结构尺寸进行优化设计.根据以上比较分析结果,本文将基于Kriging代理模型的序贯设计方法引进到滑动轴承的优化设计中.选取动静压轴承单位承载力下摩擦功耗最小为目标函数进行优化,并同传统复合形法优化结果进行比较.

以圆柱动静压轴承(见图 4)为研究对象,对其几何结构参数进行序贯设计优化,其中轴承重要参数为:主轴转速为10 000 r/min、腔数为4个、偏心率为0.25、直径为70 mm、深腔包角为15°、浅腔深度0.022 mm、深腔深度0.25 mm、半径间隙0.03 mm.选择以单位承载力下的摩擦功耗最小为目标函数,设计变量取长径比l/d、轴向两侧封油边宽度z1, z2和浅腔包角θq.分别建立了相同工况条件下基于序贯设计和基于OLH设计的Kriging代理模型并与复合形法进行对比分析.

$\min F\left( x \right) = \left( {\frac{{{{\bar H}_f} + {{\bar H}_p}}}{{{{\bar F}_r}}}} \right).$

Fig. 4
图 4 带深浅腔的圆柱动静压滑动轴承 Fig. 4 Cylindrical hybrid bearing with deep and shallow pocket


式中:F(x)为目标函数;Hf为无量纲摩擦功耗;Hp为无量纲泵功耗;Fr为无量纲承载力.

${\rm{s}}.\;{\rm{t}}.\left\{ \begin{array}{l}l/d \in \left[ {0.5,0.8} \right]\\{z_1},{z_2} \in \left[ {1,7} \right]\\{\theta _q} \in \left[ {44,53} \right]\end{array} \right..$

式中:l/d为长径比;z1, z2为轴向两侧封油边宽度(mm);θq为浅腔包角(°).

浅腔包角和封油面的大小决定动压效应的作用区域,区域越大,承载力越大,摩擦功耗也越大.而长径比的大小决定了轴承侧流量,进而影响摩擦功耗.轴承长径比对轴承性能有较大影响,一般为了提高运行稳定性和减少轴承尺寸,且使摩擦阻力减小、温升下降,取用较小值.但若取值过小,将造成轴承的承载能力迅速下降等不良后果.然而,这些参数对设计目标的影响不可单一而论,各参数的耦合效应难以直观地从物理意义上进行分析,基于代理模型的优化设计有助于设计空间的全局性寻优.

由不同试验设计方案对代理模型精度的影响分析,OLH设计和均匀设计的模型精度优于LHS设计方案的模型.在此选择OLH设计和Kriging序贯设计进行比较.

首先,采用OLH方法选取25个初始样本点,并采用Kriging“并行加点准则”(EI+MSE准则),建立基于序贯设计的代理模型.此过程中设计变量及每组设计变量对应的目标函数真实值变化如表 5所示.

表 5
表 5 Kriging序贯设计建模计算过程 Tab. 5 Kriging modeling process based on sequential design 样本点数序贯建模设计变量

l/d z1, z2(mm) θq(°) F(x)

25 0.8 5.0 53.000 35.250

31 0.8 5.2 52.633 35.156

37 0.8 4.2 51.714 36.215

43 0.8 6.9 53.000 34.632



表 5 Kriging序贯设计建模计算过程 Tab. 5 Kriging modeling process based on sequential design


比较Kriging模型在最优点处的预测函数值与精确数值计算结果,随着样本点数增加,代理模型在最优点处的误差不断减小,经9次迭代,当样本点数增加至43时,Kriging代理模型在最优点处的相对误差ε*达到8.7×10-4.

按照OLH试验设计准则给出样本点数同为43的试验设计方案,并基于这些点的数值计算结果建立代理模型,求得此模型的最优点.采用复合形法初始样本点为4个,在经过39次迭代共计43个样本点情况下得到最优点.与Kriging序贯设计所建立的模型最优点对比见表 6.

表 6
表 6 不同试验设计方案动静压轴承的优化结果对比 Tab. 6 Comparison of optimization results of Hybrid Bearing with different DoE methods 设计变量 l/d z1, z2(mm) θq(°) 目标函数仿真计算结果

复合形法 0.794 5.857 48.276 36.158

OLH设计 0.800 4.800 52.630 35.457

Kriging序贯设计 0.800 6.900 44.000 34.631



表 6 不同试验设计方案动静压轴承的优化结果对比 Tab. 6 Comparison of optimization results of Hybrid Bearing with different DoE methods


对比OLH试验设计、Kriging序贯设计和复合形法的最优点参数值和对应的精确数值计算结果, 在样本点数相同的情况下,Kriging序贯设计所建立的模型能够找到更好的设计结果,相对基于OLH设计的优化结果,目标函数即单位承载力下的摩擦功耗降低了2.33%.在相同的样本点数的情况下,Kriging序贯设计相对于复合形优化方法目标函数降低了4.22%.由于模型预测的时间较短,相对于样本点仿真计算的时间几乎可忽略,因此可近似认为,样本点数正比于寻优时间的长短,可直接用于衡量寻优效率.相对于传统一次成型的试验设计以及传统优化方法,序贯试验设计能提高模型的全局精度,进而直接提高模型的寻优效率.

5 结论采用并行加点准则及相应收敛条件构建基于Kriging序贯设计的代理模型,与基于传统试验设计的模型进行比对.通过测试函数和动静压轴承的工程应用,验证了Kriging序贯设计的先进性和可靠性,揭示了序贯优化设计能有效利用设计空间的空间填充特性和目标函数的函数特性,进而提高了模型的全局和局部精度,有效缩短了基于高精度数值计算的优化设计周期.


参考文献
[1] FANG K T, LIN D K J, WINKER P, et al. Uniform design: theory and application[J]. Technometrics, 2000, 42(3): 237. DOI:10.1080/00401706.2000.10486045


[2] GIUNTA A, WOJTKIEWICZ S, ELDRED M. Overview of modern design of experiments methods for computational simulations[C]//41st Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 2003: 649.


[3] ZHAO L, CHOI K, LEE I. A metamodeling method using dynamic kriging and sequential sampling[C]//13thAIAA/ISSMO Multidisciplinary Analysis Optimization Conference.Fort Worth: AIAA, 2010: 9391.


[4] KLEIJNEN J P C. Design and analysis of simulation experiments[C]//International Workshop on Simulation. Springer, Cham, 2015: 3


[5] 孙志礼, 李瑞, 闫玉涛, 等. 一种用于结构可靠性分析的Kriging学习函数[J]. 哈尔滨工业大学学报, 2017, 49(07): 146.
SUN Zhili, LI Rui, YAN Yutao, et al. A Kriging based learning function for chemical reliability analysis[J]. Journal of Harbin Institute of Technology, 2017, 49(07): 146. DOI:10.11918/j.issn.0367-6234.201604121


[6] JONES D R, SCHONLAU M, WELCH W J. Efficient global optimization of expensive black-box functions[J]. Journal of Global Optimization, 1998, 13(4): 455. DOI:10.1023/A:1008306431147


[7] 韩忠华. Kriging模型及代理优化算法研究进展[J]. 航空学报, 2016, 37(11): 3197.
HAN Zhonghua. Kriging surrogate model and its application to design optimization: a review of recent progress[J]. Acta Aeronautica ET Astronautica Sinica, 2016, 37(11): 3197. DOI:10.7527/S1000-6893.2016.0083


[8] 韩少强, 宋文萍, 韩忠华, 等. 基于梯度增强型Kriging模型的气动反设计方法[J]. 航空学报, 2017, 38(07): 138.
HAN Shaoqiang, SONG Wenping, HAN Zhonghua, et al. Aerodynamic inverse design method based on gradient-enhanced Kriging model[J]. Acta Aeronautica ET Astronautica Sinica, 2017, 38(07): 138. DOI:10.7527/S1000-6893.2016.0299


[9] 李小刚, 程锦, 刘振宇, 等. 基于双层更新Kriging模型的机械结构动态特性稳健优化设计[J]. 机械工程学报, 2014, 50(03): 165.
LI Xiaogang, CHENG Jin, LIU Zhenyu, et al. Robustoptimization for dynamic characteristics of mechanical structures based on double renewal Kriging model[J]. Journal of Mechanical Engineering, 2014, 50(03): 1653. DOI:10.3901/JME.2014.03.165


[10] 张宇, 覃刚, 张云清, 等. 基于克里金元模型的离心泵水力性能多目标优化[J]. 华中科技大学学报(自然科学版), 2015, 43(04): 54.
ZHANG Yu, ZHAO Gang, ZHANG Yunqing, et al. Kriging based multi-objective optimization for hydraulic performance of centrifugal pump[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2015, 43(04): 54. DOI:10.13245/j.hust.150411


[11] HUANG D, ALLEN T T, NOTZ W I. Sequential Kriging optimization using multiple-fidelity evaluations[J]. Structural and Multidisciplinary Optimization, 2006, 32(5): 369. DOI:10.1007/s00158-005-0587-0


[12] 郭红, 岑少起, 张绍林.圆柱、圆锥动静压滑动轴承设计[D].郑州: 郑州大学, 2013: 260
GUO Hong, CEN Shaoqi, ZHANG Shaolin. Cylindrical and conical hydrostatic sliding bearing design[D]. Zhengzhou: Zhengzhou University, 2013: 260


[13] 王丹丹.独立供油径推浮环轴承性能优化分析与软件包开发[D].郑州: 郑州大学, 2010.
WANG Dandan. Perormance optimization analysis of the journal-thrust floating ring bearing with independent oil supply and software development[D]. Zhengzhou: Zhengzhou University, 2010.


[14] ZHANG Z, DEMORY B, HENNER M, et al. Space infill study of Kriging meta-model for multi-objective optimization of an engine cooling fan[C]//ASME Turbo Expo 2014 : Turbine Technical Conference and Exposition. 2014. DOI: 10.1115/GT2014-25281


[15] LOPHAVEN S N, NIELSEN H B, S?NDERGAARD J. DACE: a matlab kriging toolbox[M]. IMM, Informatics and Mathematical Modelling, The Technical University of Denmark, 2002.


[16] FORRESTER A I J, SOBESTER A, KEANE A J. Engineering design via surrogate modelling: a practical guide[M]. Wiley, 2008.


[17] 高月华, 王希诚. 基于Kriging代理模型的多点加点序列优化方法[J]. 工程力学, 2012, 29(4): 90.
GAO Yuehua, WANG Xicheng. A sequential optimization method with mult-point sampling criterion based on Kriging surrogate model[J]. Engineering Mechanics, 2012, 29(4): 90.


[18] 黄章俊, 王成恩. 基于Kriging模型的涡轮盘优化设计方法[J]. 计算机集成制造系统, 2010, 16(5): 905.
HUANG Zhangjun, WANG Chengen. Turbine discs optimization design based on Kriging model[J]. Computer Integrated Manufacturing Systems, 2010, 16(5): 905. DOI:10.13196/j.cims.2010.05.11.huangzhj.007



相关话题/设计 优化 方案 机械工程学院 结构

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 316H堆焊UNS N10003合金参数优化、组织和硬度的研究
    316H堆焊UNSN10003合金参数优化、组织和硬度的研究杨飞1,2,黎超文2,李志军2,蒋力2,叶祥熙2,刘芳1(1.上海理工大学材料科学与工程学院,上海200093;2.中国科学院上海应用物理研究所,上海201800)摘要:研究异种合金焊接可以降低熔盐堆结构材料的成本并确保其安全性,本文采用钨 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 前驱体温度对激光化学气相沉积YBa2Cu3O 7δ超导薄膜结构及性能的影响
    前驱体温度对激光化学气相沉积YBa2Cu3O7δ超导薄膜结构及性能的影响张琼1,赵培1,吴慰1,戴武斌1,GOTOTakashi2,徐源来3(1.等离子体化学与新材料湖北省重点实验室(武汉工程大学),武汉430205;2.东北大学金属材料研究所,沈阳160001;3.绿色化工过程省部共建教育部重点实 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 泡沫金属三明治结构压印-粘接复合接头剥离性能分析
    泡沫金属三明治结构压印-粘接复合接头剥离性能分析张杰,何晓聪,雷蕾,初明明,刘可欣,黄炎宁(昆明理工大学机电工程学院,昆明650500)摘要:为研究三明治结构压印-粘接复合接头的抗剥离性能,选取AA5052铝合金板以及泡沫镍、泡沫铜以及泡沫铁镍进行压印-粘接复合连接,对接头进行拉伸剪切试验,采用扫描 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 静不定结构的热应力分析
    静不定结构的热应力分析杨益航1,王德志2,林高用2(1.福建省功能材料及应用重点实验室(厦门理工学院),福建厦门361024;2.中南大学材料科学与工程学院,长沙410083)摘要:分析了受约束机构内材料热应力存在的危害,给出静不定结构内热应力的定量描述方式:弹性变形状态时,热应力仅与线膨胀系数、弹 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 铀铌合金结构及性能调控研究进展
    铀铌合金结构及性能调控研究进展蒙大桥,陈向林,陈冬,赵雅文,黄河,法涛(中国工程物理研究院材料研究所,四川绵阳621900)摘要:铀铌合金作为一种重要的核工程材料,因其较高的密度、优异的耐蚀性能和良好的力学性能等特点,被广泛应用于核工业领域。铀铌合金受成分及热处理工艺影响显著,表现出复杂的相转变和组 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 核壳结构纳米镁基复合储氢材料研究进展
    核壳结构纳米镁基复合储氢材料研究进展张秋雨1,2,3,邹建新1,2,3,任莉1,2,3,马哲文1,2,3,朱文1,2,3,丁文江1,2,3(1.上海交通大学轻合金精密成型国家工程研究中心,上海200240;2.上海交通大学材料科学与工程学院,上海200240;3.上海交通大学氢科学中心,上海2002 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 电冲击处理(EST)对TC11钛合金微结构和力学性能的影响研究
    电冲击处理(EST)对TC11钛合金微结构和力学性能的影响研究谢乐春1,2,刘畅1,2,华林1,2(1.现代汽车零部件技术湖北省重点实验室(武汉理工大学),武汉430070;2.汽车零部件技术湖北省协同创新中心(武汉理工大学),武汉430070)摘要:为寻求优化钛合金组织和力学性能的新思路,本文采用 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 高氮过饱和奥氏体表面改性层结构研究进展
    高氮过饱和奥氏体表面改性层结构研究进展车宏龙,王克胜,梁健,雷明凯(大连理工大学材料科学与工程学院表面工程实验室,辽宁大连116024)摘要:本文针对氮过饱和奥氏体改性层结构本质的研究现状,归纳了近年来国内外相关研究,重点阐述和总结了高氮过饱和奥氏体层的相结构、微结构,以及有序化结构3个方面的相关研 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • ECAP形变对高纯铝微结构及冲击层裂损伤的影响
    ECAP形变对高纯铝微结构及冲击层裂损伤的影响佘其海,李超,钟政烨(材料先进技术教育部重点实验室(西南交通大学),成都610031)摘要:在高应变率冲击载荷下,金属材料的主要失效方式之一是层裂损伤。为探讨微结构对层裂损伤的影响,本文利用等径角挤压(EqualChannelAngularPressin ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 响应面分析法优化不锈钢激光切割工艺参数
    响应面分析法优化不锈钢激光切割工艺参数李永亮1,2,3,王敬1,3,梁强1,3(1.重庆工商大学机械工程学院,重庆400069;2.重庆工商大学工程训练中心,重庆400069;3.制造装备机构设计与控制重庆市重点实验室(重庆工商大学),重庆400069)摘要:为了获得良好的不锈钢激光切割质量,确定合 ...
    本站小编 哈尔滨工业大学 2020-12-05