《1 引言》

1 引言

在工程结构的可靠性分析中, 一次二阶矩法 (包括中心点法、JC法、映射变换法和实用分析法) 以及二次二阶矩法得到了广泛应用 [1,2]。一次 (二次) 二阶矩法主要利用随机变量的均值和方差信息以及分布概型来计算结构的失效概率, 属于低阶矩法。使用低阶矩法的前提是:采用的随机变量的分布概型是正确的, 随机变量的有关统计参数是准确的。然而, 在实际工程中, 获得的样本数量往往很小, 在一定的置信度下, 2个或者几个不同的分布概型都有可能被接受, 随机变量的统计参数可能与真实值差别较大。此外, 目前采用的随机变量分布概型几乎都是理想的数学模型, 现实中的问题由于受多种因素影响可能并不严格符合理想的数学模型, 而是服从某种未知的分布或几种分布组合而成的分布 [3,4]。为了全面反映和利用随机变量的统计信息, 利用随机变量的高阶矩对工程结构进行可靠性分析是一条新思路 [5]。四阶矩分析法用到了随机变量的前四阶矩信息, 对结构功能函数非线性程度不高的结构可靠性分析问题, 比较适用 [6]。以{1, x, \:, xn}为基的多项式逼近法在多项式的阶数不是很高的情况下, 可以较好地逼近功能函数的概率密度函数 [4]。当阶数较高时容易产生龙格现象, 计算结构失效概率时误差较大。因此, 在实际工程中, 如何利用已有的统计信息, 对结构的失效做出符合实际情况的概率测度, 在结构可靠性研究中很有必要。作者基于数值逼近原理, 充分考虑随机变量的高阶矩信息, 利用切比雪夫正交多项式来逼近计算功能函数的概率密度函数, 进而可以利用数值积分计算结构的可靠度, 并给出了计算实例。

《2 随机变量概率密度函数的正交多项式逼近》

2 随机变量概率密度函数的正交多项式逼近

根据工程结构可靠度的定义 [1,2], 结构功能函数可以用Z= (X1, X2, \:, Xn) 表示, X1, X2, \:, Xn表示结构的基本随机变量。假定结构的抗力随机变量为R和荷载效应随机变量为S是两个独立的连续随机变量, 其相应的概率密度函数为fR (r) , fS (s) , 则结构的失效概率pf为:

pf=Ρ(Ζ<0)=r<sfR(r)fS(s)drds(1)

在实际问题中, 有时功能函数是非线性的, 比较复杂, 对fR (r) 或fS (s) 积分难以得到显式表达式, 因此直接通过数值积分计算结构的失效概率, 在实际工程中难以实现。正交多项式逼近法是采用正交函数族做基, 利用功能函数的高阶矩信息, 通过计算得到功能函数概率密度函数的近似表达式, 进而可以利用数值积分计算结构的失效概率。其具体推导过程如下:

f (x) , g (x) 是闭区间[a, b]上关于x的连续函数, η (x) 是[a, b]上的权函数, 积分

(f(x),g(x))=abη(x)f(x)g(x)dx(2)

称为函数f (x) 与g (x) 在[a, b]上的内积。满足内积定义的函数空间C[a, b]称为内积空间。并且在C[a, b]上首项系数an≠0的n次多项式满足

abω(x)φi(x)φj(x)dx={0,ijAk>0,i=j(i,j=0,1,\:)(3)

称{φn (x) }是C[a, b]上带权ω (x) 的正交多项式 [7]。设f (x) ∈C[a, b], φ0 (x) , φ1 (x) , …, φn (x) 为空间C[a, b]中的一个线性无关系, a0, a1, \:, an为任意实数, 则S(x)=i=0nω(x)aiφi(x)C[a, b]中是的一个子集。根据函数最佳平方逼近原理 [7], 存在最佳平方逼近函数S* (x) , 使得

f(x)-S*22=infSf(x)-S*22=infSCabη(x)[S(x)-f(x)]2dx(4)

即求Ι=abη(x)[i=0nω(x)aiφi(x)-f(x)]2dx的最小值。

由于{φi (x) }为C[a, b]上的一个无关系, 利用多元函数求极值的方法, 可得方程:

Ιαk=2abη(x)[i=0nω(x)aiφi(x)-f(x)]φk(x)dx=0(k=0,1,2,\:,)

ab[i=0nω(x)aiφi(x)]φk(x)dx=abf(x)φk(x)dx(k=0,1,2,\:,)(5)

系数ai满足方程组

ω(φ0,φ0)a0+ω(φ0,φ1)a1+\:+ω(φ0,φn)an=(φ0,f)ω(φ1,φ0)a0+ω(φ1,φ1)a1+\:+ω(φ1,φn)an=(φ1,f)ω(φn,φ0)a0+ω(φn,φ1)a1+\:+ω(φn,φn)an=(φ0,f)}(6)

这是关于a0, a1, \:, an的线性方程组。因为φ0 (x) , φ1 (x) , \:, φn (x) 线性无关, 故系数行列式G (φ0, φ1, \:, φn) ≠0, 方程组 (6) 有唯一解。

