• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1999년도 봄 학술발표회 논문집
• /
• pp.169-176
• /
• 1999

An improved formulation for spatial stability md free vibration of thin-walled curved beams with variable curvature and non-symmetric cross sections are presented based on the displacement field considering the second order terms of finite semitangential rotations. By introducing Vlasov's assumptions, the total potential energy is derived from the principle of linearized virtual work for a continuum. In this formulation, all displacement parameters and the warping function are defined at the centroid axis so that the coupled terms of bending and torsion are added to the elastic strain energy. Also, the potential energy due to initial stress resultants is consistently derived corresponding to the semitangential rotation and moment. The cubic Hermitian polynomials are utilized as shape functions for development of the curved thin-walled beam element having eight degrees of freedom. In order to illustrate the accuracy and practical usefulness of this study, . numerical solutions for free vibration of arches are presented and compared with resells of other researchers and solutions analyzed by the ABAQUS's shell element.

• International Journal of Ocean System Engineering
• /
• 제1권2호
• /
• pp.76-88
• /
• 2011

The present study deals with the hydroelastic vibration analysis of structures in contact with fluid via coupled fluid-structure interaction (FSI) embedded with a finite element method (FEM) such that a structure displacement formulation is coupled with a fluid pressure-displacement formulation. For the preliminary study and validation of FEM based coupled FSI analysis, hydroelastic vibration characteristics of a rectangular plate in contact with fluid are first compared with the elastic vibration in terms of boundary condition and mode frequency. Numerical results from coupled FSI analysis have been shown to be rational and accurate, compared to energy method based theoretical solutions and experimental results. The effect of free surface on the vibration mode is numerically studied by changing the submerged depth of a rectangular plate. As a practical application, the hull structural vibration of 4,000 twenty-foot equivalent units (TEU) container ship is considered. Hydroelastic results of the ship hull structure are compared with those obtained from the elastic condition.

• Structural Engineering and Mechanics
• /
• 제6권8호
• /
• pp.921-937
• /
• 1998

A new four-node degenerated shell element with drilling degrees of freedom (DOF) is proposed. Allman-type displacement approximation is incorporated into the formulation of degenerated shell elements. The approximation improves in-plane performance and eliminates singularities of system matrices resulted from DOF deficiency. Transverse shear locking is circumvented by introducing assumed covariant shear strains. Two kinds of penalty energy are considered in the formulation for the purpose of suppressing spurious modes and representing true drilling rotations. The proposed element can be applied to almost all kinds of shell problems including composite laminated shell structures and folded shell structures. Numerical examples show that the element is of good accuracy and of reasonably fast convergence rate.

• 전산구조공학
• /
• 제10권2호
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.

• Structural Engineering and Mechanics
• /
• 제41권2호
• /
• pp.157-181
• /
• 2012

The paper presents a review of the application of the newly proposed mixed finite element model for seismic simulation of different types of composite frame structures. To evaluate the performance of the element, a comparison with displacement-based and force-based models is conducted. The study revealed that the mixed model is superior to the others in terms of both speed of convergence and numerical stability, and is therefore considered the most practical approach for modeling of composite structures. In this model, the element is derived using independent force and displacement shape functions. The nonlinear response of the frame element is based on the section discretization into fibers with uniaxial material models. The interfacial behavior is modeled using an inelastic interface element. Numerical examples to clarify the advantages of the model are presented for the following structural applications: anchored reinforcing bar problems, composite steel-concrete girders with deformable shear connectors, beam on elastic foundation elements, R/C girders strengthened with FRP sheets, R/C beam-columns with bond-slip, and prestressed concrete girders. These studies confirmed that the model represents a major advancement over existing elements in simulating the inelastic behavior of composite structures.

• Structural Engineering and Mechanics
• /
• 제60권1호
• /
• pp.1-19
• /
• 2016

In this study, the supersonic panel flutter of doubly curved composite sandwich panels with variable thickness is considered under aerothermoelastic loading. Considering different radii of curvatures of the face sheets in this paper, the thickness of the core is a function of plane coordinates (x,y), which is unique. For the first time in the current model, the continuity conditions of the transverse shear stress, transverse normal stress and transverse normal stress gradient at the layer interfaces, as well as the conditions of zero transverse shear stresses on the upper and lower surfaces of the sandwich panel are satisfied. The formulation is based on an enhanced higher order sandwich panel theory and the vertical displacement component of the face sheets is assumed as a quadratic one, while a cubic pattern is used for the in-plane displacement components of the face sheets and the all displacement components of the core. The formulation is based on the von $K{\acute{a}}rm{\acute{a}}n$ nonlinear approximation, the one-dimensional Fourier equation of the heat conduction along the thickness direction, and the first-order piston theory. The equations of motion and boundary conditions are derived using the Hamilton principle and the results are validated by the latest results published in the literature.

