• 대한토목학회논문집
• /
• 제26권2A호
• /
• pp.383-390
• /
• 2006

본 연구에서는 급수전개를 이용한 추계론적 유한요소해석법의 개선을 위한 등가몬테카를로 추계장함수를 제안하고 1차 Taylor전개를 이용한 추계론적 유한요소해석법인 가중적분법에 적용하였다. 일반적으로 1차 Taylor전개를 이용하는 수치해석법에서의 응답변화도는 고려하고 있는 추계장의 분산계수에 대하여 선형거동을 보인다. 그러나 몬테카를로 해석의 경우 추계장 분산계수에 대하여 비선형 거동을 나타낸다. 이는 급수전개법의 1차 Taylor전개에 따른 선형특성에 기인한다. 따라서, 가중적분법에서 사용되는 Taylor전개된 변위벡터와 몬테카를로 해석에서의 변위벡터를 비교하고 이들 두 변위벡터 사이에 상호 불일치 하는 점을 고찰하여 몬테카를로 해석에서의 변위벡터와 등가의 변위벡터를 구성하고 이를 가중적분법에 적용하였다. 제안한 등가몬테카를로 추계장은 본래의 추계장 함수에 대한 고차함수로 주어진다. 평면구조에 대한 수치해석을 통하여 제안한 등가몬테카를로 추계장을 이용한 정식화의 타당성을 고찰하였다 새로운 정식화는 기존의 l차 가중적분법을 위한 정식화 과정과 유사하게 수행할 수 있었다.

• Structural Engineering and Mechanics
• /
• 제86권3호
• /
• pp.349-359
• /
• 2023

This paper presents a technique for determining the optimal number of elements in stochastic finite element analysis based on reliability analysis. Using the change-of-variable perturbation stochastic finite element approach, the probability density function of the dynamic responses of stochastic structures is explicitly determined. This method combines the perturbation stochastic finite element method with the change-of-variable technique into a united model. To further examine the relationships between the random fields, discretization of the random field parameters, such as the variance function and the scale of fluctuation, is also performed. Accordingly, the reliability index is calculated based on the explicit probability density function of responses with Gaussian or non-Gaussian random fields in any number of elements corresponding to the random field discretization. The numerical examples illustrate the effectiveness of the proposed method for a one-dimensional cantilever reinforced concrete column and a two-dimensional steel plate shear wall. The benefit of this method is that the probability density function of responses can be obtained explicitly without the use simulation techniques. Any type of random variable with any statistical distribution can be incorporated into the calculations, regardless of the restrictions imposed by the type of statistical distribution of random variables. Consequently, this method can be utilized as a suitable guideline for the efficient implementation of stochastic finite element analysis of structures, regardless of the statistical distribution of random variables.

• 전산구조공학
• /
• 제8권1호
• /
• pp.127-136
• /
• 1995

본 연구는 가중적분법을 이용한 추계론적 유한요소해석에 관한 것으로 구조계 내에 존재하는 재료상수와 기하학적 상수의 임의성을 해석에 고려하여 추계론적 해석을 수행하였으며 대상 구조로는 평판구조를 택하였다. 재료와 기하학적 해석인자의 임의성을 포함한 요소강성행렬의 유도를 위해서 임의장을 가장하였으며 임의장의 평균은 0이고 표준편차 값은 0.1을 사용하였다. 이러한 임의장의 특성은 auto-correlation 함수에 의해서 표현되었으며 이 함수는 반응변화도를 얻는 과정에 사용되었다. 본 연구에서는 평판의 두께에 대한 임의성을 고려하기 위해서 새로운 auto-correlation 함수가 유도되었다. 유도된 새로운 auto-correlation 함수는 재료탄성계수의 임의장 특성을 나타내는 기존의 함수와 임의장 분산 계수의 함수로 나타났다. 수치해석결과는 몬테카를로 시뮬레이션 결과와 비교되었으며 상호 잘 일치하는 좋은 결과를 나타내었고 이들 결과는 제시된 이론적인 수렴치와도 잘 일치하였다. 평판두께에 대한 해석의 경우 역시 Lawrence의 결과는 물론 몬테카를로 시뮬레이션과 제시된 이론치와도 잘 일치하였다.

• Structural Engineering and Mechanics
• /
• 제4권6호
• /
• pp.703-715
• /
• 1996

In stochastic analysis, the randomness of the structural parameters is taken into consideration and the response variability is obtained in addition to the conventional (mean) response. In the present paper the structural response variability of plate structure is calculated using the weighted integral method and is compared with the results obtained by different methods. The stochastic field is assumed to be normally distributed and to have the homogeneity. The decomposition of strain-displacement matrix enabled us to extend the formulation to the stochastic analysis with the quadratic elements in the weighted integral method. A new auto-correlation function is derived considering the uncertainty of plate thickness. The results obtained in the numerical examples by two different methods, i.e., weighted integral method and Monte Carlo simulation, are in a close agreement. In the case of the variable plate thickness, the obtained results are in good agreement with those of Lawrence and Monte Carlo simulation.

• Steel and Composite Structures
• /
• 제11권2호
• /
• pp.115-130
• /
• 2011

