• Steel and Composite Structures
• /
• v.47 no.5
• /
• pp.551-568
• /
• 2023

In this research, a new two-dimensional (2D) and quasi three-dimensional (quasi-3D) higher order shear deformation theory is devised to address the bending problem of functionally graded plates resting on an elastic foundation. The displacement field of the suggested theories takes into account a parabolic transverse shear deformation shape function and satisfies shear stress free boundary conditions on the plate surfaces. It is expressed as a combination of trigonometric and exponential shear shape functions. The Pasternak mathematical model is considered for the elastic foundation. The material properties vary constantly across the FG plate thickness using different distributions as power-law, exponential and Mori-Tanaka model. By using the virtual works principle and Navier's technique, the governing equations of FG plates exposed to sinusoidal and evenly distributed loads are developed. The effects of material composition, geometrical parameters, stretching effect and foundation parameters on deflection, axial displacements and stresses are discussed in detail in this work. The obtained results are compared with those reported in earlier works to show the precision and simplicity of the current formulations. A very good agreement is found between the predicted results and the available solutions of other higher order theories. Future mechanical analyses of three-dimensionally FG plate structures can use the study's findings as benchmarks.

• Steel and Composite Structures
• /
• v.44 no.5
• /
• pp.721-740
• /
• 2022

In this paper, an accurate kinematic model has been developed to study the mechanical response of functionally graded (FG) sandwich beams, mainly covering the bending, buckling and free vibration problems. The studied structure with homogeneous hardcore and softcore is considered to be simply supported in the edges. The present model uses a new refined shear deformation beam theory (RSDBT) in which the displacement field is improved over the other existing high-order shear deformation beam theories (HSDBTs). The present model provides good accuracy and considers a nonlinear transverse shear deformation shape function, since it is constructed with only two unknown variables as the Euler-Bernoulli beam theory but complies with the shear stress-free boundary conditions on the upper and lower surfaces of the beam without employing shear correction factors. The sandwich beams are composed of two FG skins and a homogeneous core wherein the material properties of the skins are assumed to vary gradually and continuously in the thickness direction according to the power-law distribution of volume fraction of the constituents. The governing equations are drawn by implementing Hamilton's principle and solved by means of the Navier's technique. Numerical computations in the non-dimensional terms of transverse displacement, stresses, critical buckling load and natural frequencies obtained by using the proposed model are compared with those predicted by other beam theories to confirm the performance of the proposed theory and to verify the accuracy of the kinematic model.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.151-158
• /
• 2009

To implement the insects' flapping flight for developing flapping MAVs(micro air vehicles), the unsteady flow characteristics of the insects' forward flight is investigated. In this paper, two-dimensional FSI(Fluid-Structure Interaction) simulations are conducted to examine realistic flow features of insects' flapping flight and to examine the flexibility effects of the insect's wing. The unsteady incompressible Navier-Stokes equations with an artificial compressibility method are implemented as the fluid module while the dynamic finite element equations using a direct integration method are employed as the solid module. In order to exchange physical information to each module, the common refinement method is employed as the data transfer method. Also, a simple and efficient dynamic grid deformation technique based on Delaunay graph mapping is used to deform computational grids. Compared to the earlier researches of two-dimensional rigid wing simulations, key physical phenomena and flow patterns such as vortex pairing and vortex staying can still be observed. For example, lift is mainly generated during downstroke motion by high effective angle of attack caused by translation and lagging motion. A large amount of thrust is generated abruptly at the end of upstroke motion. However, the quantitative aspect of flow field is somewhat different. A flexible wing generates more thrust but less lift than a rigid wing. This is because the net force acting on wing surface is split into two directions due to structural flexibility. As a consequence, thrust and propulsive efficiency was enhanced considerably compared to a rigid wing. From these numerical simulations, it is seen that the wing flexibility yields a significant impact on aerodynamic characteristics.

• Wind and Structures
• /
• v.22 no.4
• /
• pp.429-454
• /
• 2016

