• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.8
• /
• pp.745-752
• /
• 2009

In robust design, the mean and variance of design performance are frequently used to measure the design performance and its robustness under uncertainties. In this paper, we present the Gauss-type quadrature formula as a rigorous method for mean and variance estimation involving arbitrary input distributions and further extend its use to robust design optimization. One dimensional Gauss-type quadrature formula are constructed from the input probability distributions and utilized in the construction of multidimensional quadrature formula such as the tensor product quadrature (TPQ) formula and the univariate dimension reduction (UDR) method. To improve the efficiency of using it for robust design optimization, a semi-analytic design sensitivity analysis with respect to the statistical moments is proposed. The proposed approach is applied to a simple bench mark problems and robust topology optimization of structures considering various types of uncertainty.

• Journal of applied mathematics & informatics
• /
• v.14 no.1_2
• /
• pp.473-484
• /
• 2004

In this paper we give an expression for the kernel and remainder term of Gauss-Radau and Gauss-Lobatto quadratures and study on an error bound with end points of multiplicity γ.

• Communications of the Korean Mathematical Society
• /
• v.12 no.3
• /
• pp.797-812
• /
• 1997

For analytic functions we give an expression for the kernel $K_n$ of the remainder terms for the Gauss-Radau and the Gauss-Lobatto rules with end points of multiplicity r and prove the convergence of the kernel we obtained. The error bound are obtained for the type $$\mid$R_n(f)$\mid$ \leq \frac{1}{\pi}l(\Gamma) max_{z \in \Gamma} $\mid$K_n(z)$\mid$ max_{z \in \Gamma} $\mid$f(z)$\mid$$, where $l(\Gamma)$ denotes the length of contour $\Gamma$.

• Structural Engineering and Mechanics
• /
• v.76 no.4
• /
• pp.529-540
• /
• 2020

Concurrent topology optimization of macrostructure and microstructure has attracted significant interest due to its high structural performance. However, most of the existing works are carried out under deterministic conditions, the obtained design may be vulnerable or even cause catastrophic failure when the load position exists uncertainty. Therefore, it is necessary to take load position uncertainty into consideration in structural design. This paper presents a computational method for robust concurrent topology optimization with consideration of load position uncertainty. The weighted sum of the mean and standard deviation of the structural compliance is defined as the objective function with constraints are imposed to both macro- and micro-scale structure volume fractions. The Bivariate Dimension Reduction method and Gauss-type quadrature (BDRGQ) are used to quantify and propagate load uncertainty to calculate the objective function. The effective properties of microstructure are evaluated by the numerical homogenization method. To release the computation burden, the decoupled sensitivity analysis method is proposed for microscale design variables. The bi-directional evolutionary structural optimization (BESO) method is used to obtain the black-and-white designs. Several 2D and 3D examples are presented to validate the effectiveness of the proposed robust concurrent topology optimization method.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1053-1067
• /
• 2009

The objective of this study is to develop an efficient and accurate quadratic finite element model based on Streamline Upwind/Petrov Galerkin (SU/PG) scheme for analyzing and predicting two dimensional flow features in complex natural rivers. For a development of model, quadratic tin, quadrilateral and mixed elements as well as linear tin, quadrilateral and mixed elements were used in the model. Also, this model was developed through reinforcement of Gauss Quadrature which was necessary to integral of governing equation. Several tests for bottom-rising channel and U-type channel were performed for the purpose of validation and verification of the developed model. Such results showed that solutions of second order elements are better accurate and improved than those of linear elements. Results obtained by the developed model and RMA-2 model are compared, and the results for the developed model were better accurate than those of RMA-2 model. In the future if the developed model is applied in natural rivers, it can provide better accurate results than those of existing model.

• Structural Engineering and Mechanics
• /
• v.66 no.3
• /
• pp.329-341
• /
• 2018

This paper considers the smooth receding contact problem between a homogeneous half-plane and a composite laminate composed of an inhomogeneously coated elastic layer. The inhomogeneity of the elastic modulus of the coating is approximated by an exponential function along the thickness dimension. The three-component structure is pressed together by either a concentrated force or uniform pressures applied at the top surface of the composite laminate. Both semianalytical and finite element analysis are performed to solve for the extent of contact and the contact pressure. In the semianalytical formulation, Fourier integral transformation of governing equations and boundary conditions leads to a singular integral equation of Cauchy-type, which can be numerically integrated by Gauss-Chebyshev quadrature to a desired degree of accuracy. In the finite element modeling, the functionally graded coating is divided into homogeneous sublayers and the shear modulus of each sublayer is assigned at its lower boundary following the predefined exponential variation. In postprocessing, the stresses of any node belonging to sublayer interfaces are averaged over its surrounding elements. The results obtained from the semianalytical analysis are successfully validated against literature results and those of the finite element modeling. Extensive parametric studies suggest the practicability of optimizing the receding contact peak stress and the extent of contact in multilayered structures by the introduction of functionally graded coatings.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.173-195
• /
• 2013

Meshfree methods are known to have the capability to overcome the strict regularization requirements and numerical instabilities that encumber the finite element method (FEM) in large deformation problems. They are also more naturally suited for problems involving material perforation and fragmentation. To take advantage of the high efficiency of FEM and high accuracy of meshfree methods, a coupled finite element (FE) and reproducing kernel (RK, one of the meshfree approximations) formulation is described in this paper. The coupling of FE and RK approximation is implemented in an evolutionary fashion, where the extent and location of the evolution is dependent on a triggering criteria provided by the material constitutive laws. To enhance computational efficiency, Gauss quadrature is applied to integrate both FE and RK domains so that no state variable transfer is required when mesh conversion is performed. To control the hourglassing that might occur with 1-point integrated hexahedral grids, viscous type hourglass control is implemented. Meanwhile, the FEM version of the K&C concrete (KCC) model was modified to make it applicable in both FE and RK formulations. Results using this code and the KCC model are shown for the modeling of concrete responses under quasi-static, blast and impact loadings. These analyses demonstrate that fragmentation phenomena of the sort commonly observed under blast and impact loadings of concrete structures was able to be realistically captured by the coupled formulation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.109-122
• /
• 1983

An accurate thick beam element (TB4) which includes the effects of the shear deformation and rotatory inertia has been degenerated from the three dimensional continuum by employing the Timoshenko beam assumptions. The proposed TB4 element has four nodes and two degrees of freedom at each node, totally eight degrees of freedom. The transverse deflection W and plane rotation ${\theta}$ with the cubic interpolation functions are selected as nodal variables. The element characteristics are formulated by discretizing the beam equations of motion, using the Galerkin weighted residual method, and are numerically integrated by the reduced shear integration technique, using the three-point Gauss quadrature with the various shear coefficients. Several numerical examples are analyzed to demonstrate the accuracy and the monotonic convergence behavior of the proposed TB4 beam element. The result indicates that the TB4 element shows the more excellent performance and the monotonic convergence behavior than the other existing Timoshenko beam type elements for the whole range of the beam aspect ratios, in both static and free vibration analyses.