构件损伤状态判定方法.pdfVIP

附录G构件损伤状态判定方法

G.0.1构件损伤状态应采用蒙特卡洛模拟方法确定，蒙特卡洛模拟流程见图

F.0.2。

模拟开始

工程需求参数

矩阵扩充

i=0

i=i+1

取第i组工程需求参数

生成0~1之间的随机数

根据随机数和构件易损

性曲线确定损伤状态

修复费用修复时间

否

i1000

是

模拟结束

图G.0.1蒙特卡洛模拟流程

G.0.2每次时程分析可得到一组工程需求参数，将各次时程分析的工程需求参

数组装为矩阵，每列表示一个特定的工程需求参数的取值，每行表示一次时程

分析的结果。

38

G.0.3扩充后的工程需求参数与分析得到的工程需求参数应具有相同的联合分

布，扩充后的工程需求参数矩阵与时程分析得到的工程需求参数矩阵应具有相

同的均值与方差。

G.0.4工程需求参数的分布取联合对数正态分布。对工程需求参数矩阵的值应

取对数，新矩阵中参数分布应满足联合正态分布。

G.0.5工程需求参数的对数分布可由下式确定

Z=LU+M（G.0.5）

Y

式中

T

Z——Z=[Z,Z…,Z]，Z为第i个工程需求参数的取对数变量；

12ni

T

U——U=[U,U…,U]，U为第i个独立正态分布的变量；

12mi

L——∑矩阵Cholesky分解得到的下三角矩阵，∑为取对数后工程需求

YYYY

参数矩阵Y的协方差矩阵；

M——矩阵Y的均值矩阵。

Y

G.0.6当矩阵∑满秩时，G.0.5中m应等于工程需求参数的个数n；当矩阵

YY

∑不满秩时，m应等于∑的秩。L为n×m的矩阵，当i小于j时，l取0。

YYYYij

G.0.7独立正态分布变量的取值可由计算机的伪随机数生成。

G.0.8工程需求参数的分布由其取对数的分布取指数确定。

G.0.9利用式（G.0.5）计算得到一组工程需求参数的对数取值，取指数后得

到工程需求参数的取值，多次利用式（G.0.5），将计算结果组装为矩阵，即可

得到扩充后的工程需求参数分布的取值。

G.0.10构件各