24.12.2014 Views

Earthquake Engineering Research - HKU Libraries - The University ...

Earthquake Engineering Research - HKU Libraries - The University ...

Earthquake Engineering Research - HKU Libraries - The University ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

133<br />

1st Generation<br />

dGenerate individuals randomly !<br />

generations<br />

[ Evaluation of fitness )<br />

I Preserve the elite solutions j<br />

Selection<br />

O<br />

D<br />

parent<br />

child<br />

.0 a<br />

..-Parent-]<br />

/<br />

0<br />

•Parent .2<br />

F<br />

Parent-3<br />

\.<br />

0"<br />

gQ-^arent _2<br />

a-'<br />

%r'DD<br />

0-Sn<br />

.©"'<br />

..•Parent- !<br />

Parent-3<br />

Q<br />

/ ••.<br />

f*B\<br />

/%£ )D<br />

••' Q<br />

•'* n<br />

(£"""<br />

Parent- 1<br />

••*"'"<br />

Type-1 Type-2 Type-3<br />

Fzg. 4 Outline of Crossover in Real Coded GA<br />

..-••"" Parent-<br />

Crossover and Mutation J modified Newton methods (DFP, BFGS, etc.) can not be<br />

applicable. In this study, we adopt GA (Genetic algorithms)<br />

Optimum solution<br />

which attracts attention as practical global optimization<br />

methods.<br />

Fig.3 Flowchart of Genetic Algorithm<br />

Traditionally, binary codes have been used to represent<br />

the design variables including real number parameters, and<br />

the crossover operation is performed by the exchange of a<br />

part of the binary bit. However, some researchers are using<br />

the GA with real coded crossover operators (Tsutsui et al,1997) and reported the advantage of their<br />

methods. Among them, the unimodal normal distribution crossover (UNDX) is reported to show good<br />

performance in optimizing problems (Kita, et al, 1999).<br />

Basically, the outline of real coded GA which is shown in Fig.3 is same as that of binary coded GA,<br />

except crossover operator and mutation. <strong>The</strong>re are several models proposed for real coded crossover<br />

operator with Gaussian noise. <strong>The</strong> concepts of the crossover operator which are used in this paper are<br />

illustrated in Fig.4 (Yoshida, 2000). In the case of Type-1, the children are derived by following<br />

equation.<br />

g_ A+P 2 Pi-P* here, a * N(0,l) (2)<br />

V<br />

2 2<br />

c, p l and p 2 are vectors in design parameter space, c is a child made from parents p l and p 2 . <strong>The</strong><br />

parents are randomly selected in proportion to their fitness, which is the inverse of safety factor in this<br />

slope stability study. <strong>The</strong> children are generated in one-dimensional subspace of the parameter space in<br />

Type-1. In Type-2, n dimensional perturbation is given to Type-1 children, consequently Type-2<br />

children are generated in n dimensional space, when there are n design variables. <strong>The</strong> magnitude of the<br />

perturbation is defined based on the distance from Parent-3 to the line with Parent-1 and 2. In Type-3,<br />

the children are generated in 2-dimensional subspace made by three parents as shown in Fig.4. In all<br />

types, the generated children are in the neighbourhood of the gravity center of their parents.<br />

Perturbation v is added to each design variable as mutation. In this study, Cauchy distribution is<br />

used.<br />

NUMERICAL EXAMPLES<br />

Models and Characteristics of Dynamic Response<br />

Two types of slope model are used for the numerical calculations. FEM model of them are shown<br />

in Fig.5. <strong>The</strong> height and width of the slope are indicated in Fig.6. <strong>The</strong> number of node and element of<br />

Type-1 model is 8848 and 7557, while those of Type-2 are 11336 and 9883. <strong>The</strong> sides and base of the

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

Saved successfully!

Ooh no, something went wrong!