• 한국해양학회지
• /
• v.30 no.1
• /
• pp.13-26
• /
• 1995

In this paper three turbulence models including two-equation model by Blumberg and Mellor (1987), one-equation model with mixing length formula of Blackadar's (1962), and zero-equation model of Prandtl's (1925) were compared in homogeneous, unstratified channel flows. Steady flows which a steep-sided trapezoidal trench with uniform discharge, tidal flow and steady wind-driven flow in finite channels are considered in detail. Steady flows in a trench and tidal flows in a finite channel were reproduced fairly accurately and there was virtually no difference among results of three turbulence models. However, In case of steady wind-driven flow only two-equation model reproduced the important features of experimental data. the other two models underestimated the surface velocity. In tidal and wind-driven flows with negligibly small adjective and diffusive effects, the two-equation model gives rise to parabolic profile of eddy viscosity with maximum at the mid0depth, and the one and zero equation model based on Blackadar formula linear profile with maximum at the surface.

• Structural Engineering and Mechanics
• /
• v.62 no.6
• /
• pp.675-693
• /
• 2017

In this paper, based on the first-order shear deformation plate theory, buckling analysis of piezoelectric coupled transversely isotropic rectangular plates is investigated. By assuming the transverse distribution of electric potential to be a combination of a parabolic and a linear function of thickness coordinate, the equilibrium equations for buckling analysis of plate with surface bonded piezoelectric layers are established. The Maxwell's equation and all boundary conditions including the conditions on the top and bottom surfaces of the plate for closed and open circuited are satisfied. The analytical solution is obtained for Levy type of boundary conditions. The accurate buckling load of laminated plate is presented for both open and closed circuit conditions. From the numerical results it is found that, the critical buckling load for open circuit is more than that of closed circuit in all boundary and loading conditions. Furthermore, the critical buckling loads and the buckling mode number increase by increasing the thickness of piezoelectric layers for both open and closed circuit conditions.

• Communications of the Korean Mathematical Society
• /
• v.10 no.1
• /
• pp.39-48
• /
• 1995

The ADI method was introduced by Peaceman and Rachford [6] in 1955, to solve the discretized boundary value problems for elliptic and parabolic PDEs. The finite difference discretization of the model elliptic problem $$ (1) -\Delta u = f, \Omega = [0, 1] \times [0, 1] $$ $$ u = 0 on \delta \Omega $$ with 5-point centered finite difference discretization, with n +2 mesh-points in the x - direction and m + 2 points in the y direction, leads to the solution of a linear system of equations of the form $$ (2) Au = b $$ where A is a matrix of dimension $N = n \times m$. Without loss of generality and for the sake of simplicity, we will assume for the remainder of this paper that m = n, so that $N = n^2$.

• Structural Engineering and Mechanics
• /
• v.13 no.6
• /
• pp.713-730
• /
• 2002

This paper presents a finite element approach for determining the natural frequencies for planar inclined arches of various shapes vibrating in three-dimensional space. The profile of inclined arches, represented by undeformed centriodal axis of cross-section, is defined by the equation of plane curves expressed in the rectangular coordinates which are : circular, parabolic, sine, elliptic, and catenary shapes. In free vibration state, the arch is slightly displaced from its undeformed position. The linear relationship between curvature-torsion and axial strain is expressed in terms of the displacements in three-dimensional space. The finite element discretization along the span length is used rather than the total are length. Numerical results for arches of various shapes are given and they are in good agreement with those reported in literature. The natural frequency parameters and mode shapes are reported as functions of two nondimensional parameters: the span to cord length ratio (e) and the rise to cord length ratio (f).

• Smart Structures and Systems
• /
• v.19 no.4
• /
• pp.403-413
• /
• 2017

The present article encompasses a nonlinear finite element (FE) and genetic algorithm (GA) based optimal vibration energy harvesting from nonprismatic piezo-laminated cantilever beams. Three cases of cross section profiles (such as linear, parabolic and cubic) are modelled to analyse the geometric nonlinear effects on the output responses such as displacement, voltage, and power. The simultaneous effects of taper ratios (such as breadth and height taper) on the output power are also studied. The FE based nonlinear dynamic equation of motion has been solved by an implicit integration method (i.e., Newmark method in conjunction with the Newton-Raphson method). Besides this, a real coded GA based constrained optimization scheme has also been proposed to determine the best set of design variables for optimal harvesting of power within the safe limits of beam stress and PZT breakdown voltage.

