• Structural Engineering and Mechanics
• /
• v.3 no.2
• /
• pp.121-133
• /
• 1995

It is well known that two-dimensional simplified third-order theories satisfy the layer interface continuity of transverse shear strains, thus these theories violate the continuity of transverse shear stresses when two consecutive layers differ either in fibre orientation or material. The third-order theories considered herein involve four/or five dependent unknowns in the displacement field and satisfy the condition of vanishing of transverse shear stresses at the bounding planes of the plate. The objective of this investigation is to examine (i) the flexural response prediction accuracy of these third-order theories compared to exact elasticity solution (ii) the effect of layer interface continuity conditions on the flexural response. To investigate the effect of layer interface continuity conditions, three-dimensional elasticity solutions are developed by enforcing the continuity of different combinations of transverse stresses and/or strains at the layer interfaces. Three dimensional twenty node solid finite element (having three translational displacements as degrees of freedom) without the imposition of any of the conditions on the transverse stresses and strains is also employed for the flexural analysis of the laminated plates for the purposes of comparison with the above theories. These shear deformation theories and elasticity approaches in terms of accuracy, adequacy and applicability are examined through extensive numerical examples.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.65-92
• /
• 2019

We survey a recently developed immersed finite element method (IFEM) for the interface problems. The IFEM uses structured grids such as uniform grids, even if the interface is a smooth curve. Instead of fitting the curved interface, the bases are modified so that they satisfy the jump conditions along the interface. The early versions of IFEM [1, 2] were suboptimal in convergence order [3]. Later, the consistency terms were added to the bilinear forms [4, 5], thus the scheme became optimal and the error estimates were proven. For elasticity problems with interfaces, we modify the Crouzeix-Raviart based element to satisfy the traction conditions along the interface [6], but the consistency terms are not needed. To satisfy the Korn's inequality, we add the stabilizing terms to the bilinear form. The optimal error estimate was shown for a triangular grid. Lastly, we describe the multigrid algorithms for the discretized system arising from IFEM. The prolongation operators are designed so that the prolongated function satisfy the flux continuity condition along the interface. The W-cycle convergence was proved, and the number of V-cycle is independent of the mesh size.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.1
• /
• pp.65-76
• /
• 1994

A higher-oder laminated plate theory including the effect of transverse shear deformation is developed to calculate the gross response and the detailed stress distribution. The theory satisfies the continuity condition of transverse shear stress, and accounts for parabolic variation of the transverse shear stresses through the thickness of each layer. Exact closed-ply laminates are obtained and the results are compared with three-dimensional elasticity solutions and a simple higher-order theory solutions. The results of the present work exhibit acceptable accuracy when compared to the three-dimensional elasticity solutions.

• Interaction and multiscale mechanics
• /
• v.2 no.2
• /
• pp.129-151
• /
• 2009

This paper presents a meshfree procedure using a convex generalized meshfree (GMF) approximation for the large deformation analysis of particle-reinforced rubber compounds on microscopic level. The convex GMF approximation possesses the weak-Kronecker-delta property that guarantees the continuity of displacement across the material interface in the rubber compounds. The convex approximation also ensures the positive mass in the discrete system and is less sensitive to the meshfree nodal support size and integration order effects. In this study, the convex approximation is generated in the GMF method by choosing the positive and monotonic increasing basis function. In order to impose the periodic boundary condition in the unit cell method for the microscopic analysis, a singular kernel is introduced on the periodic boundary nodes in the construction of GMF approximation. The periodic boundary condition is solved by the transformation method in both explicit and implicit analyses. To simulate the interface de-bonding phenomena in the rubber compound, the cohesive interface element method is employed in corporation with meshfree method in this study. Several numerical examples are presented to demonstrate the effectiveness of the proposed numerical procedure in the large deformation analysis.

• Computational Structural Engineering
• /
• v.10 no.3
• /
• pp.185-193
• /
• 1997

The order of the stress singularity that occurs at the termination of an interface between materials exhibiting bilinear stress-strain response under plane strain conditions has been calculated, The governing equation of elasticity together with traction-free boundary condition and interface continuity condition defines a two-point boundary value problem. The stress components near the free edge are assumed to be proportional to r/sup s-1/, with solutions existing only for certain values of s. Finding these values entails the solution of an eigenvalue problem. Because it has been impossible to integrate the differential equations analytically, the integration has been performed numerically with a shooting method coupled with a Newton improvement scheme.

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

• Journal of Mechanical Science and Technology
• /
• v.14 no.9
• /
• pp.928-934
• /
• 2000

This paper describes a method developed for the simulation of ring pack lubrication characteristic in an internal combustion engine. In general, the quantity of oil supply for piston ring lubrication may be insufficient in filling the entire volume formed at the interference between the piston ring and the cylinder liner. Thus the oil starvation condition should be considered in analyzing piston ring lubrication. In order to reasonably estimate the amount of oil left over on the cylinder liner, the flow rate at the posterior portion of the interface should be calculated with an adequate boundary condition that confirms flow continuity condition. In this analysis, oil starvation and open-end boundary conditions are considered at the inlet and outlet of the piston rings. The lubrication characteristic of each piston ring is obtained by an iterative method with sequential steps. It is revealed that piston rings are operated under oil starvation in most operating cycles and the result under these conditions are quite different from that with the fully-flooded assumption.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2230-2241
• /
• 2006

This paper presents a sole application of boundary element method to the conjugate heat transfer problem of thermally developing laminar flow in a thick walled pipe when the fluid velocities are fully developed. Due to the coupled mechanism of heat conduction in the solid region and heat convection in the fluid region, two separate solutions in the solid and fluid regions are sought to match the solid-fluid interface continuity condition. In this method, the dual reciprocity boundary element method (DRBEM) with the axial direction marching scheme is used to solve the heat convection problem and the conventional boundary element method (BEM) of axisymmetric model is applied to solve the heat conduction problem. An iterative and numerically stable BEM solution algorithm is presented, which uses the coupled interface conditions explicitly instead of uncoupled conditions. Both the local convective heat transfer coefficient at solid-fluid interface and the local mean fluid temperature are initially guessed and updated as the unknown interface thermal conditions in the iterative solution procedure. Two examples imposing uniform temperature and heat flux boundary conditions are tested in thermally developing region and compared with analytic solutions where available. The benchmark test results are shown to be in good agreement with the analytic solutions for both examples with different boundary conditions.