• Structural Engineering and Mechanics
• /
• v.75 no.6
• /
• pp.747-769
• /
• 2020

This paper provides a semi-analytical approach to investigate the variations of 3D displacement components, electric potential, stresses, electric displacements and transverse vibration frequencies in laminated piezoelectric composite plates based on the scaled boundary finite element method (SBFEM) and the precise integration algorithm (PIA). The proposed approach can analyze the static and dynamic responses of multilayered piezoelectric plates with any number of laminae, various geometrical shapes, boundary conditions, thickness-to-length ratios and stacking sequences. Only a longitudinal surface of the plate is discretized into 2D elements, which helps to improve the computational efficiency. Comparing with plate theories and other numerical methods, only three displacement components and the electric potential are set as the basic unknown variables and can be represented analytically through the transverse direction. The whole derivation is built upon the three dimensional key equations of elasticity for the piezoelectric materials and no assumptions on the plate kinematics have been taken. By virtue of the equilibrium equations, the constitutive relations and the introduced set of scaled boundary coordinates, three-dimensional governing partial differential equations are converted into the second order ordinary differential matrix equation. Furthermore, aided by the introduced internal nodal force, a first order ordinary differential equation is obtained with its general solution in the form of a matrix exponent. To further improve the accuracy of the matrix exponent in the SBFEM, the PIA is employed to make sure any desired accuracy of the mechanical and electric variables. By virtue of the kinetic energy technique, the global mass matrix of the composite plates constituted by piezoelectric laminae is constructed for the first time based on the SBFEM. Finally, comparisons with the exact solutions and available results are made to confirm the accuracy and effectiveness of the developed methodology. What's more, the effect of boundary conditions, thickness-to-length ratios and stacking sequences of laminae on the distributions of natural frequencies, mechanical and electric fields in laminated piezoelectric composite plates is evaluated.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

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

This paper presents a numerical simulation method for the flow around advancing ships in regular waves by using a rectilinear grid system. Because the grid lines do not consist with body surface in the rectilinear grid system, the body geometries are defined by the interaction points of those grid lines and the body surface. For the satisfaction of body boundary conditions, no-slip and divergence free conditions are imposed on the body surface and body boundary cells, respectively. Meanwhile, free surface is defined with the modified marker density method. The pressure on the free surface is determined to make the pressure gradient terms of the governing equations continuous, and the velocity around the free surface is calculated with the pressure on the free surface. To validate the present numerical method, a vortex induced vibration (VIV) phenomenon and flows around an advancing Wigley III ship model in various regular waves are simulated, and the results are compared with existing and corresponding research data. Also, to check the applicability to practical ship model, flows around KRISO Container Ship (KCS) model advancing in calm water are numerically simulated. On the simulations, the trim and the sinkage are set free to compare the running attitude with some other experimental data. Moreover, flows around the KCS model in regular waves are also simulated.

• Journal of Ocean Engineering and Technology
• /
• v.29 no.6
• /
• pp.427-435
• /
• 2015

This paper provides the numerical results of a simulation of the flow around a propeller working beneath the free surface. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes (URANS) equations, where the wave-making problem is solved using a volume-of-fluid (VOF) method. The numerical analysis focuses on the propeller wake structure affected by the free surface, where we consider another free surface boundary condition that treats the free surface as a rigid wall surface. The propeller wake under the effect of these two free surface conditions shows a reduction in the magnitude of the longitudinal and vertical flow velocities, and its vortical structures strongly interact with the free surface. The thrust and torque coefficient under the free surface effect decrease about 3.7% and 3.1%, respectively. Finally, the present numerical results show a reasonable agreement with the available experimental data.

• Explosives and Blasting
• /
• v.22 no.3
• /
• pp.13-26
• /
• 2004

This work reports results of our study on the dynamic response of buried pipelines depending on their boundary conditions. We have studied behavior of the buried pipelines both along the axial and the transverse direction with a boundary condition of free ends. The buried pipelines are modeled as beams on elastic foundation while the seismic wave as a ground displacement in the form of a sinusoidal wave. The natural frequency, its mode, and the effect of parameters have been interpreted in terms of the free vibration. In order to investigate the response on the earthquake, the resulting frequency and the mode shape obtained from the free vibration have been utilized to derive the mathematical formula for the farced vibration. We have also completed the computer program to simulate the time-displacement graphs of the pipe lines with free ends for both cases of vibrations.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.5
• /
• pp.564-571
• /
• 2008

