• Steel and Composite Structures
• /
• v.39 no.3
• /
• pp.291-306
• /
• 2021

Based on Reissner's mixed variational theorem (RMVT), the authors develop a semi-analytical finite element (FE) method for a three-dimensional (3D) bending analysis of nonhomogeneous orthotropic, complete and incomplete toroidal shells subjected to uniformly-distributed loads. In this formulation, the toroidal shell is divided into several finite annular prisms (FAPs) with quadrilateral cross-sections, where trigonometric functions and serendipity polynomials are used to interpolate the circumferential direction and meridian-radial surface variations in the primary field variables of each individual prism, respectively. The material properties of the toroidal shell are considered to be nonhomogeneous orthotropic over the meridianradial surface, such that homogeneous isotropic toroidal shells, laminated cross-ply toroidal shells, and single- and bi-directional functionally graded toroidal shells can be included as special cases in this work. Implementation of the current FAP methods shows that their solutions converge rapidly, and the convergent FAP solutions closely agree with the 3D elasticity solutions available in the literature.

• Interaction and multiscale mechanics
• /
• v.2 no.3
• /
• pp.277-293
• /
• 2009

The Enriched Free Mesh Method (EFMM) is a patch-wise procedure in which both a displacement field on an element and a stress/strain field on a cluster of elements connected to a node can be defined. On the other hand, the Superconvergent Patch Recovery (SPR) is known to be an efficient post-processing procedure of the finite element method to estimate the error norm at a node. In this paper, we discuss the relationship between solutions of the EFMM and those of the SPR through several convergence studies. In addition, in order to solve the demerit of the smoothing effect on the fracture mechanics fields, we implement a singular stress field to a local patch in the EFMM, and its effectiveness is investigated.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.11a
• /
• pp.173-176
• /
• 2000

Finite element model based on the Naghdi's shell theory in the general tensor-based form is formulated in the present study. Partial mixed variational functional for assumed strain is formulated in order to avoid the severe locking troubles known as transverse shear and membrane locking. The proposed assumed strain element in general tensor Naghdi's shell model provides very accurate solutions for thin shells in benchmark problems. In additions, linear elastic constitutive equations are given in the general curvilinear coordinate system including anisotropic layered structures. Thus laminated composited shell structures are easily analyzed in the present formulation.

• Journal of the Korean Mathematical Society
• /
• v.57 no.5
• /
• pp.1239-1266
• /
• 2020

In this paper, some novel discrete formulations for stabilizing the mixed finite element method Q1-Q0 (bilinear velocity and constant pressure approximations) are introduced and discussed for the generalized Stokes problem. These are based on stabilizing discontinuous pressure approximations via local jump operators. The developing idea consists in a reduction of terms in the local jump formulation, introduced earlier, in such a way that stability and convergence properties are preserved. The computer implementation aspects and numerical evaluation of these stabilized discrete formulations are also considered. For illustrating the numerical performance of the proposed approaches and comparing the three versions of the local jump methods alongside with the global jump setting, some obtained results for two test generalized Stokes problems are presented. Numerical tests confirm the stability and accuracy characteristics of the resulting approximations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1327-1334
• /
• 2003

Flow characteristics of the fluid laden with many particles in the two-dimensional channel are investigated using the Navier-Stokes equations coupled with the equation of motion of particles by direct numerical simulation. A four-step fractional step method with Crank-Nicolson scheme and ALE technique is used for P2P1 mixed finite element method. The motion and distribution of particles in the fluid is virtually described as a result of direct numerical simulation and the increase of viscosity is compared with theoretical equations. The effect of channel height on the relative viscosity and the tubular pinch effect are discussed.

• Composites Research
• /
• v.30 no.6
• /
• pp.343-349
• /
• 2017

A multifield variational formulation is developed for the finite element (FE) cross-sectional analysis of composite beams. The cross-sectional warping displacements and sectional stresses are considered to be the primary variables through the application of Reissner's partially mixed principle. The warping displacements are modeled using generic FE shape functions with nonlinear distribution over the beam section. A generalized Timoshenko level stiffness matrix is derived which incorporates the effects of elastic couplings, transverse shear, and Poisson's deformations. The accuracy of the present analysis is validated for the stiffness constants and elastostatic responses of composite box beams which correlate well with the experimental data and other state-of-the-art approaches.

• Structural Engineering and Mechanics
• /
• v.4 no.6
• /
• pp.645-653
• /
• 1996

In order to avoid pathological mesh dependency in finite element modelling of strain localization, an isotropic elasto-plastic model with a yield function depending on the Laplacian of the equivalent plastic strain is implemented in a 4-node quadrilateral finite element with one integration point based on a mixed formulation derived from Hu-Washizu principle. The evaluation of the Laplacian is based on a least square polynomial approximation of the equivalent plastic strain around each integration point. This non local approach allows to satisfy exactly the consistency condition at each integration point. Some examples are treated to illustrate the effectiveness of the method.

• Journal of the Korean Mathematical Society
• /
• v.59 no.3
• /
• pp.519-548
• /
• 2022

In this paper, a stabilized-penalized collocated finite volume (SPCFV) scheme is developed and studied for the stationary generalized Navier-Stokes equations with mixed Dirichlet-traction boundary conditions modelling an incompressible biological fluid flow. This method is based on the lowest order approximation (piecewise constants) for both velocity and pressure unknowns. The stabilization-penalization is performed by adding discrete pressure terms to the approximate formulation. These simultaneously involve discrete jump pressures through the interior volume-boundaries and discrete pressures of volumes on the domain boundary. Stability, existence and uniqueness of discrete solutions are established. Moreover, a convergence analysis of the nonlinear solver is also provided. Numerical results from model tests are performed to demonstrate the stability, optimal convergence in the usual L² and discrete H¹ norms as well as robustness of the proposed scheme with respect to the choice of the given traction vector.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.11
• /
• pp.1348-1358
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation fur flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among seven filling patterns at each tetrahedral control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. The collapse of a water dam and the filling of a fluidity spiral have been analyzed. The numerical results have been in good agreement with the experimental results and the efficiency of the adaptive grid refinement and filling pattern techniques have been verified.