• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.3
• /
• pp.267-275
• /
• 2009

This paper presents the p-convergent coupling element on the basis of the ESSE(equivalent single layer shell element) and the PLLE(partial-linear layerwise element) to analyze laminated composite plates. The ESSE is formulated by the degenerated shell theory, on the other hand, the assumption of the PLLE is piecewise linear variation of the in-plane displacement and a constant value of lateral displacement across the thickness. The proposed finite element model is based on p-convergence approach. The integrals of Legendre polynomials and Gauss-Lobatto technique are chosen to interpolate displacement fields and to implement numerical quadrature, respectively. This study has been focused on the verification of p-convergent element. For this purpose, various finite element multiple models associated with the combination of ESSE and PLLE elements are tested to show numerical stability. The simple examples such as a cantilever beam subjected vertical load and a plate with tension are adopted to evaluate the performance of proposed element.

• 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.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• Structural Engineering and Mechanics
• /
• v.86 no.3
• /
• pp.291-309
• /
• 2023

This paper addresses the finite element modeling of functionally graded porous (FGP) beams for free vibration and buckling behaviour cases. The formulated finite element is based on simple and efficient higher order shear deformation theory. The key feature of this formulation is that it deals with Euler-Bernoulli beam theory with only three unknowns without requiring any shear correction factor. In fact, the presented two-noded beam element has three degrees of freedom per node, and the discrete model guarantees the interelement continuity by using both C⁰ and C¹ continuities for the displacement field and its first derivative shape functions, respectively. The weak form of the governing equations is obtained from the Hamilton principle of FGP beams to generate the elementary stiffness, geometric, and mass matrices. By deploying the isoparametric coordinate system, the derived elementary matrices are computed using the Gauss quadrature rule. To overcome the shear-locking phenomenon, the reduced integration technique is used for the shear strain energy. Furthermore, the effect of porosity distribution patterns on the free vibration and buckling behaviours of porous functionally graded beams in various parameters is investigated. The obtained results extend and improve those predicted previously by alternative existing theories, in which significant parameters such as material distribution, geometrical configuration, boundary conditions, and porosity distributions are considered and discussed in detailed numerical comparisons. Determining the impacts of these parameters on natural frequencies and critical buckling loads play an essential role in the manufacturing process of such materials and their related mechanical modeling in aerospace, nuclear, civil, and other structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.1A
• /
• pp.11-18
• /
• 2012

This paper deals with the numerical determination of the stress intensity factors of cracked aluminum plates under the mixed mode of $K_I$ and $K_{II}$ in glass-epoxy fiber reinforced composites. For the stress intensity factors, two different models are reviewed such as VCCT and two-step extension method. The p-convergent partial layerwise model is adopted to determine the fracture parameters in terms of energy release rates and stress intensity factors. The p-convergent approach is based on the concept of subparametric element. In assumed displacement field, strain-displacement relations and 3-D constitutive equations of a layer are obtained by combination of 2-D and 1-D higher-order shape functions. In the elements, Lobatto shape functions and Gauss-Lobatto technique are employed to interpolate displacement fields and to implement numerical quadrature. Using the models and techniques considered, effects of composite laminate configuration according to inclined angles and adhesive properties on the performance of bonded composite patch are investigated. In addition to these, the out-of-plane bending effect has been investigated across the thickness of patch repaired laminate plates due to the change of neutral axis. The present model provides accuracy and simplicity in terms of stress intensity factors, stress distribution, number of degrees of freedom, and energy release rates as compared with previous works in literatures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.109-122
• /
• 1983

An accurate thick beam element (TB4) which includes the effects of the shear deformation and rotatory inertia has been degenerated from the three dimensional continuum by employing the Timoshenko beam assumptions. The proposed TB4 element has four nodes and two degrees of freedom at each node, totally eight degrees of freedom. The transverse deflection W and plane rotation ${\theta}$ with the cubic interpolation functions are selected as nodal variables. The element characteristics are formulated by discretizing the beam equations of motion, using the Galerkin weighted residual method, and are numerically integrated by the reduced shear integration technique, using the three-point Gauss quadrature with the various shear coefficients. Several numerical examples are analyzed to demonstrate the accuracy and the monotonic convergence behavior of the proposed TB4 beam element. The result indicates that the TB4 element shows the more excellent performance and the monotonic convergence behavior than the other existing Timoshenko beam type elements for the whole range of the beam aspect ratios, in both static and free vibration analyses.

• Journal of Radiation Protection and Research
• /
• v.17 no.1
• /
• pp.43-56
• /
• 1992

As one of the methods to ameliorate the ray effects which are the nature of anomalous computational effects due to the discretization of the angular variable in discrete ordinates approximations, a computational program, named TWODET (TWO dimensional Discrete Element Transport), has developed in 2 dimensional cartesian coordinates system using the discrete elements method, in which the discrete angle quadratures are steered by the spatially dependent angular fluxes. The results of the TWODET calculation with K-2, L-3 discrete angular quadratures, in the problem of a centrally located, isotropically emitting flat source in an absorbing square, are shown to be more accurate than that of the DOT 4.3 calculation with S-10 full symmetry angular quadratures, in remedy of the ray effect at the edge flux distributions of the square. But the computing time of the TWODET is about 4 times more than that of the DOT 4.3. In the problem of vacuum boundaries just outside of the source region in an absorbing square, the results of the TWODET calculation are shown severely anomalous ray effects, due to the sudden discontinuity between the source and the vacuum, like as the results of the DOT 4.3 calculation. In the probelm of an external source in an absorbing square in which a highly absorbing medium is added, the results of the TWODET calculation with K-3, L-4 show a good ones like as, somewhat more than, that of the DOT 4.3 calculation with S-10.