In this paper, a shear deformation plate theory based on neutral surface position is developed for free vibration analysis of functionally graded material (FGM) plates. The material properties of the FGM plates are assumed to vary through the thickness of the plate by a simple power-law distribution in terms of the volume fractions of the constituents. During manufacture, defects such as porosities can appear. It is therefore necessary to consider the vibration behavior of FG plates having porosities in this investigation. The proposed theory is based on assumption that the in-plane and transverse displacements consist of bending and shear components, in which the bending components do not contribute toward shear forces and, likewise, the shear components do not contribute toward bending moments. The neutral surface position for a functionally graded plate which its material properties vary in the thickness direction is determined. The equation of motion for FG rectangular plates is obtained through Hamilton's principle. The closed form solutions are obtained by using Navier technique, and then fundamental frequencies are found by solving the results of eigenvalue problems. Numerical results are presented and the influences of the volume fraction index and porosity volume fraction on frequencies of FGM plates are clearly discussed.

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

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.11
• /
• pp.895-903
• /
• 2007

CMBT(Curved Moving Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the curved solid wall of moving obstacle in a flow field. In our research CMBT was used to perform LBM simulation of a flow over a moving circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of CMBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a fixed circular cylinder in a channel and the results were compared against the solution of Navier-Stokes equation with deforming mesh technique. The simulations were performed in a moderate range of reynolds number at each moving cylinder to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is ar Re=250 and the result is the same as the case of fixed cylinder. As the cylinder approaching to one wall, the 2nd vortex is developed by interacting with the wall boundary-layer vorticity. As the velocity ratio increase the third vortex are generated by interacting with the 2nd vortexes developed on the upper and lower wall boundary layer. The resultant $C_d$ decrease as reynolds number increasing and the Cd approached to a value when Re>1000.

• Steel and Composite Structures
• /
• v.34 no.4
• /
• pp.511-524
• /
• 2020

In this research, a simple four-variable trigonometric integral shear deformation model is proposed for the static behavior of advanced functionally graded (AFG) ceramic-metal plates supported by a two-parameter elastic foundation and subjected to a nonlinear hygro-thermo-mechanical load. The elastic properties, including both the thermal expansion and moisture coefficients of the plate, are also supposed to be varied within thickness direction by following a power law distribution in terms of volume fractions of the components of the material. The interest of the current theory is seen in its kinematics that use only four independent unknowns, while first-order plate theory and other higher-order plate theories require at least five unknowns. The "in-plane displacement field" of the proposed theory utilizes cosine functions in terms of thickness coordinates to calculate out-of-plane shear deformations. The vertical displacement includes flexural and shear components. The elastic foundation is introduced in mathematical modeling as a two-parameter Winkler-Pasternak foundation. The virtual displacement principle is applied to obtain the basic equations and a Navier solution technique is used to determine an analytical solution. The numerical results predicted by the proposed formulation are compared with results already published in the literature to demonstrate the accuracy and efficiency of the proposed theory. The influences of "moisture concentration", temperature, stiffness of foundation, shear deformation, geometric ratios and volume fraction variation on the mechanical behavior of AFG plates are examined and discussed in detail.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Structural Engineering and Mechanics
• /
• v.74 no.3
• /
• pp.365-380
• /
• 2020

The focus of this paper is to develop an analytical approach based on an efficient shear deformation theory with stretching effect for bending stress analysis of cross-ply laminated composite plates subjected to transverse parabolic load and line load by using a new kinematic model, in which the axial displacements involve an undetermined integral component in order to reduce the number of unknowns and a sinusoidal function in terms of the thickness coordinate to include the effect of transverse shear deformation. The present theory contains only five unknowns and satisfies the zero shear stress conditions on the top and bottom surfaces of the plate without using any shear correction factors. The governing differential equations and its boundary conditions are derived by employing the static version of principle of virtual work. Closed-form solutions for simply supported cross-ply laminated plates are obtained applying Navier's solution technique, and the numerical case studies are compared with the theoretical results to verify the utility of the proposed model. Lastly, it can be seen that the present outlined theory is more accurate and useful than some higher-order shear deformation theories developed previously to study the static flexure of laminated composite plates.