• Smart Structures and Systems
• /
• 제14권4호
• /
• pp.699-718
• /
• 2014

The indirect displacement estimation using acceleration and strain (IDEAS) method is extended to various types of beam structures beyond the previous validation on the prismatic or near-prismatic beams. By fusing different types of responses, the IDEAS method is able to estimate displacements containing pseudo-static components with high frequency noise to be significantly reduced. However, the concerns to the IDEAS method come from possible disagreement of the assumed sinusoidal mode shapes to the actual mode shapes, which allows the IDEAS method to be valid only for simply-supported prismatic beams and limits its applicability to real world problems. In this paper, the extension of the IDEAS method to the general types of beams is investigated by the mathematical formulation of the modal mapping matrix only for the monitored substructure, so-called monitoring span. The formulation particularly considers continuous and wide beams to extend the IDEAS method to general beam structures that reflect many real bridges. Numerical simulations using four types of beams with various irregularities are presented to show the effectiveness and accuracy of the IDEAS method in estimating displacements.

• Journal of Mechanical Science and Technology
• /
• 제19권spc1호
• /
• pp.283-291
• /
• 2005

Accurate seismic analyses of large deformable moving structures are still unsolved problems in the field of earthquake engineering. In order to analyze these problems, the nonlinear finite element method formulated by the absolute nodal coordinate approach is noticed. Because, this formulation has several advantages over the standard procedures on mass matrix, elastic forces and damping forces in the case of large displacement problems. But, it has not been fully studied to build frame structure models by using beam elements in the absolute nodal coordinate formulation. In this paper, we propose the connecting method of the beam elements formulated by the absolute nodal coordinate. The coordinate transformation matrix of this element is introduced into the frame structure. This beam element has the characteristic that the mass matrix and bending stiffiness matrix are constant even if in the case of large displacement problems, and this characteristic is being kept after the transformation. In order to verify the proposed method, we show the numerical simulation results of frame structures for a vibration problem and a large displacement problem.

• Computers and Concrete
• /
• 제2권6호
• /
• pp.499-516
• /
• 2005

A composite model is used to represent the heterogeneity of plain concrete consisting of coarse aggregates, mortar matrix and the mortar-aggregate interface. The composite elements of plain concrete are modeled using triangular finite element units which have six interface nodes along the sides. Fracture is captured through a constitutive single branch softening-fracture law at the interface nodes, which bounds the elastic domain inside each triangular unit. The inelastic displacement at an interface node represents the crack opening or sliding displacement and is conjugate to the internodal force. The path-dependent softening behaviour is developed within a quasi-prescribed displacement control formulation. The crack profile is restricted to the interface boundaries of the defined mesh. No re-meshing is carried out. Solutions to the rate formulation are obtained using a mathematical programming procedure in the form of a linear complementary problem. An event by event solution strategy is adopted to eliminate solutions with simultaneous formation of softening zones in symmetric problems. The composite plain concrete model is compared to experimental results for the tensile crack growth in a Brazilian test and three-point bending tests on different sized specimens. The model is also used to simulate wedge-type shear-compression failure directly under the loading platen of a Brazilian test.

• Structural Engineering and Mechanics
• /
• 제76권4호
• /
• pp.435-449
• /
• 2020

This work presents the formulation of the isogeometric collocation method to solve the strong form equation of a unified and integrated approach of Reissner Mindlin plate theory (UI-RM). In this plate theory model, the total displacement is expressed in terms of bending and shear displacements. Rotations, curvatures, and shear strains are represented as the first, the second, and the third derivatives of the bending displacement, respectively. The proposed formulation is free from shear locking in the Kirchhoff limit and is equally applicable to thin and thick plates. The displacement field is approximated using the B-splines functions, and the strong form equation of the fourth-order is solved using the collocation approach. The convergence properties and accuracy are demonstrated with square plate problems of thin and thick plates with different boundary conditions. Two approaches are used for convergence tests, e.g., increasing the polynomial degree (NELT = 1×1 with p = 4, 5, 6, 7) and increasing the number of element (NELT = 1×1, 2×2, 3×3, 4×4 with p = 4) with the number of control variable (NCV) is used as a comparable equivalent variable. Compared with DKMQ element of a 64×64 mesh as the reference for all L/h, the problem analysis with isogeometric collocation on UI-RM plate theory exhibits satisfying results.