DETERMINATION OF INTERFACIAL MECHANICAL PARAMETERS FOR AN Al/EPOXY/Al_2O_3 SYSTEM BY USING PEEL TEST SIMULATIONS

Peel test measurements and simulations of the interfacial mechanical parameters for the Al/Epoxy/Al2O3 system are performed in the present investigation. A series of Al film thicknesses between 20 and 250 microns and three peel angles of 90, 135 and 180 degrees are considered. Two types of epoxy adhesives are adopted to obtain both strong and weak interface adhesions. A finite element model with cohesive zone elements is used to identify the interfacial parameters and simulate the peel test process. By simulating and recording normal stress near the crack tip, the separation strength is obtained, Furthermore, the cohesive energy is identified by comparing the simulated steady-state peel force and the experimental result. It is found from the research that both the cohesive energy and the separation strength can be taken as the intrinsic interfacial parameters which are dependent on the thickness of the adhesive layer and independent of the film thickness and peel angle.     

Tough and tunable adhesion of hydrogels: experiments and models

As polymer networks infiltrated with water, hydro-gels are major constituents of animal and plant bodies and have diverse engineering applications. While natural hydro-gels can robustly adhere to other biological materials, such as bonding of tendons and cartilage on bones and adhe-sive plaques of mussels, it is challenging to achieve such tough adhesions between synthetic hydrogels and engineer-ing materials. Recent experiments show that chemically anchoring long-chain polymer networks of tough synthetic hydrogels on solid surfaces create adhesions tougher than their natural counterparts, but the underlying mechanism has not been well understood. It is also challenging to tune sys-tematically the adhesion of hydrogels on solids. Here, we provide a quantitative understanding of the mechanism for tough adhesions of hydrogels on solid materials via a com-bination of experiments, theory, and numerical simulations. Using a coupled cohesive-zone and Mullins-effect model val-idated by experiments, we reveal the interplays of intrinsic work of adhesion, interfacial strength, and energy dissipation in bulk hydrogels in order to achieve tough adhesions. We fur-ther show that hydrogel adhesion can be systematically tuned by tailoring the hydrogel geometry and silanization time of solid substrates, corresponding to the control of energy dis-sipation zone and intrinsic work of adhesion, respectively. The current work further provides a theoretical foundation for rational design of future biocompatible and underwater adhesives.     

CRACK PROPAGATION IN POLYCRYSTALLINE ELASTIC-VISCOPLASTIC MATERIALS USING COHESIVE ZONE MODELS

Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated firstly in an elastic-viscoplastic polycrystals. The grain interior's material properties (e.g. strain rate sensitivity) characterize the competitions between plastic and cohesive energy dissipation mechanisms. The higher the strain rate sensitivity is, the larger amount of the external work is transformed into plastic dissipation energy than into cohesive energy, which delays the cohesive zone rupturing. With the strain rate sensitivity decreased, the material property tends to approach the elastic-plastic responses. In this case, the plastic dissipation energy decreases and the cohesive dissipation energy increases which accelerates the cohesive zones debonding. Increasing the cohesive strength or the critical separation displacement will reduce the stress triaxiality at grain interiors and grain boundaries. Enhancing the cohesive zones ductility can improve the matrix materials resistance to void damage.     

Inverse analysis determining interfacial properties between metal film and ceramic substrate with an adhesive layer

In the present study, peel tests and inverse analysis were performed to determine the interracial mechanical parameters for the metal film/ceramic system with an epoxy interface layer between film and ceramic. Al films with a series of thicknesses between 20 and 250 μm and three peel angles of 90°, 135° and 180° were considered. A finite element model with the cohesive zone elements was used to simulate the peeling process. The finite element results were taken as the training data of a neural network in the inverse analysis. The interracial cohesive energy and the separation strength can be determined based on the inverse analysis and peel experimental result     

“Equilibrium”and“non-equilibrium”turbulence

We review the concept of ‘‘equilibrium' in turbulence. It generally means a property of the energy spectrum, it can also be understood in terms of a scalar property, the Taylor–Kolmogorov formula relating the dissipation rate to the total energy and integral length scale. The implications of equilibrium and strong departure from equilibrium for turbulence modeling are stressed.     

SHEAR BEAM MODEL FOR INTERFACE FAILURE UNDER ANTIPLANE SHEAR (Ⅰ)—FUNDAMENTAL BEHAVIOR