This paper is concerned with the free vibration analysis of an annular plate with boundary conditions of simply supported-free, clamped-free and free-free, respectively. Exact solutions for the natural frequency and mode of an annular plate can be obtained by solving the differential equation but other methods such as the Rayleigh-Ritz method and the finite element method can be also used. In this research, we applied the Independent Coordinate Coupling Method(ICCM) to the annular plate and prove that the ICCM can accurately predict the natural frequency and mode shape of the annular plate. The numerical results show that the ICCM can be used effectively for the free vibration problem of plate with a hole compared to the Rayleigh-Ritz method and the finite element method.

• Advances in nano research
• /
• v.10 no.2
• /
• pp.165-174
• /
• 2021

This paper based on Sanders theory aims to investigate the vibration of SWCNTs considering the clamped-simply supported, clamped-free, clamped-clamped and simply supported-simply supported end conditions. After developing the governing equation of the objective system, the Rayleigh-Ritz technique is implemented for the purpose of obtaining the frequency equation in the eigen form. In addition, the applicability of this model for the analysis of vibration of CNTs is examined with the effect of length and ratio of height-to-radius. A detailed description of different types of SWCNTs with different indices is provided in the theoretical methodology. The effect of extended length is stimulated with increasing the radii and the model is effective because it also predicts the effect of thickness on vibration of SWCNTs. For different boundary conditions, the present results are verified with earlier literature.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.35 no.2
• /
• pp.41-48
• /
• 2023

In this study, simplified linear numerical method that can simulate wave generation and transformation by a moving bottom is introduced. Numerical analysis is conducted in wave number domain after continuity equation, linear dynamic and kinematic free surface boundary conditions and linear kinematic bottom boundary condition are Fourier transformed, and the results are expressed in space domain by an inverse Fourier transform. In the wavenumber domain, the dynamic free water surface boundary condition and the kinematic free water surface boundary condition are numerically calculated, and the velocity potential in the mean water level (z = 0) satisfies the continuity equation and the kinematic bottom boundary condition. Wave generation and transformation are investigated when the triangular and rectangular shape of bottoms move periodically. The results of the simplified numerical method are compared with the results of previous analytical solutions and agree well with them. Stability of numerical results according to the calculation time interval (Δt) and the calculation wave number interval (Δk) was also investigated. It was found that the numerical results were appropriate when Δt ≤ T(period)/1000 and Δk ≤ π/100.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.6
• /
• pp.2267-2275
• /
• 2013

Existing conventional model for analysis of shallow water flow just assumed the internal boundary condition as free-slip, which resulted in the wrong prediction about the velocity, vorticity, water level, shear stress distribution, and time variation of drag and lift force around a structure. In this study, a finite element model that can predict flow characteristics around the structure accurately was developed and internal boundary conditions were generalized as partial slip condition using slip length concept. Laminar flow characteristics behind circular cylinder were analyzed by varying the internal boundary conditions. The simulation results of (1) time variations of longitudinal and transverse velocities, and vorticity; (2) wake length; (3) vortex shedding phenomena by slip length; (4) and mass conservation showed that the vortex shedding had never observed and laminar flow like creeping motion was occurred under free-slip condition. Assignment of partial slip condition changed the velocity distribution on the cylinder surface and influenced the magnitude of the shear stress and the occurrence of vorticity so that the period of vortex shedding was reduced compared with the case of no slip condition. The maximum mass conservation error occurred in the case of no slip condition, which had the value of 0.73%, and there was 0.21 % reduction in the maximum mass conservation error by changing the internal boundary condition from no slip to partial slip condition.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.11b
• /
• pp.1260-1265
• /
• 2001

It is important to model the mechanical structure precisely and reasonably in predicting the dynamic characteristics, controlling the vibration, and designing the structural dynamics. In the finite element modeling, the errors can be contained from the physical parameters, the approximation of the boundary conditions, and the element modeling, From the dynamic test. more precise dynamic characteristics can be obtained. Model updating using parameter modification is appropriate when the design parameter is used to analyze the input parameter like finite element method. Finite element analysis for free-free-free-free(FFFF) and clamped-free-free-free(CFFF) plate with uniform area and shape change are carried out as model updating examples, Mass and stiffness matrices are updated by comparing test and analytical modal frequencies. The result shows that the updated frequencies become closer to the test frequencies.