CN104636563B - 高速压力机上横梁可靠性设计方法 - Google Patents
高速压力机上横梁可靠性设计方法 Download PDFInfo
- Publication number
- CN104636563B CN104636563B CN201510079808.8A CN201510079808A CN104636563B CN 104636563 B CN104636563 B CN 104636563B CN 201510079808 A CN201510079808 A CN 201510079808A CN 104636563 B CN104636563 B CN 104636563B
- Authority
- CN
- China
- Prior art keywords
- reliability
- design
- constraint
- mrow
- entablature
- 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
- 238000013461 design Methods 0.000 title claims abstract description 92
- 238000013459 approach Methods 0.000 title claims abstract description 11
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 238000013401 experimental design Methods 0.000 claims abstract description 7
- 230000002068 genetic effect Effects 0.000 claims description 22
- 238000004422 calculation algorithm Methods 0.000 claims description 19
- 238000000034 method Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000005314 correlation function Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000000205 computational method Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000007476 Maximum Likelihood Methods 0.000 claims description 2
- 230000006978 adaptation Effects 0.000 claims description 2
- 230000002457 bidirectional effect Effects 0.000 claims description 2
- 238000005315 distribution function Methods 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 230000010429 evolutionary process Effects 0.000 claims 1
- 230000008520 organization Effects 0.000 claims 1
- 238000002921 genetic algorithm search Methods 0.000 abstract 1
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000003754 machining Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种高速压力机上横梁可靠性设计方法。包括以下步骤:根据实际上横梁设计中可靠性要求选择设计变量,建立以区间描述不确定性因素的高速压力机上横梁可靠性设计模型;在实验设计中采用LHS获取拟合Kriging所需样本点,通过协同仿真技术获取各样本点所对应的目标函数和约束函数值,并以此构建Kriging模型;基于均布区间优势度计算可靠性设计模型中的可靠性约束值;采用基于区间约束违反度的双层嵌套遗传算法搜寻符合可靠性要求的最优设计方案。本发明根据高速压力机上横梁实际可靠性需求,在可靠性设计中采用均布区间优势度计算可靠性指标值,可便捷地获得符合可靠性要求的高速压力机上横梁设计方案。
Description
技术领域
本发明属于高速压力机设计领域,尤其涉及一种高速压力机上横梁可靠性设计方法。
技术背景
高速压力机的性能直接影响着冲压产品的精度、使用性能和生产效率,而上横梁作为高速压力机的重要组成部分,其刚度和强度对压力机的整体工作性能和加工精度具有直接而重要的影响。为了提高压力机的加工精度和工作效率,降低生产成本和能源消耗,应尽可能地在提高上横梁刚度与强度的同时减轻上横梁的重量,因此,需要对其进行以高刚度轻量化为目标的可靠性设计。
高速压力机实际设计制造中,上横梁的受力情况和材料属性存在一定的不确定性波动,使得其性能具有一定的变异性,若仍然按照确定性问题求解思路进行优化,不仅最终的结果不一定最优,甚至可能存在较大的偏差,无法满足实际精度需求。因此,为了保证设计结果的可靠性,在设计过程必须充分考虑这些客观存在的不确定性,对高速压力机上横梁进行可靠性设计,才能获得真正可靠的最优设计方案。
基于概率可靠性的优化设计是处理不确定性的有效途径之一,在方法和应用上都已有较为充分的研究。概率可靠性分析需要大量的样本数据以得到关于不确定量的精确概率分布信息,然而工程中往往只能得到非常有限的样本数据,并且,概率可靠性对随机参数的分布信息可能是敏感的,即概率模型参数的小误差可引起结构可靠性计算的较大误差。在实际决策中,不确定量的概率分布虽不易准确得到,但其变差的界限则易于确定,因此,可以利用区间数学理论对不确定量的边界进行描述,在不确定量的整个范围内确定结构的可靠性,这实际上得到了更为可靠的结构***。在基于区间的非概率可靠性设计中,一个变量取值为区间数反映了该变量取值的不确定性,因而区间数的大小比较应该能够反映出变量取值的不确定性,一般情况下不应该将区间数的大小关系绝对化,而应该给出一个区间数大于、等于或小于另一个区间数的程度。故在区间描述不确定变量的上横梁可靠性设计模型求解中,需要用一种更为通用、客观、全面的区间数比较方法来计算约束的可靠性指标值。
发明内容
为解决实际工程中高速压力机上横梁高刚度轻量化设计中存在的多重不确定性问题,本发明的目的在于提供了一种高速压力机上横梁可靠性设计方法,根据实际上横梁设计中的可靠性要求选择设计变量,建立以区间描述不确定性因素的高速压力机上横梁可靠性设计模型,基于均布区间优势度计算可靠性设计模型中的可靠性约束值,采用双层嵌套的遗传算法和高精度的Kriging代理模型相结合直接寻找可靠性设计模型最优解。该方法能在存在多重不确定性的情况下获得高可靠性和高精度的压力机上横梁设计方案。
本发明是通过以下技术方案实现的:一种高速压力机上横梁可靠性设计方法,包括以下步骤:
(1)建立以区间描述不确定性因素的高速压力机上横梁可靠性设计模型:根据实际设计需求,确定高速压力机上横梁可靠性设计中的优化目标和约束条件、设计变量及其取值范围、设计中需考虑的不确定性因素及其波动区间,建立如下基于区间变量的上横梁可靠性设计模型:
x=(x1,x2,…,xn)∈Rn;
U=(U1,U2,…,Uq)∈Iq;
其中,x为n维设计向量,U为q维区间向量,f(x,U)为上横梁可靠性设计的目标函数,gi(x,U)为第i个需考虑可靠性的力学性能指标,Bi、Rgi和Rsi分别为其对应的允许变化区间、实际可靠性和给定的可靠性约束值,分别为Bi的下界和上界;hj(x,U)为第j个无需考虑可靠性的力学性能指标,Cj为其对应的允许变化区间,分别为Cj的下界和上界;
(2)采用拉丁超立方采样法(Latin Hypercube Sampling,LHS)在输入变量空间内进行实验设计,获取拟合样本点;
在实验设计中,根据设计向量x和区间向量U的波动范围,在由x和U组成的输入变量空间内采用LHS进行抽样,S个输入变量、N次试验运行的拉丁超立方实验设计表示为取值范围为[0,1]的N×S阶矩阵,获得具有空间均布性和投影均匀性的样本点群,再将其反归一化到x和U组成的输入变量空间中,完成对设计向量x和区间向量U的初始采样;
(3)建立参数化模型,通过协同仿真技术得到样本点对应的目标和约束函数中力学性能指标的响应值;
利用三维建模软件,以设计向量x为独立控制参数,建立高速压力机上横梁参数化三维模型;通过接口技术实现建模软件和有限元分析软件间参数的实时双向传递;通过协同仿真,调用动态更新的参数化三维模型进行有限元分析计算,得到各样本所对应的目标和约束函数中力学性能指标的响应值;
(4)利用完整的输入-输出样本点数据,以上横梁设计变量和不确定性因素为输入参数,以上横梁力学性能指标的响应值为输出参数,建立Kriging响应面模型;
Kriging模型近似表达为一个随机分布函数和一个多项式之和,如下式所示:
y(x)=f(x)β+z(x)
式中,y(x)为未知的Kriging模型,f(x)为已知的关于x的函数,提供了设计空间内的全局近似模拟,β为回归函数待定系数,其值通过已知的响应值估计得到;z(x)为一随机过程,是在全局模拟的基础上创建的期望为0、方差为σ2的局部偏差,其协方差矩阵cov[z(xi),z(xj)]表示为
cov[z(xi),z(xj)]=σ2R[R(xi,xj)]
式中,R为相关矩阵;R(xi,xj)表示任意两个样本点xi,xj的相关函数,选择高斯函数作为相关函数,其表达式为:
其中,n'为样本点的个数,根据无偏条件以及方差最小条件,结合拉格朗日乘子法和极大似然估计法,求得参数β、R和θk的值,进而得到所需要的Kriging模型;
(5)建立基于均布区间优势度的可靠性计算准则;
根据区间数学理论,区间A=[aL,aR]相对于区间B=[bL,bR]的优势度P(A>B)的计算方法:
(a)当aL≥bR时,P(A>B)=1;
(b)当bL≤aL≤bR≤aR时,
(c)当aL≤bL≤bR≤aR时,
(d)当aL≤bL≤aR≤bR时,
(e)当bL≤aL≤aR≤bR时,
(f)当aL≤aR≤bL≤bR时,P(A>B)=0;
利用上述区间优势度计算方法计算高速压力机上横梁各设计约束性能的区间可靠性指标Rgi[gi(x,U)≤Bi];
(6)采用双层嵌套的遗传算法求解上横梁可靠性设计模型,对外层遗传优化当前代种群中的所有个体,利用内层单目标遗传算法和步骤4中建立的Kriging模型计算出其所对应的目标函数和约束函数中力学性能指标区间值的上下界fR(x),fL(x), 并求出其中目标函数和非可靠性约束函数区间值的中点及半径fC(x),fW(x), 再结合步骤5中均布区间优势度的可靠性计算准则得到可靠性约束值Rgi[gi(x,U)≤Bi];其中,上标R、L、C、W分别表示区间上界、区间下界、区间中点和区间半径;
对可靠性约束Rgi[gi(x,U)≤Bi]≥Rsi而言,其约束违反度的计算方式为:
(a)若Rgi[gi(x,U)≤Bi]≥Rsi,则其约束违反度Vi(x)=0;
(b)若Rgi[gi(x,U)≤Bi]<Rsi,则其约束违反度为Vi(x)=Rsi-Rgi[gi(x,U)≤Bi];
对非可靠性约束hj(x,U)≤Cj而言,其约束违反度的计算方式为:
(c)当时,约束违反度Vj(x)=<0,0>;
(d)当时,若则Vj(x)=<0,0>,若则
(e)当时,约束违反度为
由此可以计算出当前代种群所有个体的总约束违反度p为上横梁可靠性设计模型中总的约束个数,则VT(x)=0的解为可行解,否则为不可行解;
将fC(x),fW(x),VT(x)的计算结果由内层优化传递到外层优化各样本点之后,基于区间约束违反度的优于关系准则对外层优化种群中的所有个体进行优劣排序,确定其优劣序位,从而计算获得当前代种群中所有个体的适应度,确定设计向量x1与x2优劣关系的方式为:
(a)若x1为可行解,x2为不可行解,则始终有x1优于x2;
(b)若x1与x2均为可行解,则以目标函数区间值判断两者的相对优劣,当fC(x1)<fC(x2)时,或fC(x1)=fC(x2)且fW(x1)<fW(x2)时,x1优于x2;
(c)若x1与x2均为不可行解,则根据约束违反度来判断其优劣,若VT(x1)<VT(x2),则x1优于x2,否则,x2优于x1;
若外层遗传算法进化代数达到给定最大值或者达到收敛性要求,则终止外层遗传算法进化过程,输出具有最大适应度值的个体作为最优个体,将其所对应的设计向量作为最优设计向量,得到满足可靠性要求的高速压力机上横梁设计方案;否则,生成新一代种群个体,进化代数加1,继续外层遗传进化过程。
本发明具有的有益效果是:
1)采用三维建模软件和有限元软件协同仿真,对于大型复杂的装配体结构能够方便地建立参数化模型,实现参数的双向传递和模型的动态更新。
2)考虑实际工程问题中客观性存在的不确定性因素,采用区间数的形式进行描述,建立更为客观和真实的非概率可靠性优化模型,将均布区间优势度应用于可靠性指标的计算,在不确定变量概率分布未知的情况下能有效获得各样本点的约束可靠性指标值。
3)构建具有良好全局统计性的Kriging模型,避免可靠性设计求解过程中反复调用有限元分析软件计算给定上横梁设计方案在不确定性因素作用下的力学性能指标值,可在提高求解效率的同时保证良好的计算精度和稳健性。
4)采用基于约束违反度和均布区间优势度的双层嵌套遗传算法进行上横梁区间可靠性设计模型的直接求解,避免了现有基于区间的可靠性设计方法需先将其转换为确定性模型再进行求解时不确定性信息的丢失和模型转换过程中参数选择的主观随意性。
附图说明
图1是高速压力机上横梁可靠性设计流程图;
图2是高速压力机上横梁1/4模型三维图;
图3是高速压力机上横梁所受约束和载荷的示意图;
图4是高速压力机上横梁横截面图和各主要尺寸平面图。
具体实施方式
以下结合附图和实施例对本发明作进一步说明。
如图1所示,本发明一种高速压力机上横梁可靠性设计方法,包括以下步骤:
(1)建立以区间描述不确定性因素的高速压力机上横梁可靠性设计模型:
选择某型号高速压力机上横梁为研究对象,如图2所示为其1/4模型,根据实际设计需求,以图4中上横梁的截面尺寸h1、h2、l1、l2和l3为尺寸设计变量,以图3所示外力P1、P2、P3和密度ρ为区间描述的不确定性因素,以变形量为优化设计目标,上横梁的最大许用应力为[60,61]MPa,可靠度要求为0.98,对上横梁进行高刚度轻量化可靠性设计。
定义功能函数δ(x,U)为上横梁的最大等效应力,则其可靠度约束为Rg[δ(x,U)≤[60,61]MPa]≥Rs=0.98,即
Rg[δ(x,U)=δ(h1,h2,l1,l2,l3,P1,P2,P3,ρ)≤[60,61]MPa]≥Rs=0.98
定义函数w(x,ρ)为上横梁的重量,作为设计模型的约束函数,即
w(x,ρ)=w(h1,h2,l1,l2,l3,ρ)≤5500kg
则基于区间的上横梁可靠性设计模型可以表示为:
s.t.Rg[δ(x,U)=δ(h1,h2,l1,l2,l3,P1,P2,P3,ρ)≤[60,61]MPa]≥Rs=0.98;
w(x,ρ)=w(h1,h2,l1,l2,l3,ρ)≤5500kg;
200mm≤h1≤260mm,250mm≤h2≤300mm;
80mm≤l1≤120mm,20mm≤l2≤60mm,330mm≤l3≤400mm;
P1=[P1 L,P1 R]=[2.44×105kN,2.56×105kN];
P2=[P2 L,P2 R]=[4.90×105kN,5.10×105kN];
P3=[P3 L,P3 R]=[7.80×105kN,8.20×105kN];
ρ=[ρL,ρR]=[7290kg·m-2,7310kg·m-2]。
其中,设计向量x=(h1,h2,l1,l2,l3),区间描述的不确定性向量U=(P1,P2,P3,ρ),d(x,U)为表征上横梁最大变形量的目标函数,Rg[δ(x,U)≤[60,61]MPa]为表征上横梁最大等效应力的可靠性约束,可靠度要求为0.98,w(x,ρ)为表征上横梁重量的非可靠性约束。
(2)采用LHS在输入变量空间内获取拟合样本点。
在设计向量x和不确定因素向量U取值范围已确定的情况下,在由x和U组成的输入变量空间内,采用基于最大最小优化准则的拉丁超立方采样法生成取值范围是[0,1]的矩阵,获得具有空间均布性的70个样本点,并将其反归一化到输入向量空间中去,完成对设计向量和不确定性区间向量的初始采样;
(3)建立高速压力机上横梁的参数化三维模型,通过协同仿真技术得到样本点对应的目标和约束函数中各力学性能指标的响应值;
利用三维建模软件,以设计向量x为独立控制参数,建立目标型号高速压力机上横梁参数化模型;在有限元软件中添加不确定因素向量U为二次输入参数;通过接口技术实现三维建模软件和有限元分析软件间参数的双向传递;通过协同仿真,调用三维参数化模型进行有限元分析计算,得到70个样本所对应目标和约束函数中力学性能指标的响应值;
(4)利用样本点数据建立目标和约束函数中各力学性能指标的Kriging模型;
根据包含完整输入-输出信息的样本点数据建立目标和约束函数中力学性能指标的Kriging模型,在拟合过程中,选定高斯函数为基函数,选定二阶回归函数进行拟合,并对拟合结果进行验证,保证其拟合精度和泛化能力满足实际需求,建立的Kriging模型如下:
(5)将可靠性设计模型、完整样本点数据和拟合的Kriging模型代入到双层嵌套的遗传算法中,给定内外层遗传算法的最大进化代数分别为200和400、内外层遗传算法的种群规模分别为100和200、内外层遗传算法的交叉概率分别为0.95和0.90、内外层遗传算法的变异概率分别为0.01和0.05。当外层遗传算法达到终止条件,输出结果方案为h1=229.7mm,h2=264.2mm,l1=119.6mm,l2=55.1mm,l3=337.5mm,将其结果与初始方案进行对比,其目标函数区间中点值且在满足可靠性要求的情况下降低了重量,符合高速压力机上横梁高刚度轻量化可靠性设计要求。
Claims (3)
1.一种高速压力机上横梁可靠性设计方法,其特征在于,该方法包括以下步骤:
(1)建立以区间描述不确定性因素的高速压力机上横梁可靠性设计模型:根据实际设计需求,确定高速压力机上横梁可靠性设计中的优化目标和约束条件、设计变量及其取值范围、设计中需考虑的不确定性因素及其波动区间,建立如下基于区间变量的上横梁可靠性设计模型:
<mrow>
<munder>
<mi>min</mi>
<mi>x</mi>
</munder>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>U</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>U</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>b</mi>
<mi>i</mi>
<mi>L</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>b</mi>
<mi>i</mi>
<mi>R</mi>
</msubsup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>&GreaterEqual;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>U</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>C</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>c</mi>
<mi>j</mi>
<mi>L</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>c</mi>
<mi>j</mi>
<mi>R</mi>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
x=(x1,x2,…,xn)∈Rn;
U=(U1,U2,…,Uq)∈Iq;
其中,x为n维设计向量,U为q维区间向量,f(x,U)为上横梁可靠性设计的目标函数,gi(x,U)为第i个需考虑可靠性的力学性能指标,Bi、Rgi和Rsi分别为其对应的允许变化区间、实际可靠性和给定的可靠性约束值,分别为Bi的下界和上界;hj(x,U)为第j个无需考虑可靠性的力学性能指标,Cj为其对应的允许变化区间,分别为Cj的下界和上界;
(2)采用拉丁超立方采样法LHS在输入变量空间内进行实验设计,获取拟合样本点:在实验设计中,根据设计向量x和区间向量U的波动范围,在由x和U组成的输入变量空间内采用LHS进行抽样,S个输入变量、N次试验运行的拉丁超立方实验设计表示为取值范围为[0,1]的N×S阶矩阵,获得具有空间均布性和投影均匀性的样本点群,再将其反归一化到x和U组成的输入变量空间中,完成对设计向量x和区间向量U的初始采样;
(3)建立参数化模型,通过协同仿真技术得到样本点对应的目标和约束函数中力学性能指标的响应值:利用三维建模软件,以设计向量x为独立控制参数,建立高速压力机上横梁参数化三维模型;通过接口技术实现建模软件和有限元分析软件间参数的实时双向传递;通过协同仿真,调用动态更新的参数化三维模型进行有限元分析计算,得到各样本所对应的目标和约束函数中力学性能指标的响应值;
(4)利用完整的输入-输出样本点数据,以上横梁设计变量和不确定性因素为输入参数,以上横梁力学性能指标的响应值为输出参数,建立Kriging响应面模型;
Kriging模型近似表达为一个随机分布函数和一个多项式之和,如下式所示:
y(x)=f(x)β+z(x)
式中,y(x)为未知的Kriging模型,f(x)为已知的关于x的函数,提供了设计空间内的全局近似模拟,β为回归函数待定系数,其值通过已知的响应值估计得到;z(x)为一随机过程,是在全局模拟的基础上创建的期望为0、方差为σ2的局部偏差,其协方差矩阵cov[z(xi),z(xj)]表示为
cov[z(xi),z(xj)]=σ2R[R(xi,xj)]
式中,R为相关矩阵;R(xi,xj)表示任意两个样本点xi,xj的相关函数,选择高斯函数作为相关函数,其表达式为:
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mi>i</mi>
</msup>
<mo>,</mo>
<msup>
<mi>x</mi>
<mi>j</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&theta;</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>x</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>k</mi>
<mi>j</mi>
</msubsup>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
1
其中,根据无偏条件以及方差最小条件,结合拉格朗日乘子法和极大似然估计法,求得参数β、R和θk的值,进而得到所需要的Kriging模型;
(5)建立基于均布区间优势度的可靠性计算准则;
根据区间数学理论,区间A=[aL,aR]相对于区间B=[bL,bR]的优势度P(A>B)的计算方法:
(a)当aL≥bR时,P(A>B)=1;
(b)当bL≤aL≤bR≤aR时,
(c)当aL≤bL≤bR≤aR时,
(d)当aL≤bL≤aR≤bR时,
(e)当bL≤aL≤aR≤bR时,
(f)当aL≤aR≤bL≤bR时,P(A>B)=0;
利用上述区间优势度计算方法计算高速压力机上横梁各设计约束性能的区间可靠性指标Rgi[gi(x,U)≤Bi];
(6)采用双层嵌套的遗传算法求解上横梁可靠性设计模型,对外层遗传优化当前代种群中的所有个体,利用内层单目标遗传算法和步骤(4)中建立的Kriging模型计算出其所对应的目标函数和约束函数中力学性能指标区间值的上下界fR(x),fL(x), 并求出其中目标函数和非可靠性约束函数区间值的中点及半径fC(x),fW(x), 再结合步骤(5)中均布区间优势度的可靠性计算准则得到可靠性约束值Rgi[gi(x,U)≤Bi];其中,上标R、L、C、W分别表示区间上界、区间下界、区间中点和区间半径;
对可靠性约束Rgi[gi(x,U)≤Bi]≥Rsi而言,其约束违反度的计算方式为:
(a)若Rgi[gi(x,U)≤Bi]≥Rsi,则其约束违反度Vi(x)=0;
(b)若Rgi[gi(x,U)≤Bi]<Rsi,则其约束违反度为Vi(x)=Rsi-Rgi[gi(x,U)≤Bi];
对非可靠性约束hj(x,U)≤Cj而言,其约束违反度的计算方式为:
(c)当时,约束违反度Vj(x)=<0,0>;
(d)当时,若则Vj(x)=<0,0>,若则
(e)当时,约束违反度为
由此可以计算出当前代种群所有个体的总约束违反度p为上横梁可靠性设计模型中总的约束个数,则VT(x)=0的解为可行解,否则为不可行解;
将fC(x),fW(x),VT(x)的计算结果由内层优化传递到外层优化各样本点之后,基于区间约束违反度的优于关系准则对外层优化种群中的所有个体进行优劣排序,确定其优劣序位,从而计算获得当前代种群中所有个体的适应度,确定设计向量x1与x2优劣关系的方式为:
(a)若x1为可行解,x2为不可行解,则始终有x1优于x2;
(b)若x1与x2均为可行解,则以目标函数区间值判断两者的相对优劣,当fC(x1)<fC(x2)时,或fC(x1)=fC(x2)且fW(x1)<fW(x2)时,x1优于x2;
(c)若x1与x2均为不可行解,则根据约束违反度来判断其优劣,若VT(x1)<VT(x2),则x1优于x2,否则,x2优于x1;
若外层遗传算法进化代数达到给定最大值或者达到收敛性要求,则终止外层遗传算法进化过程,输出具有最大适应度值的个体作为最优个体,将其所对应的设计向量作为最优设计向量,得到满足可靠性要求的高速压力机上横梁设计方案;否则,生成新一代种群个体,进化代数加1,继续外层遗传进化过程。
2.根据权利要求1所述的一种高速压力机上横梁可靠性设计方法,其特征在于:所述步骤(5)中,利用均布区间优势度计算各设计向量所对应上横梁结构方案的可靠性指标值,从而计算出可靠性设计模型中各约束条件的违反度,并由此计算出总约束违反度。
3.根据权利要求1所述的一种高速压力机上横梁可靠性设计方法,其特征在于:所述步骤(6)中,利用基于约束违反度和均布区间优势度的双层嵌套遗传算法实现了压力机上横梁可靠性设计模型的直接求解,其中,内层遗传算法利用Kriging模型计算出可靠性设计目标和约束中各上横梁力学性能指标值,外层遗传算法在利用均布区间优势度计算各约束条件的违反度和设计方案总约束违反度的基础上,确定设计方案是否可行,对可行方案根据其目标函数区间响应值进行优劣排序,对不可行方案,根据其约束违反度进行优劣排序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510079808.8A CN104636563B (zh) | 2015-02-14 | 2015-02-14 | 高速压力机上横梁可靠性设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510079808.8A CN104636563B (zh) | 2015-02-14 | 2015-02-14 | 高速压力机上横梁可靠性设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104636563A CN104636563A (zh) | 2015-05-20 |
CN104636563B true CN104636563B (zh) | 2017-11-10 |
Family
ID=53215304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510079808.8A Active CN104636563B (zh) | 2015-02-14 | 2015-02-14 | 高速压力机上横梁可靠性设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104636563B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104915478A (zh) * | 2015-05-24 | 2015-09-16 | 浙江大学 | 基于多参数不确定性分析的产品设计模型等效简化方法 |
CN105930562A (zh) * | 2016-04-13 | 2016-09-07 | 浙江大学 | 一种非概率条件下的结构性能优化设计方法 |
CN105956368B (zh) * | 2016-04-21 | 2018-07-03 | 北京航空航天大学 | 一种复合材料层合结构的非概率动力可靠性评估方法 |
CN106096127B (zh) * | 2016-06-07 | 2019-11-19 | 浙江大学 | 含区间参数不确定性结构的稳健性优化设计方法 |
CN107766601B (zh) * | 2016-08-22 | 2021-03-02 | 南京理工大学 | 制退机节制杆不确定性尺寸优化及优选方法 |
CN107016173B (zh) * | 2017-03-22 | 2020-04-17 | 浙江大学 | 一种高速压力机底座动态特性的可靠性设计方法 |
CN107391804B (zh) * | 2017-06-27 | 2020-12-22 | 广东省建筑设计研究院有限公司 | 基于综合造价法的高层建筑结构抗震性能优化方法 |
CN107563046A (zh) * | 2017-08-30 | 2018-01-09 | 国家电网公司 | 基于大坝失效功能函数模型的失效风险率计算方法及装置 |
CN108427832B (zh) * | 2018-02-11 | 2020-04-14 | 浙江大学 | 一种基于区间约束三维违反矢量的结构性能稳健优化方法 |
CN109992848B (zh) * | 2019-03-14 | 2020-11-13 | 浙江大学 | 一种基于负理想解贴近距离的压力机上横梁稳健优化设计方法 |
CN110847894B (zh) * | 2019-10-30 | 2022-11-04 | 中国石油天然气股份有限公司 | 一种井下节流气井流压的确定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867083A (zh) * | 2012-08-30 | 2013-01-09 | 浙江大学 | 考虑不确定性的压力机滑块机构高刚度轻量化设计方法 |
CN103699720A (zh) * | 2013-12-09 | 2014-04-02 | 浙江大学 | 基于区间约束违反度的高速压力机滑块机构尺寸优化方法 |
CN103942375A (zh) * | 2014-04-11 | 2014-07-23 | 浙江大学 | 基于区间的高速压力机滑块尺寸稳健设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8001515B2 (en) * | 2007-12-21 | 2011-08-16 | National Semiconductor Corporation | Simultaneous optimization of analog design parameters using a cost function of responses |
-
2015
- 2015-02-14 CN CN201510079808.8A patent/CN104636563B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867083A (zh) * | 2012-08-30 | 2013-01-09 | 浙江大学 | 考虑不确定性的压力机滑块机构高刚度轻量化设计方法 |
CN103699720A (zh) * | 2013-12-09 | 2014-04-02 | 浙江大学 | 基于区间约束违反度的高速压力机滑块机构尺寸优化方法 |
CN103942375A (zh) * | 2014-04-11 | 2014-07-23 | 浙江大学 | 基于区间的高速压力机滑块尺寸稳健设计方法 |
Non-Patent Citations (1)
Title |
---|
"基于双层更新Kriging模型的机械结构动态特性稳健优化设计";李小刚等;《机械工程学报》;20141231;第50卷(第3期);165-173 * |
Also Published As
Publication number | Publication date |
---|---|
CN104636563A (zh) | 2015-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104636563B (zh) | 高速压力机上横梁可靠性设计方法 | |
CN102867083B (zh) | 考虑不确定性的压力机滑块机构高刚度轻量化设计方法 | |
CN111353656B (zh) | 一种基于生产计划的钢铁企业氧气负荷预测方法 | |
CN106096127A (zh) | 含区间参数不确定性结构的稳健性优化设计方法 | |
CN103761423B (zh) | 一种基于pso-elm的热轧板材组织-性能预测方法 | |
Izadyar et al. | Intelligent forecasting of residential heating demand for the District Heating System based on the monthly overall natural gas consumption | |
CN103942375A (zh) | 基于区间的高速压力机滑块尺寸稳健设计方法 | |
CN109063234A (zh) | 一种考虑多类型不确定性的高速压力机施力部件可靠性设计方法 | |
CN105930562A (zh) | 一种非概率条件下的结构性能优化设计方法 | |
CN102509240A (zh) | 一种基于多指标多层次的电网投资效益评估方法 | |
CN108664682A (zh) | 一种变压器顶层油温的预测方法及其*** | |
CN103942422B (zh) | 一种基于粒度计算的冶金企业转炉煤气柜位长期预测方法 | |
CN100445901C (zh) | 基于AR(p)模型的过程工业生产过程动态成本控制方法 | |
CN101863088A (zh) | 一种橡胶混炼过程中门尼粘度的预报方法 | |
CN104504280A (zh) | 考虑规划需求的电动汽车充电桩集群管理***通信方式综合评价方法 | |
CN103985000A (zh) | 基于函数型非参数回归的中长期典型日负荷曲线预测方法 | |
CN104408643A (zh) | 一种行为感知偏好获取***及其获取方法 | |
Geng et al. | A new fuzzy process capability estimation method based on kernel function and FAHP | |
CN103853939A (zh) | 一种基于社会经济因素影响的电力***月度负荷的组合预测方法 | |
Chen et al. | Combining fuzzy iteration model with dynamic programming to solve multiobjective multistage decision making problems | |
CN107016173A (zh) | 考虑概率和区间不确定性的高速压力机底座动态特性的可靠性设计方法 | |
CN101739614A (zh) | 层级组合预测通信业务的方法 | |
Singh Sivam et al. | An investigation of hybrid models FEA coupled with AHP-ELECTRE, RSM-GA, and ANN-GA into the process parameter optimization of high-quality deep-drawn cylindrical copper cups | |
CN110880044B (zh) | 一种基于马尔科夫链的负荷预测方法 | |
CN108932562A (zh) | 综合能源***综合效益评价模型建立方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |