Determination of stress-intensity factors from photoelastic data with applications to surface-flaw problems

ATaylor-series correction to the maximum inplane shear stress was studied as a means of extending the data zone in photoelastic determination of stress-intensity factors beyond the singular region of a two-degree-of-freedom analysis. Convergence properties were obtained by comparing with several complete two-dimensional solutions. Experiments were performed on two kinds of three-dimensional problems, plates containing surface flows in both bending and extension. Results were analyzed by both a two-degree-of-freedom and aTaylor-series correction method (TSCM). Results were compared to theories of F. W. Smith and A. S. Kobayashi and R. C. Shah. It was concluded that:

1. The TSCM program converges rapidly to accurateK _I values and will accommodate the scatter inherent in experimental data if the series is properly truncated.
2. The TSCM program is essentially equivalent to the two-parameter representation when only the crack-surface effects dominate.
3. When effects other than crack surfaces are important, TSCM requires more terms but still predictsK _I with reasonable accuracy.

     

Oscillation modes of laval nozzle flow

Experimental investigations of Laval nozzle flow show for relatively low supply to exit pressure ratios, which correspond to shock wave positions close to the nozzle throat, three different, oscillatory instabilities.

1. Shock pattern oscillations where the root of a λ-like shock front remains nearly in constant position, but where the proportion between the normal part and the oblique part of the shock changes periodically.
2. Shock wave and separation bubble oscillations where the motion of the shock wave is accompanied by displacements of the separation bubble.
3. Flow rate oscillations where the shock waves leave periodically through the nozzle throat in upstream direction.

     

Some pending problems in polymer melt rheology,as seen from the point of view of Doi's slip-link model

An exposition is given of results, as obtained with the aid of Doi's sliplink model, being considered as the most simple version of the famous “reptation model”. It turns out that this model which exhibits three distinct phases of relaxation (an extremely fast phase, an equilibration phase and a slow disengagement phase) is capable of explaining several peculiar features of polymer melt rheology:

1. The molecular mass dependence of the breadth of the rubber plateau in the storage modulus, of the zero-shear viscosity and of the normal-stress coefficients.
2. The molecular mass independence of the equilibrium (shear and tensile) compliances for monodisperse polymers (semi-quantitative prediction).
3. The seemingly contradictory sensitivity of these compliances for the breadth of the molecular mass distribution.
4. The critical value of the shear stress at which melt fracture occurs in capillary flow.
5. An equilibration phase in tensile experiments on unvolcanized rubber.

In this evaluation optical (flow birefringence) measurements are preferentially used.     

Marangoni convection at alcohol aqueous solutions-air interfaces

For aqueousn-heptanol solutions and in a nearly two-dimensional flow, two strikingfeatures have been detected:

1. a shift of the minimum of the surface tension
2. a discrepancy between the observed Marangoni flow velocities and the expected ones from static surface tension values.

A qualitative explanation is given.     

Flow past a suddenly cooled horizontal plate

The problem of unsteady laminar, incompressible free convection above a horizontal semi-infinite flat plate is studied theoretically. It is assumed that for timet<0 the plate is hotter than its surroundings and at timet=0 the plate is suddenly cooled to the same temperature of its surroundings. Three solutions of the momentum and energy equations are obtained, namely

1. an analytical solution which is valid for small time,
2. an asymptotic analytical solution which is valid for large time, and
3. a numerical solution which matches these two limiting analytical solutions.

It is found that the numerical solution matches the small and large time solutions accurately. Finally, the variation of the velocity, temperature, skin friction and heat transfer on the plate with time are discussed.     

A logrithmic bound on the location of the poles of the scattering matrix

It is shown that the complex poles z of the scattering matrix satisfy the inequality: Im z≧a+b log ¦z¦, b>0, in three instances of classical scattering in three space dimensions described by the wave equation u_{t t}?c²Δu+qu=0.

1. c and q smooth with c=1 and q=0 for ¦x¦>p, all rays going to infinity, and the energy form positive definite.
2. c=1 and q=0 outside of a convex body on which u=0.
3. c=1, q bounded and measurable, q=0 for ¦x¦>p, and the energy form not necessarily positive definite.

     

The initial-value problem for elastodynamics of incompressible bodies

We prove short-time well-posedness of the Cauchy problem for incompressible strongly elliptic hyperelastic materials. Our method consists in:

1. Reformulating the classical equations in order to solve for the pressure gradient (The pressure is the Lagrange multiplier corresponding to the constraint of incompressibility.) This formulation uses both spatial and material variables.
2. Solving the reformulated equations by using techniques which are common for symmetric hyperbolic systems. These are:

1. Using energy estimates to bound the growth of various Sobolev norms of solutions.
2. Finding the solution as the limit of a sequence of solutions of linearized problems.

Our equations differ from hyperbolic systems, however, in that the pressure gradient is a spatially non-local function of the position and velocity variables.     

Guidance properties of nonlinear planar waveguides

We discuss the propagation of electromagnetic waves through a stratified dielectric. The ability of such a device to support guided waves depends upon the way in which the refractive index varies across the layers. In the present discussion, we show how nonlinear effects and appropriate stratification can be used to obtain any one of the following behaviours:

1. guidance occurs only at low power.
2. guidance occurs only at high power.
3. guidance occurs at all powers.
4. there is no guidance.

The situation (i) is obtained by using materials with a defocusing dielectric response, whereas the situation (ii) is obtained for suitable configurations of self-focusing materials. The situations (iii) and (iv) can be obtained by using either defocusing or self-focusing materials. By seeking solutions of a particular form, we reduce the problem to the study of solutions in the Sobolev space H ¹(?) of a second-order differential equation. The discussion of defocusing nonlinearities is based in the study of the global behaviour of the branch of solutions bifurcating from a simple eigenvalue. For self-focusing nonlinearities we use a variational approach.     

Variable isochromatic/isopachic-fringe visibility in photoholoelasticity

This paper presents a new “hybrid” method whereby the ratio of the isochromatic-fringe visibility/isopachic-fringe visibility may be easily and continuously varied. This simple procedure merely combines a conventional polariscope with a holographic system. A variable beam splitter permits an incoherent superposition of the reconstruction of a doubly exposed hologram with real-time isochromatics, either dark or light field. By varying the ratios of the above two, in the image plane, numerous interesting results may be obtained including:

1. Isochromatics only, without errors in position
2. Isochromatics-isopachic fringes identical to those obtained through classical interferometry
3. Isochromatic-isopachic fringes whereby the amplitude modulation between the two may be minimized
4. Continuously variable isopachic/isochromatic-fringe visibility.

     

Basic equations for the microscopic theory of plasticity of polycrystalline materials with given texture

An expression for the yield stress of anisotropic materials is applied to the anisotropic strength of hard rolled copper foils whose crystallographic texture is known. We assume that this crystallographic texture is the only cause of the anisotropic plastic behaviour of the material. The model used for the yield stress is also used to deduce:

1. Stress-strain relations for isotropic polycrystalline materials;
2. A formula for the fully plastic strain tensor, applied to anisotropic hard rolled copper foils.

For the anisotropic copper foils considered the calculated curves of the yield stress and of the strain tensor as a function of the angle x between rolling and tensile direction agree qualitatively with the measured values. However, the theory is not complete, since the yield stress and the plastic strain tensor are both a function of a parameter Q, the fraction of the number of available crystallographic slip planes on which the maximum shear stress has reached the critical value τ_a. We assume that for “fully” plastic deformation a certain critical fraction Q _e of the total number of slip planes has to be active. The fraction Q _e is called the critical active quantity. With the parameter Q _e we adjust the calculated curves to the measured ones. The dependence of Q _e on the properties of the material (e.g. the crystallographic texture) is discussed in Appendix I.     

Pressure drop in alternating curved tubes

Pressure drop measurements in the laminar and turbulent regions for water flowing through an alternating curved circular tube (x=h sin 2πz/λ) are presented. Using the minimum radius of curvature of this curved tube in place of that of the toroidally curved one in calculating the Dean number (ND=Re(D/2R _c )², it is found that the resulting Dean number can help in characterizing this flow. Also, the ratio between the height and length of the tube waves which represents the degree of waveness affects significantly the pressure drop and the transition Dean number. The following correlations have been found:

1. For laminar flow: $$F_w \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} = F_s \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} + 0.03,\operatorname{Re}< 2000.$$
2. For turbulent flow: $$F_w \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} = F_s \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} + 0.005,2000< \operatorname{Re}< 15000.$$
3. The transition Dean number: $$ND_{crit} = 5.012 \times 10^3 \left( {\frac{D}{{2R}}} \right)^{2.1} ,0.0111< {D \mathord{\left/ {\vphantom {D {2R_c }}} \right. \kern-\nulldelimiterspace} {2R_c }}< 0.71.$$

     

