• Structural Engineering and Mechanics
• /
• v.48 no.1
• /
• pp.41-56
• /
• 2013

Initiation and growth of delamination is a great concern of designers of composite structures. Interface elements with de-cohesive constitutive law in the content of continuum damage mechanics can be used to predict initiation and growth of delamination in single and mixed mode conditions. In this paper, an interface element based on the cohesive zone method has been developed to simulate delaminatoin growth of post-buckled laminate under fatigue loading. The model was programmed as the user element and user material by the "User Programmable Features" in ANSYS finite element software. The interface element is a three-dimensional 20 node brick with small thickness. Because of mixed-mode condition of stress field at the delamination-front of post-buckled laminates, a mixed-mode bilinear constitutive law has been used as user material in this model. The constitutive law of interface element has been verified by modelling of a single element. A composite laminate with initial delamination under quasi-static compressive Loading available from literature has been remodeled with the present approach. Moreover, it will be shown that, the closer the delamination to the free surface of laminate, the slower the delamination growth under compressive fatigue loading. The effects of laminate configuration on delamination growth are also investigated.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.27-42
• /
• 2018

Behavior of soil is usually described with continuum type of failure models such as Mohr-Coulomb or Drucker-Prager model. The main advantage of these models is in a relatively simple and efficient way of predicting the main tendencies and overall behavior of soil in failure analysis of interest for engineering practice. However, the main shortcoming of these models is that they are not able to capture post-peak behavior of soil nor the corresponding failure modes under extreme loading. In this paper we will significantly improve on this state-of-the-art. In particular, we propose the use of a discrete beam lattice model to provide a sharp prediction of inelastic response and failure mechanisms in coupled soil-foundation systems. In the discrete beam lattice model used in this paper, soil is meshed with one-dimensional Timoshenko beam finite elements with embedded strong discontinuities in axial and transverse direction capable of representing crack propagation in mode I and mode II. Mode I relates to crack opening, and mode II relates to crack sliding. To take into account material heterogeneities, we determine fracture limits for each Timoshenko beam with Gaussian random distribution. We compare the results obtained using the discrete beam lattice model against those obtained using the modified three-surface elasto-plastic cap model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.3
• /
• pp.53-61
• /
• 1984

A nonlinear formulation of a beam finite element(NB6) on the total Lagrangian mode for the geometrically nonlinear analysis of two-dimensional elastic framed structures is presented. The NB6 beam element has been degenerated from the three-dimensional continuum by introducing the deep beam assumptions and consists of three reference nodes and three relative nodes. The element characteristics are derived by discretizing the beam equations of motion using the Galerkin weighted residual method and are reduced-integrated repeatedly for each loading step by the Newton-Raphson iteration techpique. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed nonlinear NB6 beam element.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2002.02a
• /
• pp.274-284
• /
• 2002

Deep drawing process for rectangular drawn section is different with that for axisymmetric circular one. Therefore deep drawing process for rectangular drawn section requires several intermediate steps to generate the final configuration without any significant defect. In this study, finite element analysis for multi-stage deep drawing process for high precision rectangular cases is carried out especially for an extreme aspect ratio. The analysis is performed using rigid-plastic finite element method with an explicit time integration scheme of the commercial program, LS-DYNA3D. The sheet blank is modeled using eight-node continuum brick elements. The results of analysis show that the irregular contact condition between blank and die affects the occurrence of failure, and the difference of aspect ratio in the drawn section leads to non-uniform metal flow, which may cause failure. A series of experiments for multi-stage deep drawing process for the rectangular cases are conducted, and the deformation configuration and the thickness distribution of the drawn rectangular cases are investigated by comparing with the results of the numerical analysis. The numerical analysis with an explicit time integration scheme shows good agreement with the experimental observation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.427-432
• /
• 2007

In general, the response of bulk material is independent of its size when it comes to considering classical elasticity theory. Because the surface to bulk ratio of the large solids is very small, the influence of surface can be negligible. But the surface effect plays important role as the surface to bulk ratio becomes larger, that is, the contribution of the surface effect must be considered in nano-size elements such as thin film or beam structure. Molecular dynamics computation has been a conventional way to analyze these ultra-thin structures but this method is limited to simulate on the order of $10^6-10^8$ atoms for a few nanoseconds, and besides, very time consuming. Analysis of structures in submicro to micro range(thin-film, wire etc.) is difficult with classical molecular dynamics due to the restriction of computing resources and time. Therefore, in this paper, the continuum-based method is considered to simulate the overall physical and mechanical properties of the structures in nano-scale, especially, for the thin-film.

