碳纤维布补强梁裂尖奇异性的光学焦散线实验研究

通过光学焦散线方法对碳纤维布加固含裂纹梁进行三点弯荷载下的静荷载和冲击实验,研究了其断裂特性。在静态实验研究中,记录了不同载荷状态下裂纹的焦散斑光学图像,计算了与之对应的裂纹尖端的应力强度因子,分析了试件在粘贴加固不同长度的碳纤维布情况下焦散斑和应力强度因子演化规律。在动态实验中,记录了裂纹扩展的过程和在此过程中不同时刻的焦散斑和裂纹扩展的路线。通过各瞬间的焦散斑图像,分析裂纹的扩展长度及速度,并得到各时刻的动态应力强度因子。     

Dynamic caustics and its applications

For the study of elastodynamic problems of propagating cracks it is necessary to evaluate the dynamic stress intensity factor K^d_I which depends on the form of expressions for the stress components existing at the running crack tip at any instant of the propagation of the crack and the corresponding dynamic mechanical and optical properties of the material of the specimen under identical loading conditions. In this paper the distortion of the form of the corresponding reflected caustic from the lateral faces of a dynamically loaded transparent and optically inert specimen containing a transverse crack running under constant velocity was studied on the basis of complex potential elasticity theory and the influence of this form on the value of the dynamic stress intensity factor was given. The method was applied to the study of a propagating Mode I crack in a PMMA specimen under various propagation velocities and the corresponding dynamic stress intensity factor K^d_I evaluated. Also, crack propagation behaviour of notched composites in dynamic loading modes are reviewed and evaluated. A relatively large data base using metal-epoxy particulates, rubber-toughened poly(methyl methacrylate), and Sandwich plates are given. In all cases, a combination of high-speed photography and the optical method of dynamic caustics has been used. Results on the dynamic crack propagation mode, fracture toughness and crack propagation velocities of several rubber-modified composite models are presented. The composite models studied include specimens with one and/or two ‘complex’ two-stage inclusions, i.e. PMMA round inclusions surrounded by concentric rubber rings, one and/or press-fifting inclusions without rubber interface, all under dynamic loading. In all cases both qualitative and quantitative results were obtained. Also, results on crack propagation mode, crack propagation velocity, stress intensity factors and on the influence of the sandwich phases on crack propagation mode are presented.     

Optical method of caustics applied in viscoelastic fracture analysis

The optical method of caustics is developed here to study the fracture of viscoelastic materials. By adopting a distribution of viscoelastic stress fields near the crack tip, the method of caustics is used to determine the viscoelastic fracture parameters from the caustic patterns near the crack tip. Two viscoelastic materials are studied. These are PMMA and ternary composites of HDPE/POE-g-MA/CaCO₃. The transmitted and reflective methods of caustics are performed separately to investigate viscoelastic fracture behaviors. The stress intensity factors (SIFs) versus time is determined by a series of shadow spot patterns combined with viscoelastic parameters evaluated by creep tests. In order to understand the viscoelastic fracture mechanisms of HDPE/POE-g-MA/CaCO₃ composites, their fracture surfaces are observed by a Scanning Electron Microscope (SEM). The results indicate that the method of caustics can be used to characterize the fracture behaviors of viscoelastic materials and further to optimize the design of polymer composites.     

Linear least squares approach for evaluating crack tip fracture parameters using isochromatic and isoclinic data from digital photoelasticity

In this work, digital photoelasticity technique is used to estimate the crack tip fracture parameters for different crack configurations. Conventionally, only isochromatic data surrounding the crack tip is used for SIF estimation, but with the advent of digital photoelasticity, pixel-wise availability of both isoclinic and isochromatic data could be exploited for SIF estimation in a novel way. A linear least square approach is proposed to estimate the mixed-mode crack tip fracture parameters by solving the multi-parameter stress field equation. The stress intensity factor (SIF) is extracted from those estimated fracture parameters. The isochromatic and isoclinic data around the crack tip is estimated using the ten‐step phase shifting technique. To get the unwrapped data, the adaptive quality guided phase unwrapping algorithm (AQGPU) has been used. The mixed mode fracture parameters, especially SIF are estimated for specimen configurations like single edge notch (SEN), center crack and straight crack ahead of inclusion using the proposed algorithm. The experimental SIF values estimated using the proposed method are compared with analytical/finite element analysis (FEA) results, and are found to be in good agreement.     

