李伟平+窦现东+王振兴+柳超
文章编号:16742974(2014)05003207
收稿日期:20130719
基金项目:国家高新技术研究发展计划(863计划)资助项目(2012AA111802)
作者简介:李伟平(1971-),男,湖南邵阳人,湖南大学副教授,博士
通讯联系人,Email: lwpzlbb@yeah.cn
摘 要: 提出基于误差反向传播神经网络(Back Propagation Neural Network, BPNN)的高维模型表示(high dimensional model representation, HDMR)方法,即BPNNHDMR方法.BPNNHDMR方法的显著优势是将BP神经网络的非线性函数逼近能力与高维模型的层级结构理论相结合来构建近似模型,并能够揭示输入变量之间固有的线性或非线性关系及其相关性,将构造模型复杂度由指数级增长降阶为多项式级,有效地解决了高维建模问题.通过测试和对比验证了该算法的效率和建模能力,并将该算法应用于矿用自卸车安全驾驶室翻车保护装置(RollOver Protective Structure, ROPS)的优化设计.通过优化结果验证了所提方法的可行性和有效性.
关键词:近似模型;高维模型;误差反向传播神经网络;非线性;结构优化
中图分类号:U463.4 文献标识码:A
BPNNHDMR Nonlinear Metamodeling Technique
and Its Application
LI Weiping, DOU Xiandong, WANG Zhenxing, LIU Chao
(State Key Laboratory of Advanced Designed and Manufacture for Vehicle Body,
Hunan Univ, Changsha, Hunan 410082,China)
Abstract:This paper proposed a new highdimension model representation (HDMR) based on back propagation neural network (BPNN), which is called BPNNHDMR. The most remarkable advantage of this method lies in its ability to integrate the nonlinear function approximation capability of BP neural network and the hierarchy structure theory of high dimensional model to build an approximation model. Moreover, this method can reveal the inherent linearity or nonlinearity relationship as well as correlation with respect to input variables. The problem of modeling high dimension model is effectively tackled by reducing the computation cost from exponential growing to polynomial. Testing and comparative analysis confirm the efficiency and capability of BPNNHDMR for high dimension nonlinear problems. Furthermore, the algorithm was applied to optimize the ROPS of Mining Dump Truck's Safety Cab. The optimized results verify the feasibility and effectiveness of the method proposed.
Key words: metamodel; high dimensional model representation(HDMR); back propagation neural network(BPNN); nonlinearity; structural optimization
工程优化问题中,对于基于真实模型的嵌套优化,每次计算目标函数值都要调用费时的仿真计算模型,其计算代价不可小视.而利用近似模型可以有效地解决这一问题[1],即通过对近似模型的求优近似得到真实模型的优化值.
目前广泛应用的一些近似方法,如响应面法、Kriging插值、人工神经网络等,在处理较低维问题时有很好的效果.而对于工程中复杂的高维非线性模型,随着维数和非线性程度的增加,构造近似模型所需的样本点数量和计算花费呈指数增长,使解决此类问题的计算效率大大降低[2] .针对这一问题,Sobol证明了可积函数可以分解为不同维数函数的叠加理论[3].该理论表明,对于任意一个可以积分的函数,在其积分空间内存在唯一的、可以扩展的高维模型(High Dimensional Model Representation,HDMR).这一模型是精确的,有确定的阶数,并包含一个层级结构.高维模型可以将计算时间随非线性程度和维数增加按指数增长的隐函数,转化为可以忽略高阶耦合项的多项式函数,并揭示了每个设计变量对近似函数的贡献量,大大减少了计算时间.同时反映了输入变量之间固有的线性或非线性关系及其相关性,在近似高维非线性问题时非常有效.由此,一系列不同特性的高维模型开始发展起来,被研究和应用于不同的领域.其中,Shan等提出了基于径向基的高维模型(RBFHDMR)[4],汤龙等提出了基于Kriging的高维模型(KrigingHDMR)[5].
本文采用误差反向传播神经网络(Back Propagation Neural Network, BPNN)与中心切面高维模型(CutHDMR[6])方法相结合的BPNNHDMR.BP神经网络是一种多层前馈型神经网络,其神经元的传递是S型函数,它可以实现从输入到输出的任意非线性映射.对于非线性问题,在精度表达上具有一定优势.CutHDMR方法是用过指定点的特定的直线、平面和超平面上的信息来建立模型,计算效率高,方便易行.
1 高维模型(HDMR)基本理论
设待求问题的设计变量可行域为An(n维实数空间),那么多元函数f(x)∈R与输入变量x∈An之间的映射关系可以用HDMR[6-7]来表示为:
f(x)=f0+∑ni=1fi(xi)+∑1≤i
其中f0为函数在中心点的函数值,后面依次为不同阶耦合项对近似函数的贡献量.
为了方便计算,本文引入CutHDMR.与其他类型的高维模型相比,CutHDMR用少量简单的算术运算来表达计算花费高昂的真实模型,并达到了其他类型高维模型相似的精度,计算效率高.CutHDMR展开式的各分项表达请参考文献[6].
2 BPNNHDMR
2.1 BP神经网络
BP网络是一种多层前馈型神经网络,它由一个输入层、一个输出层和至少一层隐含层组成[8].该网络的主要特点是信号向前传递,误差反向传播.在前向传递中,输入信号从输入层经隐含层逐层处理,直至输出层.每一层的神经元状态只影响下一层神经元状态.如果输出层得不到期望输出,则转入反向传播,根据预测误差调整网络权值和阈值,从而使BP神经网络预测输出不断逼近期望输出.因此,BP网络可以看作是解决函数逼近的工具.
本文中,隐层神经元的传递函数采用双曲正切S型函数h,它类似于一个平滑的阶梯函数:
h(x)=(ex-e-x)/(ex+e-x) (2)
考虑一个分层神经网络(x,W),其中x是输入矢量,W是可变权系数矢量.在这里要训练(x,W)来近似函数f:[a,b]R→R,其中f([a,b])是一个有界集.在区间[a, b]中随机选取xk作为神经网络的输入,(xk,Wk)为网络的输出,f(xk)为在区间f([a,b])中的期望输出值.神经网络的任务是调整网络中的可变权系数来减小误差Ek,Ek定义为:
Ek=12[xk,W(k)]-f(xk)2(3)
假设wi(k)为W(k)中的任意元素.调整wi(k)对误差Ek的影响直接取决于偏导数(Ek)/[wi(k)].BP算法的作用在于如何评价(Ek)/[wi(k)].详细内容可参考文献[9].这样,wi(k)就向着减小误差Ek的方向调整:
wi(k+1)=wi(k)-μ(Ek)/[wi(k)],
wi(k)∈W(k)(4)
其中μ为常量,是指定的更新率.可以看出,如果μ足够小时,则:
{[xk,W(k+1)]-f(xk)}2<
{[xk,W(k)]-f(xk)}2(5)
其中假设Ek>0.由式(5)可知,对于同样的输入参数x1,x2,…,xk,更新后的神经网络的输出更接近于真实函数输出值.
HechtNielsen[9]研究了分层神经网络近似非线性函数的能力.他在文献[9]中表明,在确定的条件下,对于任意ε>0,存在一个三层神经网络(包含一个输入层,一个输出层和一个非线性隐含层)可以在均方误差精确度ε范围内近似函数f.因此,BP神经网络可以应用于工程实际中的大多数情况.
2.2 BPNNHDMR的构建
对大多数工程问题而言,非耦合项和低阶耦合项对响应函数影响较大.为此,BPNNHDMR方法只考虑到一阶耦合项,表达式如下:
f(x)≈f0+∑ni=1i(xi)+∑1≤i 式中带顶标^的项表示BPNN近似项,∑ni=1i(xi),∑1≤i 高维模型的一般构建流程[4-5]如下: 1)选取各设计变量中心位置的点x0=[x10,x20,…,xn0]作为中心点,计算得到f0. 2)在每个变量xi([x10,x20,…,xi,…,xn0])的取值区间上为非耦合项fi(xi)=f([x10,…,xi,…,xn0]T)-f0布点.选取xi取值区间上下界的两个端点,并计算这两点的函数值,用BPNN构建近似函数i(xi). 3)检验i(xi)的线性关系,如果(x0)-f0f0≤10-3,则认为i(xi)是线性的,程序终止;否则,继续采样和xi上下界三个点构建一个新的i(xi),看i(xi)是否满足给定的精度要求,如果满足,则程序终止;否则,继续采样构建一个新的近似函数i(xi),直至满足精度要求.循环执行第2步和第3步直到所有的非耦合项构造完毕. 4)判断模型中是否存在一阶耦合项.创建新检验点(xi,xj,xij0)=[x10,x20,…,xi,…,xj,…,xn0],不失一般性,选取构造非耦合项时用的采样点分量的上下界中的一个作为新样本的第i维分量.在精确度准则所允许的误差范围内,若f0+∑ni=1fi(xi)=f0+∑ni=1i(xi),就认为不存在一阶耦合项,程序终止;否则转入第5步. 5)识别相互耦合的变量.构造新样本点[x10,…,xei,…,xej,…,xn0]T(1≤i 6)将上面构建的各阶BPNN近似函数代入式(6)就得到了高维模型近似响应函数. 在用BP神经网络构建近似模型过程中,一般隐含层结点数取2P+1[10] (P为输入结点数),由于本文样本数据准确,对于非耦合项和一阶耦合项我们选择的隐结点数分别为3和5.根据文献[11],只有当学习率为0<η<43N+1时,神经网络算法是收敛的,这里η为学习率,N为神经网络训练样本点数.同时,本文中神经网络最大训练步数设定为100次,训练目标为0.000 04,单隐含层.其他更多BPNN具体参数设置和程序设计请参看文献[10]和[12]. 本文高维模型构建时的精确度准则和收敛准则都是通过相对误差来定义的.精确度准则主要用于输入变量之间耦合性的识别,所允许的相对误差一般不超过10-4.由于本文方法具有较高精确度,所以本文针对测试函数收敛准则所允许的相对误差均取0.001;而对于工程问题,考虑到计算效率和一般工程要求的5%的近似精度,设置的相对误差均取0.01.通过以上设置,在样本点数一定时,由于本文方法将BP神经网络的非线性函数逼近能力与高维模型的层级结构理论相结合,所以本文近似模型方法在近似非线性问题时较传统方法具有更高的精度,建模效率也有所提高. 3 数值算例 3.1 评价指标 为了测试近似模型近似效果,采用3种比较常用的评价指标,R2(R square),相对平均绝对误差(relative average absolute error, RAAE)和相对最大绝对误差(relative maximum absolute error, RMAE).这些指标反映了BPNNHDMR在新样本点上的预测能力,其具体表达式请参考文献[4].其中,R2是从整体上反映近似模型的精度,其值越接近1越好;RAAE也是从整体上反映近似模型的精度,其值越小越好;RMAE是一个局部指标,描述了设计空间中某个局部区域的误差,其值越小越好.
3.2 函数算例
首先选择一个高维非线性测试函数:
f(x)=∑9i=1[(x2i)(x2i+1+1)+(x2i+1)(x2i+1)],
xi∈[0,1](7)
采用相同数量的训练样本点(计算费用相同),分别采用BPNN,BPNNHDMR和Kriging 3种方法进行建模并比较它们的精度,计算结果如表1所示(表中数据是计算100次的平均值).
表1 BPNN,BPNNHDMR和Kriging方法比较
Tab.1 Comparison between BPNN, BPNNHDMR
and Kriging algorithm
方法
R2
RAAE
RMAE
BPNN
0.886 0
0.221 6
0.225 1
BPNNHDMR
0.993 9
0.058 2
0.049 5
Kriging
0.326 4
2.975 8
0.836 6
分析表中数据,对于同一测试函数,3种近似方法中,BPNNHDMR方法的R2值最接近1,RAAE和RMAE也都最小.由对比可知,对于高维非线性问题,基于同样数量的一组训练样本,采用BPNNHDMR方法得到的近似模型精度更高.
其次,增加高维非线性测试函数式(7)的维数
f(x)=∑d-1i=1[(x2i)(x2i+1+1)+(x2i+1)(x2i+1)],
0≤xi≤1(8)
来测试BPNNHDMR方法的效率.式(8)中,维数d分别取d=10,30,50等不同值,假设每一维的训练样本点数为9(经测试,基本可以满足精度要求),表2中列出了各阶HDMR计算费用的比较.
表2 各阶HDMR建模费用的比较
Tab.2 Comparison of modeling cost among variousorder
HDMRs for the study problem
维数
BPNNHDMR
高维模型二阶项全部
展开(多项式级增长)
全因子设计sd
(指数级增长)
10
172
2 961
3.486 8×109
30
532
28 081
4.239 1×1028
50
892
78 801
5.153 8×1047
4 工程应用
为了验证BPNNHDMR在处理工程实际问题时的可行性,本文以某矿用自卸车全驾驶室式翻车保护装置(ROPS)的优化设计为例.根据标准ISO 3471:2008[13]对ROPS进行非线性有限元分析,确保ROPS满足标准ISO 3164:1992[14].通过仿真分析得到ROPS的变形量,为进一步优化打下基础.安全驾驶室ROPS的有限元模型如图1所示.
图1 ROPS加载分析模型
Fig.1 Loads and constraints of the ROPS model
ROPS框架采用壳单元模拟,焊接使用刚性壳单元模拟为一个载荷变换器[15].焊接连接处网格局部加密,计算时假设焊缝与母材材料特性相同(实际中焊缝一般不先破坏).单元大小设定为20 mm,其中四边形单元54 951个,三角形单元1 389个,单元类型分别为S4R和S3R.本文ROPS材料为进口钢板A710,屈服极限为690 MPa,断裂极限为792 MPa,材料模型选择金属塑性材料模型,定义材料的应力应变曲线.依据标准ISO 3471:2008进行载荷和约束的施加.其中,约束ROPS底部与车架连接处3个平动方向的自由度(UX,UY,UZ),FC为侧向力加载,FV为垂向力加载,FL为纵向力加载.标准挠曲极限量(Deflectionlimiting volume ,DLV)与ROPS的侧向间距为200 mm,垂向间距为110 mm,纵向间距为320 mm.分析时考虑几何非线性和材料非线性.
4.1 优化变量的选取
本文优化的ROPS主要由一些不同厚度的矩形管和加强筋焊接而成.考虑对称性和加工可行性,对称的选取ROPS不同部位的矩形管厚度作为设
计变量.各设计变量选取见图1,取值范围如表3所示.
表3 各优化变量取值范围
Tab.3 Range of optimization variablesmm
参数
x1
x2
x3
x4
x5
x6
x7
x8
x9
x10
取值范围
[8,12]
[6,14]
[8,12]
[6,14]
[6,10]
[6,10]
[5,7]
[5,7]
[6,10]
[6,10]
4.2 目标函数的建立
通过非线性有限元分析可知,由于要考虑ROPS侧向吸能问题,ROPS侧向变形最接近DLV,因此侧向变形为危险工况.根据企业要求,本文选取侧向变形作为一个优化目标,以垂向和纵向变形为约束.因此以ROPS侧向变形最小和ROPS质量最小建立多目标优化问题,ROPS的多目标优化模型可以描述为:
min XUC(X)min Xm(X)s.t.UV(X)
式(9)中UC为ROPS的侧向变形量;UV为ROPS的垂直变形量;UL为ROPS的纵向变形量;B1,B2分别为ROPS垂向和纵向的最大允许变形量;m为ROPS的总质量.
4.3 近似模型的构建
本文采用BPNNHDMR构建ROPS变形量和质量响应近似模型.对每个设计变量在其取值范围内等间距取点,取点数量根据构建近似模型时是否达到预先设定的精度决定,通过ABAQUS计算获得样本点处的真实响应值.首先构建一阶项,然后进行精度判断.在设计变量空间,采用拉丁超立方实验设计方法,选择10个采样点.分别对真实模型和近似模型进行计算,并计算两者的相对误差.如果相对误差满足给定精度要求(设定为5%),则构建近似模型成功,若精度不满足要求,则继续构建高阶项.
通过计算,本文所求近似模型只需计算到一阶项,就可以得到比较精确的近似结果,可以省略高阶耦合项.构建过程中只需56个样本点即得到满足工程要求的近似模型.
为了客观、全面地反映所建近似模型在设计域内的精确程度,在设计域内采用拉丁超立方实验设计方法随机生成10个测试样本点,用它们分别对BPNN,BPNNHDMR和Kriging模型进行精度测试,并以侧向变形响应精度为例,测试结果如表4所示.
表4 采用BPNN,BPNNHDMR和Kriging
方法所得侧向变形近似模型的比较
Tab.4 Comparison of the lateral deformation between
BPNN, BPNNHDMR and Kriging algorithm
方法
样本点数
R2
RAAE
BPNN
56
0.603 3
0.298 9
BPNNHDMR
56
0.984 0
0.085 1
Kriging
56
0.238 3
0.716 7
由表4中数据可知,BPNNHDMR解决此类工程问题具有高效准确的优点.同时,可以看出BPNN和Kriging方法在构建高维非线性有限元模型时的局限性,本文需要增加采样点数量才能用这两种方法获得比较精确的近似模型,这样就增加了真实有限元模型的计算次数,而真实非线性有限元模型一次计算要数小时,导致计算花费大大增加.
4.4 优化过程及结果分析
本文采用基于Pareto概念的多目标优化遗传算法,该方法是求解多目标问题非劣最优解的有效途径之一 [16].在基于Pareto最优概念的遗传算法中,NSGAⅡ[17-18]是最有效的.因此,本文采用NSGAⅡ,在近似模型的基础上对ROPS进行多目标优化.
初始种群设为100,最大迭代次数设为200,交叉概率0.9,变异概率0.1.收敛规则为:达到最大迭代次数作为终止条件.得到Pareto最优解集如图2所示.
图2中星号表示多目标优化非劣解.在优化解集中取有代表性的10组解,多目标优化的Pareto最优解集见表5.ROPS优化前侧向变形和质量如表6所示.对比表6与表5中第4组和第6组数据可以看出,在质量相当的情况下,优化后的侧向变形大约减小了26.7 mm,在侧向变形相当的情况下,优化后的质量大约减小了75.1 kg.因此,通过多目标优化,ROPS的质量和变形情况明显改善.同时,根据得到的Pareto最优解集,可以根据设计者经验和需求,高效率地实现ROPS各矩形管厚度的选取,以满足不同性能需要.
把优化结果代入有限元计算模型验证,计算结果如表7所示.
由ABAQUS计算验证可知,所选优化解都符合标准ISO 3164:1992的要求.这也证明了本文所提方法在工程优化实际应用中的可行性和有效性.同时,由于模型简化掉了一些蒙皮和附属部件,以及真实情况下车架也有一部分吸能作用,所以本文的分析结果是偏向安全的.
侧向变形/mm
图2 多目标优化的Pareto最优解集
Fig.2 Feasible Pareto optimal solutions
表5 双目标优化的Pareto最优解
Tab.5 The Pareto optimal solutions of twoobjective uncertain optimization
Pareto解
x1
/mm
x2
/mm
x3
/mm
x4
/mm
x5
/mm
x6
/mm
x7
/mm
x8
/mm
x9
/mm
x10
/mm
质量
/kg
侧向变形
/mm
1
12.0
6.0
12.0
13.2
10.0
6.7
9.9
5.8
7.0
9.8
1 508.5
88.0
2
11.9
6.0
11.5
10.3
10.0
6.6
9.8
5.8
6.8
9.9
1 453.1
94.5
3
11.8
6.0
10.7
10.6
10.0
6.5
9.6
5.9
6.6
9.9
1 421.2
100.9
4
12.0
8.4
8.5
10.1
10.0
6.9
7.8
6.4
5.0
9.9
1 347.9
116.4
5
11.0
8.5
8.0
9.6
9.3
6.4
7.5
6.4
5.0
9.5
1 284.3
136.5
6
9.9
6.0
9.4
8.6
9.5
7.0
9.0
6.0
5.0
7.5
1 269.8
143.9
7
9.6
6.2
8.8
6.3
9.3
7.7
6.4
5.3
5.4
7.1
1 215.1
152.6
8
10.3
6.0
8.9
6.8
7.6
6.5
6.9
5.3
5.1
8.3
1 188.2
164.2
9
10.0
6.1
8.7
6.0
7.5
6.5
7.0
5.6
5.2
8.6
1 168.5
171.5
10
10.0
6.0
8.0
6.3
6.7
6.5
7.2
5.3
5.0
8.3
1 125.0
186.3
表6 ROPS侧向变形优化前数据
Tab.6 Finite element analysis results of ROPS before optimized
变量
x1
/mm
x2
/mm
x3
/mm
x4
/mm
x5
/mm
x6
/mm
x7
/mm
x8
/mm
x9
/mm
x10
/mm
质量
/kg
侧向变形
/mm
初始值
10
10
10
10
8
8
8
6
6
8
1 344.9
143.1
表7 优化后的响应值ABAQUS计算验证
Tab.7 Check the optimization results
by ABAQUS calculation mm
优化后最大
侧向变形量
ABAQUS
计算的侧向
最大变形
ABAQUS
计算的垂向
最大变形
ABAQUS
计算的纵向
最大变形
解一
88.0
92.22
35.70
33.64
解二
94.5
98.99
34.77
35.85
解三
100.9
102.9
36.47
39.74
解四
116.4
112.8
46.09
51.62
解五
136.5
134.1
60.10
68.26
解六
143.9
142.9
62.63
74.29
解七
152.6
159.0
73.41
88.32
解八
164.2
167.1
69.32
82.35
解九
171.5
174.0
70.30
85.04
解十
186.3
192.1
78.05
100.6
5 结 论
本文提出的BPNNHDMR建模方法,很好地利用了BP神经网络的非线性函数逼近能力与高维模型的层级结构理论,并能够反映输入变量之间固有的线性关系和耦合性.尤其对于高维问题,它可以将构造模型计算花费由维数的指数级增长降解为多项式级,有效地解决了高维建模问题,且与传统算法相比,具有更高的精度.通过数值算例和工程优化问题的对比,验证了BPNNHDMR近似方法的精度和效率.另外,该算法对于高度非线性问题精确的数学建模还需进一步发展和完善,BP神经网络的计算效率也有待进一步提高.
同时,本文所采用的结构分析与优化方法在工程领域的其他方面也可以广泛应用,具有一定的理论和工程实际意义.
参考文献
[1] 穆雪峰. 多学科设计优化代理模型技术的研究和应用[D].南京:南京航空航天大学航空宇航学院,2004.
MU Xuefeng. The research and application of multidiscipline design optimization surrogate model technology[D]. Nanjing: College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, 2004.(In Chinese)
[2] 赵子衡,韩旭,姜潮.基于近似模型的非线性区间数优化方法及其应用[J].计算力学学报,2010,27(3):451-456.
ZHAO Ziheng, HAN Xu, JIANG Chao. Approximation model based nonlinear interval number optimization method and its application [J]. Chinese Journal of Computational Mechanics, 2010, 27(3):451-456. (In Chinese)
[3] SOBOL I M. Sensitivity estimates for nonlinear mathematical models [J]. Mathematical Modeling & Computational Experiment, 1993, 1(4): 407-414.
[4] SHAN Songqing,WANG G Gary. Metamodeling for high dimensional simulationbased design problems[J]. Journal of Mechanical Design, 2010, 132:051009.
[5] 汤龙,李光耀,王琥. KrigingHDMR非线性近似模型方法[J].力学学报,2011,43(4):780-784.
TANG Long, LI Guangyao, WANG Hu. KrigingHDMR metamodeling technique for nonlinear problems [J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(4):780-784. (In Chinese)
[6] RABITZ H, AL1S F. General foundations of highdimensional model representations[J]. Journal of Mathematical Chemistry,1999,25(2/3):197-233.
[7] LI G, WANG S W, ROSENTHAL C, et al. High dimensional model representations generated from low dimensional data samples. I. MpCutHDMR[J]. Journal of Mathematical Chemistry, 2001,30(1): 1-30.
[8] CHEN Fuchang .Backpropagation neural networks for nonlinear self tuning adaptive control[J]. Control Systems Magazine, IEEE April, 1990, 10(3): 44-48.
[9] HECHTNIELSEN R.Theory of the backpropagation neural network[C]//Proc IEEE 2002 Intl Conf Neural Networks.Califormia:International Toint Conference in Nemal Networks, 2002.
[10]谢庆生,尹健,罗延科,等.机械工程中的神经网络方法[M].北京:机械工业出版社,2003: 39-51.
XIE Qingsheng, YIN Jian, LUO Yanke,et al. Neural network method in mechanical engineering [M].Beijing:China Machine Press, 2003:39-51. (In Chinese)
[11]曾喆昭,李仁发.高阶多通带滤波器优化设计研究[J].电子学报,2002,30(1):87-89.
ZENG Zhezhao, LI Renfa. Study on the optimum design of the highorder multibandpass filters [J].Acta Electronica Sinica, 2002, 30(1):87-89. (In Chinese)
[12]张德丰. MATLAB 神经网络应用设计[M]. 北京:机械工业出版社,2012:49-75.
ZHANG Defeng. MATLAB neural network application design [M]. Beijing:China Machine Press, 2012:49-75. (In Chinese)
[13]ISO 3471:2008 Earthmoving machineryrollover protective structureslaboratory tests and performance requirements[S]. British: Standards Policy and Strategy Committee,2009.
[14]ISO 3164: 1992Earthmoving machinery —Laboratory evaluation of rollover and fallingobject protective structuresspecifications for deflectionlimiting volume[S]. British: The Authority of the Standards Board, 1992.
[15]周传月. MSC.Fatigue疲劳分析应用与实例[M]. 北京:科学出版社,2005:86-91.
ZHOU Chuanyue. Fatigue analysis of application and examples with MSC.Fatigue[M]. Beijing:Science Press, 2005:86-91. (In Chinese)
[16]谢涛,陈火旺,康立山.多目标优化的演化算法[J].计算机学报,2003,26(8):997-1003.
XIE Tao, CHEN Huowang, KANG Lishan. Evolutionary algorithms of multiobjective optimization problems [J].Chinese Journal of Computers, 2003, 26(8):997-1003. (In Chinese)
[17]DEB K, PRATAP A, AGARWAL S, et al. A fast elitist nondominated sorting genetic algorithm for multiobjective optimization: NSGAII,Kan GAL Report No.200001[R]. Kanpur: Indian Istitute of Technology, 2000.
[18]李伟平,王世东,周兵,等. 基于响应面法和NSGAⅡ算法的麦弗逊悬架优化[J]. 湖南大学学报:自然科学版, 2011,38(6):27-32.
LI Weiping, WANG Shidong, ZHOU Bin, et al. Macpherson suspension parameter optimization based on response surface method and NSGAⅡ algorithm [J]. Journal of Hunan University:Natural Sciences, 2011,38(6):27-32. (In Chinese)
同时,本文所采用的结构分析与优化方法在工程领域的其他方面也可以广泛应用,具有一定的理论和工程实际意义.
参考文献
[1] 穆雪峰. 多学科设计优化代理模型技术的研究和应用[D].南京:南京航空航天大学航空宇航学院,2004.
MU Xuefeng. The research and application of multidiscipline design optimization surrogate model technology[D]. Nanjing: College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, 2004.(In Chinese)
[2] 赵子衡,韩旭,姜潮.基于近似模型的非线性区间数优化方法及其应用[J].计算力学学报,2010,27(3):451-456.
ZHAO Ziheng, HAN Xu, JIANG Chao. Approximation model based nonlinear interval number optimization method and its application [J]. Chinese Journal of Computational Mechanics, 2010, 27(3):451-456. (In Chinese)
[3] SOBOL I M. Sensitivity estimates for nonlinear mathematical models [J]. Mathematical Modeling & Computational Experiment, 1993, 1(4): 407-414.
[4] SHAN Songqing,WANG G Gary. Metamodeling for high dimensional simulationbased design problems[J]. Journal of Mechanical Design, 2010, 132:051009.
[5] 汤龙,李光耀,王琥. KrigingHDMR非线性近似模型方法[J].力学学报,2011,43(4):780-784.
TANG Long, LI Guangyao, WANG Hu. KrigingHDMR metamodeling technique for nonlinear problems [J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(4):780-784. (In Chinese)
[6] RABITZ H, AL1S F. General foundations of highdimensional model representations[J]. Journal of Mathematical Chemistry,1999,25(2/3):197-233.
[7] LI G, WANG S W, ROSENTHAL C, et al. High dimensional model representations generated from low dimensional data samples. I. MpCutHDMR[J]. Journal of Mathematical Chemistry, 2001,30(1): 1-30.
[8] CHEN Fuchang .Backpropagation neural networks for nonlinear self tuning adaptive control[J]. Control Systems Magazine, IEEE April, 1990, 10(3): 44-48.
[9] HECHTNIELSEN R.Theory of the backpropagation neural network[C]//Proc IEEE 2002 Intl Conf Neural Networks.Califormia:International Toint Conference in Nemal Networks, 2002.
[10]谢庆生,尹健,罗延科,等.机械工程中的神经网络方法[M].北京:机械工业出版社,2003: 39-51.
XIE Qingsheng, YIN Jian, LUO Yanke,et al. Neural network method in mechanical engineering [M].Beijing:China Machine Press, 2003:39-51. (In Chinese)
[11]曾喆昭,李仁发.高阶多通带滤波器优化设计研究[J].电子学报,2002,30(1):87-89.
ZENG Zhezhao, LI Renfa. Study on the optimum design of the highorder multibandpass filters [J].Acta Electronica Sinica, 2002, 30(1):87-89. (In Chinese)
[12]张德丰. MATLAB 神经网络应用设计[M]. 北京:机械工业出版社,2012:49-75.
ZHANG Defeng. MATLAB neural network application design [M]. Beijing:China Machine Press, 2012:49-75. (In Chinese)
[13]ISO 3471:2008 Earthmoving machineryrollover protective structureslaboratory tests and performance requirements[S]. British: Standards Policy and Strategy Committee,2009.
[14]ISO 3164: 1992Earthmoving machinery —Laboratory evaluation of rollover and fallingobject protective structuresspecifications for deflectionlimiting volume[S]. British: The Authority of the Standards Board, 1992.
[15]周传月. MSC.Fatigue疲劳分析应用与实例[M]. 北京:科学出版社,2005:86-91.
ZHOU Chuanyue. Fatigue analysis of application and examples with MSC.Fatigue[M]. Beijing:Science Press, 2005:86-91. (In Chinese)
[16]谢涛,陈火旺,康立山.多目标优化的演化算法[J].计算机学报,2003,26(8):997-1003.
XIE Tao, CHEN Huowang, KANG Lishan. Evolutionary algorithms of multiobjective optimization problems [J].Chinese Journal of Computers, 2003, 26(8):997-1003. (In Chinese)
[17]DEB K, PRATAP A, AGARWAL S, et al. A fast elitist nondominated sorting genetic algorithm for multiobjective optimization: NSGAII,Kan GAL Report No.200001[R]. Kanpur: Indian Istitute of Technology, 2000.
[18]李伟平,王世东,周兵,等. 基于响应面法和NSGAⅡ算法的麦弗逊悬架优化[J]. 湖南大学学报:自然科学版, 2011,38(6):27-32.
LI Weiping, WANG Shidong, ZHOU Bin, et al. Macpherson suspension parameter optimization based on response surface method and NSGAⅡ algorithm [J]. Journal of Hunan University:Natural Sciences, 2011,38(6):27-32. (In Chinese)
同时,本文所采用的结构分析与优化方法在工程领域的其他方面也可以广泛应用,具有一定的理论和工程实际意义.
参考文献
[1] 穆雪峰. 多学科设计优化代理模型技术的研究和应用[D].南京:南京航空航天大学航空宇航学院,2004.
MU Xuefeng. The research and application of multidiscipline design optimization surrogate model technology[D]. Nanjing: College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, 2004.(In Chinese)
[2] 赵子衡,韩旭,姜潮.基于近似模型的非线性区间数优化方法及其应用[J].计算力学学报,2010,27(3):451-456.
ZHAO Ziheng, HAN Xu, JIANG Chao. Approximation model based nonlinear interval number optimization method and its application [J]. Chinese Journal of Computational Mechanics, 2010, 27(3):451-456. (In Chinese)
[3] SOBOL I M. Sensitivity estimates for nonlinear mathematical models [J]. Mathematical Modeling & Computational Experiment, 1993, 1(4): 407-414.
[4] SHAN Songqing,WANG G Gary. Metamodeling for high dimensional simulationbased design problems[J]. Journal of Mechanical Design, 2010, 132:051009.
[5] 汤龙,李光耀,王琥. KrigingHDMR非线性近似模型方法[J].力学学报,2011,43(4):780-784.
TANG Long, LI Guangyao, WANG Hu. KrigingHDMR metamodeling technique for nonlinear problems [J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(4):780-784. (In Chinese)
[6] RABITZ H, AL1S F. General foundations of highdimensional model representations[J]. Journal of Mathematical Chemistry,1999,25(2/3):197-233.
[7] LI G, WANG S W, ROSENTHAL C, et al. High dimensional model representations generated from low dimensional data samples. I. MpCutHDMR[J]. Journal of Mathematical Chemistry, 2001,30(1): 1-30.
[8] CHEN Fuchang .Backpropagation neural networks for nonlinear self tuning adaptive control[J]. Control Systems Magazine, IEEE April, 1990, 10(3): 44-48.
[9] HECHTNIELSEN R.Theory of the backpropagation neural network[C]//Proc IEEE 2002 Intl Conf Neural Networks.Califormia:International Toint Conference in Nemal Networks, 2002.
[10]谢庆生,尹健,罗延科,等.机械工程中的神经网络方法[M].北京:机械工业出版社,2003: 39-51.
XIE Qingsheng, YIN Jian, LUO Yanke,et al. Neural network method in mechanical engineering [M].Beijing:China Machine Press, 2003:39-51. (In Chinese)
[11]曾喆昭,李仁发.高阶多通带滤波器优化设计研究[J].电子学报,2002,30(1):87-89.
ZENG Zhezhao, LI Renfa. Study on the optimum design of the highorder multibandpass filters [J].Acta Electronica Sinica, 2002, 30(1):87-89. (In Chinese)
[12]张德丰. MATLAB 神经网络应用设计[M]. 北京:机械工业出版社,2012:49-75.
ZHANG Defeng. MATLAB neural network application design [M]. Beijing:China Machine Press, 2012:49-75. (In Chinese)
[13]ISO 3471:2008 Earthmoving machineryrollover protective structureslaboratory tests and performance requirements[S]. British: Standards Policy and Strategy Committee,2009.
[14]ISO 3164: 1992Earthmoving machinery —Laboratory evaluation of rollover and fallingobject protective structuresspecifications for deflectionlimiting volume[S]. British: The Authority of the Standards Board, 1992.
[15]周传月. MSC.Fatigue疲劳分析应用与实例[M]. 北京:科学出版社,2005:86-91.
ZHOU Chuanyue. Fatigue analysis of application and examples with MSC.Fatigue[M]. Beijing:Science Press, 2005:86-91. (In Chinese)
[16]谢涛,陈火旺,康立山.多目标优化的演化算法[J].计算机学报,2003,26(8):997-1003.
XIE Tao, CHEN Huowang, KANG Lishan. Evolutionary algorithms of multiobjective optimization problems [J].Chinese Journal of Computers, 2003, 26(8):997-1003. (In Chinese)
[17]DEB K, PRATAP A, AGARWAL S, et al. A fast elitist nondominated sorting genetic algorithm for multiobjective optimization: NSGAII,Kan GAL Report No.200001[R]. Kanpur: Indian Istitute of Technology, 2000.
[18]李伟平,王世东,周兵,等. 基于响应面法和NSGAⅡ算法的麦弗逊悬架优化[J]. 湖南大学学报:自然科学版, 2011,38(6):27-32.
LI Weiping, WANG Shidong, ZHOU Bin, et al. Macpherson suspension parameter optimization based on response surface method and NSGAⅡ algorithm [J]. Journal of Hunan University:Natural Sciences, 2011,38(6):27-32. (In Chinese)