根据函数族{φn (x) }的正交性, 方程组 (6) 左边系数矩阵除对角线上元素外, 其余都为0, 所以得到方程组 (6) 的解, 即:

ak=(φk,f)ω(φk,φk)(7)

f (x) 的n次最小二次方逼近多项式为

Sk=k=0n(φk,f)ω(φk,φk)φk(x)=k=0nakφk(x)(8)

考虑f (x) ∈C[-1, 1], 利用切比雪夫正交多项式{Tk (x) }逼近f (x) , 其中权函数ω(x)=11-x2, 并且满足如下关系:

Τ0(x)=1,Τ1(x)=x,Τn+1(x)=2xΤn(x)-Τn-1(x).-11Τm(x)Τn(x)dx1-x2={0,mnπ2,m=n0π,m=n=0(9)

式 (6) 的右边可写为

(φk,f)=(Τk,f)=-11Τk(x)f(x)dx(10)

其中Tk (x) 为一般多项式, 设其形式为Bk0+Bk1x+…+Bkkxk, 代入式 (10) , 有

-11Τk(x)f(x)dx=-11(Bk0+Bk1x+\:+Bkkxk)f(x)dx=Bk0-11f(x)dx+Bk1-11xf(x)dx++Bkk-11xkf(x)dx

μ=-11xf(x)dxμk=-11xkf(x)dx}(11)

μ, μk分别为随机变量x的均值和k阶原点矩。

把式 (11) 代入式 (10) , 则有

(φk,f)=(Τk,f)=-11Τk(x)f(x)dx=Bk0μ0+Bk1μ1+\:+Bkkμk=i=0kBkiμi(12)

式 (6) 的左边可写为:

ω(φk,φk)=ω(Τk,Τk)=-11Τk(x)Τk(x)1-x2dx={π,k=0π2,k0(13)

所以

α0=(Τ0,f)ω(Τ0,Τ0)=B00μ0π,a1=(Τ1,f)ω(Τ1,Τ1)=B10μ0+B11μ1π/2,\:,ak=(Τk,f)ω(Τk,Τk)=i=0kBkiμiπ/2,\:(14)

得到概率密度函数f (x) 的多项式逼近表达式fn (x) , 即f(x)fn(x)=k=0nω(x)akΤk(x)

一般情况下, f (x) ∈C[a, b]。因为切比雪夫多项式定义在区间[-1, 1]上, 所以求解在[a, b]上的逼近多项式, 需做区间变换

x=b-a2t+b+a2(-1t1)

xt之间为线性关系, 当x在区间[a, b]上变化时, 对应的t在区间[-1, 1]上变化, 则

Τ0(t)=Τ0(2x-(b+a)b-a)=1Τ1(t)=t=2x-(b+a)b-a,Τ2(t)=3t2-12=3(2x-(b+a)b-a)2-12

最后, 得到f (x) 的多项式逼近表达式fn (x) , 即

f(x)fn(x)=k=0nω(t)akΤk(t)=k=0nakω(2x-(b+a)b-a)Τk(2x-(b+a)b-a)(15)

在结构可靠性分析中, 以切比雪夫正交多项式做基, 利用功能函数的高阶矩 (即f (x) 的前k阶原点矩) , 求解方程组 (6) , 得到各项系数ak, 即可求得功能函数概率密度函数的逼近表达式, 对结构的可靠性进行分析。

《3 结构可靠性分析的高阶矩法》

3 结构可靠性分析的高阶矩法

基于数值逼近的结构可靠性高阶矩法是以切比雪夫正交多项式{Tk (x) }做基, 利用功能函数f (x) 的高阶原点矩, 求解线性方程组, 得到功能函数概率密度函数的逼近表达式fn (x) , 然后估算结构的失效概率。具体计算步骤如下:

1) 分析工程结构的实际情况, 建立功能函数Z=g (X) = (X1, X2, …, Xn) 。

2) 根据结构随机变量的统计特征或分布概型, 对g (X) 进行抽样。

3) 考虑正交多项式性质, 区间做适当变换, 计算功能函数g (X) 的各阶样本矩。

4) 以切比雪夫正交多项式{Tk (x) }做基, 利用得到的各阶样本矩, 求解功能函数概率密度函数f (x) 的近似多项式fn (x) 。

5) 计算结构的失效概率, 对工程结构进行可靠性分析。

在计算失效概率时 [1,2], 引入标准化随机变量z, 即 (μz=0, σz=1) :

z=(Ζ-μz)/σz

f (x) 的近似多项式fn (x) 代入式 (1) , 得

pf=Ρ(Ζ<0)=Ρ(z<-μz/σz)=--μzσzfn(x)(15)

在现实情况中, 随机变量不可能无限制地趋近于-∞, 因此, 对上式积分时, 积分下界可根据工程结构的实际情况估计一个界值, 使误差在可接受的范围之内。

《4 结构可靠性高阶矩法的检验》

4 结构可靠性高阶矩法的检验

《4.1 经典分布函数的数值检验》

