• Steel and Composite Structures
• /
• v.23 no.1
• /
• pp.41-51
• /
• 2017

A $C^0$ FE model developed based on an efficient higher order zigzag theory is used for hygrothermal analysis of laminated composite plates. The $C^0$ FE model satisfies the inter-laminar shear stress continuity at the interfaces and zero transverse shear stress conditions at plate top and bottom. In this model the first derivatives of transverse displacement have been treated as independent variables to circumvent the problem of $C^1$ continuity associated with the above plate theory. In the present theory the above mentioned $C^0$ continuity of the present element is compensated in the stiffness matrix formulation by using penalty parameter approach. In order to avoid stress oscillations observed in the displacement based finite element, the stress field derived from temperature/moisture fields (initial strains) must be consistent with total strain field. Special steps are introduced by field consistent approach (e.g., sampling at gauss points) to compensate this problem. A nine noded $C^0$ continuous isoparametric element is used in the proposed FE model. Comparison of present numerical results with other existing solutions shows that the proposed FE model is efficient, accurate and free of locking.

• International Journal of Aerospace System Engineering
• /
• v.3 no.1
• /
• pp.10-16
• /
• 2016

A trigonometric Layerwise higher order shear deformation theory (TLHSDT) is developed and implemented for free vibration and buckling analysis of laminated composite and sandwich plates by analytical and finite element formulation. The present model assumes parabolic variation of out-plane stresses through the depth of the plate and also accomplish the zero transverse shear stresses over the surface of the plate. Thus a need of shear correction factor is obviated. The present zigzag model able to meet the transverse shear stress continuity and zigzag form of in-plane displacement continuity at the plate interfaces. Hence, botheration of shear correction coefficient is neglected. In the case of analytical method, the governing differential equation and boundary conditions are obtained from the principle of virtual work. For the finite element formulation, an efficient eight noded $C^0$ continuous isoparametric serendipity element is established and employed to examine the dynamic analysis. Like FSDT, the considered mathematical model possesses similar number of variables and which decides the present models computationally more effective. Several numerical predictions are carried out and results are compared with those of other existing numerical approaches.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.5
• /
• pp.275-282
• /
• 2018

Using isogeometric analysis(IGA) gives more accurate results for higher order mode in eigenvalue problem than using the finite element method(FEM). This is because the FEM has $C^0$ continuity between elements, whereas IGA guarantee $C^{P-1}$ between elements for p-th order basis functions. In this paper, a mode based reduced model is constructed by using IGA and dynamic behavior analysis is performed using this advantage. Craig-Bampton(CB) method is applied to construct the reduced model. Several numerical examples were performed to compare the eigenvalue analysis results for various order of element basis function by applying the IGA and FEM to simple rod analysis. We have confirmed that numerical error increases in the higher order mode as the continuity between elements decreases in the IGA by allowing internal knots multiplicity. The accuracy of the solution can be improved by using the IGA with high inter-element continuity when high-frequency external force acts on the reduced model for dynamic behavior analysis.

• Structural Engineering and Mechanics
• /
• v.14 no.3
• /
• pp.245-262
• /
• 2002

The free vibration analysis of stiffened laminated composite plates has been performed using the layered (zigzag) finite element method based on the first order shear deformation theory. The layers of the laminated plate is modeled using nine-node isoparametric degenerated flat shell element. The stiffeners are modeled as three-node isoparametric beam elements based on Timoshenko beam theory. Bilinear in-plane displacement constraints are used to maintain the inter-layer continuity. A special lumping technique is used in deriving the lumped mass matrices. The natural frequencies are extracted using the subspace iteration method. Numerical results are presented for unstiffened laminated plates, stiffened isotropic plates, stiffened symmetric angle-ply laminates, stiffened skew-symmetric angle-ply laminates and stiffened skew-symmetric cross-ply laminates. The effects of fiber orientations (ply angles), number of layers, stiffener depths and degrees of orthotropy are examined.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.183-189
• /
• 2014

In this paper, a polyhedral element is presented to solve three-dimensional problems by developing shape functions based on Wachspress coordinates and moving least square approximation. A subdivision of polyhedrons into tetrahedral domains is performed for the construction of shape functions of polyhedral elements, and numerical integration of the weak form is carried out consistently over the tetrahedral domains. The weight functions for moving least square approximation are defined by solving Laplace equation with boundary values based on Wachspress coordinates on polyhedral element faces. Polyhedral elements presented in this paper have similar properties to conventional finite element regarding the continuity, the completeness, the node-element connectivity and the inter-element compatibility. Numerical examples show the effectiveness of the present method for solving three-dimensional problems using polyhedral elements.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.3
• /
• pp.9-18
• /
• 1985

