• Steel and Composite Structures
• /
• 제23권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
• /
• 제3권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.

• 한국전산구조공학회논문집
• /
• 제31권5호
• /
• pp.275-282
• /
• 2018

등기하 해석법을 이용한 고유치 해석은 유한요소를 이용한 결과보다 고차 모드에서 더 정확한 결과를 주는 것으로 알려져 있다. 이는 유한요소법이 차수에 상관없이 요소 간에 $C^0$ 연속성을 보이는 것과 다르게 등기하 해석법은 p차 요소에 대해서 $C^{p-1}$의 연속성을 보장하기 때문이다. 본 논문에서는 이러한 장점을 이용하여 등기하 해석법을 이용하여 모드 기반의 축소 모델을 구성하고 동적 거동 해석을 수행하였다. 축소 모델 구성을 위해 Craig-Bampton(CB) 기법을 적용하였다. 수치 예제를 통해 간단한 봉 요소에 대해 등기하 해석법과 유한요소 해석법을 적용하여 요소의 차수에 따른 고유치 해석 결과를 비교분석하였다. 등기하 해석법에 중첩 노트를 허용하여 요소 간 연속성을 조절하고, 요소 간 연속성이 줄어듦에 따라 고차 모드에서의 수치 오차가 커짐을 확인하였다. 동적 거동 해석을 위한 축소 모델에 높은 차수의 외력이 주어지는 경우 요소간 연속성이 높은 등기하해석법을 사용하면, 해의 정확도를 높일 수 있다.

• Structural Engineering and Mechanics
• /
• 제14권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.

• 한국전산구조공학회논문집
• /
• 제27권3호
• /
• pp.183-189
• /
• 2014

본 연구에서는 다면체 요소의 개발을 위하여 Wachspress 좌표계와 이동최소자승 근사를 기반으로하는 형상함수와 수치적분 방법을 제시하고 있다. 사면체 요소를 사면체 영역으로 분할하여 형상함수가 구성이 되고 이 영역을 사용한 일관성있는 수치적분이 수행되게 된다. 다면체 요소 면에서 Wachspress 좌표계를 사용하고 요소 내부에서 라플라스 방정식을 적용하여 이동최소자승 근사의 가중함수를 정의하게 된다. 본 연구에서 개발되는 다면체 요소의 형상함수와 수치적분 방법은 일반적인 유한요소와 유사한 특성을 갖게 되는데 수치 예제를 통하여 유효성을 보여주었다.

• 대한조선학회지
• /
• 제22권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
• /
• 제6권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
• /
• 제90권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
• /
• 제90권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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• 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.