• Structural Engineering and Mechanics
• /
• 제7권2호
• /
• pp.203-216
• /
• 1999

A previously proposed finite element formulation method is refined and modified to generate a new type of elements. The method is based on selecting a set of general solution modes for element formulation. The constant strain modes and higher order modes are selected and the formulation method is designed to ensure that the element will pass the basic single element test, which in turn ensures the passage of the basic patch test. If the element is to pass the higher order patch test also, the element stiffness matrix is in general asymmetric. The element stiffness matrix depends only on a nodal displacement matrix and a nodal force matrix. A symmetric stiffness matrix can be obtained by either modifying the nodal displacement matrix or the nodal force matrix. It is shown that both modifications lead to the same new element, which is demonstrated through numerical examples to be more robust than an assumed stress hybrid element in plane stress application. The method of formulation can also be used to arrive at the conforming displacement and hybrid stress formulations. The convergence of the latter two is explained from the point of view of the proposed method.

• Structural Engineering and Mechanics
• /
• 제51권4호
• /
• pp.601-625
• /
• 2014

Based on continuum mechanics and the principle of virtual displacements, incremental total Lagrangian formulation (T.L.) and incremental updated Lagrangian formulation (U.L.) were presented. Both T.L. and U.L. considered the large displacement stiffness matrix, which was modified to be symmetrical matrix. According to the incremental updated Lagrangian formulation, small strain, large displacement, finite rotation of three dimensional Timoshenko fiber beam element tangent stiffness matrix was developed. Considering large displacement and finite rotation, a new type of tangent stiffness matrix of the beam element was developed. According to the basic assumption of plane section, the displacement field of an arbitrary fiber was presented in terms of nodal displacement of centroid of cross-area. In addition, shear deformation effect was taken account. Furthermore, a nonlinear finite element method program has been developed and several examples were tested to demonstrate the accuracy and generality of the three dimensional beam element.

• Structural Engineering and Mechanics
• /
• 제42권1호
• /
• pp.39-53
• /
• 2012

This paper presents an alternative way to derive the exact element stiffness matrix for a beam on Winkler foundation and the fixed-end force vector due to a linearly distributed load. The element flexibility matrix is derived first and forms the core of the exact element stiffness matrix. The governing differential compatibility of the problem is derived using the virtual force principle and solved to obtain the exact moment interpolation functions. The matrix virtual force equation is employed to obtain the exact element flexibility matrix using the exact moment interpolation functions. The so-called "natural" element stiffness matrix is obtained by inverting the exact element flexibility matrix. Two numerical examples are used to verify the accuracy and the efficiency of the natural beam element on Winkler foundation.

• 전산구조공학
• /
• 제9권3호
• /
• pp.169-177
• /
• 1996

구조물에 있어서 위상학적 최적화 문제는 최적화를 구하는 과정에서 구조체가 변화함으로 인한 어려움 때문에 최적화 분야에서 가장 어려운 문제로 간주되어 왔다. 종래의 방법으로는 일반적으로 구조요소 사이즈가 영으로 접근할 때 강성 매트릭스의 singularity를 발생시킴으로써 최적의 해를 얻지 못하고 도중에 계산이 종료되어 버린다. 본 연구에 있어서는 이러한 문제점들을 해결하기 위한 비선형 프로그래밍 formulation을 제안하는 것을 목적으로 한다. 이 formulation의 주된 특성은 요소 사이즈가 영이 되는 것을 허용한다. 평형방정식을 등제약조건으로 간주함으로써 강성 매트릭스의 singularity를 피할 수 있다. 이 formulation을 하중을 받는 구조물에 있어서 응력과 변위의 제약조건하에서 중량을 최소화할때의 유한요소의 두께를 구하는 디자인 문제에 적용하여, 이 formulation이 위상학적 최적화에 있어서의 효과를 입증하였다.

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

The focus of this study is on the problem of the design of structure of undetermined topology. This problem has been regarded as being the most challenging of structural optimization problems, because of the difficulty of allowing topology to change. Conventional approaches break down when element sizes approach to zero, due to stiffness matrix singularity. In this study, a novel nonlinear Programming formulation of the topology Problem is developed and examined. Its main feature is the ability to account for topology variation through zero element sizes. Stiffness matrix singularity is avoided by embedding the equilibrium equations as equality constraints in the optimization problem. Although the formulation is general, two dimensional plane elasticity examples are presented. The design problem is to find minimum weight of a plane structure of fixed geometry but variable topology, subject to constraints on stress and displacement. Variables are thicknesses of finite elements, and are permitted to assume zero sizes. The examples demonstrate that the formulation is effective for finding at least a locally minimal weight.

