28.05.2014 Views

r - The Hong Kong Polytechnic University

r - The Hong Kong Polytechnic University

r - The Hong Kong Polytechnic University

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Sigmoid function<br />

Another kernel, that might seem appealing, is the hyperbolic tangent kernel<br />

K xc , = tanh v( xc T<br />

) + c<br />

(12)<br />

( )<br />

i<br />

(<br />

i<br />

)<br />

Supposing y = 1+ exp( − 2 ( v( xc ) ))<br />

+ c<br />

i<br />

i<br />

z<br />

m<br />

⎛ 2<br />

( − α + α )<br />

i i<br />

1<br />

i=<br />

1<br />

i<br />

⎞<br />

, then<br />

= ∑ ⎜ − ⎟+<br />

b<br />

y<br />

(13)<br />

⎝<br />

∂y<br />

i<br />

∂x<br />

j<br />

⎠<br />

( ( ))<br />

( j ) T<br />

2 exp 2 ( ) 2 ( j<br />

vc v c vc )<br />

( y 1)<br />

i i i i<br />

=− − xc + =− −<br />

(14)<br />

m m m<br />

∂z<br />

∂z<br />

∂y<br />

1 1 1<br />

= = − − − −<br />

∂x ∂y<br />

∂x y<br />

y y<br />

i<br />

*<br />

( )( ( j ) *<br />

) ( j<br />

∑ 2∑ − α + α 2 vc ( y 1) = 4 vc<br />

)<br />

2<br />

( α α )( )<br />

i i i i ∑ − +<br />

i<br />

i i<br />

2 (15)<br />

j i= 1 i j i= 1 i i=<br />

1<br />

i i<br />

2<br />

m<br />

m<br />

∂ z ∂ ⎛<br />

( j) *<br />

1 1 ⎞<br />

( j)<br />

*<br />

∂ 1 1<br />

= ⎜4 ∑vc<br />

( − α + α )( − ) 4<br />

2<br />

( ) ( )<br />

i<br />

i i ⎟ = ∑vc<br />

− α + α −<br />

i<br />

i i<br />

2<br />

∂x ∂x ∂x j k k ⎝ i= 1 y y<br />

i i ⎠ i=<br />

1<br />

∂x y y<br />

k i i<br />

2 1 ∂y<br />

2 1<br />

= − = − − + − −<br />

SVM-based reliability analysis method with UDM<br />

m<br />

( j) 2 ( j) ( k) *<br />

4 ∑vc ( ) 8 v c c ( )( y 1)<br />

i 3 2 i i i i 3 2 i<br />

i= 1 y y ∂x i i k i=<br />

1<br />

y y<br />

i i<br />

m<br />

*<br />

i<br />

( − α + α<br />

i i<br />

) ∑ ( α α )<br />

Once the first-order and second-order derivatives of the structural response have been calculated, FORM or<br />

SORM can be used to get the failure probability.<br />

<strong>The</strong> calculation steps of SVM-based reliability method can be described as follows:<br />

1) Select the random variables, specify their probabilistic characters, and define the limit state function g ( x ) = 0;<br />

2) Generate the samples from the distribution functions of the basic variables with UDM and compute the<br />

corresponding value of the performance function g(x) and divide the samples into two groups, namely the<br />

training set and the validation set;<br />

3) Establish the SVM model of the performance function;<br />

4) Calculate the support vectors according to the training set and test SVM according to the validation set. If the<br />

tolerance defined before is met, the next step will be executed. Otherwise, repeat steps 3 through 4 until the<br />

error bellows the tolerance;<br />

5) Choose the reliability method, such as FORM, SORM or MCSM and compute failure probability according<br />

to the functional approximation by SVM.<br />

6) Calculate the reliability index β and failure probability P<br />

f<br />

.<br />

APPLICATIONS<br />

Example 1<br />

<strong>The</strong> limit state function is defined as in [17]<br />

g ( x ) = exp(0.2× x + 1.4) −x<br />

1 2<br />

(17)<br />

where x and x are assumed to be independent and have a standard normal distribution with zero mean and<br />

1 2<br />

unit standard deviation.<br />

2<br />

Select the uniform design table U<br />

25 ( 25 ). Take 20 samples which are uniformly distributed within[ − 3, 3]<br />

and compose the experimental designs with 2 factors and 20 levels. <strong>The</strong> results are shown in table 1. It can be<br />

seen that the results obtained by UDM-SVM RSM are in good agreement with reference(Kim et al. 1997; Zheng<br />

et al. 2000; Elhewy et al. 2006). It suggests that even in the case of small sample sizes, SVM is also more<br />

accurate fitting of the limit state surface. Moreover, for the same number of training samples taken (eg 25),<br />

UDM selected samples obtained better results than the orthogonal design.<br />

(16)<br />

-402-

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!