• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1992.10a
• /
• pp.117-123
• /
• 1992

In recent years, the finite element method has become one of the most popular numerical technique for obtaining solutions of engineering science problems. However, there exist various uncertainties in modeling the problems, such as the dimensions(geometry shape), the material properties, boundary conditions, etc. The consideration for the uncertainties inherent in the problems can be made by understanding the influences of uncertain parameters[1]. Determining the influences of uncertainties as statistical quantities using the standard finite element method requires enormous computing time, while the probabilistic finite element method is realized as an efficient scheme[2,3] yielding statistical solution with just a few direct computations. In this paper, a formulation of the probabilistic fluid-structure interaction problem accounting for the first order perturbation of geometric shape is derived, and especially probabilistical acoustic pressure scattering from the structure with surrounding fluid is focused on. In Section 2, governing equations for the fluid-structure problems are given. In Section 3, a finite element formulation, based on the functional, is presented. First order perturbation of geometric shape with randomness is incorporated into the finite element formulation in conjunction with discretization of the random fields in Section 4 and 5. Finally, the proposed formulation is applied to a acoustic pressure scattering problem from an infinitely long cylindrical shell structure with randomness of radial perturbation.

• Structural Engineering and Mechanics
• /
• v.8 no.6
• /
• pp.623-637
• /
• 1999

This paper presents a four-noded quadrilateral $C^0$ strain plate element for the analysis of thick laminated composite plates. The element formulation is based on: 1) the third-order shear deformation theory; 2) assumed strain element formulation; and 3) interrelated edge displacements and rotations along element boundaries. Unlike the existing displacement-type composite plate elements based on the third-order theory, which rely on the $C^1$-continuity formulation, the present plate element is of $C^0$-continuity, and its element stiffness matrix is evaluated explicitly. Because of the third-order expansion of the in-plane displacements through the thickness, the resulting theory and hence elements do not need shear correction factors. The explicit element stiffness matrix makes the present element more computationally efficient than the composite plate elements using numerical integration for the analysis of thick layered composite plates.

• Advances in aircraft and spacecraft science
• /
• v.1 no.2
• /
• pp.125-143
• /
• 2014

This paper provides a new technique for solving the static analysis of arbitrarily shaped composite plates by using Strong Formulation Finite Element Method (SFEM). Several papers in literature by the authors have presented the proposed technique as an extension of the classic Generalized Differential Quadrature (GDQ) procedure. The present methodology joins the high accuracy of the strong formulation with the versatility of the well-known Finite Element Method (FEM). The continuity conditions among the elements is carried out by the compatibility or continuity conditions. The mapping technique is used to transform both the governing differential equations and the compatibility conditions between two adjacent sub-domains into the regular master element in the computational space. The numerical implementation of the global algebraic system obtained by the technique at issue is easy and straightforward. The main novelty of this paper is the application of the stress and strain recovery once the displacement parameters are evaluated. Computer investigations concerning a large number of composite plates have been carried out. SFEM results are compared with those presented in literature and a perfect agreement is observed.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Computers and Concrete
• /
• v.6 no.4
• /
• pp.319-341
• /
• 2009

A new formulation based on lumped plasticity and inelastic hinges is presented in this paper for nonlinear analysis of Reinforced Concrete (RC) frame structures. Inelastic hinge behaviour is described using the principles of Continuum Damage Mechanics (CDM). Member formulation contains provisions to model stiffness degradation due to cracking of concrete and yielding of reinforcing steel. Depending on its nature, cracking is classified as concentrated or distributed. Concentrated cracking is accounted through a damage variable and its growth is defined based on strain energy principles. Presence of distributed flexural cracks in a member is taken care of by modelling it as non-prismatic. Plasticity theory supported by effective stress concept of CDM is applied to describe the post-yield response. Nonlinear quasi-static analysis is carried out on a RC column and a wide two-storey RC frame to verify the formulation. The column is subjected to constant axial load and monotonic lateral load while the frame is subjected to only lateral load. Computed results are compared with those due to experiments or other numerical methods to validate the performance of the formulation and also to highlight the contribution of distributed cracking on global response.

• Journal of applied mathematics & informatics
• /
• v.10 no.1_2
• /
• pp.175-191
• /
• 2002

The equations of a polar fluid of hydromagnetic fluctuating through a porous medium axe cast into matrix form using the state space and Laplace transform techniques the resulting formulation is applied to a variety of problems. The solution to a problem of an electrically conducting polar fluid in the presence of a transverse magnetic field and to a problem for the flow between two parallel fixed plates is obtained. The inversion of the Laplace transforms is carried out using a numerical approach. Numerical results for the velocity, angular velocity distribution and the induced magnetic field are given and illustrated graphically for each problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.9-10
• /
• 2006

A numerical method for simulating tree surface flows including the surface tension is presented. Numerical scheme is based an a fractional-step method with a finite volume formulation and the interface between liquid and gas is tracked by Volume of Fluid (VOF) method. Piecewise Linear Interface Calculation (PLIC) method is used to reconstruct the interface and the surface tension is considered using a Continuum Surface Force (CSF) model. Several free surface flow phenomena were simulated to show its effectiveness to find such phenomena.

• Structural Engineering and Mechanics
• /
• v.45 no.6
• /
• pp.853-865
• /
• 2013

Inflatable panels made of modern and new textile materials can be inflated at high pressure to have a high mechanical strength. This paper is based on the finite element method as a general solution to determine the characteristics of deformed inflatable panels at high pressure in various end and loading conditions. Proposed method is based on the construction of weak form of formulation and application of Reduced Integration Element method (RIE) to solve the numerical problem of shear locking. The numerical results are validated as an outcome of comparison with other published results.

• Journal of the Korean Mathematical Society
• /
• v.38 no.4
• /
• pp.723-761
• /
• 2001

The objects of this paper are to formulated a model for the transport of a chain of radioactive waste products in a fractured porous medium, to devise an effective and efficient numerical method for approximating the solution of the model, and to demonstrated the convergence of the numerical method. The formulation begins from a model in an unfractured (single porosity) medium, passes through a double porosity model in a fractured medium, and ends with a modified single porosity model that takes the relevant time scales of the flow and the nuclear decay.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.71-78
• /
• 2018

This paper presents the assembling of multiple patches based on the single patch isogeometric formulation for the shear deformable shell element given in the previous study. The geometrically exact shell formulation has been accomplished with the shell theory based formulation and the generalized curvilinear coordinate system directly derived from the given NURBS geometry. For the knot elements matching across adjacent surfaces, the zero-th and first parametric continuity conditions are considered and the corresponding coupling constraints are implemented by a master-slave formulation between adjacent patches. The constraints are then enforced by a substitution method for condensation of the slave variables, thereby reducing the model size. Through numerical investigations, the important features of the first parametric continuity condition are confirmed. The performance of the multi-patch shell models is also examined comparing the rate of convergence of response coefficients for the zero and first order continuity conditions and continuity in coupling boundary between two patches is confirmed.