• Structural Engineering and Mechanics
• /
• 제8권6호
• /
• pp.623-637
• /
• 1999

This paper presents a four-noded quadrilateral $C^0$ strain plate element for the analysis of thick laminated composite plates. The element formulation is based on: 1) the third-order shear deformation theory; 2) assumed strain element formulation; and 3) interrelated edge displacements and rotations along element boundaries. Unlike the existing displacement-type composite plate elements based on the third-order theory, which rely on the $C^1$-continuity formulation, the present plate element is of $C^0$-continuity, and its element stiffness matrix is evaluated explicitly. Because of the third-order expansion of the in-plane displacements through the thickness, the resulting theory and hence elements do not need shear correction factors. The explicit element stiffness matrix makes the present element more computationally efficient than the composite plate elements using numerical integration for the analysis of thick layered composite plates.

• Structural Engineering and Mechanics
• /
• 제24권6호
• /
• pp.647-664
• /
• 2006

In this study, an improved finite element formulation with a scheme of solution for the large deflection analysis of inextensible prismatic and nonprismatic slender beams is developed. For this purpose, a three-noded Lagrangian beam-element with two dependent degrees of freedom per node (i.e., the vertical displacement, y, and the actual slope, $dy/ds=sin{\theta}$, where s is the curved coordinate along the deflected beam) is used to derive the element stiffness matrix. The element stiffness matrix in the global xy-coordinate system is achieved by means of coordinate transformation of a highly nonlinear ($6{\times}6$) element matrix in the local sy-coordinate. Because of bending with large curvature, highly nonlinear expressions are developed within the global stiffness matrix. To achieve the solution after specifying the proper loading and boundary conditions, an iterative quasi-linearization technique with successive corrections are employed considering these nonlinear expressions to remain constant during all iterations of the solution. In order to verify the validity and the accuracy of this study, the vertical and the horizontal displacements of prismatic and nonprismatic beams subjected to various cases of loading and boundary conditions are evaluated and compared with analytic solutions and numerical results by available references and the results by ADINA, and excellent agreements were achieved. The main advantage of the present technique is that the solution is directly obtained, i.e., non-incremental approach, using few iterations (3 to 6 iterations) and without the need to split the stiffness matrix into elastic and geometric matrices.

• Structural Engineering and Mechanics
• /
• 제1권1호
• /
• pp.47-60
• /
• 1993

The eigen-equation of a wave traveling over repetitive structure is derived directly form the stiffness matrix formulation, in a form which can be used for the case of the cross stiffness submatrix $K_{ab}$ being singular. The weighted adjoint symplectic orthonormality relation is proved first. Then the general method of solution is derived, which can be used either to find all the eigensolutions, or to find the main eigensolutions for large scale problems.

• 한국CDE학회논문집
• /
• 제15권6호
• /
• pp.411-417
• /
• 2010

In this paper, the boom of a floating crane is modeled as a 3-dimensional elastic beam in order to analyze the dynamic response of the crane and its cargo. The boom is divided into more than two elements based on finite element formulation, and deformation of each element is expressed in terms of shape matrix and nodal coordinates. The equations of motion for the elastic boom consist of a mass matrix, a stiffness matrix, and a quadratic velocity vector that contains the gyroscopic and Coriolis forces. The size and complicity of the matrices increase in proportion with the number of elements. Therefore, it is not possible to derive the equations of motion explicitly for different number of elements. To overcome this difficulty, matrices for one 3-dimensional element are expressed with elementary sub-matrices. In particular, the quadratic velocity vector is derived as a product of a shape matrix and a 3-dimensional rotation matrix. By using the derived matrices, the equations of motion for the multi-element boom are automatically constructed. To verify the implementation of the elastic boom based on finite element formulation, we simulated a simple vibration of the elastic boom and compared the average deformation with the analytic solution. Finally, heave motion of the floating crane and surge motion of the cargo are presented as application examples of the elastic boom.

• Advances in aircraft and spacecraft science
• /
• 제1권3호
• /
• pp.291-310
• /
• 2014

