Search

A new approach of multi-scale simulation via coupling of boundary element and finite element is proposed to predict fracture toughness and crack propagation of a single layer graphene sheet. In this simulation a molecular dynamics-based nonlinear beam element is developed for the atomistic model near the crack tip, whereas a special two-dimensional boundary element is employed for the continuum model in the remaining field of the cracked specimen. The material and section properties required in the nonlinear beam element are estimated through the equivalence between the potential energy of molecular dynamics and the elastic strain energy of continuum mechanics. With the estimated properties of beam element, the material properties of boundary element are further estimated by applying loads on the specimen of graphene, which is formed by a hexagonal lattice of carbon atoms. Coupling the nonlinear beam elements by boundary elements with the local-global concept of multi-scale modeling, a vast of computational time can be saved. The accuracy of near tip stresses obtained in our simulation remedies the inaccuracy of linear elastic fracture mechanics, and can be used for the prediction of atomic bond-breaking, which leads to crack propagation. The associated critical load can then be applied in the continuum model for the cracked specimen to predict mode I and mode II fracture toughness of graphene. The obtained values are then verified by the published results measured or predicted by the other methods. By varying the size of cracks and orientation of applied loads, several interesting phenomena have been observed and discussed in this paper. [ABSTRACT FROM AUTHOR]

In the cohesive surface model cohesive tractions are transmitted across a two-dimensional surface, which is embedded in a three-dimensional continuum. The relevant kinematic quantities are the local crack opening displacement and the crack sliding displacement, but there is no kinematic quantity that represents the stretching of the fracture plane. As a consequence, in-plane stresses are absent, and fracture phenomena as splitting cracks in concrete and masonry, or crazing in polymers, which are governed by stress triaxiality, cannot be represented properly. In this paper we extend the cohesive surface model to include in-plane kinematic quantities. Since the full strain tensor is now available, a three-dimensional stress state can be computed in a straightforward manner. The cohesive band model is regarded as a subgrid scale fracture model, which has a small, yet finite thickness at the subgrid scale, but can be considered as having a zero thickness in the discretisation method that is used at the macroscopic scale. The standard cohesive surface formulation is obtained when the cohesive band width goes to zero. In principle, any discretisation method that can capture a discontinuity can be used, but partition-of-unity based finite element methods and isogeometric finite element analysis seem to have an advantage since they can naturally incorporate the continuum mechanics. When using interface finite elements, traction oscillations that can occur prior to the opening of a cohesive crack, persist for the cohesive band model. Example calculations show that Poisson contraction influences the results, since there is a coupling between the crack opening and the in-plane normal strain in the cohesive band. This coupling holds promise for capturing a variety of fracture phenomena, such as delamination buckling and splitting cracks, that are difficult, if not impossible, to describe within a conventional cohesive surface model. [ABSTRACT FROM AUTHOR]

The modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations of complex crack topologies including branching. This drawback can be overcome by a diffusive crack modeling based on the introduction of a crack phase field as proposed in Miehe et al. (Comput Methods Appl Mech Eng 19:2765-2778, ; Int J Numer Meth Eng 83:1273-1311, ), Hofacker and Miehe (Int J Numer Meth Eng, ). In this work, we summarize basic ingredients of a thermodynamically consistent, variational-based model of diffusive crack propagation under quasi-static and dynamic conditions. It is shown that all coupled field equations, in particular the balance of momentum and the gradient-type evolution equation for the crack phase field, follow as the Euler equations of a mixed rate-type variational principle that includes the fracture driving force as the mixed field variable. This principle makes the proposed formulation extremely compact and provides a perfect basis for the finite element implementation. We then introduce a local history field that contains a maximum energetic crack source obtained in the deformation history. It drives the evolution of the crack phase field. This allows for the construction of an extremely robust operator split scheme that updates in a typical time step the history field, the crack phase field and finally the displacement field. We demonstrate the performance of the phase field formulation of fracture by means of representative numerical examples, which show the evolution of complex crack patterns under dynamic loading. [ABSTRACT FROM AUTHOR]

