• 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.

• Structural Engineering and Mechanics
• /
• 제51권4호
• /
• pp.601-625
• /
• 2014

Based on continuum mechanics and the principle of virtual displacements, incremental total Lagrangian formulation (T.L.) and incremental updated Lagrangian formulation (U.L.) were presented. Both T.L. and U.L. considered the large displacement stiffness matrix, which was modified to be symmetrical matrix. According to the incremental updated Lagrangian formulation, small strain, large displacement, finite rotation of three dimensional Timoshenko fiber beam element tangent stiffness matrix was developed. Considering large displacement and finite rotation, a new type of tangent stiffness matrix of the beam element was developed. According to the basic assumption of plane section, the displacement field of an arbitrary fiber was presented in terms of nodal displacement of centroid of cross-area. In addition, shear deformation effect was taken account. Furthermore, a nonlinear finite element method program has been developed and several examples were tested to demonstrate the accuracy and generality of the three dimensional beam element.

• Nuclear Engineering and Technology
• /
• 제49권1호
• /
• pp.29-42
• /
• 2017

The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.

• Advances in Computational Design
• /
• 제5권3호
• /
• pp.277-289
• /
• 2020

This article describes the technology for constructing of a multiply-connected three-dimensional area's finite element representation. Representation of finite-element configuration of an area is described by a discrete set that consist of the number of nodes and elements of the finite-element grid, that are orderly set of nodes' coordinates and numbers of finite elements. Corresponding theorems are given, to prove the correctness of the solution method. The adequacy of multiply-connected area topology's finite element model is shown. The merging of subareas is based on the criterion of boundary nodes' coincidence by establishing a simple hierarchy of volumes, surfaces, lines and points. Renumbering nodes is carried out by the frontal method, where nodes located on the outer edges of the structure are used as the initial front.

• 소성∙가공
• /
• 제12권3호
• /
• pp.194-201
• /
• 2003

An automatic quadrilateral mesh generator for large deformation finite element analysis such as metal forming simulation was developed. The NURBS interpolation method is used for modeling arbitrary 2-D free surface. This mesh generation technique is the modified paving algorithm, which is an advancing front technique with element-by-element resolving method for paving boundary intersection problem. The mesh density for higher analysis accuracy and less analysis time can be easily controlled with high-density points, maximum and minimum element size. A couple of application to large deformation finite element analysis is given as an example, which shows versatility and applicability of the proposed approach and the developed mesh generator for large deformation finite element analysis.

• 소음진동
• /
• 제8권2호
• /
• pp.267-273
• /
• 1998

This study aims at performing sensitivity analysis of piezoelectric smart structure for minimizing radiated noise from the structure, The structure consists of a flat plate on which disk shaped piezoelectric actuator is mounted, and finite element modeling is used for the structure. The finite element modeling uses a combination of three dimensional piezoelectric, flat shell and transition elements so thus it can take into account the coupling effects of the piezoelectric device precisely and it can also reduce the degrees of freedom of the finite element model. Electric potential on the piezoelectric actuator is taken as a design variable and total radiated power of the structure is chosen as an objective function. The objective function can be represented as Rayleigh's integral equation and is a function of normal displacements of the structure. For the convenience of computation, all degrees of freedom of the finite element equation is condensed out except the normal displacements of the structure. To perform the design sensitivity analysis, the derivative of the objective function with respect to the normal displacements is found, and the derivative of the norma displacements with respect to the design variable is calculated from the finite element equation by using so called the adjoint variable method. The analysis results are compared with those of the finite difference method, and shows a good agreement. This sensitivity analysis is faster and more accurate than the finite difference method. Once the sensitivity analysis program is used for gradient-based optimizations, one could achieve a better convergence rate than non-derivative methods for optimal design of piezoelectric smart structures.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2007년도 춘계학술대회논문집
• /
• pp.859-865
• /
• 2007

This paper presents a method to analyze the unbalance response of a high speed polygon mirror scanner motor supported by sintered bearing and flexible supporting structures by using the finite element method and the mode superposition method. The appropriate finite element equations for polygon mirror are described by rotating annular sector element using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. The rotating components except for the polygon mirror are modeled by Timoshenko beam element including the gyroscopic effect. The flexible supporting structures are modeled by using a 4-node tetrahedron element and 4-node shell element with rotational degrees of freedom. Finite element equations of each component of the polygon mirror scanner motor and the flexible supporting structures are consistently derived by satisfying the geometric compatibility in the internal boundary between each component. The rigid link constraints are also imposed at the interface area between sleeve and sintered bearing to describe the physical motion at this interface. A global matrix equation obtained by assembling the finite element equations of each substructure is transformed to a state-space matrix-vector equation, and both damped natural frequencies and modal damping ratios are calculated by solving the associated eigenvalue problem by using the restarted Arnoldi iteration method. Unbalance responses in time and frequency domain are performed by superposing the eigenvalues and eigenvectors from the free vibration analysis. The validity of the proposed method is verified by comparing the simulated unbalance response with the experimental results. This research also shows that the flexibility of supporting structures plays an important role in determining the unbalance response of the polygon mirror scanner motor.

• Structural Engineering and Mechanics
• /
• 제16권6호
• /
• pp.655-674
• /
• 2003

Investigated in this study are the natural vibration characteristics of the coaxial cylindrical shells coupled with a fluid. Theoretical method is developed to find the natural frequencies of the shell using the finite Fourier series expansion, and their results are compared with those of finite element method to verify the validation of the method developed. The effect of the fluid-filled annulus and the boundary conditions on the modal characteristics of the coaxial shells is investigated using a finite element modeling.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집A
• /
• pp.458-463
• /
• 2000

The mechanical structures generally have discontinuous parts such as the cracks, notches and holes owing to various reasons. In this paper, in order to analyze effectively these singularity problems using the finite element method, a mixed analysis method which an analytical solution and finite element solutions are simultaneously used is newly proposed. As the analytical solution is used in the singularity region and the finite element solutions are used in the remaining regions except this singular zone, this analysis method reasonably provides for the numerical solution of a singularity problem. Through various numerical examples, it is shown that the proposed analysis method is very convenient and gives comparatively accurate solution.

• 한국광학회:학술대회논문집
• /
• 한국광학회 1989년도 제4회 파동 및 레이저 학술발표회 4th Conference on Waves and lasers 논문집 - 한국광학회
• /
• pp.157-160
• /
• 1989

The most serious difficulty in using the finite element method is the appearance of the so-called spurious, nonphysical modes. We have proposed the finite element formulation of the variational expression in the three-component magnetic field based on Galerkin's method. In this approach, the divergence relation H is satisfied and spurious modes does not appear and finite-element solutions agree with the exact solutions.