• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.8
• /
• pp.1822-1831
• /
• 1995

A general formulation of the design optimization problem with the random parameters is presented here. The formulation is based on the stochastic finite element second-order perturbation method ; it takes into full account of the stress and displacement constraints together with the rates of change of the random variables. A method of direct differentiation for calculating the sensitivity coefficients in regard to the governing equation and the second-order perturbed equation is derived. A gradient-based nonlinear programming technique is used to solve the problem. The numerical results are specifically noted, where the stiffness parameter and external load are treated as random variables.

• Structural Engineering and Mechanics
• /
• v.11 no.4
• /
• pp.373-392
• /
• 2001

The main idea behind the paper is to present two alternative methods of homogenization of the heat conduction problem in composite materials, where the heat conductivity coefficients are assumed to be random variables. These two methods are the Monte-Carlo simulation (MCS) technique and the second order perturbation second probabilistic moment method, with its computational implementation known as the Stochastic Finite Element Method (SFEM). From the mathematical point of view, the deterministic homogenization method, being extended to probabilistic spaces, is based on the effective modules approach. Numerical results obtained in the paper allow to compare MCS against the SFEM and, on the other hand, to verify the sensitivity of effective heat conductivity probabilistic moments to the reinforcement ratio. These computational studies are provided in the range of up to fourth order probabilistic moments of effective conductivity coefficient and compared with probabilistic characteristics of the Voigt-Reuss bounds.

• Steel and Composite Structures
• /
• v.26 no.1
• /
• pp.1-15
• /
• 2018

This study reported, the effect of random variation in system properties on bending response of single wall carbon nanotube reinforced composite (SWCNTRC) plates subjected to transverse uniform loading is examined. System parameters such as the SWCNT armchair, material properties, plate thickness and volume fraction of SWCNT are modelled as basic random variables. The basic formulation is based on higher order shear deformation theory to model the system behaviour of the SWCNTRC composite plate. A C0 finite element method in conjunction with the first order perturbation technique procedure developed earlier by the authors for the plate subjected to lateral loading is employed to obtain the mean and variance of the transverse deflection of the plate. The performance of the stochastic SWCNTRC composite model is demonstrated through a comparison of mean transverse central deflection with those results available in the literature and standard deviation of the deflection with an independent First Order perturbation Technique (FOPT), Second Order perturbation Technique (SOPT) and Monte Carlo simulation.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.32 no.4
• /
• pp.432-446
• /
• 1996

A stochastic Hamilton variational principle(SHVP) is formulated for dynamic problems of linear continuum. The SHVP allows incorporation of probabilistic distributions into the finite element analysis. The formulation is simplified by transformation of correlated random variables to a set of uncorrelated random variables through a standard eigenproblem. A procedure based on the Fourier analysis and synthesis is presented for eliminating secularities from the perturbation approach. In addition to, a method to analyse stochastic design sensitivity for structural dynamics is present. A combination of the adjoint variable approach and the second order perturbation method is used in the finite element codes. An alternative form of the constraint functional that holds for all times is introduced to consider the time response of dynamic sensitivity. The algorithms developed can readily be adapted to existing deterministic finite element codes. The numerical results for stochastic analysis by proceeding approach of cantilever, 2D-frame and 3D-frame illustrates in this paper.

• Structural Engineering and Mechanics
• /
• v.20 no.3
• /
• pp.293-312
• /
• 2005

In this paper, free vibration of rotating beam with random properties is studied. The cross-sectional area, elasticity modulus, moment of inertia, shear modulus and density are modeled as random fields and the rotational speed as a random variable. To study uncertainty, stochastic finite element method based on second order perturbation method is applied. To discretize random fields, the three methods of midpoint, interpolation and local average are applied and compared. The effects of rotational speed, setting angle, random property variances, discretization scheme, number of elements, correlation of random fields, correlation function form and correlation length on "Coefficient of Variation" (C.O.V.) of first mode eigenvalue are investigated completely. To determine the significant random properties on the variation of first mode eigenvalue the sensitivity analysis is performed. The results are studied for both Timoshenko and Bernoulli-Euler rotating beam. It is shown that the C.O.V. of first mode eigenvalue of Timoshenko and Bernoulli-Euler rotating beams are approximately identical. Also, compared to uncorrelated random fields, the correlated case has larger C.O.V. value. Another important result is, where correlation length is small, the convergence rate is lower and more number of elements are necessary for convergence of final response.

• Steel and Composite Structures
• /
• v.44 no.5
• /
• pp.651-676
• /
• 2022

Most of the experimental, theoretical, and numerical studies on the stability of functionally graded composites are deterministic, while there are full of complex interactions of variables with an inherently probabilistic nature, this paper presents a non-intrusive framework to investigate the stochastic nonlinear buckling behaviors of porous functionally graded cylindrical shells exposed to inevitable source-uncertainties. Euler-Lagrange equations are theoretically derived based on the three variable refined shear deformation theory. Closed-form solutions for the shell buckling loads are achieved by solving the deterministic eigenvalue problems. The analytical results are verified with numerical results obtained from finite element analyses that are conducted in the commercial software ABAQUS. The non-intrusive framework is completed by integrating the Monte Carlo simulation with the verified closed-form solutions. The convergence studies are performed to determine the effective pseudorandom draws of the simulation. The accuracy and efficiency of the framework are verified with statistical results that are obtained from the first and second-order perturbation techniques. Eleven cases of individual and compound uncertainties are investigated. Sensitivity analyses are conducted to figure out the five cases that have profound perturbative effects on the shell buckling loads. Complete probability distributions of the first three critical buckling loads are completely presented for each profound uncertainty case. The effects of the shell thickness, volume fraction index, and stochasticity degree on the shell buckling load under compound uncertainties are studied. There is a high probability that the shell has non-unique buckling modes in stochastic environments, which should be known for reliable analysis and design of engineering structures.