Considering the randomness of material parameters in the laminated composite plate, a scheme of stochastic finite element method to analyze the displacement response variability is suggested. In the formulation we adopted the concept of the weighted integral where the random variable is defined as integration of stochastic field function multiplied by a deterministic function over a finite element. In general the elastic modulus of composite materials has distinct value along an individual axis. Accordingly, we need to assume 5 material parameters as random. The correlations between these random parameters are modeled by means of correlation functions, and the degree of correlation is defined in terms of correlation coefficients. For the verification of the proposed scheme, we employ an independent analysis of Monte Carlo simulation with which statistical results can be obtained. Comparison is made between the proposed scheme and Monte Carlo simulation.

• Wind and Structures
• /
• 제16권5호
• /
• pp.433-455
• /
• 2013

A phase delay spectrum model towards the representation of spatial coherence of stochastic wind fields is proposed. Different from the classical coherence functions used in the spectral representation methods, the model is derived from the comprehensive description of coherence of fluctuating wind speeds and from the thorough analysis of physical accounts of random factors affecting phase delay, building up a consistent mapping between the simulated fluctuating wind speeds and the basic random variables. It thus includes complete probabilistic information of spatial stochastic wind fields. This treatment prompts a ready and succinct scheme for the simulation of fluctuating wind speeds, and provides a new perspective to the accurate assessment of dynamic reliability of wind-induced structures. Numerical investigations and comparative studies indicate that the developed model is of rationality and of applicability which matches well with the measured data at spatial points of wind fields, whereby the phase spectra at defined datum mark and objective point are feasibly obtained using the numerical scheme associated with the starting-time of phase evolution. In conjunction with the stochastic Fourier amplitude spectrum that we developed previously, the time history of fluctuating wind speeds at any spatial points of wind fields can be readily simulated.

• Journal of Mechanical Science and Technology
• /
• 제16권7호
• /
• pp.959-965
• /
• 2002

Effective hydraulic conductivity of a statistically anisotropic heterogeneous medium is obtained for steady two-dimensional flows employing stochastic analysis. Flow equations are solved up to second order and the effective conductivity is obtained in a semi-analytic form depending only on the spatial correlation function and the anisotropy ratio of the hydraulic conductivity field, hence becoming a true intrinsic property independent of the flow field. Results are obtained using a statistically anisotropic Gaussian correlation function where the anisotropy is defined as the ratio of integral scales normal and parallel to the mean flow direction. Second order results indicate that the effective conductivity of an anisotropic medium is greater than that of an isotropic one when the anisotropy ratio is less than one and vice versa. It is also found that the effective conductivity has upper and lower bounds of the arithmetic and the harmonic mean conductivities.

• 한국전산구조공학회논문집
• /
• 제23권6호
• /
• pp.727-737
• /
• 2010

Together with the Young's modulus the Poisson's ratio is another independent material parameter that governs the behavior of a structural system. Therefore, it is meaningful to evaluate separately the influence of the parameter on the random response of the structural system. To this end, a formulation dealing with the spatial randomness in the Poisson's ratio in laminated composite plates is proposed. The main idea of the paper is to transform the fraction form of the constitutive coefficients into the expanded form in an ascending order of the stochastic field function. To validate the adequacy of the formulation, a square plate is chosen and the computation results are compared with those obtained using conventional Monte Carlo simulation. It is observed that the results show good agreement with those by the Monte Carlo simulation(MCS).

• Wind and Structures
• /
• 제19권4호
• /
• pp.421-441
• /
• 2014

With the assistance of typhoon field data at aerial elevation level observed by meteorological satellites and wind velocity and direction records nearby the ground gathered in Guangzhou Weather Station between 1985 and 2001, some key wind field parameters under typhoon climate in Guangzhou region were calibrated based on Monte-Carlo stochastic algorithm and Meng's typhoon numerical model. By using Peak Over Threshold method (POT) and Generalized Pareto Distribution (GPD), Wind field characteristics during typhoons for various return periods in several typical engineering fields were predicted, showing that some distribution rules in relation to gradient height of atmosphere boundary layer, power-law component of wind profile, gust factor and extreme wind velocity at 1-3s time interval are obviously different from corresponding items in Chinese wind load Codes. In order to evaluate the influence of typhoon field parameters on long-span flexible bridges, 1:100 reduced-scale wind field of type B terrain was reillustrated under typhoon and normal conditions utilizing passive turbulence generators in TJ-3 wind tunnel, and wind-induced performance tests of aero-elastic model of long-span Guangzhou Xinguang arch bridge were carried out as well. Furthermore, aerodynamic admittance function about lattice cross section in mid-span arch lib under the condition of higher turbulence intensity of typhoon field was identified via using high-frequency force-measured balance. Based on identified aerodynamic admittance expressions, Wind-induced stochastic vibration of Xinguang arch bridge under typhoon and normal climates was calculated and compared, considering structural geometrical non-linearity, stochastic wind attack angle effects, etc. Thus, the aerodynamic response characteristics under typhoon and normal conditions can be illustrated and checked, which are of satisfactory response results for different oncoming wind velocities with resemblance to those wind tunnel testing data under the two types of climate modes.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1993년도 가을 학술발표회논문집
• /
• pp.207-215
• /
• 1993

With the aid of finite element method, this paper deals with the problem of structural response variability of transmission tower subjected to the spatial variability of material properties, Young's modulus herein. The spatial variability of material property are modeled as two-dimensional stochastic field which has an isotropic auto-correlation function. Response variability has been computed based on two numerical techniques, such as the Neumann expansion method in conjunction with the Monte Carlo simulation method. The results by these numerical methods are compared with those by the deterministic approach.