An advanced model for the linear flutter analysis is introduced in this paper. Higher-order beam structural models are developed by using the Carrera Unified Formulation, which allows for the straightforward implementation of arbitrarily rich displacement fields without the need of a-priori kinematic assumptions. The strong form of the principle of virtual displacements is used to obtain the equations of motion and the natural boundary conditions for beams in free vibration. An exact dynamic stiffness matrix is then developed by relating the amplitudes of harmonically varying loads to those of the responses. The resulting dynamic stiffness matrix is used with particular reference to the Wittrick-Williams algorithm to carry out free vibration analyses. According to the doublet lattice method, the natural mode shapes are subsequently used as generalized motions for the generation of the unsteady aerodynamic generalized forces. Finally, the g-method is used to conduct flutter analyses of both isotropic and laminated composite lifting surfaces. The obtained results perfectly match those from 1D and 2D finite elements and those from experimental analyses. It can be stated that refined beam models are compulsory to deal with the flutter analysis of wing models whereas classical and lower-order models (up to the second-order) are not able to detect those flutter conditions that are characterized by bending-torsion couplings.

• Structural Engineering and Mechanics
• /
• 제15권2호
• /
• pp.199-223
• /
• 2003

The formulation of a non-linear shear deformable shell element is presented for the solution of stability problems of stiffened plates and shells. The formulation of the geometrical stiffness presented here is exactly defined on the midsurface and is efficient for analyzing stability problems of thick plates and shells by incorporating bending moment and transverse shear resultant force. As a result of the explicit integration of the tangent stiffness matrix, this formulation is computationally very efficient in incremental nonlinear analysis. The element is free of both membrane and shear locking behaviour by using the assumed strain method such that the element performs very well in the thin shells. By using six degrees of freedom per node, the present element can model stiffened plate and shell structures. The formulation includes large displacement effects and elasto-plastic material behaviour. The material is assumed to be isotropic and elasto-plastic obeying Von Mises's yield condition and its associated flow rules. The results showed good agreement with references and computational efficiency.

• Earthquakes and Structures
• /
• 제3권3_4호
• /
• pp.383-411
• /
• 2012

The objective of the present paper is to review and implement the most recent developments in the Spectral Element Method (SEM), as well as improve aspects of its implementation in the study of wave propagation by numerical simulation in elastic unbounded domains. The classical formulation of the method is reviewed, and the construction of the mass matrix, stiffness matrix and the external force vector is expressed in terms of matrix operations that are familiar to earthquake engineers. To account for the radiation condition at the external boundaries of the domain, a new absorbing boundary condition, based on the Perfectly Matched Layer (PML) is proposed and implemented. The new formulation, referred to as the Multi-Axial Perfectly Matched Layer (M-PML), results from generalizing the classical Perfectly Matched Layer to a medium in which damping profiles are specified in more than one direction.

• Steel and Composite Structures
• /
• 제24권1호
• /
• pp.79-91
• /
• 2017

Free vibrations of steel-concrete composite beams are analyzed by using the dynamic stiffness approach. The coupled equations of motion of the composite beams are derived with help of the Hamilton's principle. The effects of the shear deformation and rotary inertia of the two beams as well as the transverse and axial deformations of the stud connectors are included in the formulation. The dynamic stiffness matrix is developed on the basis of the exact general solutions of the homogeneous governing differential equations of the composite beams. The use of the dynamic stiffness method to determine the natural frequencies and mode shapes of a particular steel-concrete composite beam with various boundary conditions is demonstrated. The accuracy and effectiveness of the present model and formulation are validated by comparison of the present results with the available solutions in literature.

• 한국공간구조학회지
• /
• 제2권1호
• /
• pp.1-7
• /
• 2002

• 대한기계학회논문집
• /
• 제18권4호
• /
• pp.887-902
• /
• 1994

Within the framework of anisotropic thermoelasticity, the problem of an interlaminar crack in a laminated fiber-reinforced composite subjected to a uniform heat flow is investigated. Under a state of generalized plane deformation, dissimilar anisotropic half-spaces with different fiber orientations are considered to be bound together by a matrix interlayer containing the crack. The interlayer models the matrix-rich interlaminar region of the fibrous composite laminate. Based on the flexibility/stiffness matrix approach, formulation of the current crack problem results in having to solve two sets of singular integral equations for temperature and thermal stress analyses. Numerical results are obtained, illustrating the parametric effects of laminate stacking sequence, relative crack size, crack location, crack surface partial insulation, and fiber volume fraction on the values of mixed mode thermal stress intensity factors.