• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.33-45
• /
• 1996

A general stiffener element which includes transverse shear deformation is formulated using the p-version finite element method. Hierarchic C/sup o/-shape functions, derived from Integrals of Legendre polynomials, are used to define the assembled stiffness matrix of the stiffener with respect to the local reference frame is transformed to the plate reference system by applying the appropriate transformation matrices in order to insure compatibility of displacements at the junction of the stiffener and plate. The transformation matrices which account for the orientation and the eccentricity effects of the stiffener with respect to the plate reference axes are used to find local behavior at the junction of the stiffener and the relative contributions of the plate and stiffener to the strength of the composite system. The results obtained by the p-version finite element method are comared with the results in literatures, especially those by the h-version finite element analysis program, MICROFEAP-II.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.1
• /
• pp.17-26
• /
• 1990

A solution strategy for geometric non-linear analysis for the plane frame structures subjected to conservative and non-conservative forces is presented. By making efficient combination of the load incremental method and the displacement incremental method, this strategy can find postbuckling configuration such as snap-through and turning-back phenomena which cannot be easily found by the conventional load and displacement incrementation scheme. In the case of the analysis of the framed structure subjected to circulatory non-conservative forces, the total tangent stiffness matrix becomes unsymmetric and when calculating the incremental load and unbalanced load vector components, the direction change of the non-conservative forces is considered. Several example problems to demonstrate the feasibility of the present strategy, over ranges of deformation that are well beyond those likly to occur in practical framed structures, are given and discussed.

• 한국공간정보시스템학회:학술대회논문집
• /
• 2004.05a
• /
• pp.75-82
• /
• 2004

Many researcher's efforts have made a significant advancement of space frame structure with various portion, and it becomes the most outsanding one of space structures. However, with the characteristics of thin and long term of spacing, the unstable behavior of space structure is shown by initial imperfection, erection procedure or joint, especially space frame structure represents more. This kind of unstable problem could not be set up clearly and there is a huge difference between theory and experiment. Moreover, the discrete structure such as space frame has more complex solution, this it is not easy to derive the formulation of design about space structure. In this space frame structure, the character of rise-span ratio or load mode is represented by the instability of space frame structure with initial imperfection, and snap-through or bifurcation might be the main phenomenon. Therefore, in this study, space frame structure which has a lot of aesthetic effect and profitable for large space covering single layer is dealt. And because that the unstable behavior due to variation of inner force resistance in the elastic range is very important collapse mechanism, I would like to investigate unstable character as a nonlinear behavior with a geometric nonlinear. In order to study the instability. I derive tangent stiffness matrix using finite element method and with displacement incremental method perform nonlinear analysis of unit space structure, star dome and 3-ring star dome considering rise-span $ratio(\mu}$ and load $ratio(R_L)$ for analyzing unstable phenomenon.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.159-168
• /
• 2005

An improved stability design method for beam-columns of plane frames is proposed based on system buckling analysis and second-order elastic analysis. For this, the tangent stiffness matrix of beam-column elements is first derived using stability functions and a procedure for evaluating effective buckling lengths is reviewed using elastic system buckling analysis. And then the second-order analysis procedure is presented considering $P-\Delta$ effects and is compared with the closed-form solution through numerical examples. Design examples showing the validity of the proposed method we presented and their numerical results are compared with those obtained from the conventional stability design methods. Finally some useful conclusions are drawn.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• Journal of the Earthquake Engineering Society of Korea
• /
• v.16 no.4
• /
• pp.37-52
• /
• 2012

In this study, a new 4-node general shell element with 6 DOFs per node is presented. Drilling rotational degrees of freedom are introduced by the variational principle with an independent rotation field. In formulation of the element, substitute transverse shear strain fields are used to avoid shear locking, while four nonconforming modes are applied in the in-plane displacement fields as a remedy for membrane locking. In addition, a direct modification method for nonconforming modes is employed in the numerical implementation of nonconforming modes to represent constant strain states. A 9-points integration rule is adopted for volume integration in the computation of the element stiffness matrix. With the combined use of these techniques, the developed shell element has no spurious zero energy modes, and can represent a constant strain state. Several numerical tests are carried out to evaluate the performance of the new element developed. The test results show that the behavior of the elements is satisfactory.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.24 no.3
• /
• pp.73-83
• /
• 1982

study aims to derive a rational method for the analysis of the farm silo supported on an elastic foundation in which it is assumed that the reaction pressure of the soil at a point is proportional to the deflection at that point. In order to investigate the effects of an elastic foundation on the behaviour of the structures on it, the analysis of the farm silo resting on an elastic foundation was compared with the solution that the ground support may be assumed uniform (which was obtained from part I of this paper). To calculate the deformation of an elastic foundation, Boussinesq's solution which allows an interaction of the various parts of ground was adopted. In this case, the foundation was treated as a superparametric element additionally. In the evaluation of an element stiffness matrix, Gauss quadrature' was used. In above numerical integration, 3-point rule for the farm silo wall and the footing was introduced and 2-point rule for the evaluation of a reaction between the footing and the elastic foundation was adopted. The stresses of a farm silo on an elastic foundation were smaller than those which the distribution of contact pressure between the footing and the soil is assumed uniformly. Since the differences of stresses were remarkable in PS structures than RC structures, it is desirable that designers take into account the effect of an elastic foundation for the case of PS structures. It can be noted that while the effect of an elastic foundation was more conspicuously observed in near of the ground, the value of stresses at far from the soil was little affected by an supported soil.

• Journal of Korean Society of Steel Construction
• /
• v.18 no.2
• /
• pp.225-237
• /
• 2006

Based on the study conducted by Kim et al. (205a, b), an improved stability design method for evaluating the effective buckling lengths of beam-column members is proposed herein, using system elastic/inelastic buckling analysis and second-order elastic analysis. For this purpose, the stress-strain relationship of a column is inversely formulated from the reference load-carrying capacity proposed in design codes, so as to derive the tangent modulus of a column as a function of the slenderness ratio. The tangent stiffness matrix of a beam-column element is formulated using the so-called "stability functions," and elastic/inelastic buckling analysis Effective buckling lengths are then evaluated by extending the basic concept of a single simply-supported column to the individual members as one component of a whole frame structure. Through numerical examples of several structural systems and loading conditions, the possibilities of enhancement in stability design for frame structures are addressed by comparing their numerical results obtained when the present design method is used with those obtained when conventional stability design methods are used.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.11
• /
• pp.1157-1169
• /
• 2008

This paper presents spectral element formulation which approximates Lamb wave propagation by PZT transducers bonded on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by a piezoelectric (PZT) layer rigidly bonded on a base plate. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Euler-Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with the electro-mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are obtained through equations of motions converted into frequency domain. Detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through numerical examples.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.3
• /
• pp.41-46
• /
• 2001

The very large vessels like VLCC and container ship have been built recently and those vessels have smaller structural strength in comparison with the other convectional skips. As a result the fatigue destruction of upper deck occurs a frequently due to the springing phenomenon at the encountering frequencies. In this study, the hydrodynamic loads are calculated by three-dimensional source distribution method with the translating and pulsating Green function. A ship is longitudinally divided into 23 sections and the added mass, damping and hydrodynamic force of each section is calculated. focusing only on the vertical motion. Stiffness matrix is calculated by the Euler beam theory. The calculation is carried out for Esso Osaka.