• Structural Engineering and Mechanics
• /
• v.4 no.6
• /
• pp.655-678
• /
• 1996

The paper reports on alternative hybrid/mixed formulations being developed by the Structural Analysis Research Group of Institute Superior T$\acute{e}$cnico. These formulations open the scope and increase the power of the finite element method by allowing different fields to be independently approximated, within certain consistency criteria, and by enhancing the use of a wide range of approximation functions. They have been applied to the analysis of 2-D problems, laminar structures and solids, using different constitutive relations, both in quasi-static and dynamic regimes. The fundamental properties of the formulations are identified and assessed and their performance is illustrated using simple, linear applications.

• Structural Engineering and Mechanics
• /
• v.10 no.2
• /
• pp.125-138
• /
• 2000

The paper analyzes the problem of torsion in an adhesive non-tubular bonded single-lap joint. The joint considered consists of two thin rectangular section beams bonded together along a side surface. Assuming the materials involved to be governed by linear elastic laws, equilibrium and compatibility equations were used to arrive at an integro-differential relation whose solution makes it possible to determine torsional moment section by section in the bonded joint between the two beams. This is then used to determine the predominant stress and strain field at the beam-adhesive interface (stress field along the direction perpendicular to the interface plane, equivalent to the applied torsional moment and the corresponding strain field) and the joint's elastic strain (absolute and relative rotations of the bonded beam cross sections). All the relations presented were obtained in closed form. Results obtained theoretically are compared with those given by a three dimensional finite element numerical model. Theoretical and numerical analysis agree satisfactorily.

• Structural Engineering and Mechanics
• /
• v.81 no.3
• /
• pp.335-348
• /
• 2022

It is well known that after finding the displacement in the structural mechanics, strain and stress can be obtained in the straight-forward process. The main purpose of this paper is to unify the displacement functions for solving the solid body. By performing mathematical operations, three sets of these key relationships are found in this paper. All of them are written in the Cartesian Coordinates and in terms of a simple function. Both analytical and numerical approaches are utilized to validate the correctness of the presented formulations. Since all required conditions for the bodies with self-equilibrated loadings are satisfied accurately, the authors' relations can solve these kinds of problems. This fact is studied in-depth by solving some numerical examples. It is found that a very simple function can be used for each formulation instead of ten different and complex displacement potentials defined by previous studies.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.1
• /
• pp.95-110
• /
• 2010

A numerical model to simulate bond-slip behavior of composite beam bridges is introduced in this paper. Assuming a linear bond stress-slip relation along the interface between the slab and girder, the slip behavior is implemented into a finite element formulation. Adopting the introduced model, the slip behavior can be taken account even in a beam element which is composed of both end nodes only. Governing equation of the slip behavior, based on the linear partial interaction theory, can be determined from the force equilibrium and a constant curvature distribution across the section of a composite beam. Since the governing equation for the slip behavior requires the moment values at both end nodes, the piecewise linear distribution of the constant bending moment in an element is assumed. Analysis results by the model are compared with numerical results and experimental values, and load-displacement relations of composite beams were then evaluated to verify the validity of the proposed model.

• Structural Engineering and Mechanics
• /
• v.72 no.5
• /
• pp.617-629
• /
• 2019

Based on the finite element method of traditional straight Euler-Bernoulli beams and the coupled relations between linear displacement and angular displacement of a pre-twisted Euler-Bernoulli beam, the shape functions and stiffness matrix are deduced. Firstly, the stiffness of pre-twisted Euler-Bernoulli beam is developed based on the traditional straight Euler-Bernoulli beam. Then, a new finite element model is proposed based on the displacement general solution of a pre-twisted Euler-Bernoulli beam. Finally, comparison analyses are made among the proposed Euler-Bernoulli model, the new numerical model based on displacement general solution and the ANSYS solution by Beam188 element based on infinite approach. The results show that developed numerical models are available for the pre-twisted Euler-Bernoulli beam, and which provide more accurate finite element model for the numerical analysis. The effects of pre-twisted angle and flexural stiffness ratio on the mechanical property are investigated.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.40 no.3
• /
• pp.103-112
• /
• 1998

