共查询到20条相似文献,搜索用时 0 毫秒
1.
Lin Chen Bin Li René de Borst 《International journal for numerical methods in engineering》2019,120(4):433-446
The analysis of (dynamic) fracture often requires multiple changes to the discretisation during crack propagation. The state vector from the previous time step must then be transferred to provide the initial values of the next time step. A novel methodology based on a least-squares fit is proposed for this mapping. The energy balance is taken as a constraint in the mapping, which results in a complete energy preservation. Apart from capturing the physics better, this also has advantages for numerical stability. To further improve the accuracy, Powell-Sabin B-splines, which are based on triangles, have been used for the discretisation. Since continuity of the displacement field holds at crack tips for Powell-Sabin B-splines, the stresses at and around crack tips are captured much more accurately than when using elements with a standard Lagrangian interpolation, or with NURBS and T-splines. The versatility and accuracy of the approach to simulate dynamic crack propagation are assessed in two case studies, featuring mode-I and mixed-mode crack propagation. 相似文献
2.
Hyunil Baek Chulmin Kweon Kyoungsoo Park 《International journal for numerical methods in engineering》2020,121(24):5719-5741
In this study, a computational framework is proposed to investigate multiscale dynamic fracture phenomena in materials with microstructures. The micro- and macro-scales of a composite material are integrated by introducing an adaptive microstructure representation. Then, the far and local fields are simultaneously computed using the equation of motion, which satisfies the boundary conditions between the two fields. Cohesive surface elements are dynamically inserted where and when needed, and the Park-Paulino-Roesler cohesive model is employed to approximate nonlinear fracture processes in a local field. A topology-based data structure is utilized to efficiently handle adjacency information during mesh modification events. The efficiency and validity of the proposed computational framework are demonstrated by checking the energy balances and comparing the results of the proposed computation with direct computations. Furthermore, the effects of microstructural properties, such as interfacial bonding strength and unit cell arrangement, on the dynamic fracture behavior are investigated. The computational results demonstrate that local crack patterns depend on the combination of microstructural properties such as unit cell arrangement and interfacial bonding strength; therefore, the microstructure of a material should be carefully considered for dynamic cohesive fracture investigations. 相似文献
3.
The present work studies a self‐similar high‐speed expanding crack of mode‐I in a ductile material with a modified cohesive zone model. Compared with existing Dugdale models for moving crack, the new features of the present model are that the normal stress parallel to crack faces is included in the yielding condition in the cohesive zone and traction force in the cohesive zone can be non‐uniform. For a ductile material defined by von Mises criterion without hardening, the present model confirms that the normal stress parallel to crack face increases with increasing crack speed and can be even larger than the normal traction in the cohesive zone, which justifies the necessity of including the normal stress parallel to the crack faces in the yielding condition at high crack speed. In addition, strain hardening effect is examined based on a non‐uniform traction distribution in the cohesive zone. 相似文献
4.
Marco Musto Giulio Alfano 《International journal for numerical methods in engineering》2015,103(5):313-341
This paper presents a novel formulation of a hereditary cohesive zone model able to effectively capture rate‐dependent crack propagation along a defined interface, over a wide range of applied loading rates and with a single set of seven input parameters only, as testified by the remarkable agreement with experimental results in the case of a double cantilever beam made of steel adherends bonded along a rubber interface. The formulation relies on the assumption that the measured fracture energy is the sum of a rate‐independent ‘rupture’ energy, related to the rupture of primary bonds at the atomic or molecular level, and of additional dissipation caused by other rate‐dependent dissipative mechanisms present in the material and occurring simultaneously to rupture. The first contribution is accounted for by introducing a damage‐type internal variable, whose evolution follows a rate‐independent law for consistency with the assumption of rate independence of the rupture energy. To account for the additional dissipation, a fractional‐calculus‐based linear viscoelastic model is used, because for many polymers, it is known to capture the material response within an extremely wide range of strain rates much more effectively than classic models based on an exponential kernel. To the authors' knowledge, this is the first application of fractional viscoelasticity to the simulation of fracture. © 2015 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons Ltd. 相似文献
5.
T. Rabczuk T. Belytschko 《International journal for numerical methods in engineering》2004,61(13):2316-2343
A new approach for modelling discrete cracks in meshfree methods is described. In this method, the crack can be arbitrarily oriented, but its growth is represented discretely by activation of crack surfaces at individual particles, so no representation of the crack's topology is needed. The crack is modelled by a local enrichment of the test and trial functions with a sign function (a variant of the Heaviside step function), so that the discontinuities are along the direction of the crack. The discontinuity consists of cylindrical planes centred at the particles in three dimensions, lines centred at the particles in two dimensions. The model is applied to several 2D problems and compared to experimental data. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
6.
Low temperature cracking is the major distress observed in asphalt pavements in the northern US and Canada. In the past years
fracture mechanics concepts were introduced to investigate the fracture properties of asphalt mixtures at low temperatures.
In this paper the cohesive zone model (CZM) is used to describe the fracture behavior of asphalt mixtures at low temperatures
and the interface element is used to numerically simulate the material response under monotonic loading. The simulation is
calibrated with the experimental results from a newly proposed semi circular bend (SCB) test. A parametric analysis of the
input material properties indicates that the tensile strength has a significant effect on the peak load in the SCB configuration,
the modulus has a strong effect on the calculated stiffness of the SCB specimen, and the fracture energy influences the post-peak
behavior of the asphalt mixtures. The calibrated numerical model was applied to simulate the low temperature cracking in a
simplified asphalt pavement and to study the influence of these material parameters on the performance of asphalt pavements. 相似文献
7.
M. F. S. F. De MOURA N. DOURADO J. J. L. MORAIS F. A. M. PEREIRA 《Fatigue & Fracture of Engineering Materials & Structures》2011,34(3):149-158
The objective of this work is to verify numerically the adequacy of the ENF and the ELS tests to determine the fracture toughness under mode II loading of cortical bovine bone tissue. A data‐reduction scheme based on the specimen compliance and the equivalent crack concept is proposed to overcome the difficulties inherent to crack monitoring during its growth. A cohesive damage model was used to simulate damage initiation and growth, thus assessing the efficacy of the proposed data‐reduction scheme. The influences of the initial crack length, local strength and toughness on the measured fracture energy were analysed, taking into account the specimen length restriction. Some limitations related to spurious influence on the fracture process zone of the central loading in the ENF test, and clamping conditions in the ELS test were identified. However, it was verified that a judicious selection of the geometry allows, in both cases, a rigorous estimation of bone toughness in mode II. 相似文献
8.
Leonardo Snozzi Jean‐François Molinari 《International journal for numerical methods in engineering》2013,93(5):510-526
We present a model that combines interface debonding and frictional contact. The onset of fracture is explicitly modeled using the well‐known cohesive approach. Whereas the debonding process is controlled by a new extrinsic traction separation law, which accounts for mode mixity, and yields two separate values for energy dissipation in mode I and mode II loading, the impenetrability condition is enforced with a contact algorithm. We resort to the classical law of unilateral contact and Coulomb friction. The contact algorithm is coupled together to the cohesive approach in order to have a continuous transition from crack nucleation to the pure frictional state after complete decohesion. We validate our model by simulating a shear test on a masonry wallette and by reproducing an experimental test on a masonry wall loaded in compression and shear. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
9.
Andy Ziccarelli;Amit Kanvinde;Gregory Deierlein; 《Fatigue & Fracture of Engineering Materials & Structures》2024,47(2):319-343
A procedure is outlined for calibrating a finite element model to simulate ductile cracking that employs a continuum damage criterion for ductile fracture initiation, termed the Stress Weighted Ductile Fracture Model, with an Adaptive Cohesive Zone method to simulate crack propagation. The proposed procedure is implemented and validated with data from 44 coupon-scale tests of an A913 structural steel. The fracture model calibration is performed in two stages to determine parameters to evaluate: (1) the continuum damage mechanics criterion and (2) ductile crack initiation and propagation. Sensitivity of the calibration to parameter uncertainty is evaluated using synthetic data, followed by calibration and validation against experimental test data. The results indicate that the computational model can accurately simulate the overall response of experimental specimens, as well as other key aspects of observed behavior, including crack tunneling and crack face closure. Guidelines are provided for practical calibration of the model components. 相似文献
10.
Martin Fagerstrm Ragnar Larsson 《International journal for numerical methods in engineering》2006,66(6):911-948
A general finite element approach for the modelling of fracture is presented for the geometrically non‐linear case. The kinematical representation is based on a strong discontinuity formulation in line with the concept of partition of unity for finite elements. Thus, the deformation map is defined in terms of one continuous and one discontinuous portion, considered as mutually independent, giving rise to a weak formulation of the equilibrium consisting of two coupled equations. In addition, two different fracture criteria are considered. Firstly, a principle stress criterion in terms of the material Mandel stress in conjunction with a material cohesive zone law, relating the cohesive Mandel traction to a material displacement ‘jump’ associated with the direct discontinuity. Secondly, a criterion of Griffith type is formulated in terms of the material‐crack‐driving force (MCDF) with the crack propagation direction determined by the direction of the force, corresponding to the direction of maximum energy release. Apart from the material modelling, the numerical treatment and aspects of computational implementation of the proposed approach is also thoroughly discussed and the paper is concluded with a few numerical examples illustrating the capabilities of the proposed approach and the connection between the two fracture criteria. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
11.
D. A. Cendón A. R. Torabi M. Elices 《Fatigue & Fracture of Engineering Materials & Structures》2015,38(5):563-573
This article explores the capability of the Cohesive Zone Model in predicting the critical load of blunt notched specimens made of coarse‐grained polycrystalline graphite, a brittle material that has gained the attention of researchers because of its favourable properties for protection against thermal loads. To that aim, 39 different tests on U‐notched and V‐notched specimens made of this material, with loading modes raging from mode I to mixed mode I/II, have been modelled by using the Cohesive Zone Model. The model has been implemented through the embedded crack approach, avoiding thus the necessity of defining the crack trajectory prior to the simulation because it is automatically generated once the maximum principal stress overcomes the tensile strength of the material. The numerical predictions obtained show good agreement with the experimental results. 相似文献
12.
Julien Leclerc Ling Wu Van Dung Nguyen Ludovic Noels 《International journal for numerical methods in engineering》2018,113(3):374-410
Modelling the entire ductile fracture process remains a challenge. On the one hand, continuous damage models succeed in capturing the initial diffuse damage stage but are not able to represent discontinuities or cracks. On the other hand, discontinuous methods, as the cohesive zones, which model the crack propagation behaviour, are suited to represent the localised damaging process. However, they are unable to represent diffuse damage. Moreover, most of the cohesive models do not capture triaxiality effect. In this paper, the advantages of the two approaches are combined in a single damage to crack transition framework. In a small deformation setting, a nonlocal elastic damage model is associated with a cohesive model in a discontinuous Galerkin finite element framework. A cohesive band model is used to naturally introduce a triaxiality‐dependent behaviour inside the cohesive law. Practically, a numerical thickness is introduced to recover a 3D state, mandatory to incorporate the in‐plane stretch effects. This thickness is evaluated to ensure the energy consistency of the method and is not a new numerical parameter. The traction‐separation law is then built from the underlying damage model. The method is numerically shown to capture the stress triaxiality effect on the crack initiation and propagation. 相似文献
13.
S. Mostofizadeh M. Fagerström R. Larsson 《International journal for numerical methods in engineering》2013,96(2):63-86
In this paper, a method to analyse and predict crack propagation in thin‐walled structures subjected to large plastic deformations when loaded at high strain rates—such as impact and/or blast—has been proposed. To represent the crack propagation independently of the finite element discretisation, an extended finite element method based shell formulation has been employed. More precisely, an underlying 7‐parameter shell model formulation with extensible directors has been extended by locally introducing an additional displacement field, representing the displacement discontinuity independently of the mesh. Of special concern in the paper has been to find a proper balance between, level of detail and accuracy when representing the physics of the problem and, on the other hand, computational efficiency and robustness. To promote computational efficiency, an explicit time step scheme has been employed, which however has been discovered to generate unphysical oscillations in the response upon crack propagation. Therefore, special focus has been placed to investigate these oscillations as well as to find proper remedies. The paper is concluded with three numerical examples to verify and validate the proposed model.Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
14.
Ragnar Larsson Martin Fagerstrm 《International journal for numerical methods in engineering》2005,62(13):1763-1788
A theoretical and computational framework which covers both linear and non‐linear fracture behaviour is presented. As a basis for the formulation, we use the material forces concept due to the close relation between on one hand the Eshelby energy–momentum tensor and on the other hand material defects like cracks and material inhomogeneities. By separating the discontinuous displacement from the continuous counterpart in line with the eXtended finite element method (XFEM), we are able to formulate the weak equilibrium in two coupled problems representing the total deformation. However, in contrast to standard XFEM, where the direct motion discontinuity is used to model the crack, we rather formulate an inverse motion discontinuity to model crack development. The resulting formulation thus couples the continuous direct motion to the inverse discontinuous motion, which may be used to simulate linear as well as non‐linear fracture in one and the same formulation. In fact, the linear fracture formulation can be retrieved from the non‐linear cohesive zone formulation simply by confining the cohesive zone to the crack tip. These features are clarified in the two numerical examples which conclude the paper. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
15.
A mesoscopic fracture model of equilibrium slit cracks in brittle solids, including inter-surface cohesive tractions acting near the crack tip, is analyzed and the effects of the cohesive tractions on the in-plane stress fields, crack-opening displacement profiles, and crack driving forces examined quantitatively for linear and penny-shaped cracks. The (numerical) analysis method is described in detail, along with results for four different cohesive forces. The assumed distribution of cohesive tractions were found to suppress the in-plane stress field adjacent to cracks in a homogeneous, isotropic medium when uniformly loaded in mode-I, and the suppression was a function of crack length. The crack-opening displacement profile was also perturbed and a new regime identified between the near-field Barenblatt zone and the far-field continuum zone. The extent of this `cohesive zone' was quantified by use of an interpolating function fit to the calculated profiles and found to be independent of crack size for a given cohesive tractions distribution. Furthermore, the crack-opening displacement at the edge of the cohesive zone was found to be independent of crack size, implying that despite significant perturbations to the stress field, the crack driving force at unstable equilibrium remains unchanged with crack size. 相似文献
16.
Jeong‐Hoon Song Ted Belytschko 《International journal for numerical methods in engineering》2009,77(3):360-385
A new method for modeling discrete cracks based on the extended finite element method is described. In the method, the growth of the actual crack is tracked and approximated with contiguous discrete crack segments that lie on finite element nodes and span only two adjacent elements. The method can deal with complicated fracture patterns because it needs no explicit representation of the topology of the actual crack path. A set of effective rules for injection of crack segments is presented so that fracture behavior beginning from arbitrary crack nucleations to macroscopic crack propagation is seamlessly modeled. The effectiveness of the method is demonstrated with several dynamic fracture problems that involve complicated crack patterns such as fragmentation and crack branching. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
17.
I. S. Nijin Anuradha Banerjee 《Fatigue & Fracture of Engineering Materials & Structures》2020,43(5):893-906
The mechanistic aspects of process of initiation of a mode‐I fatigue crack in an aluminium alloy (AA 2219‐T87) are studied in detail, both computationally as well as experimentally. Simulations are carried out under plane strain conditions with fatigue process zone modelled as stress‐state–dependent cohesive elements along the expected mode‐I failure path. An irreversible damage parameter that accounts for the progressive microstructural damage due to fatigue is employed to degrade cohesive properties. The simulations predict the location of initiation of the fatigue crack to be subsurface where the triaxiality and the opening tensile stresses are higher in comparison with that at the notch surface. Examination of the fracture surface profile of fracture test specimens near notch tip reveals a few types of regions and existence of a mesoscopic length scale that is the distance of the location of highest roughness from the notch root. A discussion is developed on the physical significance of the experimentally observed length scale. 相似文献
18.
Crack Growth Across a Strength Mismatched Bimaterial Interface 总被引:3,自引:0,他引:3
Crack growth across an interface between materials with different strength is examined by a cohesive zone model. The two materials
have identical elastic properties but different fracture process properties, or different yield stresses, which is modeled
by different cohesive stresses. The fracture criteria is a critical crack opening displacement. Load is represented by a stress
intensity factor defining a remote square root singular stress field. The results show that the ratio between the cohesive
stresses of the two materials primarily determines the behavior of the critical stress intensity factor. When the crack approaches
a material with a higher cohesive stress the crack tip is shielded, but if the crack approaches a material with smaller critical
crack opening displacement the maximum level of shielding is determined by the ratio between the critical crack opening displacements.
When a crack approaches a material with a lower cohesive stress it is exposed to an amplified load.
This revised version was published online in August 2006 with corrections to the Cover Date. 相似文献
19.
A truly predictive dynamic fracture model would require detailed information about the local fracture process. For instance, recent experimental evidence has conclusively demonstrated that rapid crack propagation (RCP) in brittle polymers such as PMMA is accompanied by the nucleation, growth, interaction and coalescence of microcracks and the advancing macroscopic fracture. Some insights into this phenomenon are offered and a novel computational model of this aspect of dynamic fracture is proposed. The procedure is based upon local material (cohesive) strength considerations. Here, the initial development of such procedures is presented. Application is at this stage restricted to single crack problems so that the effects of various geometrical, but more importantly, cohesive parameters on predicted fractures may be examined. Extension of the cohesive computational procedures of this work to multiple crack problems is proposed to be straightforward and without the need for extensive reconstructing of the computational procedures outlined. 相似文献
20.
Ragnar Larsson Niklas Jansson 《International journal for numerical methods in engineering》2002,54(4):473-497
The contribution of this paper concerns the fracture modelling of an interface with a fixed internal material surface in the context of geometrically non‐linear kinematics. Typical applications are composite laminates and adhesive/frictional joints in general. In the model development, a key feature is the concept of regularized strong discontinuity, which provides a regular deformation gradient within the interface. The deformation gradient within the interface is formulated in a multiplicative structure with a continuous part and a discontinuous part, whereby the interface deformation is interpreted as a transformation between the material damaged configuration and the actual spatial configuration. In analogy with the continuum formulation of hyper‐inelasticity, a constitutive framework is defined for the relation between the induced material traction and the displacement jump vector, which are defined on the material damaged interface configuration. Within this framework, a simple, but yet still representative, model for the delamination problem is proposed on the basis of a damage–plasticity coupling for the interface. The model is calibrated analytically in the large deformation context with respect to energy dissipation in mode I so that a predefined amount of fracture energy is dissipated. The paper is concluded with a couple of numerical examples that display the properties of the interface. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献