The wave propagation in a beam on a random elastic foundation

This paper presents a study of wave propagation in an infinite beam on a random Winkler foundation. The spatial variation of the foundation spring constant is modelled as a random field and the influence of the correlation length is studied. As it is impossible to determine the general stochastic Green’s function, the configurational average of the Green’s function and its correlation function are considered. These functions are found through the transformation of the stochastic equation of motion into the Dyson equation for the mean or coherent field and the Bethe–Salpeter equation for the field correlation, as used in the study of wave propagation in random media. The approximate solutions of the Dyson and the Bethe–Salpeter equations are validated by means of a Monte Carlo simulation and compared with the results of a classical Neumann expansion method. Although both methods only involve the second order statistics of the random field, the approximation of the Dyson and the Bethe–Salpeter equations gives better results than the Neumann expansion, at the expense of a larger computational effort. Furthermore, the results show that a small spatial variation of the spring constant has an influence on the response if the correlation length and the wavelength have a similar order of magnitude, while the waves in the beam cannot resolve the spatial variation in the case where the correlation length is much smaller than the wavelength.     

Finite element modelling of multiple cohesive discrete crack propagation in reinforced concrete beams

This paper presents a finite element (FE) model for fully automatic simulation of multiple discrete crack propagation in reinforced concrete (RC) beams. The discrete cracks are modelled based on the cohesive/fictitious crack concept using nonlinear interface elements with a bilinear tensile softening constitutive law. The model comprises an energy-based crack propagation criterion, a simple remeshing procedure to accommodate crack propagations, two state variable mapping methods to transfer structural responses from one FE mesh to another, and a local arc-length algorithm to solve system equations characterised by material softening. The bond-slip behaviour between reinforcing bars and surrounding concrete is modelled by a tension-softening element. An example RC beam with well-documented test data is simulated. The model is found capable of automatically modelling multiple crack propagation. The predicted cracking process and distributed crack pattern are in close agreement with experimental observations. The load-deflection relations are accurately predicted up to a point when compressive cracking becomes dominant. The effects of bond-slip modelling and the efficiency and effectiveness of the numerical algorithms, together with the limitations of the current model, are also discussed.     

Cohesive zone modeling of interface fracture in elastic bi-materials

Some basic issues regarding the cohesive zone modeling of interface fracture between two dissimilar elastic materials are studied. The dependence of the cohesive energy density on the phase angle is first discussed under small scale cohesive zone conditions. It is then shown that in general the stress singularities in tension and shear cannot be simultaneously removed at the cohesive zone tip if a single cohesive zone length is adopted for both tensile and shear fracture modes. Finally, the energy dissipation at the tip of a prescribed cohesive zone is examined using a bilinear cohesive zone model under the uncoupled tension/shear conditions.     

Modelling crack propagation in reinforced concrete using a hybrid finite element–scaled boundary finite element method

A previously developed hybrid finite element–scaled boundary finite element method (FEM–SBFEM) is extended to model multiple cohesive crack propagation in reinforced concrete. This hybrid method can efficiently extract accurate stress intensity factors from the semi-analytical solutions of SBFEM and is also flexible in remeshing multiple cracks. Crack propagation in the concrete bulk is modelled by automatically inserted cohesive interface elements with nonlinear softening laws. The concrete–reinforcement interaction is also modelled by cohesive interface elements. The bond shear stress–slip relation of CEB-FIP Model Code 90 and an empirical confining stress–crack opening relation are used to characterise slip and split failure at the concrete–reinforcement interface, respectively. Three RC beams were simulated. The numerical results agreed well with both experimental and numerical results available in the literature. Parametric studies demonstrated the importance of modelling both slip and split failure mechanisms at the concrete–reinforcement interface.     

Dynamics steady crack propagation of antiplane shear in an elastic perfectly plastic material

Abstract

A theoretical analysis of steady‐state crack growth in an elastic ideally‐plastic material under small‐scale yielding conditions has been carried out for anti‐plane shear. Asymptotic expansion method is used to construct the solutions for the region near the crack line. Exact solutions for the distribution of strain on the crack line within the primary active plastic zone is obtained. It is shown that the solution reduces to the correct asymptotic form as the crack speed approaches zero (quasi‐static) for any point on the crack line. The results are used to discuss the applicability of quasi‐static solutions to moving steady‐state situations. It is found that if the crack propagation speed is less than 0.1 of the shear wave speed, the quasi‐static solutions can be accurately approximated for the steady state solutions.     

Limited weld residual stress measurements in fatigue crack propagation: Part II. FEM-based fatigue crack propagation with complete residual stress fields

Using a limited set of residual stress measurements acquired by neutron diffraction and an equilibrium‐based, weighted least square algorithm to reconstruct the complete residual stress tensor field from the measured residual stress data, the effect of weld residual stress on fatigue crack propagation is investigated for 2024‐T351 aluminium alloy plate joined by friction stir welding. Through incorporation of the least squares, complete equilibrated residual stress field into a finite element model of the Friction Stir Weld (FSW) region, progressive crack growth along a direction perpendicular to the welding line is simulated as part of the analysis. Both the residual stress redistribution and the stress intensity factor due to the residual stress field, K_res, are calculated during the crack extension process. Results show that (a) incorporation of the complete, self‐equilibrated residual stress field into a finite element (FE) model of the specimen provides a robust, hybrid approach for assessing the importance of residual stress on fatigue crack propagation, (b) the calculated stress‐intensity factor due to the residual stress field, K_res, has the same trend as measured experimentally by the ‘cut‐compliance method’ and (c) the da/dN results are readily explained with reference to the effect of the residual stress field on the applied stress intensity factor.     

