• Structural Engineering and Mechanics
• /
• v.52 no.1
• /
• pp.115-136
• /
• 2014

In this study, an extended Kantorovich method, employing multi-term displacement functions, is applied to analyze the vibration problem of symmetrically laminated plates with arbitrary boundary conditions. The vibration behaviors of laminated plates are determined based on the variational principle of total energy minimization and the iterative Kantorovich method. The out-of-plane displacement is represented in the form of a series of a sum of products of functions in x and y directions. With a known function in the x or y directions, the formulation for the variation of total potential energy is transformed to a set of governing equations and a set of boundary conditions. The equations and boundary conditions are then numerically solved for the natural frequency and vibration mode shape. The solutions are verified with available solutions from the literature and solutions from the Ritz and finite element analysis. In most cases, the natural frequencies compare very well with the reference solutions. The vibration mode shapes are also very well modeled using the multi-term assumed displacement function in the terms of a power series. With the method used in this study, it is possible to solve the angle-ply plate problem, where the Kantorovich method with single-term displacement function is ineffective.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.8 s.251
• /
• pp.993-1000
• /
• 2006

Operating experience of steam generators has shown that cracks of various morphology frequently occur in the steam generator tubes. These cracked tubes can stay in service if it is proved that the tubes have sufficient safety margin to preclude the risk of burst and leak. Therefore, integrity assessment using exact limit load solutions is very important for safe operation of the steam generators. This paper provides global and local limit load solutions for surface cracks in the steam generator tubes. Such solutions are developed based on three-dimensional (3-D) finite element analyses assuming elastic-perfectly plastic material behavior. For the crack location, both axial and circumferential surface cracks, and for each case, both external and internal cracks are considered. The resulting global and local limit load solutions are given in polynomial forms, and thus can be simply used in practical integrity assessment of the steam generator tubes.

• Steel and Composite Structures
• /
• v.6 no.5
• /
• pp.401-415
• /
• 2006

Based on a non-linear model taking into account flexural-torsional couplings, analytical solutions are derived for lateral buckling of simply supported I beams under some representative load cases. A closed form is established for lateral buckling moments. It accounts for bending distribution, load height application and pre-buckling deflections. Coefficients $C_1$ and $C_2$ affected to these parameters are then derived. Regard to well known linear stability solutions, these coefficients are not constant but depend on another coefficient $k_1$ that represents the pre-buckling deflection effects. In numerical simulations, shell elements are used in mesh process. The buckling loads are achieved from solutions of eigenvalue problem and by bifurcations observed on non linear equilibrium paths. It is proved that both the buckling loads derived from linear stability and eigenvalue problem lead to poor results, especially for I sections with large flanges for which the behaviour is predominated by pre-buckling deflection and the coefficient $k_1$ is large. The proposed solutions are in good agreement with numerical bifurcations observed on non linear equilibrium paths.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.7
• /
• pp.585-596
• /
• 2013

This study introduces a hybrid approach combining numerical results with pre-developed analytical calculations for the sound radiation from the modal vibration of a thick, finite length cylinder with various boundary conditions. Structural vibrations of the cylinder are numerically investigated with the finite element method, and distributions of vibratory displacements on the cylinder surface are idealized as simple mathematical expressions based on the numerical results. Sound radiations from the normal vibration of the cylinder are calculated based on idealized modal displacements using a previously introduced theoretical solution. The results are confirmed with numerical analyses using the boundary element method. Based on these results, it can be concluded that the solutions suggested in this study have good accuracies in calculating the vibro-acoustic properties of a thick, finite cylinder with various boundary conditions. Also, the sound radiation characteristics of many practical components such as brake drums and motor housings are expected to be investigated using the procedure proposed in this study.

• Coupled systems mechanics
• /
• v.7 no.6
• /
• pp.649-668
• /
• 2018

In this paper, we present a numerical model for fluid-structure interaction between structure built of porous media and acoustic fluid, which provides both pore pressure inside porous media and hydrodynamic pressures and hydrodynamic forces exerted on the upstream face of the structure in an unified manner and simplifies fluid-structure interaction problems. The first original feature of the proposed model concerns the structure built of saturated porous medium whose response is obtained with coupled discrete beam lattice model, which is based on Voronoi cell representation with cohesive links as linear elastic Timoshenko beam finite elements. The motion of the pore fluid is governed by Darcy's law, and the coupling between the solid phase and the pore fluid is introduced in the model through Biot's porous media theory. The pore pressure field is discretized with CST (Constant Strain Triangle) finite elements, which coincide with Delaunay triangles. By exploiting Hammer quadrature rule for numerical integration on CST elements, and duality property between Voronoi diagram and Delaunay triangulation, the numerical implementation of the coupling results with an additional pore pressure degree of freedom placed at each node of a Timoshenko beam finite element. The second original point of the model concerns the motion of the outside fluid which is modeled with mixed displacement/pressure based formulation. The chosen finite element representations of the structure response and the outside fluid motion ensures for the structure and fluid finite elements to be connected directly at the common nodes at the fluid-structure interface, because they share both the displacement and the pressure degrees of freedom. Numerical simulations presented in this paper show an excellent agreement between the numerically obtained results and the analytical solutions.

