• 대한수학회지
• /
• 제56권4호
• /
• pp.1083-1112
• /
• 2019

In this paper, we present and analyze a cell-centered collocated finite volume scheme for incompressible flows to compute solutions simultaneous to Stokes and Darcy equations by applying a pressure jump stabilization term to avoid locking. We prove that the new stabilized FV formulation satisfies a discrete inf-sup condition and error estimates for both problems. Finally, we present some numerical examples confirming this analysis.

• 한국전산구조공학회논문집
• /
• 제13권1호
• /
• pp.37-49
• /
• 2000

본 논문은 구조물의 동역학 및 열탄성 연성문제 해석을 위한 통합된 유한요소법을 개발하는데 초점을 두고있다. 첫째로, 열전도 방정식에 열변위라는 물리량을 도입하여 동역학의 운동 방정식과 유사하도록 유도한 후, 변분법과 일반좌표계를 이용하여 시간영역에서 정식화하였다. 둘째로, 두 방정식에 라플라스 변환을 동시에 도입하고, 공간변수만을 갖는 형상함수와 가중잔여법을 적용하여 유한요소식을 변환영역에서 표현하였다. 연성된 방정식을 문제의 특성에 따라서 분류하였고 정식화 과정을 검증하였다. 또한 수치해석 알고리듬이 갖는 수치 역 변환의 정성적인 경향에 대하여 검토하였다.

• Structural Engineering and Mechanics
• /
• 제33권6호
• /
• pp.675-696
• /
• 2009

In this paper, a unified solution method is presented for the classical beam theory. In Strength of Materials approach, the geometry, material properties and load system are known and related with the unknowns of forces, moments, slopes and deformations by applying a classical differential analysis in addition to equilibrium, constitutive, and kinematic laws. All these relations are expressed in a unified formulation for the classical beam theory. In the special case of simple beams, a system of four linear ordinary differential equations of first order represents the general mechanical behaviour of a straight beam. These equations are solved using the numerical differential quadrature method (DQM). The application of DQM has the advantages of mathematical consistency and conceptual simplicity. The numerical procedure is simple and gives clear understanding. This systematic way of obtaining influence line, bending moment, shear force diagrams and deformed shape for the beams with geometric and load discontinuities has been discussed in this paper. Buckling loads and natural frequencies of any beam prismatic or non-prismatic with any type of support conditions can be evaluated with ease.

• Steel and Composite Structures
• /
• 제34권2호
• /
• pp.241-260
• /
• 2020

This article presented a comprehensive model to study static buckling stability and associated mode-shapes of higher shear deformation theories of sandwich laminated composite beam under the compression of varying axial load function. Four higher order shear deformation beam theories are considered in formulation and analysis. So, the model can consider the influence of both thick and thin beams without needing to shear correction factor. The compression force can be described through axial direction by uniform constant, linear and parabolic distribution functions. The Hamilton's principle is exploited to derive equilibrium governing equations of unified sandwich laminated beams. The governing equilibrium differential equations are transformed to algebraic system of equations by using numerical differential quadrature method (DQM). The system of equations is solved as an eigenvalue problem to get critical buckling loads and their corresponding mode-shapes. The stability of DQM in determining of buckling loads of sandwich structure is performed. The validation studies are achieved and the obtained results are matched with those. Parametric studies are presented to figure out effects of in-plane load type, sandwich thickness, fiber orientation and boundary conditions on buckling loads and mode-shapes. The present model is important in designing process of aircraft, naval structural components, and naval structural when non-uniform in-plane compressive loading is dominated.

• Smart Structures and Systems
• /
• 제16권5호
• /
• pp.781-806
• /
• 2015

A unified formulation of finite layer methods (FLMs), based on the Reissner mixed variational theorem (RMVT), is developed for the three-dimensional (3D) coupled electro-elastic analysis of simply-supported, functionally graded piezoelectric material (FGPM) plates with open- and closed-circuit surface conditions and under electro-mechanical loads. In this formulation, the material properties of the plate are assumed to obey an exponent-law varying exponentially through the thickness coordinate, and the plate is divided into a number of finite rectangular layers, in which the trigonometric functions and Lagrange polynomials are used to interpolate the in- and out-of-plane variations of the primary field variables of each individual layer, respectively, such as the elastic displacement, transverse shear and normal stress, electric potential, and normal electric displacement components. The relevant orders used for expanding these variables in the thickness coordinate can be freely chosen as the linear, quadratic and cubic orders. Four different mechanical/electrical loading conditions applied on the top and bottom surfaces of the plate are considered, and the corresponding coupled electro-elastic analysis of the loaded FGPM plates is undertaken. The accuracy and convergence rate of the RMVT-based FLMs are assessed by comparing their solutions with the exact 3D piezoelectricity ones available in the literature.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2006년도 추계학술대회논문집
• /
• pp.934-937
• /
• 2006