4.1 经典分布函数的数值检验

以工程结构随机变量最常服从的正态分布和对数正态分布为例, 检验高阶矩法拟合各种经典分布曲线的优良性。首先根据经典概率密度函数的解析表达式, 画出其理论分布曲线, 然后用积分的方法计算出理论分布的各阶原点矩, 采用切比雪夫正交多项式为基, 做拟合曲线并进行比较, 观察后者拟合性能的优劣。

图1, 图2分别是用切比雪夫正交多项式为基对正态分布、对数正态分布 (具体表达式见表1) 的拟合曲线。从图中可以看出, 拟合曲线与两种经典分布已经基本重合。通过数值检验说明, 只要已知随机变量的前n阶原点矩, 采用切比雪夫正交多项式为基生成的概率密度函数具有良好的拟合性能。

《图1》

图1 正态密度函数的数值检验

图1 正态密度函数的数值检验  

Fig.1 Numerical test of normal density function

《图2》

图2 对数正态密度函数的数值检验

图2 对数正态密度函数的数值检验  

Fig.2 Numerical test of log-normal density function

表1 两种经典分布曲线及其表达式

Table 1 Two classical distributions and their expressions

 

《表1》


名 称
正态分布 对数正态分布

表达式
y=12π0.3e-x22(0.3)2
y=10.12πxe{-12[lnx-10.1]2}

置信区间
[-1, 1] [2, 4]

置信度β/%
99.91 99.89

 

 

《4.2 工程结构失效概率的计算》

4.2 工程结构失效概率的计算

某结构构件正截面强度的功能函数为Z= (R, S) =R-S, 其中抗力R和荷载效应S均服从正态分布, 平均值μR=100 kN·m, μS=50 kN·m, 变异系数δR=0.12, δS=0.15。试求结构构件的失效概率。

因为抗力R和荷载效应S均服从正态分布, 所以σR=μR·δR=12, σS=μS·δS=7.5。功能函数Z也服从正态分布, 并且μz=50 kN·m, σz2=200.25。-μz/σz=-3.533 3, 因此在区间[-5, -3]上进行概率密度函数的逼近, 同时将功能函数随机变量Z标准化为z。根据正交多项式逼近法, 利用样本数据的前n阶矩, 计算得到在区间[-5, -3]上标准化功能函数z的概率密度函数各阶近似表达式fn (t) 。近似表达式的各项系数见表2, 把t=x+4代入, 即可求得fn (x) 的具体表达式。

表2 切比雪夫多项式生成函数各项系数及计算结果

Table 2 The coefficients of generated probability density functions and result of pf

 

《表2》


各项系数
f7 (x) f8 (x) f9 (x) f10 (x)

a0
0.000 429 6 0.000 429 6 0.000 429 6 0.000 429 6

a1
0.000 616 3 0.000 616 3 0.000 616 3 0.000 616 3

a2
0.000 144 9 0.000 144 9 0.000 144 9 0.000 144 9

a3
-0.000 167 6 -0.000 167 6 -0.000 167 6 -0.000 167 6

a4
-0.000 236 3 -0.000 236 3 -0.000 236 3 -0.000 236 3

a5
-0.000 184 4 -0.000 184 4 -0.000 184 4 -0.000 184 4

a6
-0.000 120 5 -0.000 120 5 -0.000 120 5 -0.000 120 5

a7
-7.812×10-5 -7.812×10-5 -7.812×10-5 -7.812×10-5

a8
  -5.428×10-5 -5.428×10-5 -5.428×10-5

a9
    -4.044×10-5 -4.044×10-5

a10
      -3.166×10-5

pf
2.032×10-5 2.053×10-5 2.058×10-5 2.058×10-5

 

 

根据结构可靠性分析失效概率的定义, 利用样本数据的前n阶矩计算得到的fn (x) , 分别进行数值积分求解失效概率pf, 计算结果见表2第13行。可以看出, 当n≥7时, 计算结果已经基本趋于稳定。利用中心点法、验算点法以及蒙特卡罗模拟 [1]的计算结果分别为2.078×10-4, 2.05×10-4和2.05×10-4, 从中可以看出该方法的正确性和实用性。

《5 结论》

5 结论

1) 提出了工程结构可靠性分析的高阶矩法。该方法基于数值逼近原理, 以切比雪夫正交多项式做基, 利用功能函数的高阶矩信息, 通过求解线性方程组, 得到功能函数的逼近概率密度函数表达式fn (x) , 然后利用工程结构可靠性分析的一般式计算结构的失效概率, 对工程结构的可靠性进行分析。

2) 高阶矩分析法能充分利用结构随机变量的统计信息, 避开了复杂的数值寻优计算, 很大程度地提高了计算效率, 并且可以按任意指定的高阶矩信息, 达到满意的计算精度。通过计算, 可以依次得到多项式各项系数, 计算难度和复杂性不随矩阶数的增加而增加, 数值计算结果稳定。根据结构构件失效概率的计算, 表明可靠性分析的高阶矩数值逼近法是正确而且实用的。