Damage and fracture prediction of plastic-bonded explosive by digital image correlation processing

By digital correlation processing of Scanning electronic microscopy (SEM) images, the paper presents the deformation and damage analysis of an energetic material—the plastic-bonded explosive (PBX) on mesoscopic scale. The analysis is made by observing the deformation field resulted from the digital image correlation (DIC) processing of the images corresponding to the loading steps and comparing with the surface profiles of the composite material so as to visualize the matter damage near a preset crack. The results show that the local deformation disturbance can reveal the material damage even happened underneath the specimen surface. The strain distribution in the front of the preset crack, can be used to predict the propagating route of the microcrack initiated from the tip of the pre-crack, which is related to the splitting fracture of the granular-based composite under compressive loading.     

Non-linear dynamical analysis of crack surface perturbations and their dependence on velocity and direction

The fracture surfaces of single crystal [1 0 0] silicon specimens, fractured under three-point bending (3PB) and subjected to a high strain energy upon cracking, revealed exceptional surface perturbations, generated during the unstable propagation. While macroscopically the crack is propagating on the (1 1 1) low energy cleavage plane, microscopic examination revealed small angled deviations from and fluctuations along that plane. Furthermore, while the crack is propagating at a velocity of nearly 3000 m/s in the direction, its velocity in the direction is two orders of magnitude lower, with distinctive surface perturbations. The amplitude and complexity of the perturbations increase as the normal velocity vector changes its direction and magnitude. These perturbations were recorded with a profilometer and analyzed using non-linear dynamical analysis tools. This study provides an opportunity to interpret surface phenomena of one of the most general cases of fracture and to study the effect of major variables on the nature of the perturbations involved, such as the local crack tip velocity and the crystallographic orientations. It is shown that the surface perturbations are chaotic deterministic in nature and can be described by high order non-linear differential equations; the order of the equation varying with the variations of the local velocity and direction.     

The application of high speed photography and white light speckle to the study of dynamic fracture

A white light speckle technique is combined with high speed photography to study the dynamic fracture of a sample of PMMA. A series of single exposure speckle patterns are photographically recorded at the rate of 2 × 10⁵ frames s^-1 during fracture of the sample. The displacement field around the tip of the propagating crack may then be determined by superimposing pairs of individual frames to form a ‘double-exposure’ speckle pattern. It is shown that the method gives values for the stress intensity factor K which are consistent with those obtained from other techniques.     

钢中脆硬粒子裂纹形成机理

钢中的脆硬粒子对钢的解理脆断有直接的影响，解理断裂源大都发生在脆硬粒子上.根据微裂纹形成的热力学模型，利用钢中脆硬粒子开裂时所释放的弹性应变能、位错塞积弹性能，所产生的表面能，对脆硬粒子裂纹形成机理进行分析.模型计算表明，正应力和位错塞积力都是脆硬粒子开裂的必要条件，这与实验事实相符；同时给出脆硬粒子开裂的临界条件计算方法，计算发现，脆硬粒子临界开裂应力不仅取决于脆硬粒子尺寸及表面能，而且与晶粒直径有一定的相关关系，当晶粒直径较小时，这种关系与实验测定的材料解理断裂应力与晶粒尺寸的关系一致，说明整体失稳解 关键词： 解理断裂 裂纹形核 脆硬粒子     

The fracture mechanics of glassy polymers

The application of fracture mechanics to glassy polymers, in particular crack growth in PMMA, is discussed. Particular attention is paid to two processes which modulate the energy supply to the crack tip: viscoelastic dissipation at slow crack speeds and specimen inertia at large crack speeds. The relation between fracture energy and crack speed is reviewed, and, where possible, fracture surface observations are correlated with dynamic behavior.     

Finite strain stress fields near the tip of an interface crack between a soft incompressible elastic material and a rigid substrate

