专利名称：一种高效的拉丁超立方试验设计方法
技术领域：
本发明涉及一种高效的拉丁超立方试验设计方法，属于工程优化设计技术领域。
背景技术：
当今许多工程设计问题越来越复杂，为了提高分析精度和可信度，许多高精度分 析模型被广泛用于工程设计中。高精度分析模型在提高分析精度和可信度的同时也带来了 计算耗时的困难，虽然当今计算机软硬件技术已经有了长足的发展，然而，调用高精度分析 模型完成一次分析仍然极其耗时，例如使用CFD模型完成一次气动仿真分析需要数小时甚 至数十小时。其次，现代工程设计问题往往涉及多个相互耦合的学科。譬如，飞行器设计涉 及气动、结构、动力、隐身、控制等学科，各学科相互影响，相互制约，飞行器的性能是各学科 耦合的综合体现。由于学科之间的耦合关系，工程设计问题的系统分析表现为多学科分析。 在寻优和设计分析中，为了避免高计算量，常常使用代理模型(Metamodeling)代替这些真 实的耗时的高精度模型。构造代理模型之前，必须先在设计空间获取能够反应真实模型的试验仿真设计 点，然后通过这些设计点构造代理模型。一个可行的计算试验设计必须满足空间均勻性和 投影均勻性，拉丁超立方设计(Latin Hypercube Design, LHD) (McKay et al.，1979)就具 有很好的一维投影均勻性。拉丁超立方是一种分层抽样方法，试验点在设计空间内均勻分 布，试验次数等于水平数，适用于因素个数较多的情况。为了建立投影均勻性的量化标准， Morris和Mitchell (1995)介绍了基于准的优化拉丁超立方；Park(1994)介绍了基于 最小距离最大准则(maximin criterion)优化拉丁超立方，保证了 LHD样本点的空间均勻 性；方开泰(2002)介绍了基于中心 L2 偏差标准(Centered L2 discrepancy criterion) 现在拉丁超立方试验设计主要分为一次采样方法(one-stage sampling method) 和多次采样方法(multi-stage sampling method)。一次采样方法表示在设计空间中一次 采取足够多的样本点构造代理模型。多次采样方法就是在设计空间中逐次采取样本点。与 一次采样相比，多次采样能够使用较少的样本点找到真实目标模型的局部最优点或全局最 优点，但是多次采样所构造的代理模型在整个设计空间中精度不够，采样过程中需要建立 优化模型进行寻优得到样本点，效率较低。Dam(2007)提出一种一次采点方法，采用二维设 计中基于最小距离最大的拉丁超立方法，此方法不需要优化过程，效率较高，但此方法只限 于二维问题，且受采样点个数的限制，不具有很好的通用性。
发明内容
本发明的目的是为解决工程设计的现有技术中采样效率低，通用性差的问题， 提出一种高效的拉丁超立方试验设计(Efficient maximin Latin Hypercube Design without optimization process for n-dimensional problems,ELHD)方法。本方法不需 要优化过程即可得到满足空间均勻性和投影均勻性特点的试验仿真设计点，采样效率高， 且不受采样点个数的限制，具有很好的通用性。
本发明的目的是通过以下技术方案实现的步骤1，建立用户在实际工程中η维设计空间的优化设计问题所对应的空间超体， 预设第1个试验设计点。为了实现用户从实际工程中的η维设计空间中提取m个试验设计点的目的，将η 维空间构造成为一个#的空间超体；再将此空间超体分为相邻的m个单位超体，从一端到 另一端的空间顺序编为第1个，第2个，……，第m个，则每个单位超体中有πΓ1个单位正 方体。若将m个试验设计点分别放入某一个单位超体中的某一个单位正方体内，且每一 个单位超体只能放一个点，则使得获取的m个试验设计点更易满足投影均勻性原则。并且， 为使得速度最快，设计按从一端到另一端的顺序进行。在第1个单位超体中的任意一个单位正方体中随机设置第1个试验设计点 ^(广，&尸，…，-，…次’，/二！，…，…其中，…表示第丄个试验设计点在第丄个单位超体中
第i维的坐标值。将得到的P1存入设计点集合P = {PJ。步骤2，利用步骤1得到的第1个试验设计点，在第2个单位超体中确定第2个试 验设计点。为了更好的满足投影均勻性，第2个试验设计点P2的设计原则为各维坐 标值与第1个试验设计点P1在相应各维上的坐标值均不相同。因此只能在除去已被 乃(0丨2)，…，片V..,^)，/ = !,…，ZJ占据的各维余下坐标值所对应的单位正方体中选取。 找出满足上述各维坐标值均不相同条件的单位正方体与第1个试验设计点P1之间的距离 最大者，即 Cl1= max(Cl(P2iP1))；将该点坐标⑷”,殆2),.·.，#)，·.·,始“))，/= 1,.·，《，对)定为 P2所在位置。将P2存入设计点集合P = {P1; P2I。步骤3，利用步骤1和步骤2得到的第1个和第2个试验设计点，在第3个单位超 体中确定第3个试验设计点。为了更好的满足投影均勻性，设计第3个试验设计点P3在第3个单位超体的除 已被巧(片),斤),...,片),一,衫"))和 2树”，片),...，坊),..,片)),/ = 1，-.., ，占据的各维余下坐
标值对应的单位正方体中选取。余下单位正方体中，同一个单位正方体分别与Pi、P2距离 中的较小值，作为这个单位正方体的特征值；再找出当前单位超体即第3个单位超体中，余 下单位正方体特征值中的最大值，将此最大值对应的空间坐标值确定为第3个试验设计点 P3 Of,…,时),…Λ("))"· = 1,…,《，ki'^kl^k^。将 P3 存入设计点集合 P = (PijP2jP3I。步骤4，对于第s个点(s = 4,5,…，m-1)，利用已得到的s_l个试验设计点，在第 s个单位超体中确定第s个试验设计点。为了更好的满足投影均勻性，设计第s个试验设计点Ps在第s个单位超 体的除已被巧(k\l\k\2\-·Λ \·-,k\n)) ’ P1 沐⑴，片)，…，对)，···，#))，... ρ“ (0,^，···,^,.··,0, / = 1,···, , OO…关#关对0占据的余下坐标值对应的
单位正方体中选取。余下的单位正方体中，同一个单位正方体分别与前s-l个试验设计点 距离中的最小值，作为这个单位正方体的特征值；再找出当前单位超体即第s个单位超体 中，各单位正方体特征值的最大值，将此最大值对应的空间坐标值确定为第s个试验设计 点乃。将得到的Ps存入设计点集合P =iPl，'"Ps-I' PsI。步骤5，重复步骤4，直到(m-1)个试验设计点被确定。步骤6，在完成步骤5后，第m个单位超体中满足各维坐标值均不相同条件的单位 正方体只剩一个，故第m个试验设计点唯一确定。至此，从η维优化设计空间中提取m个试验设计点的目的达到。上述步骤1中在第1个单位超体中的任意一个单位正方体中预设第1个试验设计 点的方案，也可以改为从第m个单位超体中的任意一个单位正方体中预设第1个试验设计 点，按从第m个单位超体到第1个单位超体的方向进行试验设计点的设计，而其设计方法不 变。有益效果本发明的ELHD试验设计方法与现有技术相比，不需要优化过程即可得到满足空 间均勻性和投影均勻性特点的试验试验设计点，采样效率高；获取的试验设计点特性好，适 合应用到任何工程设计等运算量巨大的优化领域，如含有大规模有限元分析的结构优化设 计、含有高精度流体力学分析的气动优化设计等。用本方法采集的设计点可以构造出各种代理模型，且构造的代理模型精度高；适 用于任意维数的试验设计问题，有较强的通用性，尤其在解决复杂的现代工程优化设计中 的高维问题时，本方法节约时间的效果尤为明显。在基于代理模型对多种飞行器总体优化 设计时的试验设计、基于代理模型对多种流体分析模型进行优化设计时的试验设计、以及 基于代理模型对大型有限元分析系统进行优化设计时的试验设计等航空航天领域或相关 工程优化设计领域都可以有较为广泛的应用。
图1为本发明的高效的拉丁超立方试验设计方法的流程图；图2为具体实施方式
中在4X4的网格平面确定4个试验设计点空间位置的示意 图；图3为具体实施方式
中ELHD采样结果与matlab maximin-LHD采样结果对比(η =2，m = 30)；图4为具体实施方式
中ELHD采样结果与matlab maximin-LHDD采样结果对比(η =2，m = 50)；图5为具体实施方式
中ELHD采样结果与matlab maximin-LHD采样结果对比(η =2，m = 100)；图6为具体实施方式
中ELHD采样结果与matlab maximin-LHD采样结果对比(η =2，m = 200)；图7为具体实施方式
中ELHD采样结果及其在二维平面的投影(η = 3，m = 30)；图8为具体实施方式
中ELHD采样结果及其在二维平面的投影(η = 3，m = 50)；图9为具体实施方式
中ELHD采样结果及其在二维平面的投影(n = 3，m = 100)；图10为具体实施方式
中ELHD采样结果及其在二维平面的投影(n = 3，m= 200)；图11为具体实施方式
中ELHD采样结果在三维空间的投影(η = 4，m = 50)；图12为具体实施方式
中ELHD采样结果在二维平面的投影(η = 4，m = 50)；
图13为具体实施方式
中结构桁架系统的有限元模型。
具体实施例方式为了进一步说明本发明的目的和优点，下面结合附图和具体实施例对本发明作进 一步的说明。图1为本发明的高效的拉丁超立方试验设计(ELHD)方法的流程图。以二维问题为例，具体说明ELHD方法的技术方案。在二维设计空间中(η = 2)，采样问题可描述为一个棋子问题。把m个棋子(采样 点)分散放到mXm方格的棋盘(二维设计空间)上，保证任意两个棋子之间的距离最小者 的最大化。为了保证采样点的空间均勻性和投影均勻性，在ELHD方法中，按相邻的列划分 m个单位超体，即第1列、第2列、……第m列；并逐列确定采样点的位置，每一行和列都只 有一个点。ELHD方法的详细采样过程如下步骤1，第1个点在第1列中随机产生。假定第1个点随机产生在第I1行，记第1 个点为P1G^l)，I1 e {1,2,…，m}，同时样AP1存入样点集合P = {PJ。步骤2，对于第2列中的第2个点，因为第I1行已经被第1个点P1占据，所以第2 个点只能从余下的行号里选择。找出余下各行的方格与第1个点P1之间的距离中最大者， BP d2 = max ((I((I2^)jP1)), i2 e {1,2,…，ifLi^l，-mi，对应的行号用于确定第 2 个 点P2 (i2，2)，将P2存入样点集合P= {P” P2I。步骤3，对于第3列中的第3个点，因为第ipi2行已经被第1个点P1、第2个点P2 占据，所以第3个点只能从余下行里选择。找出同一方格分别与余下各行的方格与第1个 AP1、第2个点P2之间的距离的较小者，作为该方格的特征值，再取各方格特征值的最大值， 将此最大值对应的方格位置作为第3个点P3(i3，3)，将P3存入样点集合P = (P1, P2, P3I。步骤4，对于第s列中的第s个点，s e {1，2，···，πι-1}。已生成的点有{Ρ1;Ρ2，…， Ps-J,这些点已经占据掉(S-I)行，行号为Upi2,…，is-J，只能从余下的行号里选择一个 确定第S列，余下的行号有UsIl彡and is ^ I1, i2,…，is-J。找出余下行中同一 方格与先前生成的各样点之间的距离最小值，作为该方格的特征值，再取各方格特征值的 最大值，即 ds = maX(min(d((is，s)，P)))，P= (P1, P2, ...，PS_J。将此特征值最大者对应 的行号作为第s点的行号，即Ps(is，s)，并将此点存入样点集合P = IP1, P2,…，Ps^1, PJ。步骤5，重复步骤4，直到(m-1)个采样点确定。步骤6，对于第m列中的第m个点，行号只剩1个，故第m个点此时唯一确定。上述步骤2、步骤3、步骤4和步骤5中计算两点距离的方式可以为2范数、1范数、 无穷范数中的任意一种。基于ELHD方法在二维空间生成样本点，采样点个数分别取m= 30，50，100，200， 采样结果与matlab标准LHD采样结果对比如图3、图4、图5和图6所示，在空间维数η = 2采样点m = 30、50、100、200的采样结果表明ELHD方法采取的样本点比matlab标准LHD 方法采取的样本点更具有空间均勻性和投影均勻性。下面用二维问题，取采样点m = 4作为实施例说明具体实施过程，过程中采用2范 数计算两点间的距离。如图2所示，在4X 4的网格平面中，确定4个点的空间位置，以保证4个点的空间均勻性和投影均勻性特点。详细生成样本点的过程如下第(1)步在第1列中随机产生第1个点。假定第1个点产生在第3行，记为P1 (3， 1)，并将样点P1存入样点集合P = {PJ。第⑵步对于第2列中的第2个点。因为第3行已经被第IfAP1占据，所以 第2个点只能从余下的行号{1，2，4}里选择。找出余下行号各方格位置与第IfAPii 间的距离的最大者，即d2 = max(d((i2,2),P1)),i2 e {1，2，4}，对应的行号1符合要求，由 此确定第2个点P2 (1,2),将第2个点存入样点集合P = (P1, P2I。第(3)步对于第3列中的第3个点，因为已生成的点有{P1; P2}，已经占据行号 为{3，1}，只能从余下的行号{2，4}里选择一个确定第3个点。找到余下的每行各方格位 置与先前生成各个样点之间的距离最小值作为该方格的特征值，如图2所示，因为d((2， 3)，P2) < (!(口⑶^，则将仏口⑶，⑴=1.414作为方格(2,3)的特征值；d ((4，3)，P1) < d( (4，3)，P2) Jl^fd(GdhP1) = 2. 236作为方格(4,3)的特征值。再取各余下方格特 征值{1. 414，2. 236}的最大者即2. 236所对应的行号4作为第3个点的行号，即P3(4，3)， 并将此点存入样点集合P = IP1, P2, P3I。第(4)步重复第(3)步，直到3个样点确定。此实例中采样点m = 4，故此步直 接跳过。第(5)步对于第4列中的第4个点，行号只剩1个，故第4个点P4(2，4)此时唯
一确定。基于ELHD方法在η维空间生成样本点。当维数η = 3时，采样点个数分别取m = 30，50，100，200，采样结果以及在二维平面内的投影如图7-图10所示；维数η = 4时，采样 点个数m= 50，采样结果在三维空间中的投影及在二维平面内的投影结果如图11-图12所 示；对于更高维空间，ELHD方法采取的样本点在二维平面内同样具有很好的投影均勻性。采用经典的测试方法对采样结果进行测试，基于(^和Q2标准的测试结果如表1 所示。从表中可以看到，基于ELHD方法的采样结果对应的(^和。值都比标准matlab的 采样结果小，故ELHD方法采样点具有更好的空间均勻性的投影均勻性。表1. ELHD 与 matlab maximin-LHD 方法对比结果 (一 )下面通过考察采样点构造工程设计中常用的几个重要代理模型的精度也说 明ELHD方法的优越性。典型的测试函数如下Six-hump camel—back (SC)函数y = 4Xl2-2. Ix14+(1/3) χι6+χιχ2-4χ22+4χ24(1)X1,2 e [-2,2]Brain function 函数y = (χ2-5· Kxi/^JiV+G/jOxfer+lOa-a/SjOcosO^+lO (2)x1 e [-5，10]，X2 e
Peaks function 函数γ = 3(\-Χι)2β{-χ'-{χ^)2)-IOiXlZS-Xii -x52)e(-xf-xl) ~(U3)ewlf-xl)X1,2 e [-3,3]Geometric container function (GC)函数y = 0. 2/x1x2x3+4/x1+3/x3(4)X1j2j3 e
Colville-Himmelblau(CH)函数^ = 100(x2 -X12)2 +(1-x,)2 +90(x4-X32)2+(I-X3)2(5)+10.1(02 -1)2 + (x4 -1)2) + 1 9.8(X2 — 1)(JC4 -1)
, i = 1,2,3,4. Hesse function (HE)函数
-y = 25 (xr2)2+ (x2-2) 2+ (x3-l)2+ (x4-4)2+ (x5-l)2+ (x6-4)2(6)
xi e
, i = 1,2, ...，6
采用常用的函数拟合精度标准RMSE，RAAE, R2和RMAE，对利用采样点构造以上函 数代理模型精度进行测试，测试结果如表2所示。从表2中可以看出，ELHD方法构造代理 模型精度RMSE，RAAE, RMAE结果较matlab maximin-LHD方法小，且ELHD方法构造代理模型精度的R2结果也较mat lab maximin-LHD方法更接近于1，尤其对于高维问题(HE函数)，R2 结果相比matlab提高了约65%，因此ELHD方法获得了更好的采样结果。表2. ELHD与matlab LHD构造代理模型精度对比结果 本发明在上述测试标准和构造代理模型的具体测试表明，ELHD方法的采样效果较 好，具有较好的空间均勻性和投影均勻性，整个过程无需优化，采样效率高，且不受采样点 个数限制，具有很好的通用性。本方法支持N维空间的采样问题，为构造代理模型提供良好 的计算试验设计，适合应用到任何工程设计等运算量巨大的优化领域，如含有大规模有限 元分析的工程结构优化设计、含有高精度计算流体力学分析的气动优化设计等中；用本方 法采集的设计点可以构造出各种工程设计的代理模型，尤其在解决复杂的现代工程优化设 计中的高维问题时，本方法所采集的设计点构造代理模型的精度较好。在航空航天领域或 相关工程优化设计领域都可以有较为广泛的应用。( 二)工程优化设计实例以一个用于工程建筑中简化的结构桁架系统为例，如图13所示，说明拉丁超立方 试验设计方法的应用。该工程设计问题的目标是在两种载荷工况(easel和CaSe2)下，满 足给定的横截面积和刚度约束条件，使系统的总结构质量最小，结构桁架系统中材料的的 基本参数是杨氏模量为E = 7 X IO10Pa密度为P = 2. 7X103kg/m3三个杆的长度分别为Z1= 0.25λ/2, I2 = 0. 25,Z3= 0.25>/2 m两种工况作用下的力载荷为easel :FX = 7. 1187e5N, Fy = _5· 339e6Ncase2 :FX = _7· 1187e5N, Fy = _5· 339e6N根据设计参数，可将该工程优化问题数学描述为下式，即通过优化杆的横截面面 积，保证两种工况(easel和CaSe2)下4号节点的位移都小于2. Oe-4米的同时，使整个桁 架系统结构质量最小化。min mass = (l^Ai+l^A；,+^^) * P
s. t. easel and case2 :dnode4 < = 2. 0e_4m0. 001 彡 A1 彡 0· 1,0. 001 彡 A2 彡 0· 1,0. 001 彡 A3 彡 0. 1其中，mass为桁架系统的结构总质量，dnode4为4号节点的总位移，A1, A2，A3分别为 长度为I1,12,13的杆横截面面积。上述问题为一个结构有限元计算实例，属于非线性约束问题，该问题可以认为是 一个大规模有限元模型计算的简化实例。现采用本发明对该问题所在的设计空间进行采 样，然后利用采样结果进行构造代理模型，再用遗传算法对此代理模型进行寻优，将最优解 代入到有限元模型中计算得到本次优化的最优解。同时选用matlab的maximin LHD方法 进行采样，进行同样的寻优过程，比较两者的寻优结果。建立该工程设计优化数学模型如式(7)所示。min mass = (ΙΑ+Ι^+Ι^) * Ps. t. dnode4 < = 2. 0e-4m(7)0. 001 彡 A1 彡 0· 1,0. 001 彡 A2 彡 0· 1,0. 001 彡 A3 彡 0. 1因为设计变量为3个，故nv = 3。使用ELHD计算实验设计方法在初始设计空间A1, 2,3 e
内产生 m= (3+1) (3+2)/2 = 10 个样本点。步骤1 建立设计变量为横截面积对应的3维设计空间的空间超体，预设第1个试 验设计点。从3维空间中选取10个样本点，首先构造一个10X10X10 = IO3的空间超体； 将此空间超体分为相邻的10个单位超体，从一端到另一端的空间顺序编为第1个，第2 个，……，第10个，则每个单位超体中有IO2个单位正方体。若将10个试验设计点分别放入某一个单位超体中的某一个单位正方体内，且每 一个单位超体只能放一个点，则使得获取的10个试验设计点更易满足投影均勻性原则。并 且，为使得速度最快，设计按从一端到另一端的顺序进行。在第1个单位超体中的任意一个单位正方体中随机设置第1个试验设计点为 P1Gj^)，其中，8表示第1个试验设计点在第1个单位超体中第2维的坐标值。将得到的 P1存入设计点集合P= {PJ。步骤2 利用步骤1得到的第1个试验设计点P1Gj^)，在第2个单位超体中确 定第2个试验设计点。为了更好的满足投影均勻性，第2个试验设计点P2的设计原则为各维坐标值与 第1个试验设计点P1在相应各维上的坐标值均不相同。因此只能在除去已被P1IS,9)占 据的各维余下坐标值所对应的单位正方体中选取。找出满足上述各维坐标值均不相同条件 的单位正方体与第1个点P1之间的距离最大者，即Cl1 = max((KP2jP1))；将该点坐标(2，1， 1)定为P2所在位置。将&(2，1，1)存入设计点集合P= {P1;P2}。步骤3 利用步骤1和步骤2得到的第1个试验设计点P1 (1,8,9)和第2个试验设 计点P2 (2，1，1)，在第3个单位超体中确定第3个试验设计点。为了更好的满足投影均勻性，设计第3个试验设计点P3在第3个单位超体的除已 被？乂丨力，钔和&(2，1，1)占据的各维余下坐标值对应的单位正方体中选取。余下单位正 方体中，同一个单位正方体分别与Pp P2距离中的较小值，作为这个单位正方体的特征值； 再找出当前单位超体即第3个单位超体中，余下单位正方体特征值中的最大值，将此最大
11值对应的空间坐标值确定为第3个试验设计点P3 (3，10，2)。将P3存入设计点集合P= (P1,
P2，P3I ο步骤4 对于第s个试验设计点(s = 4,5,…，9)，利用已得到的s_l个试验设计 点，在第S个单位超体中确定第S个试验设计点。为了更好的满足投影均勻性，设计第s个试验设计点Ps在第s个单位超体 的除已被 P1(LSJ), P2(2，l，l)，…P“ (dCi，."，d,."，0,i = 1， η,
Af *片〗占据的余下坐标值对应的单位正方体中选取。余下的单位正方 体中，同一个单位正方体分别与前s-1个试验设计点距离中的最小值，作为这个单位正 方体的特征值；再找出当前单位超体即第s个单位超体中，各单位正方体特征值的最大 值，将此最大值对应的空间坐标值确定为第s个试验设计点Ps(、(1)^2),…，<°，…乂’， kln *片)*kf本…丰O将得到的Ps存入设计点集合P = {P” P2，P3,…Pp1，PJ。步骤5 重复步骤4，直到9个试验设计点被确定。此步骤的前9个试验设计点为 P1 (1，8，9)，P2 (2，1，1)，P3 (3，10，2)，P4 (4，2，8)，P5 (5，5，3)，P6 (6，9，7)，P7 (7，6，10)，P8 (8，3， 6)，P9 (9,7,4) ο步骤6，在完成步骤5后，第10个单位超体中满足各维坐标值均不相同条件的单位 正方体只剩一个，故第10个试验设计点唯一确定，即Pltl (10，4，5)。如表3所示。至此，从设计变量为横截面积对应的3维空间提取10个试验设计点的目的达到。将此10个采样点归一化后并映射到设计空间中，调用目标函数进行分析获取样 本点处的10个样本点的真实响应值，利用此样本点和响应值构造径向基代理模型。本实 施例选用的径向函数为高斯函数(Gaussian function) 0采用遗传算法对代理模型进行全 局优化，得到此次优化过程的全局近似最优点(0.0321，0. 0923,0.0322)以及真实响应值 mass为123. 6704kg，对应两种工况下节点最大位移分别为2. 0013e_4m和1. 9979e_4m，即约 束符合要求，有限元结构桁架系统的最小质量为123. 6704kg。同样用matlab的maximin LHD方法进行采样得到10个样本点，如表3所示。也 同样调用目标函数进行分析获取样本点处的10个样本点的真实响应值，利用此样本点和 响应值构造径向基代理模型。本实施例选用的径向函数为高斯函数(Gaussian function)。 采用遗传算法对代理模型进行全局优化，得到此次优化过程的全局近似最优点(0. 0443， 0. 0729,0. 0444)以及真实响应值mass为133. 8975kg，对应两种工况下节点最大位移分 别为2. 0016e-4m和1. 9994e-4m，即约束符合要求，有限元结构桁架系统的最小质量为 133.8975kg。表3. ELHD采样点与matlab的maximin LHD采样点对比 表4.利用采样点构造代理模型进行优化结果对比 本发明在上述有限元结构桁架系统的具体测试对比结果如表4，对比表明ELHD采 样方法能找到更小的结构质量最优值为123. 6704kg，Kmatlab Wmaximin LHD采样方法更 好，具有更好的空间均勻性和投影均勻性，通过构造代理模型，并对代理模型进行优化，优 化结果表明ELHD采样方法得到的代理模型精度较好，优化过程中更有利于找到全局最优 解。该采样方法适合应用到任何工程设计等运算量巨大的优化领域，如含有大规模有限元 分析的工程结构优化设计、含有高精度计算流体力学分析的气动优化设计等中；用本方法 采集的设计点可以构造出各种工程设计的代理模型，尤其在解决复杂的现代工程优化设计 中的高维问题时，本方法所采集的设计点构造代理模型的精度较好。在航空航天领域或相 关工程优化设计领域都可以有较为广泛的应用。以上所述的具体描述，对发明的目的、技术方案和有益效果进行了进一步详细说 明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范 围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明 的保护范围之内。
1权利要求
一种高效的拉丁超立方计算试验设计方法，其特征在于包括如下步骤步骤1，建立用户在实际工程中n维设计空间的优化设计问题所对应的空间超体，预设第1个试验设计点；从用户实际工程中的n维设计空间提取m个试验设计点，将n维空间构造成为一个mn的空间超体；再将此空间超体分为相邻的m个单位超体，从一端到另一端的空间顺序编为第1个，第2个，……，第m个，则每个单位超体中有mn 1个单位正方体；若将m个试验设计点分别放入某一个单位超体中的某一个单位正方体内，且每一个单位超体只能放一个点，则使得获取的m个试验设计点更易满足投影均匀性原则；并且，设计按从一端到另一端的顺序进行，避免重复计算，使得速度最快；在第1个单位超体中的任意一个单位正方体中随机生成第1个试验设计点其中，表示第1个试验设计点在第1个单位超体中第i维的坐标值；将得到的P1存入设计点集合P＝{P1}；步骤2，利用步骤1得到的第1个试验设计点，在第2个单位超体中确定第2个试验设计点；为了更好的满足投影均匀性，第2个试验设计点P2的设计原则为各维坐标值与第1个试验设计点P1在相应各维上的坐标值均不相同；因此只能在除去已被占据的各维余下坐标值所对应的单位正方体中选取；通过计算满足上述各维坐标值均不相同条件的单位正方体与第1个试验设计点P1之间的距离，找出距离最大者，即d1＝max(d(P2，P1))；将该点坐标定为P2所在位置；将P2存入设计点集合P＝{P1，P2}；步骤3，利用步骤1和步骤2得到的第1个试验设计点和第2个试验设计点，在第3个单位超体中确定第3个试验设计点；为了更好的满足投影均匀性，设计第3个试验设计点P3在第3个单位超体的除已被和占据的各维余下坐标值对应的单位正方体中选取；分别计算满足条件的第3个单位超体中余下单位正方体与P1、P2之间的距离，并将同一个单位正方体分别与P1、P2距离中的较小值，作为这个单位正方体的特征值；再找出当前单位超体即第3个单位超体中，各单位正方体特征值中的最大值，将此最大值对应的空间坐标值确定为第3个试验设计点将P3存入设计点集合P＝{P1，P2，P3}；步骤4，对于第s个试验设计点(s＝4，5，…，m 1)，利用已得到的s 1个试验设计点，在第s个单位超体中确定第s个试验设计点；为了更好的满足投影均匀性，设计第s个试验设计点Ps在第s个单位超体的除已被占据的余下坐标值对应的单位正方体中选取；分别计算余下的单位正方体与之前确定的s 1个试验设计点之间的距离，并将同一个单位正方体分别与前s 1个试验设计点距离中的最小值，作为这个单位正方体的特征值；再找出当前单位超体即第s个单位超体中，各单位正方体特征值的最大值，将此最大值对应的空间坐标值确定为第s个试验设计点将得到的Ps存入设计点集合P＝{P1，P2，P3，…Ps 1，Ps}；步骤5，重复步骤4，直到(m 1)个试验设计点被确定；步骤6，在完成步骤5后，第m个单位超体中满足各维坐标值均不相同条件的单位正方体只剩一个，故第m个试验设计点唯一确定；至此，从n维优化设计空间中提取m个试验设计点的目的达到。FSA00000230081800011.tif,FSA00000230081800012.tif,FSA00000230081800013.tif,FSA00000230081800014.tif,FSA00000230081800015.tif,FSA00000230081800016.tif,FSA00000230081800021.tif,FSA00000230081800022.tif,FSA00000230081800023.tif,FSA00000230081800024.tif,FSA00000230081800025.tif
2.根据权利要求1所述的一种高效的拉丁超立方计算试验设计方法，其特征在于所 述的步骤1中在第1个单位超体中的任意一个单位正方体中预设第1个试验设计点的方 案，也可以改为从第m个单位超体中的任意一个单位正方体中预设第1个试验设计点，按从 第m个单位超体到第1个单位超体的方向进行试验设计点的设计，而其设计方法不变。
3.根据权利要求1所述的一种高效的拉丁超立方计算试验设计方法，其特征在于所 述的步骤2、步骤3、步骤4和步骤5中计算两点距离的方式为2范数、1范数、无穷范数中的 任意一种。
全文摘要
本发明涉及一种高效的拉丁超立方计算试验设计方法，属于工程优化设计技术领域。本发明为实现用户从实际工程中的n维设计空间提取m个试验设计点的目的，通过建立空间超体，预设第1个试验设计点；利用第1个试验设计点，在第2个单位超体中确定第2个试验设计点；利用已得到的s-1个(s＝3，4，5，…，m-1)试验设计点，在第s个单位超体中确定第s个试验设计点，直到(m-1)个试验设计点被确定；第m个单位超体中满足设计条件的第m个点唯一确定。本发明的设计方法不需要优化过程即可得到满足空间均匀性和投影均匀性特点的试验试验设计点，采样效率高；获取的试验设计点特性好，适合应用到任何工程设计等运算量巨大的优化领域。
文档编号G06F17/50GK101923590SQ20101025392
公开日2010年12月22日 申请日期2010年8月16日 优先权日2010年8月16日
发明者刘莉, 朱华光, 龙腾 申请人:北京理工大学