• Steel and Composite Structures
• /
• v.20 no.1
• /
• pp.91-106
• /
• 2016

Two-dimensional elastoplastic finite element formulation is employed to investigate the load- carrying capacity degradation of reinforced concrete piers wrapped with steel plates due to occurrence of corrosion at the pier base. By comparing with experimental results, the employed finite element analysis method is verified to be accurate. After that, a series of parametric studies are conducted to investigate the effect of corrosion ratio and corrosion mode of steel plates located near the base of in-service pier P2 on load-carrying capacity of the piers. It is observed that the load-carrying capacity of the piers decreases with the increase in corrosion ratio of steel plates. There exists an obvious linear relationship between the load-carrying capacity and the corrosion ratio in the case of even corrosion mode. The degradation of load-carrying capacity resulted from the web's uneven corrosion mode is more serious than that under even corrosion mode, and the former case is more liable to occur than the latter case in actual engineering application. Finally, the failure modes of the piers under different corrosion state are discussed. It is found that the principal tensile strain of concrete and yield range of steel plates are distributed within a wide range in the case of slight corrosion, and they are concentrated on the column base when complete corrosion occurs. The findings obtained from the present study can provide a useful reference for the maintenance and strengthening of the in-service piers.

• Journal of the Korean Geosynthetics Society
• /
• v.8 no.3
• /
• pp.9-16
• /
• 2009

In generally, ground surface strengthening such as using geotextile is needed to secure trafficability of construction equipment. Several calculation methods for bearing capacity for very soft ground improved by geotextile have been proposed already but they are not established yet. In this paper, elastoplastic finite element analysis is carried out to estimate bearing capacity improvement for very soft ground reinforced by geotextile. In conclusion, by reinforcing with geotextile the bearing capacity improvement of very soft ground is increased with the increase of the settlement after yielding. And it is enough that the length of geotextile is about 3~4 times of the loaded width for the effect to the bearing capacity improvement. Finally, the biggest tensile force of geotextile occurs at an edge of the uniform load.

• Computers and Concrete
• /
• v.16 no.6
• /
• pp.933-961
• /
• 2015

Damage by high-speed impact fracture is a dominant mode of failure in several applications of concrete structures. Numerical modelling can play a crucial role in understanding and predicting complex fracture processes. The commonly used mesh-based Finite Element Method has difficulties in accurately modelling the high deformation and disintegration associated with fracture, as this often distorts the mesh. Even with careful re-meshing FEM often fails to handle extreme deformations and results in poor accuracy. Moreover, simulating the mechanism of fragmentation requires detachment of elements along their boundaries, and this needs a fine mesh to allow the natural propagation of damage/cracks. Smoothed Particle Hydrodynamics (SPH) is an alternative particle based (mesh-less) Lagrangian method that is particularly suitable for analysing fracture because of its capability to model large deformation and to track free surfaces generated due to fracturing. Here we demonstrate the capabilities of SPH for predicting brittle fracture by studying a slender concrete structure (column) under the impact of a high-speed projectile. To explore the effect of the projectile material behaviour on the fracture process, the projectile is assumed to be either perfectly-elastic or elastoplastic in two separate cases. The transient stress field and the resulting evolution of damage under impact are investigated. The nature of the collision and the constitutive behaviour are found to considerably affect the fracture process for the structure including the crack propagation rates, and the size and motion of the fragments. The progress of fracture is tracked by measuring the average damage level of the structure and the extent of energy dissipation, which depend strongly on the type of collision. The effect of fracture property (failure strain) of the concrete due to its various compositions is found to have a profound effect on the damage and fragmentation pattern of the structure.

• Earthquakes and Structures
• /
• v.10 no.5
• /
• pp.989-1012
• /
• 2016

The structural dynamic behavior and yield strength considering both ductility and strain rate effects are analyzed in this article. For the single-degree-of-freedom (SDOF) system, the relationship between the relative velocity and the strain rate response is deduced and the strain rate spectrum is presented. The ductility factor can be incorporated into the strain rate spectrum conveniently based on the constant-ductility velocity response spectrum. With the application of strain rate spectrum, it is convenient to consider the ductility and strain rate effects in engineering practice. The modal combination method, i.e., square root of the sum of the squares (SRSS) method, is employed to calculate the maximum strain rate of the elastoplastic multiple-degree-of-freedom (MDOF) system under uniform excitation. Considering the spatially varying ground motions, a new response spectrum method is developed by incorporating the ductility factor and strain rate into the conventional response spectrum method. In order to further analyze the effects of strain rate and ductility on structural dynamic behavior and yield strength, the cantilever beam (one-dimensional) and the triangular element (two-dimensional) are taken as numerical examples to calculate their seismic responses in time domain. Numerical results show that the permanent displacements with and without considering the strain rate effect are significantly different from each other. It is not only necessary in theory but also significant in engineering practice to take the ductility and strain rate effects into consideration.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.621-629
• /
• 1992

