Adiabatic decohesion in a thermoplastic craze thickening at constant or increasing rate

When a crack in a thermally non-diffusive material is impact loaded—or propagates at high speed—a cohesive process which resists slow crack extension may itself cause decohesion by adiabatic heating. By assuming that decohesion ultimately occurs by low-energy disentanglement within a melt layer of critical thickness, the fracture resistance of craze-forming crystalline polymers can be estimated quantitatively. Previous estimates used a simple, thermomechanically linear representation of craze fibril drawing. This paper presents a more physically realistic, numerical formulation, and demonstrates it for constant craze thickening rate (as imposed by an ideal full-notch tension test) and for linearly increasing thickening rate (as at the tip of an impact-loaded or rapidly propagating crack). For a linear material, the numerical formulation gives results which asymptotically approach those from analytical solutions, as craze density approaches zero. In more realistic model polymers, the enthalpy of fusion increasingly delays decohesion as impact speed increases, although the temperature distribution of an endotherm appears to have little effect. Increasing molecular weight, heuristically associated with decreasing craze density and increasing structural dimension, increases the predicted impact fracture resistance. In every case, fracture resistance passes through a minimum as impact speed increases. The conclusions encourage the use of impact fracture tests, and discourage the use of the full-notch tension test, to assess the dynamic fracture resistance of a craze-forming polymer.     

Theoretical development and experimental validation of a thermally dissipative cohesive zone model for dynamic fracture of amorphous polymers

A thermally dissipative cohesive zone model is developed for predicting the temperature increase at the tip of a crack propagating dynamically in a nominally brittle material exhibiting a cohesive-type failure such as crazing. The model assumes that fracture energy supplied to the crack tip region that is in excess of that needed for the creation of new free surfaces during crack advance is converted to heat within the cohesive zone. Bulk dissipation mechanisms, such as plasticity, are not accounted for. Several cohesive traction laws are examined, and the model is then used to make predictions of crack tip heating at various crack propagation speeds in the nominally brittle amorphous polymer PMMA, observed to fail by a crazing-type mechanism. The heating predictions are compared to experimental data where the temperature field surrounding a high speed crack in PMMA was measured. Measurements are made in real time using a multi-point high speed HgCdTe infrared radiation detector array. At the same time as temperature, simultaneous measurement of fracture energy is made by a strain gauge technique, and crack tip speed is monitored through a resistance ladder method. Material strength can be estimated through uniaxial tension tests, thus minimizing the need for parameter fitting in the stress-opening traction law. Excellent agreement between experiments and theory is found for two of the cohesive traction law temperature predictions, but only for the case where a single craze is active during the dynamic fracture of PMMA, i.e. crack tip speed up to approximately 0.2c_R. For higher speed fracture where subsurface damage becomes prominent, the line dissipation model of a cohesive zone is inadequate, and a distributed damage model is needed.     

A simple cohesive zone model that generates a mode-mixity dependent toughness

A simple, mode-mixity dependent toughness cohesive zone model (MDG_c CZM) is described. This phenomenological cohesive zone model has two elements. Mode I energy dissipation is defined by a traction–separation relationship that depends only on normal separation. Mode II (III) dissipation is generated by shear yielding and slip in the cohesive surface elements that lie in front of the region where mode I separation (softening) occurs. The nature of predictions made by analyses that use the MDG_c CZM is illustrated by considering the classic problem of an elastic layer loaded by rigid grips. This geometry, which models a thin adhesive bond with a long interfacial edge crack, is similar to that which has been used to measure the dependence of interfacial toughness on crack-tip mode-mixity. The calculated effective toughness vs. applied mode-mixity relationships all display a strong dependence on applied mode-mixity with the effective toughness increasing rapidly with the magnitude of the mode-mixity. The calculated relationships also show a pronounced asymmetry with respect to the applied mode-mixity. This dependence is similar to that observed experimentally, and calculated results for a glass/epoxy interface are in good agreement with published data that was generated using a test specimen of the same type as analyzed here.     

Bifurcations of Heteroclinic Orbits