In this study, database schema which is capable of maintaining various design-knowledge and applicable to the structural designs has been presented. Futhermore, by utilizing the database system, RECODSS (REinforced COncrete Design Support System) was developed as well. Knowledge-base that were formulated from the document is constructed by applying the concrete standard specifications of Korea from 1989. The basic drawing was prepared by $AutoCAD^TM$. As an inference in the preliminary design, the certainty function and the interpolation function were presented. The certainty function is made by assuming linear relations between each data while referencing $MYCIN^TM$ (expert systems). Nevertheless, the interpolation function was made by utilizing the linear Lagrange interpolation. For the search of the design-knowledge, fast search algorithm and full-text search algorithm were used, however, the design engine was programmed with C language. RECODSS was applied for designing a slab and a beam member. And the results were compared with the existing results. In conclusion, RECODSS was able to efficiently manage various design knowledge and to introduce stable design results.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Structural Engineering and Mechanics
• /
• v.82 no.4
• /
• pp.477-490
• /
• 2022

This article investigates the nonlinear behavior of two-directional functionally graded materials (TDFGM) doubly curved panels with porosities for the first time. An improved and effectual approach is established based on the improved first-order shear deformation shell theory (IFSDST) and von Karman's type nonlinearity. The IFSDST considers the effects of shear deformation without the need for a shear correction factor. The composition of TDFGM constitutes four different materials, and the modified power-law function is employed to vary the material properties continuously in both thickness and longitudinal directions. A nonlinear finite element method in conjunction with Hamilton's principle is used to obtain the governing equations. Then, the direct iterative method is incorporated to accomplish the numerical results using the frequency-amplitude, nonlinear central deflection relations. Finally, the influence of volume fraction grading indices, porosity distributions, porosity volume, curvature ratio, thickness ratio, and aspect ratio provides a thorough insight into the linear and nonlinear responses of the porous curved panels. Meanwhile, this study emphasizes the influence of the volume fraction gradation profiles in conjunction with the various material and geometrical parameters on the linear frequency, nonlinear frequency, and deflection of the TDFGM porous shells. The numerical analysis reveals that the frequencies and nonlinear deformations can be significantly regulated by changing the volume fraction gradation profiles in a specified direction with an appropriate combination of materials. Hence, TDFGM panels can overcome the drawbacks of the functionally graded materials with a gradation of properties in a single direction.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.22 no.6
• /
• pp.323-332
• /
• 2018

Recent two moderate earthquakes (2016 $M_w=5.4$ Gyeongju and 2017 $M_w=5.5$ Pohang) in Korea provided the unique chance of developing a set of relations to estimate instrumental seismic intensity in Korea by augmenting the time-history data from MMI seismic intensity regions above V to the insufficient data previously accumulated from the MMI regions limited up to IV. The MMI intensity regions of V and VI was identified by delineating the epicentral distance from the reference intensity statistics in distance derived by using the integrated MMI data obtained by combining the intensity survey results of KMA (Korea Meteorological Administration) and 'DYFI (Did You Feel It)' MMIs of USGS. The time-histories of the seismic stations from the MMI intensity regions above V were then preprocessed by applying the previously developed site-correction filters to be converted to a site-equivalent condition in a manner consistent with the previous study. The average values of the ground-motion parameters for the three ground motion parameters of PGA, PGV and BSPGA (Bracketed Summation of PGA per second for 30 seconds) were calculated for the MMI=V and VI and used to generate the dataset of the average values of the ground-motion parameters for the individual MMIs from I to VI. Based on this dataset, the linear regression analysis resulted in the following relations with proposed valid ranges of MMI. $MMI=2.36{\times}log_{10}(PGA(gal))+1.44$ ($I{\leq}MMI$$MMI=2.44{\times}log_{10}(PGV(kine))+4.86$ ($I{\leq}MMI$$MMI=2.59{\times}log_{10}(BSPGA(gal{\cdot}sec))-1.02$ ($I{\leq}MMI$