Application of the grating objective speckle method to composite materials

The grating objective speckle method has been applied in whole-field strain analysis of carbon-fiber-woven polyimide-composite materials. A sequence of rosette fringe patterns was obtained from one single specklegram, indicating displacement components along four different directions with a 45-deg interval. The spatial frequencies, which represent the sensitivities of the fringe intervals, were 2400 ℓ/mm forU _x ,U _y and 1697 ℓ/mm forU ₄₅ andU ₁₃₅. Normal strain components, ɛ _x , ɛ _y , ɛ₄₅, and ɛ₁₃₅, were trasformed to quantitative topographic representations by coupling to a digital computer system. Shear strain can be determined by the rosette equations.     

Full-field displacement and strain rosettes by moiré interferometry

A cross-line phase-type specimen grating was interrogated by moiré interferometry techniques to produce full-field fringe patterns of three displacement components:U _x ,U _y , andU ₄₅. ForU ₄₅, sensitivity was 1/1700 mm/fringe (1/43,000 in./fr). Closely packed fringes with much information content were obtained. Fringe patterns of strain components? _x ,? _y and? ₄₅ were produced by mechanical shearing (or optical differentiation) of the displacement patterns; a shearing distance of only 0.6 mm (0.025 in.) was used. The rosette method yields complete strain information from these three components of normal strain—which were derived from three direct derivatives of displacement. Consequently, the need for cross-derivatives of displacement, which are highly sensitive to accidental rigid-body rotations, is circumvented.     

High-magnification moiré interferometer for crack tip analysis of steels

A high-magnification moiré interferometer, particularly suitable for near-tip field analysis in cracked materials, is described. It has a submillimeter field of view, a high-resolution image sensor (1.4 million pixels), X-Y-Z translation stage and an optical fiber light delivery system. These features enable the microscope head to observe the crack tip while the specimen is loaded in a standard tensile test machine. Automated fringe pattern analysis, using temporal phase shifting and spatial phase unwrapping, enables thex ory displacement component to be measured and the corresponding in-plane strain component computed. The displacement placement accuracy is better than 40 nm, and the effective strain gage dimension is ∼ 25 μm. Furthermore, the interferometer has a built-in white light microscope that allows the observation of the specimen granular microstructure in exact registration with the displacement field. The interferometer has hence been employed to investigate the near-tip fields of a precracked stainless steel specimen under load. The influence of the grain boundaries on the measured displacement fields was relatively minor. The near-tip strain field shows a significant asymmetrical behavior despite pure mode lloading conditions.     

A photoelastic study of crack repair

Birefringent coatings have been employed to study the effectiveness of an adhesively bonded repair of a center-cracked tension panel. The repair was one sided, with photoelastic coatings applied to the opposite side. Photoelastic coatings were also applied over the patch. Analysis methods are presented to permit the stress intensity factor to be determined from the isochromatic fringe patterns recorded from both continuous andX- andY-edged coatings. The results showed that the one-sided adhesively bonded patch reduced the stress intensity factor; however, the repair did not markedly change the character of stress distributions. Fringe loops formed near the crack tips for both the cracked and repaired tension panels. The primary difference was in the size of the loops. The reduction inK _I due to repair was smaller than anticipated, but even small improvements in ΔK _I markedly enhance the life of a repaired panel. The Paris power law is used to show the relation between the reduction in ΔK _I and the improvement in the crack growth rateda/dN. Fringe patterns from the birefringent coatings applied to the patch provided a means not only to investigate the stresses in the patch but also to detect the initiation of the local debonding of the adhesive in the neighborhood of the crack. The birefringent coating on the patch is an approach for producing an optically “smart” repair.     

Experimental investigation of GH singularity field of the growing crack-tip for different ductile metals