This investigation presents a topology formulation to design optimal poroelastic acoustic foams to maximize absorbing ability. For successful formulation, a single set of equations based on Biot's theory is adopted and an appropriate material interpolation strategy is newly developed. Because there was no earlier attempt to solve poroelastic acoustic foam design problems in topology optimization setting, many challenging issues including modeling and interpolation must be addressed. First, the simulation accuracy by a proposed unified model encompassing acoustic air and poroelastic material was checked against analytical and numerical results. Then a material interpolation scheme yielding a distinct acoustic air-poroelastic material distribution was developed. Using the proposed model and interpolation scheme, the topology optimization of a two-dimensional poroelastic acoustic foam for maximizing its absorption coefficient was carried out. Numerical results show that the absorption capacity of an optimized foam layout considerably increases in comparison with a nominal foam layout.

• Steel and Composite Structures
• /
• 제4권4호
• /
• pp.281-301
• /
• 2004

The paper begins by presenting a unified variational approach to the lateral-torsional buckling (LTB) analysis of doubly symmetric prismatic and tapered thin-walled beams with open cross-sections, which accounts for the influence of the pre-buckling deflections. This approach (i) extends the kinematical assumptions usually adopted for prismatic beams, (ii) consistently uses shell membrane theory in general coordinates and (iii) adopts Trefftz's criterion to perform the bifurcation analysis. The proposed formulation is then applied to investigate the influence of the pre-buckling deflections on the LTB behaviour of prismatic and web-tapered I-section simply supported beams and cantilevers. After establishing an interesting analytical result, valid for prismatic members with shear centre loading, several elastic critical moments/loads are presented, discussed and, when possible, also compared with values reported in the literature. These numerical results, which are obtained by means of the Rayleigh-Ritz method, (i) highlight the qualitative differences existing between the LTB behaviours of simply supported beams and cantilevers and (ii) illustrate how the influence of the pre-buckling deflections on LTB is affected by a number of factors, namely ($ii_1$) the minor-to-major inertia ratio, ($ii_2$) the beam length, ($ii_3$) the location of the load point of application and ($ii_4$) the bending moment diagram shape.

• 한국전산유체공학회지
• /
• 제9권1호
• /
• pp.18-24
• /
• 2004

The unified simulation for the multiphase flow by predictor-corrector scheme based on CIP method is introduced. In this algorithm, the interface between different phases is identified by a density function and tracked by solving an advection equation. Solid body motion is modeled by the translation and angular motion. The mathematical formulation and numerical results are also described. To verify the efficiency, accuracy and capability of proposed algorithm, two dimensional incompressible cavity flow, the motion of a floating ball into water and a single rising bubble by buoyancy force are numerically simulated by the present scheme. As results, it is confirmed that the present scheme gives an efficient, stable and reasonable solution in the multiphase flow problem.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제6권1호
• /
• pp.1-15
• /
• 2002

Advection-dominated transport problems possess difficulties in the design of numerical methods for solving them. Because of the hyperbolic nature of advective transport, many characteristic numerical methods have been developed such as the classical characteristic method, the Eulerian-Lagrangian method, the transport diffusion method, the modified method of characteristics, the operator splitting method, the Eulerian-Lagrangian localized adjoint method, the characteristic mixed method, and the Eulerian-Lagrangian mixed discontinuous method. In this paper relationships among these characteristic methods are examined. In particular, we show that these sometimes diverse methods can be given a unified formulation. This paper focuses on characteristic finite element methods. Similar examination can be presented for characteristic finite difference methods.

• Structural Engineering and Mechanics
• /
• 제26권2호
• /
• pp.191-210
• /
• 2007

The three different approaches (reformulations) of evolution strategies (ESs) have been proposed in the literature as extensions of the technique for solving discrete problems. This study implements an extensive research on application, evaluation and comparison of them in discrete optimum design of steel frames. A unified formulation is first developed to explain these approaches, so that differences and similarities between their inherent search mechanisms can clearly be identified. Two examples from practical design of steel frames are studied next to measure their performances in locating the optimum. Extensive numerical experimentations are performed in both examples to facilitate a statistical analysis of their convergence characteristics. The results obtained are presented in the histograms demonstrating the distribution of the best designs located by each approach. In addition, an average improvement of the best design during the course of evolution is plotted in each case to compare their relative convergence rates.