The propagation of interlayer cracks and the resulting failure of the interface is a typical mode occurring in rock engineering and masonry structure. On the basis of the theory of elasto~plasticity and fracture mechanics, the shear beam model for the solution of interface failure was presented. The concept of `cohesive crack’ was adopted to describe the constitutive behavior of the cohesive interfacial layer. Related fundamental equations such as equilibrium equation, constitutive equations were presented. The behavior of a double shear beam bonded through cohesive layer was analytically calculated. The stable propagation of interface crack and process zone was investigated.     

THE DESIGN OF CHARACTERIZING-INTEGRAI SCHEMES AND THE APPLICATION TO THE SHALLOW WATER WAVE PROBLEMS

In this paper a new approach for designing upwind type schemes-the characterizing-integral method and its applied skills are introduced, The method is simple, convenient and eff ective. And the method isn't only limited to conservation laws unlike other methods and may be easily extended to multi-dimension problems. Furthermore, the numerical dissipation of the method can be flexibly regulated, so that it is especially suitable for solving various discontinuity problems. The paper shows us how to use this approach to simulate deformation and breaking of a nonlinear shallow water wave on a gentle slope, and to compute two-dimensional dam failure problem.     

Nonlinear progressive damage model for composite laminates used for low-velocity impact

In order to effectively describe the progressively intralaminar and interlam- inar damage for composite laminates, a three dimensional progressive damage model for composite laminates to be used for low-velocity impact is presented. Being applied to three-dimensional （3D） solid elements and cohesive elements, the nonlinear damage model can be used to analyze the dynamic performance of composite structure and its failure be- havior. For the intralaminar damage, as a function of the energy release rate, the damage model in an exponential function can describe progressive development of the damage. For the interlaminar damage, the damage evolution is described by the framework of the continuum mechanics through cohesive elements. Coding the user subroutine VUMAT of the finite element software ABAQUS/Explicit, the model is applied to an example, i.e., carbon fiber reinforced epoxy composite laminates under low-velocity impact. It is shown that the prediction of damage and deformation agrees well with the experimental results.     

The design of characterizing-integral schemes and the application to the shallow water wave problems

In this paper a new approach for designing upwind type schemes-the characterizing-integral method and its applied skills are introduced. The method is simple, convenient and eff ective. And the method isn ’t only limited to conservation laws unlike other methods and maybe easily extended to multi-dimension problems. Furthermore, the numerical dissipation of the method can be flexibly regulated, so that it is especially suitable for solving various discontinuity problems.The paper shows us now to use this approach to simulate deformation and breaking of a nonlinear shallow water wave on a gentle slope, and to compute two-dimensional dam failure problem.     

ANALYSIS ON THE COHESIVE STRESS AT HALF INFINITE CRACK TIP

The nonlinear fracture behavior of quasi-brittle materials is closely related with the cohesive force distribution of fracture process zone at crack tip. Based on fracture character of quasi-brittle materials, a mechanical analysis model of half infinite crack with cohesive stress is presented. A pair of integral equations is established according to the superposition principle of crack opening displacement in solids, and the fictitious adhesive stress is unknown function . The properties of integral equations are analyzed, and the series function expression of cohesive stress is certified. By means of the data of actual crack opening displacement, two approaches to gain the cohesive stress distribution are proposed through resolving algebra equation. They are the integral transformation method for continuous displacement of actual crack opening, and the least square method for the discrete data of crack opening displacement. The calculation examples of two approaches and associated discussions are give     

Mechanics of axisymmetric wavy surface adhesion: JKR–DMT transition solution

Recent work on the mechanics of detachment of a rigid sphere from an elastic axisymmetric wavy surface in the presence of JKR adhesion has shown that the presence of small-amplitude waviness introduces instabilities into the detachment process which dissipate mechanical energy. These instabilities result in interface toughening and strengthening; both the external work and peak force required for separation of a wavy interface are higher than those for a flat interface. In this paper, we summarize the key dimensionless parameters governing axisymmetric wavy surface adhesion in the JKR regime. We then proceed to derive a solution for the JKR–DMT adhesion transition for the axisymmetric wavy surface contact problem using a Maugis–Dugdale cohesive zone formulation. The phenomenon of interface toughening and strengthening due to the presence of surface waviness is seen to be restricted primarily to the JKR adhesion regime.     

Delamination of an elastic film from an elastic–plastic substrate during adhesive contact loading and unloading

Adhesive contact between a rigid sphere and an elastic film on an elastic–perfectly plastic substrate was examined in the context of finite element simulation results. Surface adhesion was modeled by nonlinear springs obeying a force-displacement relationship governed by the Lennard–Jones potential. A bilinear cohesive zone law with prescribed cohesive strength and work of adhesion was used to simulate crack initiation and growth at the film/substrate interface. It is shown that the unloading response consists of five sequential stages: elastic recovery, interface damage (crack) initiation, damage evolution (delamination), film elastic bending, and abrupt surface separation (jump-out), with plastic deformation in the substrate occurring only during damage initiation. Substrate plasticity produces partial closure of the cohesive zone upon full unloading (jump-out), residual tensile stresses at the front of the crack tip, and irreversible downward bending of the elastic film. Finite element simulations illustrate the effects of minimum surface separation (i.e., maximum compressive surface force), work of adhesion and cohesive strength of the film/substrate interface, substrate yield strength, and initial crack size on the evolution of the surface force, residual deflection of the elastic film, film-substrate separation (debonding), crack-tip opening displacement, and contact instabilities (jump-in and jump-out) during a full load–unload cycle. The results of this study provide insight into the interdependence of contact instabilities and interfacial damage (cracking) encountered in layered media during adhesive contact loading and unloading.     

Criteria on the artificial compliance inherent to the intrinsic cohesive zone

This Note presents criteria on the artificial compliance due to intrinsic cohesive zone models. A homogenized model is proposed for a collection of cohesive zones embedded between each mesh of a finite element-type discretization (cohesive-volumetric approach). The overall elastic behavior of this cohesive-volumetric medium is obtained as a function of the local properties and the mesh size. For an isotropic discretization, a criterion on the cohesive stiffnesses is derived: the additional compliance inherent to intrinsic cohesive zone models is bounded by lower value.     

Cohesive zone laws for void growth — II. Numerical field projection of elasto-plastic fracture processes with vapor pressure

Modeling ductile fracture processes using Gurson-type cell elements has achieved considerable success in recent years. However, incorporating the full mechanisms of void growth and coalescence in cohesive zone laws for ductile fracture still remains an open challenge. In this work, a planar field projection method, combined with equilibrium field regularization, is used to extract crack-tip cohesive zone laws of void growth in an elastic-plastic solid. To this end, a single row of void-containing cell elements is deployed directly ahead of a crack in an elastic-plastic medium subjected to a remote K-field loading; the macroscopic behavior of each cell element is governed by the Gurson porous material relation, extended to incorporate vapor pressure effects. A thin elastic strip surrounding this fracture process zone is introduced, from which the cohesive zone variables can be extracted via the planar field projection method. We show that the material's initial porosity induces a highly convex traction-separation relationship — the cohesive traction reaches the peak almost instantaneously and decreases gradually with void growth, before succumbing to rapid softening during coalescence. The profile of this numerically extracted cohesive zone law is consistent with experimentally determined cohesive zone law in Part I for multiple micro-crazing in HIPS. In the presence of vapor pressure, both the cohesive traction and energy are dramatically lowered; the shape of the cohesive zone law, however, remains highly convex, which suggests that diffusive damage is still the governing failure mechanism.     

Simulation of local instabilities during crack propagation in the ductile–brittle transition region

The crack propagation for a cohesive zone model within an elastic-plastic material under small-scale yielding conditions is simulated numerically. The resulting crack growth resistance curves show local instabilities, so-called pop-ins even for homogeneous cohesive properties if the cohesive strength lies sufficiently close to the maximum stress of the corresponding blunting solution. The formation of secondary cracks and unloading zones in front of the actual crack tip is identified as the underlying mechanism. It is found that the shape of the cohesive law has a considerable influence on the crack arrest behavior. Furthermore, requirements to the mesh resolution are derived.     

Cohesive zone laws for fatigue crack growth: Numerical field projection of the micromechanical damage process in an elasto-plastic medium

Cohesive zone failure models are widely used to simulate fatigue crack propagation under cyclic loading, but the model parameters are phenomenological and are not closely tied to the underlying micromechanics of the problem. In this paper, we will inversely extract the cohesive zone laws for fatigue crack growth in an elasto-plastic ductile solid using a field projection method (FPM), which projects the equivalent tractions and separations at the cohesive crack-tip from field information outside the process zone. In our small-scale yielding model, a single row of discrete voids is deployed directly ahead of a crack in an elasto-plastic medium subjected to cyclic mode I K-field loading. Damage accumulation under cyclic loading is captured by the growth of voids within the micro-voiding zone ahead of the crack, while the evolution of the cohesive zone law representing the micro-voiding zone is inversely extracted via the FPM. We show that the field-projected cohesive zone law captures the essential micromechanisms of fatigue crack growth in the ductile medium: from loading and unloading hysteresis caused by void growth and plastic hardening, to the softening damage locus associated with crack propagation via a void by void growth mechanism. The results demonstrate the effectiveness of the FPM in obtaining a micromechanics-based cohesive zone law in-place of phenomenological models, which opens the way for a unified treatment of fatigue crack problems.     

On the relationship between the work of adhesion and the critical shear stress for the onset of flow instabilities

The relationship between the critical shear stress for the onset of flow instabilities and the work of adhesion at the interface has been evaluated experimentally for a series of polyolefins on several steels, brasses, and coated steels. The critical shear stress was considerably affected by the chemical nature of the die. Low-surface-energy materials were found to produce smooth extrudates presumably by inducing slippage at the wall. Steel, tungsten carbide, and brass surfaces promoted sharkskin defect at shear stresses above the first critical shear stress. A linear relationship between the critical shear stress and the work of adhesion was found to be valid for values of work of adhesion smaller than approximately 30 mN/m. This indicates that slippage occurs due to a breakdown of the adhesion at the interface. For values of work of adhesion greater than approximately 30 mN/m a stronger interaction at the interface may induce a cohesive failure at the interface and subsequently, extrudate distortions.     

A comparison of cohesive zone modeling and classical fracture mechanics based on near tip stress field

A mode III crack with a cohesive zone in a power-law hardening material is studied under small scale yielding conditions. The cohesive law follows a softening path with the peak traction at the start of separation process. The stress and strain fields in the plastic zone, and the cohesive traction and separation displacement in the cohesive zone are obtained. The results show that for a modest hardening material (with a hardening exponent N = 0.3), the stress distribution in a large portion of the plastic zone is significantly altered with the introduction of the cohesive zone if the peak cohesive traction is less than two times yield stress, which implies the disparity in terms of the fracture prediction between the classical approach of elastic–plastic fracture mechanics and the cohesive zone approach. The stress distributions with and without the cohesive zone converge when the peak cohesive traction becomes infinitely large. A qualitative study on the equivalency between the cohesive zone approach and the classical linear elastic fracture mechanics indicates that smaller cracks require a higher peak cohesive traction than that for longer cracks if similar fracture initiations are to be predicted by the two approaches.     

Extraction of cohesive-zone laws from elastic far-fields of a cohesive crack tip: a field projection method

A solution method of an inverse problem is developed to extract cohesive-zone laws from elastic far-fields surrounding a crack-tip cohesive zone. The solution method is named the “field projection method (FPM).” In the process of developing the method a general form of cohesive-crack-tip fields is obtained and used for eigenfunction expansions of the plane elastic field in a complex variable representation. The closing tractions and the separation-gradients at the cohesive zone are expressed in terms of orthogonal polynomial series expansions of the general-form complex functions. The series expansion forms a set of cohesive-crack-tip eigenfunctions, which is complete and orthogonal in the sense of the interaction J-integral in the far field as well as at the cohesive-zone faces. The coefficients of the eigenfunctions in the J-orthogonal representation are extracted directly, using interaction J-integrals in the far field between the physical field of interest and auxiliary probing fields. The path-independence of the interaction J-integral enables us to identify the cohesive-zone variables, i.e. tractions and separations, and thus the cohesive-zone constitutive laws uniquely from the far-field data. A set of numerical algorithms is developed for the inversion method and the results from numerical experiments suggest that the proposed algorithms are well suited for extracting cohesive-zone laws from the far-field data. The set includes methods to find the position and size of a cohesive zone. Further included are discussions on error analysis and stability of the inversion scheme.