The displacementsu _x ,u _y inx, y directions at growing crack tip of the specimens with double edge cracks for four different alumium alloys and two coppers are measured by using moire method and optical spatial filtering technique. From experimental displacement fields, the displacement singularity fields are obtained and compared with GH theoretical field. Unknown constantsA and _yo in theoretical solution are determined from experimental data. The theoretical singularity field thus compared is given for plane-stress, mode-I, strain-hardening materials. The error in both the experimental and theoretical evaluations is within ± 10%. The experiments show that there exists dominant singularity region around a growing crack-tip. In the experiments, the strain hardening indexn amounts from 3.158 to 14. The shape of this dominant region ranges from butterfly wing to oblate or circular shape. The size and shape of GH dominant region depend on the material property, the specimen geometry and loading type. Inside GH-field, there is a 3-D deformed damage zone, where no GH singularity exists. Very near to the crack-tip, there is a fracture process zone.The project is supported by the National Science Foundation of China.     

Symmetry of Mountain Pass Solutions of Some Vector Field Equations

We prove radial symmetry (or axial symmetry) of the mountain pass solution of variational elliptic systems − AΔu(x) + ∇ F(u(x)) = 0 (or − ∇.(A(r) ∇ u(x)) + ∇ F(r,u(x)) = 0,) u(x) = (u ₁(x),...,u _N (x)), where A (or A(r)) is a symmetric positive definite matrix. The solutions are defined in a domain Ω which can be , a ball, an annulus or the exterior of a ball. The boundary conditions are either Dirichlet or Neumann (or any one which is invariant under rotation). The mountain pass solutions studied here are given by constrained minimization on the Nehari manifold. We prove symmetry using the reflection method introduced in Lopes [(1996), J. Diff. Eq. 124, 378–388; (1996), Eletron. J. Diff. Eq. 3, 1–14].     

Experimental GH singularity field around growing crack-tip

The displacement fieldsu _x ,u _y at growing crack tip of LY12-M specimens with double edge cracks are measured using moire method. The experimental singularity fields are compared with GH theoretical field [12–14]. The size and shape of the experimental GH singularity fields are obtained. The error in both the experimental and theoretical evaluations is controlled within ±10%. The experiments show that there is singularity dominant around a growing crack tip. The shape of this dominant region ranges from butterfly wing to oblate and circular. Inside GH-field, there is a 3-D deformed damage zone where no GH singularity exists. The project suppoted by National Natural Science Foundation of China     

Elastic properties of a dilute surfactant solution in the shear induced state

The first normal stress differences N ₁ of a highly dilute cationic surfactant solution are investigated in a cone-and-plate rheometer. In continuation of a previous paper (Nowak 1998), where the buildup of a shear induced structure in such a solution was attained after a reduced deformation, the N ₁ turned out to be in proportion to the square of the shear rate γ˙ reduced by a critical value γ˙ _c in a first range above γ˙ _c . At higher shear rates the N ₁ tend to lower values than predicted by this relation. Relaxation experiments were performed in the same geometry to determine the characteristic time scales of the shear induced state's decay. In the lower range above &γdot; _c the stress decay is a monoexponential process, while a second time constant has to be introduced to describe the relaxation in that range, where the N ₁ deviate from the parabolic dependence of the reduced shear rate. Received: 10 May 1999 Accepted: 15 November 2000     

The high accuracy explicit difference scheme for solving parabolic equations 3-dimension

ＩｎｔｒｏｄｕｃｔｉｏｎＦｒｏｍｐｒａｃｔｉｃａｌｐｒｏｂｌｅｍ,ｗｅｃｏｕｌｄｃｏｎｃｌｕｄｅｍａｎｙｐｒｏｂｌｅｍｓａｂｏｕｔｓｏｌｖｉｎｇｐａｒａｂｏｌｉｃｐａｒｔｉａｌｄｉｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎ．Ｎｏｗｔｈｅｒｅａｒｅｍａｎｙｎｕ...     

A correction method for measuring turbulence kinetic energy dissipation rate by PIV