Multiple random crack propagation using a boundary element formulation

This paper proposes a boundary element method (BEM) model that is used for the analysis of multiple random crack growth by considering linear elastic fracture mechanics problems and structures subjected to fatigue. The formulation presented in this paper is based on the dual boundary element method, in which singular and hyper-singular integral equations are used. This technique avoids singularities of the resulting algebraic system of equations, despite the fact that the collocation points coincide for the two opposite crack faces. In fracture mechanics analyses, the displacement correlation technique is applied to evaluate stress intensity factors. The maximum circumferential stress theory is used to evaluate the propagation angle and the effective stress intensity factor. The fatigue model uses Paris’ law to predict structural life. Examples of simple and multi-fractured structures loaded until rupture are considered. These analyses demonstrate the robustness of the proposed model. In addition, the results indicate that this formulation is accurate and can model localisation and coalescence phenomena.     

Dynamic crack propagation with cohesive elements: a methodology to address mesh dependency

In this paper, two brittle fracture problems are numerically simulated: the failure of a ceramic ring under centrifugal loading and crack branching in a PMMA strip. A three‐dimensional finite element package in which cohesive elements are dynamically inserted has been developed. The cohesive elements' strength is chosen to follow a modified weakest link Weibull distribution. The probability of introducing a weak cohesive element is set to increase with the cohesive element size. This reflects the physically based effect according to which larger elements are more likely to contain defects. The calculations illustrate how the area dependence of the Weibull model can be used to effectively address mesh dependency. On the other hand, regular Weibull distributions have failed to reduce mesh dependency for the examples shown in this paper. The ceramic ring calculations revealed that two distinct phenomena appear depending on the magnitude of the Weibull modulus. For low Weibull modulus, the fragmentation of the ring is dominated by heterogeneities. Whereas many cracks were generated, few of them could propagate to the outer surface. Monte Carlo simulations revealed that for highly heterogeneous rings, the number of small fragments was large and that few large fragments were generated. For high Weibull modulus, signifying that the ring is close to being homogeneous, the fragmentation process was very different. Monte Carlo simulations highlighted that a larger number of large fragments are generated due to crack branching. Copyright © 2003 John Wiley & Sons, Ltd.     

On the intersonic crack propagation in an orthotropic medium

Motivaded by recent theoretical studies the elastodynamic response of an orthotropic material with a semi-infinite line crack, which propagates intersonically. is revisited through an approach which differs from those used in previous studies. The near tip stress and displacement fields are obtained for Mode I and Mode II of steady state crack propagation. The strain energy release rate analysis confirms that the Mode I is physically impossible due to the order of stress singularity, which is larger then one half. For Model II the order of stress is less than one half and it is shown that a steady state intersonic propagation is allowed only for a particular crack tip velocity which is a function of the material orthotropy.     

Predictions of crack propagation using a variational multiscale approach and its application to fracture in laminated fiber reinforced composites

Predictions of crack propagation is a valuable resource for ensuring structural integrity and damage tolerance of aerospace structures. Towards that end, a variational multiscale approach to predict mixed mode in-plane cohesive crack propagation is presented here. To demonstrate applicability and to provide validation of the finite element based predictive methodology, a comparative study of the numerical results with the corresponding experimental observations of crack propagation in laminated fiber reinforced composite panels is presented.     

A rate-dependent cohesive model for simulating dynamic crack propagation in brittle materials

Numerical investigations are conducted to simulate high-speed crack propagation in pre-strained PMMA plates. In the simulations, the dynamic material separation is explicitly modeled by cohesive elements incorporating an initially rigid, linear-decaying cohesive law. Initial attempts using a rate-independent cohesive law failed to reproduce available experimental results as numerical crack velocities consistently overestimate experimental observations. As proof of concept, a phenomenological rate-dependent cohesive law, which bases itself on the physics of microcracking, is introduced to modulate the cohesive law with the macroscopic crack velocity. We then generalize this phenomenological approach by establishing a rate-dependent cohesive law, which relates the traction to the effective displacement and rate of change of effective displacement. It is shown that this new model produces numerical results in good agreement with experimental data. The analysis demonstrates that the simulation of high-speed crack propagation in brittle structures necessitates the use of rate-dependent cohesive models, which account for the complicated rate-process of dynamic fracture at the propagating crack tip.     

An adaptive multiscale method for crack propagation and crack coalescence