• Steel and Composite Structures
• /
• v.2 no.1
• /
• pp.21-34
• /
• 2002

A model of the process of local buckling in tubular steel structural elements is presented. It is assumed that this degrading phenomenon can be lumped at plastic hinges. The model is therefore based on the concept of plastic hinge combined with the methods of continuum damage mechanics. The state of this new kind of inelastic hinge is characterized by two internal variables: the plastic rotation and the damage. The model is valid if only one local buckling appears in the plastic hinge region; for instance, in the case of framed structures subjected to monotonic loadings. Based on this damage model, a new finite element that can describe the development of local buckling is proposed. The element is the assemblage of an elastic beamcolumn and two inelastic hinges at its ends. The stiffness matrix, that depends on the level of damage, the yielding function and the damage evolution law of the two hinges define the new finite element. In order to verify model and finite element, several small-scale frames were tested in laboratory under monotonic loading. A lateral load at the top of the frame was applied in a stroke-controlled mode until local buckling appears and develops in several locations of the frame and its ultimate capacity was reached. These tests were simulated with the new finite element and comparison between model and test is presented and discussed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.12 s.255
• /
• pp.1518-1525
• /
• 2006

Finite element analyses of structures made of the fiber reinforced composites require an adequate method to characterize the high anisotropic behavior induced by one or several layers of fiber cords with different spatial orientation embedded in a rubber matrix. This paper newly proposes a continuum based rebar element considering change of the orientation of the fiber during deformation of the composite. The mechanical behavior of the embedded fiber is modeled using two-node bar elements in order to consider the relative deformation and spatial orientation of the embedded fiber. For improvement of the analysis accuracy, the load-displacement curve of fiber is applied to the stiffness matrix of fiber. A finite element program is constructed based on the total Lagrangian formulation considering both geometric and material nonlinearity. Finite element analyses of the tensile test are carried out in order to evaluate the validity of the proposed method. Analysis results obtained with the proposed method provides realistic representation of the fiber reinforced rubber composite compared to results of other two models by the Halpin-Tsai equation and a rebar element in ABAQUS/Standard.

• Computers and Concrete
• /
• v.7 no.6
• /
• pp.483-496
• /
• 2010

This work deals with the determination of crack openings in 2D reinforced concrete structures using the Finite Element Method with a smeared rotating crack model or an embedded crack model. In the smeared crack model, the strong discontinuity associated with the crack is spread throughout the finite element. As is well known, the continuity of the displacement field assumed for these models is incompatible with the actual discontinuity. However, this type of model has been used extensively due to the relative computational simplicity it provides by treating cracks in a continuum framework, as well as the reportedly good predictions of reinforced concrete members' structural behavior. On the other hand, by enriching the displacement field within each finite element crossed by the crack path, the embedded crack model is able to describe the effects of actual discontinuities (cracks). This paper presents a comparative study of the abilities of these 2D models in predicting the mechanical behavior of reinforced concrete structures. Structural responses are compared with experimental results from the literature, including crack patterns, crack openings and rebar stresses predicted by both models.

• Advances in Computational Design
• /
• v.8 no.2
• /
• pp.97-112
• /
• 2023

The Finite Element (FE) modeling of Reinforced Concrete (RC) under seismic loading has a sensitive impact in terms of getting good contribution compared to experimental results. Several idealized model types for simulating the nonlinear response have been developed based on the plasticity distribution alone the model. The Continuum Models are the most used category of modeling, to understand the seismic behavior of structural elements in terms of their components, cracking patterns, hysteretic response, and failure mechanisms. However, the material modeling, contact and nonlinear analysis strategy are highly complex due to the joint operation of concrete and steel. This paper presents a numerical simulation of a chosen RC column under monotonic and cyclic loading using the FE Abaqus, to assessthe hysteretic response and failure mechanisms in the RC columns, where the perfect bonding option is used for the contact between concrete and steel. While results of the numerical study under cyclic loading compared to experimental tests might be unsuccessful due to the lack of bond-slip modeling. The monotonic loading shows a good estimation of the envelope response and deformation components. In addition, this work further demonstrates the advantage and efficiency of the damage distributions since the obtained damage distributions fit the expected results.