We present a numerical study of finite strain stress fields near the tip of an interface crack between a rigid substrate and an incompressible hyperelastic solid using the finite element method (FEM). The finite element (FE) simulations make use of a remeshing scheme to overcome mesh distortion. Analyses are carried out by assuming that the crack tip is either pinned, i.e., the elastic material is perfectly bonded (no slip) to the rigid substrate, or the crack lies on a frictionless interface. We focus on a material which hardens exponentially. To explore the effect of geometric constraint on the near tip stress fields, simulations are carried out under plane stress and plane strain conditions. For both the frictionless interface and the pinned crack under plane stress deformation, we found that the true stress field directly ahead of the crack tip is dominated by the normal opening stress and the crack face opens up smoothly. This is also true for an interface crack along a frictionless boundary in plane strain deformation. However, for a pinned interface crack under plane strain deformation, the true opening normal stress is found to be lower than the shear stress and the transverse normal stress. Also, the crack opening profile for a pinned crack under plane strain deformation is completely different from those seen in plane stress and in plane strain (frictionless interface). The crack face flips over and the tip angle is almost tangential to the interface. Our results suggest that interface friction can play a very important role in interfacial fracture of soft materials on hard substrates.     

Atomistic simulations to study crack tip behaviour in single crystal of bcc niobium and hcp zirconium

In this article, molecular dynamics based simulations were carried out to study the fracture toughness of single crystals of niobium (Nb) and zirconium (Zr). Separate set of simulations were performed with different orientations of crack plane in Nb and Zr. In each case, efforts were made to align the crack front with principal planes of corresponding crystal structure, that was bcc for Nb and hcp for Zr. Deformation in single crystal of Nb and Zr was governed either by twinning or emergence of dislocations from the crack tip and surfaces. The deformation mechanism in single crystal of Nb and Zr also helps in governing the overall toughness of the material, while deforming via twinning leads to higher change in the crack length, whereas dislocations emerging from the crack tip blunts the opening and improves the toughness.     

Finite element analysis of quasistatic crack propagation in brittle media with voids or inclusions

A three-level finite element scheme is proposed for simulation of crack propagation in heterogeneous media including randomly distributed voids or inclusions. To reduce total degrees of freedom in the view of mesh gradation, the entire domain is categorized into three regions of different-level meshes: a region of coarse-level mesh, a region of intermediate-level mesh, and a region of fine-level mesh. The region of coarse-level mesh is chosen to be far from the crack to treat the material inhomogeneities in the sense of coarse-graining through homogenization, while the region near the crack is composed of the intermediate-level mesh to model the presence of inhomogeneities in detail. Furthermore, the region very near the crack tip is refined into the fine-level mesh to capture a steep gradient of elastic field due to the crack tip singularity. Variable-node finite elements are employed to satisfy the nodal connectivity and compatibility between the neighboring different-level meshes. Local remeshing is needed for readjustment of mesh near the crack tip in accordance with crack growth, and this is automatically made according to preset values of parameters determining the propagation step size of crack, and so the entire process is fully automatic. The effectiveness of the proposed scheme is demonstrated through several numerical examples. Meanwhile, the effect of voids and inclusions on the crack propagation is discussed in terms of T-stresses, with the aid of three-level adaptive scheme.     

Elastic behaviour of an edge dislocation near a sharp crack emanating from a semi-elliptical blunt crack

The interaction between an edge dislocation and a crack emanating from a semi-elliptic hole is dealt with. Utilizing the complex variable method, closed form solutions are derived for complex potentials and stress fields. The stress intensity factor at the tip of the crack and the image force acting on the edge dislocation are also calculated. The influence of the morphology of the blunt crack and the position of the edge dislocation on the shielding effect to the crack and the image force is examined in detail. The results indicate that the shielding or anti-shielding effect to the stress intensity factor increases acutely when the dislocation approaches the tip of the crack. The effect of the morphology of the blunt crack on the stress intensity factor of the crack and the image force is very significant.     

Mechanical characterization of polytetrafluoroethylene polymer using full-field displacement method