This study is concerned with an application of the Boundary Element Method to 2-dimensional elastoplastic stress analysis on the material nonlinearities. The boundary integral formulation adopted an initial stress equation in the inelastic term. In order to determine the initial stress increment, the increment of initial elastic strain energy due to elastic increment in stressstrain curve was used as the convergence criterion during iterative process. For the validity of this procedure, the results of B.E.M. with constant elements and NISA with linear elements where compared on the thin plate with 2 edge v-notches under static tension and the thick cylinder under internal pressure. And this paper compared the results of using unmedical integral with the results of using semi-analytical integral on the plastic domain integral.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.5
• /
• pp.337-343
• /
• 2017

As a shop practice, a strength estimation method for die cast parts is suggested, in which various defects such as pores can be allowed. The equivalent porosity is evaluated by combining the stiffness data from a simple elastic test at the part level during the shop practice and the theoretical stiffness data, which are defect free. A porosity equation is derived from Eshelby's inclusion theory. Then, using the Mori-Tanaka method, the porosity value is used to draw a stress-strain curve for the porous material. In this paper, the Hollomon equation is used to capture the strain hardening effect. This stress-strain curve can be used to estimate the strength of a die cast part with porous defects. An elastoplastic theoretical solution is derived for the three-point bending of a die cast beam by using the plastic hinge method as a reference solution for a part with porous defects.

• Advances in Computational Design
• /
• v.5 no.4
• /
• pp.417-443
• /
• 2020

Tunnels have been an integral part of human civilization. Due to complexity in its design and structure, the stability of underground structures under extreme loading conditions has utmost importance. Increased terrorism and geo-political conflicts have forced the engineers and researchers to study the response of underground structures, especially tunnels under blast loading. The present study has been carried out to seek the response of tunnel structures under blast load using the finite element technique. The tunnel has been considered in quartzite rock of northern India. The Mohr-Coulomb constitutive model has been adopted for the elastoplastic behaviour of rock. The rock model surrounding the tunnel has dimensions of 30 m x 30 m x 35 m. Both unlined and lined (concrete) tunnel has been studied. Concrete Damage Plasticity model has been considered for the concrete lining. Four different parameters (i.e., tunnel diameter, liners thickness, overburden depth and mass of explosive) have been varied to observe the behaviour under different condition. To carry out blast analysis, Coupled-Eulerian-Lagrangian (CEL) modelling has been adopted for modelling of TNT (Trinitrotoluene) and enclosed air. JWL (Jones-Wilkins-Lee) model has been considered for TNT explosive modelling. The paper concludes that deformations in lined tunnels follow a logarithmic pattern while in unlined tunnels an exponential pattern has been observed. The stability of the tunnel has increased with an increase in overburden depth in both lined and unlined tunnels. Furthermore, the tunnel lining thickness also has a significant effect on the stability of the tunnel, but in smaller diameter tunnel, the increase in tunnel lining thickness has not much significance. The deformations in the rock tunnel have been decreased with an increase in the diameter of the tunnel.

• Composites Research
• /
• v.12 no.2
• /
• pp.10-25
• /
• 1999

A new evaluation method for the interfacial properties of fibrous composites based on a fragmentation technique is proposed by using the gradual multi-fiber composite, in which the inter-fiber spacing is gradually changed. The results showed that as the inter-fiber distance increased, the aspect ratio of broken fibers decreased while the interfacial shear strength between the fiber and matrix increased. When the reciprocal of the inter-fiber destance was taken for the above relations, both the aspect ratio and interfacial shear strength showed a saturated value. This means that the gradual multi-fiber composite indicates an upper bound in aspect ratio and an upper bound in interfacial shear strength. It was concluded that this fragmentation test could be a new method for composite evaluation, since reducing a difference between these two bounds is effective for composite strengthening. In addition an elastoplastic finite element analysis was carried out to relate the above results with fiber stress a distribution around fiber breaks. It was proved that the bound obtained in the gradual multi-fiber composite test is closely related to stress concentrations caused by a group of multi-fiber breaks.

• Structural Engineering and Mechanics
• /
• v.49 no.6
• /
• pp.705-726
• /
• 2014

In general, cylindrically curved plates are used in ships and offshore structures such as wind towers, spa structures, fore and aft side shell plating, and bilge circle parts in merchant vessels. In a number of studies, it has been shown that curvature increases the buckling strength of a plate under compressive loading, and the ultimate load-carrying capacity is also expected to increase. In the present paper, a series of elastic and elastoplastic large deflection analyses were performed using the commercial finite element analysis program (MSC.NASTRAN/PATRAN) in order to clarify and examine the fundamental buckling and collapse behaviors of curved plates subjected to combined axial compression and lateral pressure. On the basis of the numerical results, the effects of curvature, the magnitude of the initial deflection, the slenderness ratio, and the aspect ratio on the characteristics of the buckling and collapse behavior of the curved plates are discussed. On the basis of the calculated results, the design formula was developed to predict the buckling and ultimate strengths of curved plates subjected to combined loads in an analytical manner. The buckling strength behaviors were simulated by performing elastic large deflection analyses. The newly developed formulations were applied in order to perform verification analyses for the curved plates by comparing the numerical results, and then, the usefulness of the proposed method was demonstrated.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.