The search for traveling wave solutions of a semilinear diffusion partial differential equation can be reduced to the search for heteroclinic solutions of the ordinary differential equation ü − cu̇ + f(u) = 0, where c is a positive constant and f is a nonlinear function. A heteroclinic orbit is a solution u(t) such that u(t) → γ ₁ as t → −∞ and u(t) → γ ₂ as t → ∞ where γ ₁, γ ₂ are zeros of f. We study the existence of heteroclinic orbits under various assumptions on the nonlinear function f and their bifurcations as c is varied. Our arguments are geometric in nature and so we make only minimal smoothness assumptions. We only assume that f is continuous and that the equation has a unique solution to the initial value problem. Under these weaker smoothness conditions we reprove the classical result that for large c there is a unique positive heteroclinic orbit from 0 to 1 when f(0) = f(1) = 0 and f(u) > 0 for 0 < u < 1. When there are more zeros of f, there is the possibility of bifurcations of the heteroclinic orbit as c varies. We give a detailed analysis of the bifurcation of the heteroclinic orbits when f is zero at the five points −1 < −θ < 0 < θ < 1 and f is odd. The heteroclinic orbit that tends to 1 as t → ∞ starts at one of the three zeros, −θ, 0, θ as t → −∞. It hops back and forth among these three zeros an infinite number of times in a predictable sequence as c is varied.     

On the determination of the cohesive zone properties of an adhesive layer from the analysis of the wedge-peel test

An extensive numerical study of the mechanics of the “wedge-peel test” is performed in order to analyze the mode I steady state debonding of a sandwich structure made of two thin plastically deforming metallic plates bonded with an adhesive. The constitutive response of the metallic plates is modeled by J₂ flow theory, and the behavior of the adhesive layer is represented with a cohesive zone model characterized by a maximum separation stress and the fracture energy. A steady-state finite element code accounting for finite rotation has been developed for the analysis of this problem. Calculations performed with the steady-state formulation are shown to be much faster than simulations involving both crack initiation and propagation within a standard, non-steady-state code. The goal of this study is to relate the measurable parameters of the test to the corresponding fracture process zone characteristics for a representative range of adherent properties and test conditions. An improved beam bending model for the energy release rate is assessed by comparison with the numerical results. Two procedures are proposed for identifying the cohesive zone parameters from experimental measurements.     

Stability of an Initially, Stably Stratified Fluid Subjected to a Step Change in Temperature

The onset of convective instability in an initially quiescent, stably stratified fluid layer between two horizontal plates is analyzed with linear theory. The bottom boundary is heated suddenly from below, subjected to a step change in surface temperature. The critical time t _c to mark the onset of Rayleigh-Bénard convection is predicted by propagation theory. This theory uses the length scaled by , where α denotes thermal diffusivity. Under the normal mode analysis the dimensionless disturbance equations are obtained as a function of τ(=αt/d ²) and ζ(=Z/), where d is the fluid layer depth and Z is the vertical distance. The resulting equations are transformed to self-similar ones by using scaling and finally fixing τ as τ_c under the frame of coordinates τ and ζ. For a given γ, Pr and τ_c, the minimum value of Ra is obtained from the marginal stability curve. Here γ denotes the temperature ratio to represent the degree of stabilizing effect, Pr is the Prandtl number and Ra is the Rayleigh number. With γ=0, the minimum Ra value approaches the well-known value of 1708 as τ_c increases. However, it is inversely proportional to τ_c ^3/2 as τ_c decreases. With increasing γ, the system becomes more stable. It is interesting that in the present system, propagation theory produces the stability criteria to bound the available experimental data over the whole domain of time. Received 5 November 2001 and accepted 29 March 2002 Published online: 2 October 2002 RID="*" ID="*" This work has been supported by both SK Chemicals Co. Ltd. and LG Chemical Ltd., Seoul under the Brain Korea 21 Project of the Ministry of Education. Communicated by H.J.S. Fernando     

Asymptotic Variational Wave Equations

We investigate the equation (u _t +(f(u)) _x ) _x =f ^{′ ′}(u) (u _x )²/2 where f(u) is a given smooth function. Typically f(u)=u ²/2 or u ³/3. This equation models unidirectional and weakly nonlinear waves for the variational wave equation u _tt − c(u) (c(u)u _x ) _x =0 which models some liquid crystals with a natural sinusoidal c. The equation itself is also the Euler–Lagrange equation of a variational problem. Two natural classes of solutions can be associated with this equation. A conservative solution will preserve its energy in time, while a dissipative weak solution loses energy at the time when singularities appear. Conservative solutions are globally defined, forward and backward in time, and preserve interesting geometric features, such as the Hamiltonian structure. On the other hand, dissipative solutions appear to be more natural from the physical point of view.We establish the well-posedness of the Cauchy problem within the class of conservative solutions, for initial data having finite energy and assuming that the flux function f has a Lipschitz continuous second-order derivative. In the case where f is convex, the Cauchy problem is well posed also within the class of dissipative solutions. However, when f is not convex, we show that the dissipative solutions do not depend continuously on the initial data.     