This work presents a new multiscale technique for the efficient simulation of crack propagation and crack coalescence of macrocracks and microcracks. The fully adaptive multiscale method is able to capture localization effect mesh independently. By modeling macrocracks and microcracks, the extended finite element method offers an accurate solution and captures cracks and their propagation without changing the mesh topology. Propagating and coaliting cracks of different length scales can therefore be easily modeled and updated during the computation process. Hence, the presented method is an efficient and accurate option for modeling cracks of different length scales. This is demonstrated in several numerical examples showing the interaction of microcracks and macrocracks. Copyright © 2012 John Wiley & Sons, Ltd.     

陶瓷基复合材料基体随机开裂的损伤模拟

首先介绍了描述复合材料出现损伤时细观应力场的剪滞模型;给出了模拟陶瓷基复合材料基体随机开裂的有关公式;然后采用Monte Carlo方法模拟了基体随机开裂的过程;最后通过计算机模拟了陶瓷基复合材料基体裂纹的随机演化过程,分析了有关模型参数对基体裂纹演化的影响,并与文献模拟结果进行了对比。研究表明：最终裂纹间距主要分布在 1～2倍滑移长度之间;初始开裂应力越大,最终名义裂纹间距越小;Weibull 模量越大,最终名义裂纹间距越小;热残余应力越大,最终名义裂纹间距越大;模拟长度对裂纹演化影响很小;MonteCarlo方法可以有效地模拟陶瓷基复合材料的随机开裂过程。     

Cohesive modeling of fatigue crack retardation in polymers: Crack closure effect

A numerical model of fatigue crack growth retardation in polymers induced by artificial crack closure is proposed. The approach relies on the combination of cohesive modeling and a contact algorithm incorporated in the wake of the advancing crack to account for the effect of the introduced wedge. Numerical results are compared with existing experimental observations, showing the ability of the cohesive model to capture the key features of wedge-induced crack retardation. A study is conducted to quantify the effects of relevant parameters such as applied load levels, wedge distance to the crack tip and wedge stiffness. The model is also discussed in the context of self-healing polymers [White SR, Sottos NR, Moore J, Geubelle PH, Kessler M, Brown E, et al. Autonomic healing of polymer composites. Nature 2001;409:794-7], where the wedging effect is associated with the polymerization of the healing agent.     

Investigation of crack propagation scatter in a gear tooth’s root

This paper describes the problem of determining crack initiation location and its influence on crack propagation in a gear tooth’s root. Three different load positions on the gear tooth’s flank were considered for this investigation of crack initiation and propagation. A special test device was used for the single tooth test. It can be concluded from the measurements that a crack can be initiated at very different locations in a tooth’s root and then propagate along its own paths. A numerical investigation into a crack initiation’s position and its influences on its propagation were carried out within the framework of linear fracture mechanics. The influence of a tooth’s load position, the geometry of the tooth’s root, and the influence of non-parallel load distribution on the tooth’s flank were considered when investigating the crack initiation’s position. Results show that linear fracture mechanics can be used for determining crack propagation, if better initial conditions for crack initiation are considered.     

A model for short crack propagation in polycrystalline materials

A model for microstructurally short crack propagation in a grain structure of a polycrystalline material is developed. The crack propagation model is based on a crystal plasticity model and a microstructurally short crack propagation model in the spirit of the model by Navarro and de los Rios [A model for short fatigue crack propagation with an interpretation of the short-long crack transition. Fatigue Fract Eng Mater Struct 1987;10:169-86]. Numerical examples, where the combined crystal plasticity and crack propagation model is implemented in a model of a microstructure representing a duplex stainless steel, concludes the paper. Results showing how the misorientation of the crack- and slip-directions between two adjacent austenitic grains influences the crack propagation rate, as the crack propagates across their common grain boundary, are given.     

任意分布参数的涡轮盘裂纹扩展寿命可靠性分析

讨论了某型涡轮盘随机参数服从任意分布时的可靠性问题.在基本随机参数前四阶矩已知的情况下,以Pairs-Erdogan裂纹扩展模型为基础,应用随机摄动理论和Edgeworth 级数技术,采用疲劳寿命模型对某型涡轮盘随机参数服从任意分布时的可靠性进行分析,建立了涡轮盘疲劳寿命可靠性分析模型,并求得了涡轮盘裂纹扩展寿命的可靠度.模型计算结果与Monte-Carlo仿真结果非常接近,文章提出的方法对涡轮盘可靠性设计具有一定的参考价值.     

Elastic wave propagation in a 3‐D unbounded random heterogeneous medium coupled with a bounded medium. Application to seismic soil–structure interaction (SSSI)

We present a sub‐structuring method for the coupling between a large elastic structure, and a stratified soil half‐space exhibiting random heterogeneities over a bounded domain and impinged by incident waves. Both media are also weakly dissipative. The concept of interfaces classically used in sub‐structuring methods is extended to ‘volume interfaces’ in the proposed approach. The random dimension of the stochastic fields modelling the heterogeneities in the soil is reduced by introducing a Karhunen–Loéve expansion of these stochastic fields. The coupled overall problem is solved by Monte‐Carlo simulation techniques. A realistic example of a large industrial structure interacting with an uncertain stratified soil medium under earthquake is finally presented. This case study and others validate the presented methodology and its ability to handle complex mechanical systems. Copyright © 2002 John Wiley & Sons, Ltd.