• Earthquakes and Structures
• /
• 제3권3_4호
• /
• pp.271-295
• /
• 2012

This paper presents the solution scheme of using the continuous formulation of 1-D linear member for the dynamic analysis of structures consisting of axially loaded members. The context describes specific applications of such scheme to the verification of experimental data obtained from field test of bridges carried out by a microwave interferometer system and velocimeters. Attention is focused on analysis outlines that may be applicable to in-situ assessment for cable-stayed bridges. The derivation of the dynamic stiffness matrix of a prismatic member with distributed properties is briefly reviewed. A back calculation formula using frequencies of two arbitrary modes of vibration is next proposed to compute the tension force in cables. Derivation of the proposed formula is based on the formulation of an axially loaded flexural member. The applications of the formulation and the proposed formula are illustrated with a series of realistic examples.

• 전산구조공학
• /
• 제9권3호
• /
• pp.179-186
• /
• 1996

본 연구에서는 구조체의 위상학적 최적화를 위한 비선형 formulation(NLP)가 개발, 검토되었다. 이 NLP는 multiple-loading하에서 임의의 오브젝티브 함수, 응력, 변위 제약조건들을 쉽게 다룰 수가 있다. 또한 이 NLP는 해석과 최적화 디자인을 동시에 실시함으로써 요소 사이즈가 영으로 접근함에 따른 강성 매트릭스의 singularity를 피할 수 있다. 즉, 평형 방정식을 등제약조건으로 치환함으로써 강성 매트릭스 그 자체나 그의 역매트릭스를 구할 필요도 없어진다. 이 NLP는 multiple-loading conditon하에서 테스트되었으며, 이를 통해 이 NLP가 다양한 제약조건하에서 강력하게 작용함이 입증되었다.

• Structural Engineering and Mechanics
• /
• 제21권5호
• /
• pp.539-551
• /
• 2005

Finite elements based on isoparametric formulation are known to suffer spurious stiffness properties and corresponding stress oscillations, even when care is taken to ensure that completeness and continuity requirements are enforced. This occurs frequently when the physics of the problem requires multiple strain components to be defined. This kind of error, commonly known as locking, can be circumvented by using reduced integration techniques to evaluate the element stiffness matrices instead of the full integration that is mathematically prescribed. However, the reduced integration technique itself can have a further drawback - rank deficiency, which physically implies that spurious energy modes (e.g., hourglass modes) are introduced because of reduced integration. Such instability in an existing stiffness matrix is generally detected by means of an eigenvalue test. In this paper we show that a knowledge of the dimension of the solution space spanned by the column vectors of the strain-displacement matrix can be used to identify the instabilities arising in an element due to reduced/selective integration techniques a priori, without having to complete the element stiffness matrix formulation and then test for zero eigenvalues.

• Structural Engineering and Mechanics
• /
• 제13권4호
• /
• pp.437-453
• /
• 2002

In this paper, dynamic stiffness and flexibility for circular membranes are analytically derived using an efficient mixed-part dual boundary element method (BEM). We employ three approaches, the complex-valued BEM, the real-part and imaginary-part BEM, to determine the dynamic stiffness and flexibility. In the analytical formulation, the continuous system for a circular membrane is transformed into a discrete system with a circulant matrix. Based on the properties of the circulant, the analytical solutions for the dynamic stiffness and flexibility are derived. In deriving the stiffness and flexibility, the spurious resonance is cancelled out. Numerical aspects are discussed and emphasized. The problem of numerical instability due to division by zero is avoided by choosing additional constraints from the information of real and imaginary parts in the dual formulation. For the overdetermined system, the least squares method is considered to determine the dynamic stiffness and flexibility. A general purpose program has been developed to test several examples including circular and square cases.

• Structural Engineering and Mechanics
• /
• 제35권6호
• /
• pp.717-733
• /
• 2010

The dynamic stiffness matrix is formulated for an axially loaded slender double-beam element in which both beams are homogeneous, prismatic and of the same length by directly solving the governing differential equations of motion of the double-beam element. The Bernoulli-Euler beam theory is used to define the dynamic behaviors of the beams and the effects of the mass of springs and axial force are taken into account in the formulation. The dynamic stiffness method is used for calculation of the exact natural frequencies and mode shapes of the double-beam systems. Numerical results are given for a particular example of axially loaded double-beam system under a variety of boundary conditions, and the exact numerical solutions are shown for the natural frequencies and normal mode shapes. The effects of the axial force and boundary conditions are extensively discussed.