Flow of a dipolar fluid between parallel plates

The flow of a dipolar fluid between two parallel plates with and without heat transfer is studied. The following cases are discussed:

1. Isothermal flow due to the relative motion of the plates,
2. Isothermal flow due to a constant pressure gradient with the plates at rest,
3. Nonisothermal flow with linearly varying plate temperatures.

Case (ii) is of particular interest to the experimentalists as it shows the effect of the material constants even when there are no externally applied dipolar tractions on the plates.     

Effects of artificial cracks and poisson's ratio upon photoelastic stress-intensity determination

A series of stress-freezing photoelastic experiments were performed with multiple replications upon edge-cracked strips for three types of “cracks” in current use:

1. Rectangular slots 0.152 mm wide,
2. 1.59-mm-wide slots terminating in a 30-deg vee notch of approximately 0.025-mm root radius, and
3. Natural cracks (approximately 0.0025-mm root radius).

Stress-intensity results were compared with the Gross-Srawley analysis; in addition (1) was compared with Savin's solution. It was concluded that (2) and (3) yield the same results but (1) was slightly higher. Both (2) and (3) were about 12 percent higher than the Gross-Srawley results. This is shown to be related to a Poisson's ratio effect.     

C 1 Solutions for Semi-Implicit Systems of Differential Equations

The present note is a continuation of the author??s effort to study the existence of continuously differentiable solutions to the semi-implicit system of differential equations (1) $$f(x^{\prime}(t)) = g(t, x(t))$$ (2) $$\quad x(0) = x_0,$$ where

${\quad\Omega_g \subseteq \mathbb{R} \times\mathbb{R}^n}$ is an open set containing (0, x ₀) and ${g:\Omega_g \rightarrow\mathbb{R}^n}$ is a continuous function,

${\quad\Omega_f \subseteq \mathbb{R}^n}$ is an open set and ${f:\Omega_f\rightarrow\mathbb{R}^n}$ is a continuous function.

The transformation of (1)?C(2) into a solvable explicit system of differential equations is trivial if f is locally injective around an element ${\gamma\in \Omega_f\cap f^{-1}(g(0,x_0))}$ . In this paper, we study (1)?C(2) when such a translation is not possible because of the inherent multivalued nature of f ^?1.     

A fast and accurate graphical method for designing an air conditioning cooling coil

The aim of this paper is to develop a fast and simple accurate graphical method for designing the required cooling coil for an air conditioning system in which both the sensible and latent heat are transferred. The method has the following advantages:

1. Direct solution utilizing only Psychrometer charts.
2. Solution of combined heat and mass transfer problems occuring at pressures other than that of atmosphere.
3. Direct determination of boundary temperature at which dehumidification begins for that type of problem where a portion of the surface is in a dry condition.
4. Less effort required than trial and error method in determination of air condition leaving a counterflow coil of a given area.

The calculated values of the method showed a good agreement with the experimental results. The average deviation for the total heat is about + 15 % and for sensible heat is + 9 %, which are of positive nature and on the safe side for practical design purposes. The method is also useful for similar practical application.     

The interaction of creep and fatigue for a rotor steel

Twenty tests were performed on a 1 Cr?1 Mo?1/4 V rotor steel at 1000° F (538°C) to determine the interaction of creep and low-cycle fatigue. These tests involved five different types of strain-controlled cycling: creep at constant tensile stress; linearly varying strain at different frequencies; and hold periods at maximum compressive strain, maximum tensile strain, or both. The experimental data were then used to characterize the interaction of creep and fatigue by the:

1. Frequency-modified strain-range approach of Coffin;
2. Total time to fracture vs. the time of one cycle relation as proposed by Conway and Berling;
3. Total time to fracture vs. the number of cycles to fracture characterization of Ellis and Esztergar;
4. Summation of damage fractions obtained from tests using interspersed creep and fatigue as proposed by the Metal Properties Council;
5. Strain-range-partitioning method of Manson, Halford, and Hirschberg.

In order to properly assess the strain-range-partitioning approach, seven additional tests were performed at the NASA Lewis Research Center. Visual, ultrasonic, and acoustic-emission methods of crackinitiation determination were unsuccessful. An approximate indication of crack initiation was obtained by finding the cycle N_o where the stress-cycle curve first deviated from a constant slope. Predictive methods (based on monotonic tests) for determining the fatigue life in the creep range were examined and found deficient, though they may still be useful for preliminary comparison of materials and temperatures. The extension of the frequency-modified strain-range approach to notched members was developed and the results of notched-bar tests were shown to corroborate this approach, when crack initiation for the plain and notched bars was campared.     