Waterhammer with fluid-structure interaction

The classical theory of waterhammer is a well-known and accepted basis for the prediction of pressure surges in piping systems. In this theory the piping system is assumed not to move. In practice however piping systems move when they are loaded by severe pressure surges, which for instance occur after rapid valve closure or pump failure. The motion of the piping system induces pressure surges which are not taken into account in the classical theory.In this article the interaction between pressure surges and pipe motion is investigated. Three interaction mechanisms are distinguished: friction, Poisson and junction coupling. Numerical experiments on a single straight pipe and a liquid loading line show that interaction highly influences the extreme pressures during waterhammer occurrences.Nomenclature A _f cross-sectional discharge area - A _t cross-sectional pipe wall area - c _f pressure wave speed in fluid - c _t axial stress wave speed in pipe wall - E Young's modulus for pipe wall material - e pipe wall thickness - f Darch-Weisbach friction coefficient - g gravitational acceleration - H fluid pressure head - h elevation of pipe - K _f fluid bulk modulus - P fluid pressure - R internal radius of pipe - T _c valve closure time - t time - u _r radial displacement of pipe - u _z axial displacement of pipe - u _z axial velocity of pipe - V fluid velocity - z distance along pipe - elevation angle of pipe - length of pressure wave - Poisson's ratio - _f fluid density - _t density of pipe wall material - t _z axial pipe stress - hoop stress     

COUPLING BETWEEN LAMINAR FILM CONDENSATION AND NATURAL CONVECTION ON OPPOSITE SIDES OF A VERTICAL PLATE

Theoretical analyses which incorporate one-dimensional heat conduction along a plate and transverse heat conduction approximations are presented to predict the net heat transfer between laminar film condensation of a saturated vapour on one side of a vertical plate and boundary layer natural convection on the other side. It is assumed that countercurrent boundary layer flows are formed on the two sides. The governing boundary layer equations of this problem and their corresponding boundary conditions are all cast into dimensionless forms by using a non-similarity transformation. Thus the resulting system of equations can be solved by using the local non-similarity method for the boundary layer equations and a finite difference method for the heat conduction equation of the plate. The plate temperature and the heat flux through the plate are repetitively determined until the solutions for each side of the plate match. The predicted results show that the effect of Pr_c is not negligible for larger values of A* (thermal resistance ratio between natural convecti on side and condensing film side) and the approximation of transverse heat conduction overpredicts the plate temperature for lower values of R_t (thermal resistance ratio between plate and condensing film). However, no significant differences are observed between the two different approximations for higher values of R_t. © by 1997 John Wiley & Sons, Ltd.     

Self-similar collapse of 2D and 3D vortex filament models

In this article, a very simple toy model for a candidate blow-up solution of the Euler equation by Boratav and Pelz (vortex dodecapole) is investigated. In this model, vortex tubes are replaced with straight vortex filaments of infinitesimal thickness, and the entire motion is monitored by tracing the motion of a representative point on one vortex filament. It is demonstrated that this model permits a self-similar collapse solution which provides the time dependence of the length scale as (t _c ? t)^1/2, (t < t _c), where the collapse time t _c depends on the initial configuration. From the conservation of circulation, this time dependence implies that vorticity ω scales as (t _c ? t) ^?1, which agrees with the one observed in the direct numerical (pseudo spectral) simulations of the vortex dodecapole. Finally, possible modification of the model is considered.     

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.     

Asymptotic Solutions of Problems of One-Dimensional Time-Dependent Combustible Gas Flows in the Presence of Thermal Action

