本发明属于机械产品多学科设计优化技术领域，尤其涉及一种时变灵敏度分析的风力发电机组多学科稳健设计方法。
背景技术：
复杂机械系统工程设计优化问题通常很复杂，目标函数及约束条件通常是高度非线性的。此外，不确定性因素广泛存在于复杂耦合系统的整个生命周期中，如载荷波动、材料特性、几何尺寸、操作方式和运行环境变化，以及建立数学模型时所作假设带来的不确定性等等。这些客观存在的不确定性因素及其在学科间的传递导致最终设计方案很难达到理想状态，从而无法直接应用到工程实践中。因此，有效处理不确定因素对复杂机械系统的性能及其演变过程的影响，进而全面提高复杂机械系统的设计水平，已成为影响机械系统质量、性能、时间、成本、价格等的关键因素。许多研究机构和学者对于时不变不确定性因素条件下的mdo研究已有一些成果。
然而，大量的事实表明导致产品最终失效的主要原因之一是复杂机械系统的退化失效。这种失效主要是由于时变不确定性因素造成的，例如机械零部件磨损、强度降低、电子元器件性能衰退、环境腐蚀氧化、材料老化等。目前，对于时变不确定因素的处理上基本依然在单学科分析上。对于时变因素的单学科分析已经有了一定的研究成果。但是，在多学科设计优化中，时变不确定性因素的分析十分困难。首先，时变不确定因素表现形式多样化、相互关联；其次，在mdo中各个子系统之间存在着层次与非层次的混合耦合关系，系统各个层次模型中的时变不确定因素经过传递对系统的最终输出有着不同程度的影响，加之mdo本身就是针对大型系统的一种协调优化方法体系，这就造成了时变不确定条件下mdo的困难。
工程实际中各种不确定性因素的变化将会导致机械系统的性能发生波动，严重时可能导致失效。消除不确定因素是相当困难，减低不确定因素的影响是首选和相对容易的方法，即系统可靠性对这些不确定因素变化不敏感。已有很多学者通过稳健设计优化方法使复杂系统在多种不确定因素作用下，仍然能正常发挥其功能，并且对各种不确定因素变化不敏感。通过对时变灵敏度分析以及多学科稳健设计优化相结合，进一步提出多学科多目标时变灵敏度稳健设计优化方法。这就大大减少系统分析代价和降低计算时间。尽量降低时变不确定性因素的作用，使状态方程的可靠度随时变不确定性因素变化而变化的波动最小。
风力发电机组的设计是一个典型的时变多学科问题，需要综合考虑时变风电场特性及风电机组的结构、时变载荷、气动、控制和电磁等多个学科之间的耦合特性。
技术实现要素：
本发明的发明目的是：为了解决现有技术因忽略多学科系统中时变不确定因素的客观存在而导致的设计结果不理想等问题，本发明提出了一种时变灵敏度分析的风力发电机组多学科稳健设计方法。
本发明的技术方案是：一种时变灵敏度分析的风力发电机组多学科稳健设计方法，包括以下步骤：
a、综合考虑风力发电机组的整体性能，建立风力发电机组的mdo模型；
b、考虑时变不确定因素的影响，并结合修正公式给出时变可靠性灵敏度具体计算方法；
c、分析时变不确定因素传播，利用简化的iup给出了时变不确定因素在mdo中各学科传递的量化方法；
d、结合稳健设计，搭建时变不确定条件下mrdo优化求解框架，进而降低时变不确定因素对系统性能的影响。
进一步地，所述步骤a中建立风力发电机组mdo模型，具体为将风力发电机组设计分解为n个子学科，第i个学科级优化目标函数的数学模型表示为：
min:fi(xi)
s.t.gi(xi,y1,…,yi,…,yn)≤0
yi＝tdiscipline-i(xi,y1,…,yi-1,yi+1…,yn)
其中，xi为第i学科级设计变量，fi(·)为第i个学科的目标函数，yi为第i学科级耦合状态变量，gi(·)为系统级第i个学科的不等式约束条件，tdispline-i(·)为第i个学科的分析函数；
系统级优化目标函数的数学模型表示为：
s.t.gi(xi,y1,…,yi,…,yn)≤0
其中，x为系统级级设计变量，n表示子学科数，f为系统级目标函数，ωi为加权因子。
进一步地，所述步骤b中假设随机变量x＝[x1,x2,…,xn]t，那么可靠度对随机变量均值的灵敏度为
式中：
其中
上式对状态方程函数为线性时或者非线性不强时适用，当状态函数非线性程度较高时则不适用，上式(1)假定状态函数的方差和随机变量的均值是相互独立的，若状态函数非线性程度较高就不成立。当设计变量为正态分布时通过改进修正得到计算公式如下：
式中：
根据edgeworth级数方法和公式(2)可以得到任意分布的基本随机向量x的时变灵敏度为：
其中
上式中
hj(y)为j阶hermite多项式，递推关系如下：
计算实验表明采用经验修正公式比edgeworth级数估计可靠度所获得结果更接近montecarlo结果，edgeworth级数经验修正公式：
对可靠度指标β(t)的灵敏度通过经验公式的导数计算：
用上式替换(3)中的
进一步地，所述步骤c，
在简化的基于iup的方法中，无需对gse方程和局部偏导数进行求解计算，而是将耦合状态变量的变差值作为辅助设计变量处理。考虑设计变量误差和学科分析误差，模型表述如下：
s.t.gi(xd,x,y,p(t))+δgi(xd,x,y,p(t))≤0
其中xd为确定性设计变量，x为时变设计变量，y为耦合状态变量，p(t)为时变设计参数，△xi表示第i个设计变量的变差，其中△yj表示第j个辅助变量的变差，其中△f(·)为f(·)的变差。t(·)为学科分析模型，sdesiagn为设计变量个数，sauxiliary为耦合状态变量个数，ω1，ω2为权重因子，f*和△f*分别是[ω1，ω2]＝[0,1]及[ω1，ω2]＝[1,0]时fo的值。△ε为分析模型误差。(·)u和(·)l分别为设计向量上下限，(·)ru和(·)rl分别为稳健设计向量上下限。
进一步地，所述步骤d，
通过对时变灵敏度分析以及多学科稳健设计优化相结合，进一步提出多学科时变灵敏度稳健设计方法。其模型如下：
s.t.r(t)-[r]≥0
gi(xd,x,y,p(t))+δgi(xd,x,y,p(t))≤0
y＝tdispline(xd,x,y,p(t))
其中fr(·)为时变灵敏度分析产生的目标函数，[r]为给定的可靠度要求。
本发明的有益效果是：本发明的时变灵敏度分析的风力发电机组多学科稳健设计方法，解决了现有技术因忽略多学科系统中时变不确定因素的客观存在而导致的设计结果不理想等问题，从而更能满足产品的实际设计需求。
附图说明
图1是本发明的基于时变灵敏度分析的风力发电机组多学科稳健设计方法流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
如图1所示，为本发明的基于时变灵敏度分析的风力发电机组多学科稳健设计方法流程示意图。一种基于时变灵敏度分析的风力发电机组多学科稳健设计方法，包括以下步骤：
a、综合考虑风力发电机组的整体性能，建立风力发电机组的mdo模型；
b、考虑时变不确定因素的影响，并结合修正公式给出时变可靠性灵敏度具体计算方法；
c、分析时变不确定因素传播，利用简化的iup给出了时变不确定因素在mdo中各学科传递的量化方法；
d、结合稳健设计，搭建时变不确定条件下mrdo优化求解框架，进而降低时变不确定因素对系统性能的影响。
在步骤a中建立风力发电机组mdo模型，参考有关文献将风力发电机组设计分解叶片、齿轮箱、塔架3个子学科，其中子学科1优化目标是风力机叶片装机成本最低。风力机叶片装机成本与叶片长度、材料成本、制造水平、劳动力成本等多因素有关，根据国际可再生能源室的研究然后简化得到叶片装机总成本(单位：人民币)
min:f1＝1.4d3+13.3d2.5-26269
d为风轮直径。
子学科2以齿轮箱箱体体积和齿轮体积加权和为优化目标，
min:f2＝w1v1+w2v2
其中v1为齿轮箱箱体体积，v2为齿轮体积；w1,w2为加权系数；σh,σf分别为接触疲劳强度和弯曲疲劳强度；σh,σf分别为许用接触疲劳强度和许用弯曲疲劳强度；gba为第a个子学科的第b个不等式约束；hba为第a个子学科的第b个等式约束；l,s,h为齿轮箱长宽高，[l],[s],[h]为长宽高的极限最大值；z为行星轮系内齿轮齿数，zs为太阳轮齿数，zp为行星轮齿数；i1,i2,i3为齿轮箱1-3级传动比；n和n'为输入和输出转速。
子学科3为风力机塔架优化通常将塔架简化成三段等高锥形薄壁筒体。假设塔架从上到下三段筒体的薄壁壁厚分别为δ1,δ2,δ3，外径分别为d1,d2,d3,d4,塔架总高为h。则以塔架装机
总成本为优化目标函数，则可得到优化模型如下(假设塔架的成本约为1万元/吨)：
其中sv为工作风速下的弯曲应力，sev为极端风速下的弯曲应力，sσ为许用弯曲应力。
这时系统级优化目标函数的数学模型可表示为：
其中，x为系统级级设计变量，n表示子学科数，f为系统级目标函数，ωi为加权因子。
在步骤b中确定性设计变量xd＝[l,s,h,z,zs,zp,h]＝[xd1,xd2,…,xd7]t，其中随机变量x＝[d,δ1,δ2,δ3,d1,d2,d3,d4]＝[x1,x2,…,x8]t，时变设计参数p＝[σh,σf,sv,sev]，该例中
就可计算出时变可靠性灵敏度。
进一步地，所述步骤c，考虑设计变量误差和学科分析误差，模型表述如下：
其中△xi表示第i个设计变量的变差，(·)u和(·)l分别为设计向量上下限，(·)ru和(·)rl分别为稳健设计向量上下限。
进一步地，所述步骤d，通过对时变灵敏度分析以及多学科稳健设计优化相结合，进一步提出多学科时变灵敏度稳健设计方法。其模型如下：
s.t.r(t)-[r]≥0
其中fr(·)为时变灵敏度分析产生的目标函数，[r]为给定的可靠度要求。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。