The fidelity of the peridynamic theory in predicting fracture is investigated through a comparative study. Peridynamic predictions for fracture propagation paths and speeds are compared against various experimental observations. Furthermore, these predictions are compared to the previous predictions from extended finite elements (XFEM) and the cohesive zone model (CZM). Three different fracture experiments are modeled using peridynamics: two experimental benchmark dynamic fracture problems and one experimental crack growth study involving the impact of a matrix plate with a stiff embedded inclusion. In all cases, it is found that the peridynamic simulations capture fracture paths, including branching and microbranching that are in agreement with experimental observations. Crack speeds computed from the peridynamic simulation are on the same order as those of XFEM and CZM simulations. It is concluded that the peridynamic theory is a suitable analysis method for dynamic fracture problems involving multiple cracks with complex branching patterns. [ABSTRACT FROM AUTHOR]

In the paper, the results of numerical failure analysis of plain concrete beams loaded by impact three-point bending load are presented and discussed. The theoretical framework for the numerical analysis is continuum mechanics and irreversible thermodynamics. The spatial discretization is performed by the finite element method using update Lagrange formulation. Green–Lagrange stain tensor is used as a strain measure. To account for cracking and damage of concrete, the beam is modeled by the rate sensitive microplane model with the use of the so-called co-rotational stress tensor. Damage and cracking phenomena are modeled within the concept of smeared cracks. To assure objectivity of the analysis with respect to the size of the finite elements, crack band method is used. The contact-impact analysis is based on the mechanical interaction between two bodies—concrete beam (master) and dropping hammer (slave) falling on the mid span of the beam. The contact constrains are satisfied by Lagrange multiplier method, which is adapted for the explicit time integration scheme. To investigate the influence of loading rate on the failure mode of the beam parametric study is carried out. The numerical results are evaluated, discussed and compared with test results known from the literature. It is shown that the beam resistance and failure mode strongly depend on loading rate. For lower loading rates beam fails in bending (mode-I fracture). However, with increasing loading rate there is a transition of the failure mechanism from bending to shear. The results are in good agreement with theoretical and experimental results known from the literature. [ABSTRACT FROM AUTHOR]

A finite element model of crack propagation along a sinusoidal interface with amplitude A and wavelength λ between identical elastic materials is presented. Interface decohesion is modeled with the Xu and Needleman (J Mech Phys Solid 42(9):1397, 1994) cohesive traction–separation law. Ancillary calculations using linear elastic fracture mechanics theory were used to explain some aspects of stable and unstable crack growth that could not be directly attained from the cohesive model. For small aspect ratios of the sinusoidal interface ( A/λ ≤ 0.25), we have used the analytical Cotterell–Rice (Intl J Fract 16:155–169, 1980) approximation leading to a closed-form expression of the effective toughness, K Ic , given by $$K_{Ic}\sqrt{(1-\nu^{2})/E\phi_n}=2/\left(1+\left[1+4\pi^{2}(A/\lambda)^{2}\right]^{-1/2}\right),$$ where $$\phi_n$$ is the work of separation, E is Young’s modulus, and ν is Poisson’s ratio. For A/λ > 0.25, both the cohesive zone model and numerical J-integral estimates of crack tip stress intensity factors suggest the following linear relationship: $$K_{Ic} \sqrt{(1-\nu^{2})/E\phi_n}=0.81+1.89(A/\lambda).$$ Parametric studies show that the length of the cohesive zone does not significantly influence K Ic , although it strongly influences the behavior of the crack between the initiation of stable crack growth and the onset of unstable fracture. [ABSTRACT FROM AUTHOR]

The application of rate-dependent cohesive elements is validated in simulation of ductile fracture in aluminum round bars under dynamic loading conditions. Smooth and notched round bars made of AA6060-T6 are tested and simulated under quasi-static and dynamic loadings. The smooth round bar is modeled using finite elements that obey Gurson–Tvergaard–Needleman (GTN) formulation as the constitutive equation. Comparing with experimental results, corresponding GTN parameters and rate-dependent plasticity of the alloy are obtained. A single strain rate-dependent GTN element with the obtained parameters is examined under different values of stress triaxiality and loading rates. The resulting stress-elongation curves represent the traction separation law (TSL) for cohesive elements and the variations of the maximum traction and the energy absorbed are investigated. The notched round bars are modeled by axisymmetric continuum and cohesive elements. The undamaged bulk material is elastic-visco plastic and the cohesive elements obey the TSL defined from the single element calculations. The experiments are simulated by these models in which the cohesive elements are rate sensitive and automatically obtain the values of the total strain rate from their adjacent continuum elements to update the values of the cohesive strength during the analysis. The results of the analysis, including maximum load, time of failure and diameter reduction are validated with the experimental results. The effects of element size, rate-dependent plasticity of the material and stress triaxiality are also discussed. [ABSTRACT FROM AUTHOR]