In a development of studies [1, 2], asymptotic solutions of the Navier-Stokes equations are found for one-dimensional combustible gas flows in the presence of various forms of thermal action on a moving surface (x=x _w(t)). In the problems considered, the temperature or the heat flux q _w(t) is specified on the surface or the surface is the interface between a combustible gas and a moving heated piston or another gas (for example, in a shock tube). Use is made of the fact that, as t , in many cases the values of v _w=(dx/dt)_w and q _w are bounded. This leads to a steady-state flow in the flame zone in the coordinate system moving with its front and homogeneous uniform flow ahead of and behind it. Solutions of all these problems are given for the burnt-gas boundary layer region adjacent to the surface. The numerical calculations performed confirm the results obtained. A velocity law leading to time invariability of the flow pattern obtained with allowance for the interaction between the boundary layer and the burnt-gas homogeneous flow is found, including in the problem of the breakdown of an arbitrary discontinuity. The results are generalized to include the case of motion at an angle of incidence with an additional velocity component aligned with the surface.     

受载高聚物裂尖的损伤和银纹化

采用扫描电子显微镜(SEM)，对高聚物裂尖银纹损伤的引发和演化过程进行了原位观测．将固态高聚物本体材料视为线黏弹体，裂尖银纹区视为非线性损伤区，通过构造银纹区的损伤演化方程，给出了银纹区应力模型和银纹生长规律，数值结果与已有实验吻合良好。     

Entire Solutions with Merging Fronts to Reaction–Diffusion Equations

We deal with a reaction–diffusion equation u _t = u _xx + f(u) which has two stable constant equilibria, u = 0, 1 and a monotone increasing traveling front solution u = φ(x + ct) (c > 0) connecting those equilibria. Suppose that u = a (0 < a < 1) is an unstable equilibrium and that the equation allows monotone increasing traveling front solutions u = ψ₁(x + c ₁ t) (c ₁ < 0) and ψ₂(x + c ₂ t) (c ₂ > 0) connecting u = 0 with u = a and u = a with u = 1, respectively. We call by an entire solution a classical solution which is defined for all . We prove that there exists an entire solution such that for t≈ − ∞ it behaves as two fronts ψ₁(x + c ₁ t) and ψ₂(x + c ₂ t) on the left and right x-axes, respectively, while it converges to φ(x + ct) as t→∞. In addition, if c > − c ₁, we show the existence of an entire solution which behaves as ψ₁( − x + c ₁ t) in and φ(x + ct) in for t≈ − ∞.     

Feasibility study of three-dimensional PIV by correlating images of particles within parallel light sheet planes

One approach to obtain information about the out-of-plane velocity component from PIV recordings is to analyze the height of the peak in the correlation plane. This value depends on the portion of paired particle images, which itself depends on the out-of-plane velocity component and on other parameters. To circumvent problems with other influences (e.g. background light, amount and size of images), images from another light sheet plane parallel to the first one were also captured for peak height normalization. Our experimental results show the feasibility of an out-of-plane velocity estimation by analyzing images of particles within parallel light sheets by spatial cross-correlation.List of Symbols C particle density in the flow - d particle image diameter - f ₀, f ₁ frames containing images of particles within the first light sheet at t=t ₀ (frame f ₀) and at t=t ₀ + t (frame f ₁) - f ₂ frame containing images of particles within the second light sheet parallel to the first one at t=t ₀ + 2t - F ₁ estimator of the loss of image pairs due to in-plane motion - F ₀ estimator of the loss of image pairs due to out-plane motion - F convolution of the particle image intensity distributions - K factor containing constant parameters in the correlation plane - M imaging magnification (image size/object size) - n ₀ number of particles in the measurement volume at t=t ₀ - n _0,1 number of particle image pairs in interrogation windows of f ₀ andf ₁ - n _1,2 number of particle image pairs in interrogation windows off ₁ and f ₂ - O _z overlap of the light sheets - R _C (s) convolution of the mean intensity distributions - R _D (s) correlation which gives the image displacement - R _F (s) fluctuating noise component of the cross correlation estimator - R _0,1(s _D ) cross-correlation peak height of interrogation windows off ₀ and f ₁ - R _1,2(s _iuD) cross-correlation peak height of interrogation windows of f ₁ and f ₂ - s two-dimensional separation vector in the correlation plane - s _D mean particle image displacement in the interrogation cell - t _e light pulse duration - t _f frame-transfer time of the video camera - u three-dimensional local flow velocity vector (u,v,w) - X _i position of the center of an interrogation window in the image plane (2d) - x _i position of the center of an interrogation volume in the flow (3d) - (z ₂ — Z ₁) displacement of the light sheets in z-direction - t separation time of the light pulses - x ₀ x-extension of an interrogation volume - y ₀ y-extension of an interrogation volume - z ₀ light sheet thickness The authors would like to thank DLR for supporting Markus Raffel's and Olaf Ronneberger's visit to Caltech (Center for Quantitative Visualisation), and the Office of Naval Research through the URI grant ONR-URI-N00014-92-J-1610. Dr. Alexander Weigand's generous offer of his experimental set-up and stimulating discussions with Dr. Jerry Westerweel and Dr. Thomas Roesgen are greatly appreciated. Special thanks also to Dr. Christian Willert for his advice regarding the modifications to the DPIV software.     