The aim of this work is to estimate two important material properties of the polytetrafluoroethylene (PTFE) polymer by means of a single experimental test. The displacement fields around a crack tip are used for estimating the modulus of elasticity (or, Young's modulus) and Poisson's ratio. These parameters are evaluated by fitting linear fracture mechanic expression of displacement fields in the vicinity of the crack, for mode I, to the experimental data. Measurements of these displacements are carried out using digital image correlation (DIC) method. In this way, the experimental procedure is conducted by loading a double-edge-cracked plate specimen. In order to validate the results, two available experimental tests have been performed. The modulus of elasticity is determined by means of the tensile test, using a standard test machine. Moreover, the Poisson's ratio is obtained by measuring lateral compressive and longitudinal extensional strain using DIC method.     

Elastic behaviour of a wedge disclination dipole near a sharp crack emanating from a semi-elliptical blunt crack

The interaction between a wedge disclination dipole and a crack emanating from a semi-elliptic hole is investigated. Utilising the complex variable method, the closed form solutions are derived for complex potentials and stress fields. The stress intensity factor at the tip of the crack and the image force acting on the disclination dipole center are also calculated. The influence of the morphology of the blunt crack and the position of the disclination dipole on the shielding effect to the crack and the image force is examined in detail. The results indicate that the shielding or anti-shielding effect to the stress intensity factor increases when the wedge disclination dipole approaches the tip of the crack. The effects of the morphology of the blunt crack on the stress intensity factor of the crack and the image force are very significant.     

Modeling of thermomechanical fracture of functionally graded materials with respect to multiple crack interaction

Different aspects of thermomechanical fracture of functionally graded materials (FGMs) are considered. Among them are the crack interaction problems in a functionally graded coating on a homogeneous substrate (FGM/H). The interaction between systems of edge cracks is investigated, as well as, how this mutual interaction influences the fracture process and the formation of crack patterns. The problem is formulated with respect to singular integral equations which are referred to the boundary equation methods. The FGM properties are modeled by exponential functions. The main fracture characteristics are calculated, namely, the stress intensity factors, the angles of deviation of the cracks from their initial propagation direction and the critical stresses when the crack starts to propagate. The last two characteristics are calculated using an appropriate fracture criterion. The problem contains different parameters, such as the geometry (location and orientation of cracks, their lengths, and the width of the FGM layer) and material parameters, i.e. the inhomogeneity parameters of elastic and thermal coefficients of the functionally graded material. The influence of these parameters on the thermo-mechanical fracture of FGM/H is investigated. As examples the following real material combinations are discussed: TiC/SiC, Al₂O₃/MoSi₂, MoSi₂/SiC, ZrO₂/nickel and ZrO₂/steel.     

Multiparadigm modeling of dynamical crack propagation in silicon using a reactive force field

We report a study of dynamic cracking in a silicon single crystal in which the ReaxFF reactive force field is used for several thousand atoms near the crack tip, while more than 100,000 atoms are described with a nonreactive force field. ReaxFF is completely derived from quantum mechanical calculations of simple silicon systems without any empirical parameters. Our results reproduce experimental observations of fracture in silicon including changes in crack dynamics for different crack orientations.     

Photoelastic investigation of interfacial fracture between orthotropic and isotropic materials

Static and dynamic fracture of interfaces between orthotropic and isotropic materials were studied using photoelasticity. In this study, a bi-material specimen made of PSM-1^® and Scotchply^® 1002, a unidirectional glass fiber reinforced epoxy composite, was used. Two fiber orientations, fibers parallel to the interface (α=0°) and fibers perpendicular to the interface (α=90°) were considered. Center crack bi-material specimens having different crack lengths were loaded quasi-statically and the full-field isochromatics were recorded using a digital camera. The complex stress intensity factor corresponding to each crack length was calculated from the isochromatics and the values were compared to that obtained from boundary collocation method. Dynamic interfacial fracture was studied with an edge crack bi-material geometry for the two different fiber orientations. The isochromatics around the propagating crack were recorded using a digital high-speed camera. The fracture parameters such as crack speed, complex stress intensity factor and energy release rate were extracted from the isochromatics using the asymptotic stress field equations. The complex stress intensity factor obtained from the static experiments was in close agreement with that calculated using the boundary collocation method. The results also indicated that the fiber orientation with respect to the interface influences the fracture parameters for stationary and propagating cracks.     

Molecular dynamics simulation of brittle fracture in silicon

Brittle fracture in silicon is simulated with molecular dynamics utilizing a modified embedded atom method potential. The simulations produce propagating crack speeds that are in agreement with previous experimental results over a large range of fracture energy. The dynamic fracture toughness is found to be equal to the energy consumed by creating surfaces and lattice defects in agreement with theoretical predictions. The dynamic fracture toughness is approximately 1/3 of the static strain energy release rate, which results in a limiting crack speed of 2/3 of the Rayleigh wave speed.