• Korean Management Science Review
• /
• v.9 no.1
• /
• pp.3-15
• /
• 1992

In this paper, we study a manufacturing system of serial stages with general service times, in which the production of each stage and the coordination of stages are controlled by Kanban discipline. This Kanban discipline is modeled as a Discrete Event Dynamic System and a system of recursive equations is applied to study the dynamics of the system. The recursive relationship enables us to compare this Kanban discipline with the other blocking disciplines such as transfer blocking, service blocking, block-and-hold b, and block-and-hold K, and the Kanban is shown to be superior to the other disciplines in terms of makespan and throughput. As a special case, two stages Kanban system is modeled as $C_2/C_2/1/N$ queueing system, and a recursive algorithm is developed to calculate the system performance. In optimizing the system performance, the stochastic optimization approach of Robbins-Monro is employed via perturbation analysis, the way to estimate the stochastic partial derivative based on only one sample trajectory of the system, and the required commuting condition is verified. Then the stochastic convexity result is established to provide second-order optimality condition for this parametric optimization problem.

• Advances in aircraft and spacecraft science
• /
• v.4 no.5
• /
• pp.585-611
• /
• 2017

This article covenants with the post buckling witticism of carbon nanotube reinforced composite (CNTRC) beam supported with an elastic foundation in thermal atmospheres with arbitrary assumed random system properties. The arbitrary assumed random system properties are be modeled as uncorrelated Gaussian random input variables. Unvaryingly distributed (UD) and functionally graded (FG) distributions of the carbon nanotube are deliberated. The material belongings of CNTRC beam are presumed to be graded in the beam depth way and appraised through a micromechanical exemplary. The basic equations of a CNTRC beam are imitative constructed on a higher order shear deformation beam (HSDT) theory with von-Karman type nonlinearity. The beam is supported by two parameters Pasternak elastic foundation with Winkler cubic nonlinearity. The thermal dominance is involved in the material properties of CNTRC beam is foreseen to be temperature dependent (TD). The first and second order perturbation method (SOPT) and Monte Carlo sampling (MCS) by way of CO nonlinear finite element method (FEM) through direct iterative way are offered to observe the mean, coefficient of variation (COV) and probability distribution function (PDF) of critical post buckling load. Archetypal outcomes are presented for the volume fraction of CNTRC, slenderness ratios, boundary conditions, underpinning parameters, amplitude ratios, temperature reliant and sovereign random material properties with arbitrary system properties. The present defined tactic is corroborated with the results available in the literature and by employing MCS.

• Computers and Concrete
• /
• v.32 no.1
• /
• pp.99-117
• /
• 2023

The crack under the influence of the higher intensities of the stresses grows and the structure gets collapsed with the time when the crack length reaches to critical value. Therefore, the fracture behavior of a structure in terms of stress intensity factors (SIF) becomes important to determine the remaining fracture strength and capacity of material and structure for avoiding catastrophic failure, increasing safety and further improvement in the design. The robustness of the method has been demonstrated by comparing the numerical results with analytical and experimental results of some problems. XFEM is used to model cracks and holes in structures and predict their strength and reliability under service conditions. Further, XFEM is extended with a stochastic method for predicting the sensitivity in terms of output COVs and fracture strength in terms of mean values of stress intensity factors (SIFs) of a structure with discontinuities (cracks and holes) under tensile loading condition with input individual and combined randomness in different system parameters. In stochastic technique, the second order perturbation technique (SOPT) has been used for the predicting the fracture behavior of the structures. The stochastic/perturbation technique is also known as Taylor series expansion method and it provides the reliable results if the input randomness is less than twenty percentage. From the present numerical analysis it is observed that, the crack tip near to the hole is under the influence of the stress concentration and the variational effect of the input random parameters on the crack tip in terms of the SIFs are lesser so the COVs are the less sensitive. The COVs of mixed mode SIFs are the most sensitive for the crack angles (α=45° to 90°) for all the values of c₁ and d₁. The plate with the shorter distance between hole and crack is the most sensitive with all the crack angles but the crack tip which is much nearer to the hole has the highest sensitivity.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.35 no.4
• /
• pp.435-447
• /
• 1999

It is main objective of this approach to present a method to analyse stochastic design sensitivity for problems of structural dynamics with randomness in design parameters. A combination of the adjoint variable approach and the second order perturbation method is used in the finite element approach. An alternative form of the constant functional that holds for all times is introduced to consider the time response of dynamic sensitivity. The terminal problem of the adjoint system is solved using equivalent homogeneous equations excited by initial velocities. The numerical procedures are shown to be much more efficient when based on the fold superposition method: the generalized co-ordinates are normalized and the correlated random variables are transformed to uncorrelated variables, whereas the secularities are eliminated by the fast Fourier transform of complex valued sequences. Numerical algorithms have been worked out and proved to be accurate and efficient : they can be readily adapted to fit into the existing finite element codes whose element derivative matrices can be explicitly generated. The numerical results of two cases -2 dimensional portal frame for the comparison with reference and 3-dimensional frame structure - for the deterministic sensitivity analysis are presented.