The accuracy of turbulent kinetic energy (TKE) dissipation rate measured by PIV is studied. The critical issue for PIV-based dissipation measurements is the strong dependency on the spatial resolution, Δx, as reported by Saarenrinne and Piirto (Exp Fluids Suppl:S300–S307, 2000). When the PIV spacing is larger than the Kolmogorov scale, η, the dissipation is underestimated because the small scale fluctuations are filtered. For the case of Δx smaller than the Kolmogorov scale, the error rapidly increases due to noise. We introduce a correction method to eliminate the dominant error for the small Δx case. The correction method is validated by using a novel PIV benchmark, random Oseen vortices synthetic image test (ROST), in which quasi-turbulence is generated by randomly superposing multiple Oseen vortices. The error of the measured dissipation can be more than 1,000% of the analytical dissipation for the small Δx case, while the dissipation rate is underestimated for the large Δx case. Though the correction method does not correct the underestimate due to the low resolution, the dissipation was accurately obtained within a few percent of the true value by using the correction method for the optimal resolution of η/10 < Δx < η/2.     

Electrorheological determination of the Leslie viscosity coefficient (α2) of a main-chain liquid crystal polymer in a nematic solvent

Analysis of the electrorheological response of thermotropic solutions of a main-chain liquid crystal polymer (MCLCP) in a low molar mass nematic solvent is performed at a fixed shear rate as a function of the applied field strength. The Leslie viscosity coefficient α₂ can be obtained by least squares fits to an equation describing the balance between the viscous and electric torques, formulated via the two dimensional Leslie-Ericksen-Parodi theory. We find that the increment Δα₂ on dissolving the MCLCP increases linearly with molecular weight, consistent with earlier light scattering measurements of the increment in the twist viscosity, Δγ₁, and also with previous electrorheological measurements of the increment in the Miesowicz viscosity Δη_c. Received: 1 December 1998 Accepted: 28 April 1999     

Effects of magnetic field,porosity, and wall properties for anisotropically elastic multi-stenosis arteries on blood flow characteristics

A mathematical model for blood flow through an elastic artery with multistenosis under the effect of a magnetic field in a porous medium is presented. The considered arterial segment is simulated by an anisotropically elastic cylindrical tube filled with a viscous incompressible electrically conducting fluid representing blood. An artery with mild local narrowing in its lumen forming a stenosis is analyzed. The effects of arterial wall parameters represent viscoelastic stresses along the longitudinal and circumferential directions T _t and T _θ , respectively. The degree of anisotropy of the vessel wall γ, total mass of the vessel, and surrounding tissues M and contributions of the viscous and elastic constraints to the total tethering C and K respectively on resistance impedance, wall shear stress distribution, and radial and axial velocities are illustrated. Also, the effects of the stenosis shape m, the constant of permeability X, the Hartmann number H _α and the maximum height of the stenosis size δ on the fluid flow characteristics are investigated. The results show that the flow is appreciably influenced by surrounding connective tissues of the arterial wall motion, and the degree of anisotropy of the vessel wall plays an important role in determining the material of the artery. Further, the wall shear stress distribution increases with increasing T _t and γ while decreases with increasing T _θ , M, C, and K. Transmission of the wall shear stress distribution and resistance impedance at the wall surface through a tethered tube are substantially lower than those through a free tube, while the shearing stress distribution at the stenosis throat has inverse characteristic through totally tethered and free tubes. The trapping bolus increases in size toward the line center of the tube as the permeability constant X increases and decreases with the Hartmann number Ha increased. Finally, the trapping bolus appears, gradually in the case of non-symmetric stenosis, and disappears in the case of symmetric stenosis. The size of trapped bolus for the stream lines in a free isotropic tube (i.e., a tube initially unstressed) is smaller than those in a tethered tube.     

Flow field characteristics in the near field region of particle-laden coaxial jets