• Structural Engineering and Mechanics
• /
• v.60 no.2
• /
• pp.271-299
• /
• 2016

In this article, based on the higher-order shear deformation plate theory, buckling analysis of a rectangular plate made of functionally graded piezoelectric materials and its effective parameters are investigated. Assuming the transverse distribution of electric potential to be a combination of a parabolic and a linear function of thickness coordinate, the equilibrium equations for the buckling analysis of an FGP rectangular plate are established. In addition to the Maxwell equation, all boundary conditions including the conditions on the top and bottom surfaces of the plate for closed and open circuited are satisfied. Considering double sine solution (Navier solution) for displacement field and electric potential, an analytical solution is obtained for full simply supported boundary conditions. The accurate buckling load of FGP plate is presented for both open and closed circuit conditions. It is found that the critical buckling load for open circuit is more than that of closed circuit in all loading conditions. Furthermore, it is observed that the influence of dielectric constants on the critical buckling load is more than those of others.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.115-122
• /
• 2002

The main purpose of this paper is to determine the static optimal shapes of tapered beams with constant volume. The linear, parabolic and sinusoidal tapers with the regular polygon cross-section are considered, whose material volume and span length are always held constant. The Runge-Kutta method is used to integrate the differential equation and also Shooting method is used to calculate the unknown boundary condition. Then the static optimal shapes are determined by reading the minimum values of the deflection versus section ratio curves plotted by the deflection data. In numerical examples, the various tapered beams are analyzed and those numerical results of this study are shown in figures.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2005.04a
• /
• pp.1-4
• /
• 2005

A higher order zig-zag shell theory is developed to refine accurately predict deformation and stress of smart shell structures under the mechanical, thermal, and electric loading. The displacement fields through the thickness are constructed by superimposing linear zig-zag field to the smooth globally cubic varying field. Smooth parabolic distribution through the thickness is assumed in the transverse deflection in order to consider transverse normal deformation. The mechanical, thermal, and electric loading is applied in the sinusoidal distribution function in the in-surface direction. Thermal and electric loading is given in the linear variation through the thickness. Especially, in electric loading case, voltage is only applied in piezo-layer. The layer-dependent degrees of freedom of displacement fields are expressed in terms of reference primary degrees of freedom by applying interface continuity conditions as well as bounding surface conditions of transverse shear stresses. In order to obtain accurate transverse shear and normal stresses, integration of equilibrium equation approach is used. The numerical examples of present theory demonstrate the accuracy and efficiency of the proposed theory. The present theory is suitable for the predictions of behaviors of thick smart composite shell under mechanical, thermal, and electric loadings combined.

• Structural Engineering and Mechanics
• /
• v.70 no.5
• /
• pp.525-534
• /
• 2019

An analytical research on buckling of simply supported thin plate with variable thickness under bi-axial compression is presented in this paper. Combining the perturbation technique, Fourier series expansion and Galerkin methods, the linear governing differential equation of the plate with arbitrary thickness variation under bi-axial compression is solved and the analytical expression of the critical buckling load is obtained. Based on that, numerical analysis is carried out for the plates with different thickness variation forms and aspect ratios under different bi-axial compressions. Four different thickness variation forms including linear, parabolic, stepped and trigonometric have been considered in this paper. The calculated critical buckling loads and buckling modes are presented and compared with the published results in the tables and figures. It shows that the analytical expressions derived by the theoretical method in this paper can be effectively used for buckling analysis of simply supported thin plates with arbitrary thickness variation, especially for the stepped thickness that used in engineering widely.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.58 no.1
• /
• pp.81-90
• /
• 2016

The Agricultural Systems Application Platform (ASAP) provides bottom-up modelling and simulation environment for agricultural engineer. The purpose of this study is to expand usability of the ASAP to the second order partial differential equations: elliptic equations, parabolic equations, and hyperbolic equations. The ASAP is a general-purpose simulation tool which express natural phenomenon with capsulized independent components to simplify implementation and maintenance. To use the ASAP in continuous problems, it is necessary to solve partial differential equations. This study shows usage of the ASAP in elliptic problem, parabolic problem, and hyperbolic problem, and solves of static heat problem, heat transfer problem, and wave problem as examples. The example problems are solved with the ASAP and Finite Difference method (FDM) for verification. The ASAP shows identical results to FDM. These applications are useful to simulate the engineering problem including equilibrium, diffusion and wave problem.