• Structural Engineering and Mechanics
• /
• v.16 no.2
• /
• pp.153-176
• /
• 2003

From the equation of motion of a "bare" non-uniform beam (without any concentrated elements), an eigenfunction in term of four unknown integration constants can be obtained. When the last eigenfunction is substituted into the three compatible equations, one force-equilibrium equation, one governing equation for each attaching point of the concentrated element, and the boundary equations for the two ends of the beam, a matrix equation of the form [B]{C} = {0} is obtained. The solution of |B| = 0 (where ${\mid}{\cdot}{\mid}$ denotes a determinant) will give the "exact" natural frequencies of the "constrained" beam (carrying any number of point masses or/and concentrated springs) and the substitution of each corresponding values of {C} into the associated eigenfunction for each attaching point will determine the corresponding mode shapes. Since the order of [B] is 4n + 4, where n is the total number of point masses and concentrated springs, the "explicit" mathematical expression for the existing approach becomes lengthily intractable if n > 2. The "numerical assembly method"(NAM) introduced in this paper aims at improving the last drawback of the existing approach. The "exact"solutions in this paper refer to the numerical results obtained from the "continuum" models for the classical analytical approaches rather than from the "discretized" ones for the conventional finite element methods.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.4
• /
• pp.83-90
• /
• 1987

The applications of reduced integration technique, addition of nonconforming modes, and their coupling to the Mindlin plate elements to improve their basic behavior are reviewed and the establishment of a series of new plate elements by combined use of these schemes are presented in this paper. The element formulation is based upon quadratic Mindlin plate concept. The results obtained by new elements converged to the exact solutions very rapidly as the mesh is refined and showed reliable solutions even for severely distorted meshes. The new elements have the requisite numbers of zero eigenvalues associated with rigid body modes to avoid the spurious zero energy modes. These elements are shown to be applicable to the wide range of plate problems, giving a high accuracy for both thick and thin plates.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1053-1067
• /
• 2009

The objective of this study is to develop an efficient and accurate quadratic finite element model based on Streamline Upwind/Petrov Galerkin (SU/PG) scheme for analyzing and predicting two dimensional flow features in complex natural rivers. For a development of model, quadratic tin, quadrilateral and mixed elements as well as linear tin, quadrilateral and mixed elements were used in the model. Also, this model was developed through reinforcement of Gauss Quadrature which was necessary to integral of governing equation. Several tests for bottom-rising channel and U-type channel were performed for the purpose of validation and verification of the developed model. Such results showed that solutions of second order elements are better accurate and improved than those of linear elements. Results obtained by the developed model and RMA-2 model are compared, and the results for the developed model were better accurate than those of RMA-2 model. In the future if the developed model is applied in natural rivers, it can provide better accurate results than those of existing model.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.2
• /
• pp.123-130
• /
• 2015

This paper presents a new formulation for transient simulations of microwave propagation in heterogeneous unbounded domains. In particular, perfectly-matched-layers(PMLs) are introduced to allow for wave absorption at artificial boundaries used to truncate the infinite extent of the physical domains. The development of the electromagnetic PML targets the application to engineering mechanics problems such as structural health monitoring and inverse medium problems. To formulate the PML for plane electromagnetic waves, a complex coordinate transformation is introduced to Maxwell's equations in the frequency-domain. Then the PML-endowed partial differential equations(PDEs) for transient electromagnetic waves are recovered by the application of the inverse Fourier transform to the frequency-domain equations. A mixed finite element method is employed to solve the time-domain PDEs for electric and magnetic fields in the PML-truncated domain. Numerical results are presented for plane microwaves propagating through concrete structures, and the accuracy of solutions is investigated by a series of error analyses.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.75-82
• /
• 2004

With an aim at eliminating the numerical dispersion error arising from the numerical simulation of stress wave propagation, numerical dispersion characteristics of the wave equation based one-dimensional finite element model are analyzed and some dispersion control scheme are proposed in this paper The dispersion analyses are carried out for two types of mass matrix, namely the consistent and the lumped mass matrices. Based on the finding of the analyses, dispersion correction techniques are developed for both the implicit and explicit schemes. For the implicit scheme, either the weighting factor for the spatial derivatives of each time level or the lumping coefficient for mass matrix is adjusted to minimize the numerical dispersion. In the case of the explicit scheme an artificial dispersion term is introduced in the governing equation. The validity of the dispersion correction techniques proposed in this study is demonstrated by comparing the numerical solutions obtained using the Present techniques with the analytical ones.