The effects of solid particles on the flow structure in the near field region of a coaxial water jet are investigated non-intrusively using molecular tagging velocimetry. Glass beads of 240 μm and specific gravity SG of 2.46 are used at three volume loadings of γ_v=0.03, 0.06, and 0.09% in the central water jet with a Reynolds number of 4.1×10⁴. Measurements are acquired for four annular to central jet velocity ratios in the range 0.11≤ U _o/U _i≤1.15 at downstream distances up to six inner jet diameters and the results are analyzed for the effects of solid particles on the characteristics of flow. It is found that the addition of particles does not affect the mean fluid velocity profile in this region. The results also indicate a small and moderate enhancement of axial turbulent velocity and radial gradients of velocity fluctuations, respectively, due to the presence of particles.     

Computer simulated analyses on deformation and fracture of non-cracked and cracked specimens

A unified damage and fracture model, the combinatory work density model, which is suitable for either non-cracked body or cracked body has been suggested^[t−7]. In the present paper, the deformation and fracture of the two kinds of tensile spceimen and TPB specimen made of 40Cr steel have been simulated by using the new model together with the large elastic-plastic deformation finite element method. The results give a good picture of the whole deformation and fracture processes of the specimens in experiments; especially, the results on the TPB specimen can be used to obtain the relationship between load and displacement at the loading pointP-Δ, and between crack extension and displacement at the loading point Δa-Δ, the resistance curveJ _R -Δa and the fracture toughnessJ _1C . All the results are in remarkable agreement with those obtained by experiments. Therefore the model suggested here can be used to simulate crack initiation and propagation in non-cracked body and fracture initiation and crack stable propagation in cracked body. The project supported by National Natural Science Foundation of China     

Hydromagnetic free convection flow with induced magnetic field effects

An exact solution is presented for the hydromagnetic natural convection boundary layer flow past an infinite vertical flat plate under the influence of a transverse magnetic field with magnetic induction effects included. The transformed ordinary differential equations are solved exactly, under physically appropriate boundary conditions. Closed-form expressions are obtained for the non-dimensional velocity (u), non-dimensional induced magnetic field component (B _x ) and wall frictional shearing stress i.e. skin friction function (τ _x ) as functions of dimensionless transverse coordinate (η), Grashof free convection number (G _r ) and the Hartmann number (M). The bulk temperature in the boundary layer (Θ) is also evaluated and shown to be purely a function of M. The Rayleigh flow distribution (R) is derived and found to be a function of both Hartmann number (M) and the buoyant diffusivity parameter (ϑ ^*). The influence of Grashof number on velocity, induced magnetic field and wall shear stress profiles is computed. The response of Rayleigh flow distribution to Grashof numbers ranging from 2 to 200 is also discussed as is the influence of Hartmann number on the bulk temperature. Rayleigh flow is demonstrated to become stable with respect to the width of the boundary layer region and intensifies with greater magnetic field i.e. larger Hartman number M, for constant buoyant diffusivity parameter ϑ ^*. The induced magnetic field (B _x ), is elevated in the vicinity of the plate surface with a rise in free convection (buoyancy) parameter G _r , but is reduced over the central zone of the boundary layer regime. Applications of the study include laminar magneto-aerodynamics, materials processing and MHD propulsion thermo-fluid dynamics.     

Standing Waves for Nonlinear Schrödinger Equations with a General Nonlinearity

For elliptic equations ε²Δu − V(x) u + f(u) = 0, x ∈ R ^N , N ≧ 3, we develop a new variational approach to construct localized positive solutions which concentrate at an isolated component of positive local minimum points of V, as ε → 0, under conditions on f which we believe to be almost optimal. An erratum to this article can be found at     

Nonlinear rheology of concentrated spherical silica suspensions: 3. Concentration dependence