Dynamics as a mechanism preventing the formation of finer and finer microstructure

We study the dynamics of pattern formation in the one-dimensional partial differential equation $$u_u - (W'(u_x ))_x - u_{xxt} + u = 0{\text{ (}}u = u(x,t),{\text{ }}x \in (0,1),{\text{ }}t > 0)$$ proposed recently by Ball, Holmes, James, Pego & Swart [BHJPS] as a mathematical “cartoon” for the dynamic formation of microstructures observed in various crystalline solids. Here W is a double-well potential like 1/4((u _x )² ?1)². What makes this equation interesting and unusual is that it possesses as a Lyapunov function a free energy (consisting of kinetic energy plus a nonconvex “elastic” energy, but no interfacial energy contribution) which does not attain a minimum but favours the formation of finer and finer phase mixtures: $$E[u,u_t ] = \int\limits_0^1 {(\frac{{u_t^2 }}{2} + W(u_x ) + \frac{{u^2 }}{2})dx.}$$ Our analysis of the dynamics confirms the following surprising and striking difference between statics and dynamics, conjectured in [BHJPS] on the basis of numerical simulations of Swart & Holmes [SH]:

?While minimizing the above energy predicts infinitely fine patterns (mathematically: weak but not strong convergence of all minimizing sequences (u _nv_n) of E[u,v] in the Sobolev space W ^{1 p}(0, 1)×L²(0,1)), solutions to the evolution equation of ball et al. typically develop patterns of small but finite length scale (mathematically: strong convergence in W ₁ ^p(0,1)×L²(0,1) of all solutions (u(t),u_t(t)) with low initial energy as time t → ∞).

Moreover, in order to understand the finer details of why the dynamics fails to mimic the behaviour of minimizing sequences and how solutions select their limiting pattern, we present a detailed analysis of the evolution of a restricted class of initial data — those where the strain field u _x has a transition layer structure; our analysis includes proofs that

?at low energy, the number of phases is in fact exactly preserved, that is, there is no nucleation or coarsening

?transition layers lock in and steepen exponentially fast, converging to discontinuous stationary sharp interfaces as time t → ∞

?the limiting patterns — while not minimizing energy globally — are ‘relative minimizers’ in the weak sense of the calculus of variations, that is, minimizers among all patterns which share the same strain interface positions.

     

Glass-bonding techniques for semiconductor strain gages

Conventional organic-epoxy adhesives outgas when exposed to ultra-high vacuum and, as operating temperatures are increased, they begin to exhibit plastic behavior causing hysteresis and zero instability in the transducer. The use of an inorganic glass as the bonding material has resulted in a significant advance in transducer-fabrication technology for the following reasons:

1. The outgassing of transducers in high-vacuum applications is minimized.
2. Mechanical properties of the transducer such as hysteresis and repeatability are improved.
3. The electrical isolation of the strain gages from the metallic elements of the transducer is increased at high temperatures over that provided by epoxy. Also, the glass bond can survive and operate in severe radiation environments, wherein the epoxy adhesive will suffer either temporary or permanent loss of its dielectric strength.
4. Glass-bonding techniques are particularly useful for the extension of the temperature range of operation of silicon-strain-gage transducers.

Nispan C and 440-C stainless-steel substrates were successfully used with glass-bonded silicon strain gages to fabricate transducers for evaluation.     

Rapid interpretation of moiré photographs

The paper describes the possibilities of rapidly and conveniently obtaining results by means of the reflection moiré method. The moiré photographs are produced with the aid of a screen provided with a finely ruled grid. The contrast of the photographs is increased by optical means. By displacement of two equal moiré photographs, moiré fringes of second order for the curvature or the twist are obtained. It is shown how charts can be produced which indicate:

1. The distribution of the bending moments
2. The distribution of the bending moments for which the reinforcement should be designed
3. The distribution of the shear forces and the magnitude of the bearing reactions.

     

Interferogram processing by spatial filtering

This paper describes a technique of interferogram processing by the use of spatial-frequency filtering. Proper part of Fourier spectrum of the interferogram is blocked out by a matched mask located in the focal plane of the imaging system. Depending on the mask shape and dimensions one can obtain

1. interferogram images with delineating fringes or
2. interferogram images of one component of strain.