• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.113-123
• /
• 1990

The paper is concerned with the elasto-plastic and geometrically nonlinear analysis of shell structures using an improved degenerated shell element. In the formulation of the element stiffness, the combined use of three different techniques was made. They are; 1) an enhanced interpolation of transverse shear strains in the natural coordinate system to overcome the shear locking problem ; 2) the reduced integration technique in in-plane strains to avoid the membrane locking behavior ; and 3) selective addition of the nonconforming displacement modes to improve the element performances. This element is free of serious shear/membrane locking problems and undesirable compatible/commutable spurious kinematic deformation modes. In the formulation for plastic deformation, the concept of a layered element model is used and the material is assumed von Mises yield criterion. An incremental total Lagrangian formulation is presented which allows the calculation of arbitrarily large displacements and rotations. The resulting non-linear equilibrium equations are solved by the Netwon-Raphson method combined with load or displacement increment. The versatility and accuracy of this improved degenerated shell element are demonstrated by solving several numerical examples.

• Journal of the Korea Concrete Institute
• /
• v.28 no.3
• /
• pp.279-288
• /
• 2016

The present study simplified the moment-curvature relationship to straightforwardly determine the flexural behavior of reinforced concrete (RC) columns. For the idealized column section, moments and neutral axis depths at different stages(first flexural crack, yielding of tensile reinforcing bar, maximum strength, and 80% of the maximum strength at the descending branch) were derived on the basis of the equilibrium condition of forces and compatibility condition. Concrete strains at the extreme compression fiber beyond the maximum strength were determined using the stress-strain relationship of confined concrete, proposed by Kim et al. The lateral load-displacement curves converted from the simplified moment-curvature relationship of columns are well consistent with test results obtained from column specimens under various parameters. The moments and the corresponding neutral axis depth at different stages were formulated as a function of longitudinal reinforcement and transverse reinforcement indices and/or applied axial load index. Overall, curvature ductility of columns was significantly affected by the axial load level as well as concrete compressive strength and the amount of longitudinal and transverse reinforcing bars.

• Journal of Korean Ophthalmic Optics Society
• /
• v.6 no.2
• /
• pp.85-97
• /
• 2001

A mathematical model was proposed to analyze the damped motion of contact lens which is initially displaced from the equilibrium position. The model incorporates the differential equations and their numerical solution program, based on the formulations of restoring force arising from the capillary action in the tear-film layer between the lens and cornea. The model predicts the capillary action induced surface tension, time dependence of displacement of lens when it is released from the equilibrium position. It seems that the motion of lens is similar to the typical over-damped oscillation caused by the large viscous friction in the liquid layer between the cornea and lens. The effect of variables such as base curves, lens diameters and thickness of tear film layer were illustrated by the computer simulation of the derived program. The time required for the lens to return to the original position increases as the liquid layer thickness increases and it decreases as the diameter of lens increases. With the certain value of base curve the time interval is found to be minimum. The free vibrations of lenses were also simulated varying the parameters such as base curve, diameter, layer thickness. The resonant frequencies are inversely proportional to the liquid layer thickness and it increases as the lens diameter increases. The resonant frequency of lens has a maximum when the diameter is of certain value. If the external impulse or force of the same frequency as the natural frequency of contact lens acted on the cornea in vivo it may cause an excessive movement and thus it might cause the distortion 10 the lens or be pulled off the eye.

• Structural Engineering and Mechanics
• /
• v.35 no.6
• /
• pp.677-697
• /
• 2010

This paper focuses on geometrically non-linear static analysis of a simply supported beam made of hyperelastic material subjected to a non-follower transversal uniformly distributed load. As it is known, the line of action of follower forces is affected by the deformation of the elastic system on which they act and therefore such forces are non-conservative. The material of the beam is assumed as isotropic and hyperelastic. Two types of simply supported beams are considered which have the following boundary conditions: 1) There is a pin at left end and a roller at right end of the beam (pinned-rolled beam). 2) Both ends of the beam are supported by pins (pinned-pinned beam). In this study, finite element model of the beam is constructed by using total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element. The considered highly non-linear problem is solved by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In order to use the solution procedures of Newton-Raphson type, there is need to linearized equilibrium equations, which can be achieved through the linearization of the principle of virtual work in its continuum form. In the study, the effect of the large deflections and rotations on the displacements and the normal stress and the shear stress distributions through the thickness of the beam is investigated in detail. It is known that in the failure analysis, the most important quantities are the principal normal stresses and the maximum shear stress. Therefore these stresses are investigated in detail. The convergence studies are performed for various numbers of finite elements. The effects of the geometric non-linearity and pinned-pinned and pinned-rolled support conditions on the displacements and on the stresses are investigated. By using a twelve-node quadratic element, the free boundary conditions are satisfied and very good stress diagrams are obtained. Also, some of the results of the total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element are compared with the results of SAP2000 packet program. Numerical results show that geometrical nonlinearity plays very important role in the static responses of the beam.

• Korean Journal of Agricultural Science
• /
• v.40 no.4
• /
• pp.425-434
• /
• 2013