The numerical calculation for solving boundary-value problem related to potential flows with a free surface is carried out by application of the localized finite element method. Only forced motion of 2-D body in infinitely deep fluid is considered, although this schemes is equally applicable to any first order time-harmonic problems of similar nature. The infinite domain of the fluid is separated into the inner flow field and the outer flow field with common inter-surface boundary. The finite element method is applied to obtain the solution in the inner flow field and the Green functions are utilized to represent the solution in the outer flow field. At the inter-surface boundary, the continuity of the value of potential and the normal derivative of the potential(i.e. matching condition) is conserved. The present method has better computational efficiency than the previous LFEM and the integral equation method of Frank. This enhanced computational efficiency is presumably due to the fact that the present method gives a symmetric coefficient matrix and requires less computational time in calculating the influence coefficient matrix of Green function than the integral equation method. And the irregular frequency desen't exist because the uniqueness of the solution is assured by the such that the exact free surface condition is satisfied on the boundary of the localized finite element region(i.e. inner region). As an example of the above method, the hydrodynamic forces for the circular cylinder and the rectangular cylinders are calculated. In the computed results, the small number of singularity distribution segments($3{\sim}6$) give good result relative to Ursell's and Vugts'.

• Smart Structures and Systems
• /
• v.6 no.8
• /
• pp.921-938
• /
• 2010

Most of studies on control of beams containing piezoelectric sensors and actuators have been based on linear quadratic regulator (LQR) with state feedback or output feedback law. The aim of this study is to develop velocity-acceleration feedback law in the optimal control of smart piezoelectric beams. A new controller which is an optimal control system with velocity-acceleration feedback is presented. In finite element modeling of the beam, the variation of mechanical displacement through the thickness is modeled by a sinus model that ensures inter-laminar continuity of shear stress at the layer interfaces as well as the boundary conditions on the upper and lower surfaces of the beam. In addition to mechanical degrees of freedom, one electric potential degree of freedom is considered for each piezoelectric element layer. The efficiency of this control strategy is evaluated by applying to an aluminum cantilever beam under different loading conditions. Numerical simulations show that this new control scheme is almost as efficient as an optimal control system with state feedback. However, inclusion of the acceleration in the control algorithm increases practical value of a system due to easier and more accurate measurement of accelerations.

• Structural Engineering and Mechanics
• /
• v.90 no.3
• /
• pp.233-252
• /
• 2024

This research explores a new finite element model for the free vibration analysis of bi-directional functionally graded (BDFG) beams. The model is based on an efficient higher-order shear deformation beam theory that incorporates a trigonometric warping function for both transverse shear deformation and stress to guarantee traction-free boundary conditions without the necessity of shear correction factors. The proposed two-node beam element has three degrees of freedom per node, and the inter-element continuity is retained using both C1 and C0 continuities for kinematics variables. In addition, the mechanical properties of the (BDFG) beam vary gradually and smoothly in both the in-plane and out-of-plane beam's directions according to an exponential power-law distribution. The highly elevated performance of the developed model is shown by comparing it to conceptual frameworks and solution procedures. Detailed numerical investigations are also conducted to examine the impact of boundary conditions, the bi-directional gradient indices, and the slenderness ratio on the free vibration response of BDFG beams. The suggested finite element beam model is an excellent potential tool for the design and the mechanical behavior estimation of BDFG structures.

• Structural Engineering and Mechanics
• /
• v.90 no.3
• /
• pp.253-261
• /
• 2024

This research explores a new finite element model for the free vibration analysis of bi-directional functionally graded (BDFG) beams. The model is based on an efficient higher-order shear deformation beam theory that incorporates a trigonometric warping function for both transverse shear deformation and stress to guarantee traction-free boundary conditions without the necessity of shear correction factors. The proposed two-node beam element has three degrees of freedom per node, and the inter-element continuity is retained using both C1 and C0 continuities for kinematics variables. In addition, the mechanical properties of the (BDFG) beam vary gradually and smoothly in both the in-plane and out-of-plane beam's directions according to an exponential power-law distribution. The highly elevated performance of the developed model is shown by comparing it to conceptual frameworks and solution procedures. Detailed numerical investigations are also conducted to examine the impact of boundary conditions, the bi-directional gradient indices, and the slenderness ratio on the free vibration response of BDFG beams. The suggested finite element beam model is an excellent potential tool for the design and the mechanical behavior estimation of BDFG structures.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.577-582
• /
• 2007

A variational formulation for plane elasticity problems is derived based on an isogeometric approach. The isogeometric analysis is an emerging methodology such that the basis functions in analysis domain arc generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Thus. the solution space can be represented in terms of the same functions to represent the geometry. The coefficients of basis functions or the control variables play the role of degrees-of-freedom. Furthermore, due to h-. p-, and k-refinement schemes, the high order geometric features can be described exactly and easily without tedious re-meshing process. The isogeometric sensitivity analysis method enables us to analyze arbitrarily shaped structures without re-meshing. Also, it provides a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling. To obtain precise shape sensitivity, the normal and curvature of boundary should be taken into account in the shape sensitivity expressions. However, in conventional finite element methods, the normal information is inaccurate and the curvature is generally missing due to the use of linear interpolation functions. A continuum-based adjoint sensitivity analysis method using the isogeometric approach is derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of boundary. In isogeometric analysis, however, the geometric properties arc already embedded in the B-spline shape functions and control points. The perturbation of control points in isogeometric analysis automatically results in shape changes. Using the conventional finite clement method, the inter-element continuity of the design space is not guaranteed so that the normal vector and curvature arc not accurate enough. On tile other hand, in isogeometric analysis, these values arc continuous over the whole design space so that accurate shape sensitivity can be obtained. Through numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.