Nonlinear rheology was examined for concentrated suspensions of spherical silica particles (with radius of 40 nm) in viscous media, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture and pure ethylene glycol. The particles were randomly and isotropically dispersed in the media in the quiescent state, and their effective volume fraction φ_eff ranged from 0.36 to 0.59. For small strains, the particles exhibited linear relaxation of the Brownian stress σ_B due to their diffusion. For large step strains γ, the nonlinear relaxation modulus G(t,γ) exhibited strong damping and obeyed the time-strain separability. This damping was related to γ-insensitivity of strain-induced anisotropy in the particle distribution that resulted in decreases of σ_B/γ. The damping became stronger for larger φ_eff. This φ_eff dependence was related to a hard-core volume effect, i.e., strain-induced collision of the particles that is enhanced for larger φ_eff. Under steady/transient shear flow, the particles exhibited thinning and thickening at low and high γ˙, respectively. The thinning behavior was well described by a BKZ constitutive equation using the G(t,γ) data and attributable to decreases of a Brownian contribution, σ_B/γ˙. The thickening behavior, not described by this equation, was related to dynamic clustering of the particles and corresponding enhancement of the hydrodynamic stress at high γ˙. In this thickening regime, the viscosity growth η₊ after start-up of flow was scaled with a strain γ˙t. Specifically, critical strains γ_d and γ_s for the onset of thickening and achievement of the steadily thickened state were independent of γ˙ but decreased with increasing φ_eff. This φ_eff dependence was again related to the hard-core volume effect, flow-induced collision of the particles enhanced for larger φ_eff. Received: 26 June 1998 Accepted: 9 December 1998     

Hybrid stress analysis by digitized photoelastic data and numerical methods

A method is presented that determines photoelastic isochromatic values at the nodal points of a grid mesh which in turn is generated by a computer program that accepts digitized input. Values of σ₁ - σ₂ are computed from the digitized fringe orders. The Laplace equation is solved to separate the principal stresses at each nodal point. The method is extended to digitize isoclinics. Subsequently, σ _x - σ _y and τ _xy are calculated to be used as starting values for the solution of the pertaining partial differential equations to enhance convergence. For further accelerating the rate of convergence, superfluous boundary conditions are added from the digitized data; significant improvement is demonstrated. Estimated values of σ _x - σ _y from the digitized data are further used in conjunction with the solution of the Laplace equation to determine the state of stress without solving the boundary value problems. Paper was presented at the 1988 SEM Spring Conference on Experimental Mechanics held in Portland, OR on June 5–10.     

Photographic method to sharpen and double isochromatic fringes

A photographic method is described which gives, upon the superposition of ordinary light- and darkfield isochromatic fringe patterns, a new fringe pattern (mixed field). This mixed-field pattern has fringes at the (N/4 and 3N/4 positions. Use of the mixed-field fringe pattern coupled with the oridinary light- and dark-field patterns permit fringes to be read at theN/4,N/2, 3N/4 andN positions, and thus represents a factor of 2 increase in the number of countable fringes. The method is illustrated on two plane-stress examples, a three-dimensional example, and a case employing birefringent coatings.     

Immiscible Newtonian displacement by a viscoplastic material in a capillary plane channel

The immiscible displacement in a capillary plane channel of a Newtonian liquid by a viscoplastic one that obeys a Papanastasiou’s constitutive equation is numerically analyzed. An elliptic mesh generation technique, coupled with the Galerkin finite element method is used to determine the velocity field and the configuration of the interface between the two materials. We investigate the displacement efficiency and the flow patterns of the problem as functions of the dimensionless parameters that govern the problem: the capillary number (Ca), the viscosity ratio of the two fluids (N _η ) and the yield number, (τ′₀). The numerical results showed that for a fixed viscosity ratio, the fraction of mass attached to the wall is a decreasing function of τ′₀. We constructed maps of streamlines in the Cartesian space defined by τ′₀ and Ca for fixed viscosity ratios in order to capture the rough location of bypass and recirculating flow regimes. Higher yield number values induce bypass flow regimes, especially for high Ca. The dimensionless forms of the momentum conservation equation and the force balance at the interface were essential for the understanding of the role played by the dimensionless numbers that govern the problem.