Nowadays, the korean peninsula has endured weather changes toward to the subtropical climate. Also demographical structure has changed into multi-cultural society in which many people from subtropical areas have immigrated into Korea. Therefore, consumption and production of subtropical vegetables become important. For the analysis, we choose eight important subtropical crops. EDM (Equilibrium Displacement Model) with many parameters and elasticities is used for the forecast of consumption and required cultivation area. The simulation focuses on the changes of the number of foreign workers and immigrated women in Korea to predict the quantity of consumption and required area in Korea. The results show that we need additional land area about 581~1,065 ha for the cropping subtropical vegetables in Korea. Finally, these required area can be provided by the cities and counties in coast area in Gyeongsangnam-Do, Jeonranam-Do and Jeju-Do. Climate change will be continued in the future. Together with climate change, the change of demographical structure into multi-culture may increase consumption and production of subtropical vegetables. Forecasting of increased consumption and required cultivation area for subtropical vegetables is significant.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.2 s.45
• /
• pp.123-131
• /
• 2000

In this study, both experimental and analytical study for behavior of the existing composite steel box girder bridges, constructed along with the procedure of continuous placing slab, are conducted to establish the validity of the proposed model. The layer approach is adopted to determine the equilibrium condition in a section to consider the different material properties and concrete cracking across the sectional depth, and the beam element stiffness is constructed on the basis of the assumed displacement field formulation and the 3-points Gaussian Integration. In addition, the effects of creep and shrinkage of concrete for time-dependent behavior of the bridge are taken into consideration. Finally, both analytical and experimental results are compared.

• Steel and Composite Structures
• /
• v.39 no.3
• /
• pp.337-352
• /
• 2021

The study presented experimental and numerical investigation on the seismic performance of steel reinforced concrete (SRC) L-shaped column- reinforced concrete (RC) beam joints. Various parameters described as steel configuration form, axial compressive ratio, loading angle, and the existence of slab were examined through 4 planar joints and 7 spatial joints. The characteristics of the load-displacement response included the bearing capacity, ductility, story drift ratio, energy-dissipating capacity, and stiffness degradation were analyzed. The results showed that shear failure and flexural failure in the beam tip were observed for planar joints and spatial joint, respectively. And RC joint with slab failed with the plastic hinge in the slab and bottom of the beam. The results indicated that hysteretic curves of spatial joints with solid-web steel were plumper than those with hollow-web specimens. The capacity of planar joints was higher than that of space joints, while the opposite was true for energy-dissipation capacity and ductility. The high compression ratio contributed to the increase in capacity and initial stiffness of the joint. The elastic and elastic-plastic story deformation capacity of L-shaped column frame joints satisfied the code requirement. A design formula of joint shear resistance based on the superposition theory and equilibrium plasticity truss model was proposed for engineering application.

• Geomechanics and Engineering
• /
• v.18 no.3
• /
• pp.225-234
• /
• 2019

A numerical stepwise approach for cavity expansion problem in strain-softening rock or soil mass is investigated, which is compatible with Mohr-Coulomb and generalized Hoek-Brown failure criteria. Based on finite difference method, plastic region is divided into a finite number of concentric rings whose thicknesses are determined internally to satisfy the equilibrium and compatibility equations, the material parameters of the rock or soil mass are assumed to be the same in each ring. For the strain-softening behavior, the strength parameters are assumed to be a linear function of deviatoric plastic strain (${\gamma}p^*$) for each ring. Increments of stress and strain for each ring are calculated with the finite difference method. Assumptions of large-strain for soil mass and small-strain for rock mass are adopted, respectively. A new numerical stepwise approach for limited pressure and plastic radius are obtained. Comparisons are conducted to validate the correctness of the proposed approach with Vesic's solution (1972). The results show that the perfectly elasto-plastic model may underestimate the displacement and stresses in cavity expansion than strain-softening coefficient considered. The results of limit expansion pressure based on the generalised H-B failure criterion are less than those obtained based on the M-C failure criterion.

• Steel and Composite Structures
• /
• v.43 no.5
• /
• pp.639-660
• /
• 2022

The bending and buckling behaviours of FG-GRNC laminated sandwich plates are investigated by using novel five-variables quasi 3D higher order shear deformation plate theory by considering the modified continuum nonlocal strain gradient theory. To calculate the effective Young's modulus of the GRNC sandwich plate along the thickness direction, and Poisson's ratio and mass density, the modified Halpin-Tsai model and the rule of the mixture are employed. Based on a new field of displacement, governing equilibrium equations of the GRNC sandwich plate are solved using a developed approach of Galerkin method. A detailed parametric analysis is carried out to highlight the influences of length scale and material scale parameters, GPLs distribution pattern, the weight fraction of GPLs, geometry and size of GPLs, the geometry of the sandwich plate and the total number of layers on the stresses, deformation and critical buckling loads. Some details are studied exclusively for the first time, such as stresses and the nonlocality effect.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.1
• /
• pp.9-16
• /
• 2018

When designing a reinforced concrete member using nonlinear finite element analysis results, the bending moment at the critical section should be calculated. In this paper, a bending moment calculation method using the results of reinforced concrete finite element analysis(FEA) using continuum elements is presented and the optimum element size according to the order of the displacement function of the finite element is proposed. The bending moments calculated by integrating the stresses from the FEA are compared with the bending moments calculated using the static equilibrium conditions. In the method of integrating the stress, both the stress due to the reinforcing bar and the stress of the concrete are considered. In addition, various factors affecting the accuracy of the stresses calculated by the FEA were analyzed and the influence of the displacement function and the element size was verified. If the purpose of the analysis is to roughly observe the behavior of the members, it is appropriate to use the first order displacement function and the element size should be about 25% of the section height of the analytical model. When the bending moment of a member with high accuracy is required, it is suggested that the secondary displacement function be used and the element size be 12.5%.