• Structural Engineering and Mechanics
• /
• 제3권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
• /
• 제23권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년도 춘계 학술대회논문집
• /
• 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.

• 대한기계학회논문집
• /
• 제18권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
• /
• 제2권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.

• 전산구조공학
• /
• 제10권3호
• /
• pp.185-193
• /
• 1997

탄성 선형 경화 재료로 구성된 복합 구조물의 자유 경계면에서 나타나는 응력 특이도를 평면 변형률 상태에서 계산하였다. 자유 표면력 경계조건과 계면 연속조건을 만족해야하는 지배 탄성 방정식은 2점 경계치문제로 정의되며, 일반 고유치 문제의 해인 고유치가 응력 특이도가 될 것이다. 자유경계면 근처에서 응력 성분을 r/sup s-1/에 비례한다고 가정하여 특정한 s(고유치)를 구하는 고유치 문제를 뉴톤향상법과 사격법을 사용하여 수치적으로 해를 구하였다.

• 한국전산구조공학회논문집
• /
• 제31권2호
• /
• pp.71-78
• /
• 2018

본 연구에서는 NURBS 기반 아이소 지오메트릭 쉘 해석을 위해 다중 패치 해석 모델을 정식화하였다. 기존 연구를 통해 개발된 단일 패치로 구성된 전단 변형을 고려한 쉘 요소에 대해 일반 좌표계에서 기하학적으로 엄밀한 쉘 구조물의 아이소 지오메트릭 해석 모델을 도입하고 매개변수 연속성을 고려하여 다중 패치 모델로 확장하였다. 인접 곡면의 노트 요소가 결합 경계를 통해 조화를 이루는 경우에 대해 0차와 1차 매개변수 연속성 조건을 고려하였으며, 두 패치 간 마스터-슬레이브 관계를 정립하여 종속된 한 곡면의 자유도를 상대 곡면의 자유도로 표시하여 모델 크기를 줄이면서 두 곡면을 결합하였다. 다중 패치 쉘 예제에 대해 0차와 1차 연속성 조건을 각각 적용하여 구조해석을 수행하여 1차 연속성 조건의 주요한 특성들을 확인하였다. 또한 각 연속성 조건에 대한 해의 수렴 특성을 검토하였으며 결합 경계에서의 두 패치의 연속성을 확인하였다.

• Journal of Mechanical Science and Technology
• /
• 제14권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.

• 한국전산유체공학회지
• /
• 제4권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
• /
• 제20권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.