Existence and Stability of the Inner Boundaries of the Domain of Multiple Hysteresis of Static Aerodynamic Forces and Moments

The results of experimental investigations of the multiple static hysteresis of the aerodynamic characteristics of a rectangular high-aspect-ratio wing are presented. Schematic wing-flow structure patterns, time dependences of the coefficients c _y(t), m _z(t), and m _x(t) and their frequency spectra obtained for a fixed model are given for different boundaries of the hysteresis domain. The time dependences of aerodynamic forces and moments are analyzed at angles of incidence at which sharp changes are observed. It is shown that the static hysteresis can be described by a mathematical model used in catastrophe theory.     

Low-thrust Lambert transfer based on two-stage constant-vector thrust control method

The low-thrust Lambert transfer refers to that the spacecraft achieves the orbital transfer whose boundary conditions are represented by two sets of orbital elements at initial and final time by the low-thrust propulsion system. The modulus and direction of the low-thrust solutions in previous methods change with time, which leads to high control requirements for the engine. In this paper, to reduce the requirements of the engine, a practical two-stage constant-vector thrust control method is proposed, in which the magnitude and direction of the thrust are deemed as segmental constant value in TNH frame, where three components of the thrust are f_t, f_n, and f_h. First, the mathematical model of the two-stage constant-vector thrust is formulated, and a rapid algorithm is presented to obtain the solution based on the linearized sensitivity matrix, which describes the relationship between the constant-vector thrust and the change of the orbital elements approximately. Furthermore, two low-thrust Lambert strategies based on the two-stage constant-vector thrust are presented for cases of short-time transfer and long-time transfer. A sequence of numerical simulations demonstrated the efficiency of the proposed approaches. The proposed control strategies are solved rapidly, and they are also suitable for different types of orbits with J₂ perturbation, which are practical options for engineering applications.

     

Global Semigroup of Conservative Solutions of the Nonlinear Variational Wave Equation

We prove the existence of a global semigroup for conservative solutions of the nonlinear variational wave equation u _tt − c(u)(c(u)u _x ) _x  = 0. We allow for initial data u| _{t = 0} and u _t | _t=0 that contain measures. We assume that 0 < k^-1 \leqq c(u) \leqq k{0 < \kappa^{-1} \leqq c(u) \leqq \kappa}. Solutions of this equation may experience concentration of the energy density (u_t²+c(u)²u_x²)dx{(u_t^2+c(u)^2u_x^2){\rm d}x} into sets of measure zero. The solution is constructed by introducing new variables related to the characteristics, whereby singularities in the energy density become manageable. Furthermore, we prove that the energy may focus only on a set of times of zero measure or at points where c′(u) vanishes. A new numerical method for constructing conservative solutions is provided and illustrated with examples.     

Predicting the failure of ultrasonic spot welds by pull-out from sheet metal

A methodology for determining the cohesive fracture parameters associated with pull-out of spot welds is presented. Since failure of a spot weld by pull-out occurs by mixed-mode fracture of the base metal, the cohesive parameters for ductile fracture of an aluminum alloy were determined and then used to predict the failure of two very different spot-welded geometries. The fracture parameters (characteristic strength and toughness) associated with the shear and normal modes of ductile fracture in thin aluminum alloy coupons were determined by comparing experimental observations to numerical simulations in which a cohesive-fracture zone was embedded within a continuum representation of the sheet metal. These parameters were then used to predict the load–displacement curves for ultrasonically spot-welded joints in T-peel and lap-shear configurations. The predictions were in excellent agreement with the experimental data. The results of the present work indicate that cohesive-zone models may be very useful for design purposes, since both the strength and the energy absorbed by plastic deformation during weld pull-out can be predicted quite accurately.