Thermal residual stresses in thick graphite/epoxy composite laminates—Uniaxial approach

Thermal residual stresses have been known to be very large in laminates of continuous-fiber-reinforced polymer composites. When the thickness of the laminate is large, however, the measurement of the residual stresses raises questions on the accuracy of the conventional methods. A novel concept of layer separation is developed to measure quantitatively and precisely the tensile residual stresses in thick plates with layered distribution of residual stresses. It is applied to thick [O₂/9O₄]_13s, AS/3501-6 graphite/epoxy laminates. The test specimens were mechanically modeled into the thin strips for the application of the new concept of layer separation. The tensile residual stresses measured in the 90-deg layers of these laminates are nonuniform throughout the specimen, and vary from 55.6 MPa to 71.4 MPa. It is very interesting to compare these values with the transverse strengthF ₂ ^tu of AS/3501-6 unidirectional composites, which is 65.4 MPa.     

Least-squares finite element methods for compressible Euler equations

A method based on backward finite differencing in time and a least-squares finite element scheme for first-order systems of partial differential equations in space is applied to the Euler equations for gas dynamics. The scheme minimizes the L²-norm of the residual within each time step. The method naturally generates numerical dissipation proportional to the time step size. An implicit method employing linear elements has been implemented and proves robust. For high-order elements, computed solutions based on the L²-method may have oscillations for calculations at similar time step sizes. To overcome this difficulty, a scheme which minimizes the weighted H¹-norm of the residual is proposed and leads to a successful scheme with high-degree elements. Finally, a conservative least-squares finite element method is also developed. Numerical results for two-dimensional problems are given to demonstrate the shock resolution of the methods and compare different approaches.     

混凝土框架结构内力测量传感器研制

研制了一种可用于混凝土结构试验中测量柱底截面轴向力、剪力和弯矩的内力测量传感器。该传感器主要由加载板、基座、4根竖向测力杆和2根水平测力杆组成。使用时,柱底截面的内力通过加载板传递给6根测力杆,根据6根测力杆的结果,通过力平衡方程可以计算实际加载的轴向力、剪力和弯矩的大小。对该传感器进行了4种工况下的标定试验。结果表明,传感器的测量误差满足框架结构内力测量精度要求。将该传感器应用到单层两跨的钢筋混凝土框架推覆试验中,获得了整个试验过程中构件的内力时程。     

A stable least-squares finite element method for non-linear hyperbolic problems

A class of stable least-square finite element methods for non-linear hyperbolic problems is developed and some exploratory studies made. The methods are based on modifying the L²-norm of the. residual and a related approximation to the H¹-norm of the residual. The effect of the additional terms in these residual functionals is to introduce a dissipative effect proportional to the solution gradient. This acts to stabilize the solution for non-linear hyperbolic problems which generate shocks. Numerical results for a one-dimensional nozzle and shock tube problem demonstrate the accuracy and stability of the method. Results are for an implicit scheme and calculations for linear, quadratic and cubic elements are given.     

Molecular theory for moderately concentrated polymer solutions in shear flow

A molecular theory for the rheological properties of moderately concentrated polymer solutions is developed on the basis of a model of interacting dumbbells. The interaction is treated in a mean field approximation, leading to an effective one-particle potential and a Gaussian stationary distribution function. Various rheological functions such as birefringence, shear viscosity and first normal-stress coefficient for simple shear flow and the Trouton viscosity for simple extensional flow are calculated. Good qualitative agreement with experimental observations is found, especially at intermediate flow rates. It is predicted, for example, that the birefringence increases approximately linearly with shear rate at intermediate shear rates and that the concentration dependence of the gradient varies asc ^1/2. The typical non-Newtonian behaviour is obtained for the shear viscosity. For small concentrations the onset of shear rate dependence decreases asc ^–1/2. At intermediate shear rates an apparent power law is obtained with an exponent between – 0.5 and – 1.0, decreasing with concentration.     

Shear strains in a graphite/PEEK beam by moiré interferometry with carrier fringes

A multispan quasi-isotropic graphite-PEEK beam exhibited dramatic shear strains in the interlaminar region between plies. Shear strains in the plies themselves varied in basic accord with fiber direction. The large anomalous shear strains were developed near the center of the beam height, where shear stresses were large. High-sensitivity moiré interferometry with 2400l/mm (60,960l/in.) was used. A new technique of data extraction was developed, using carrier fringes to transform the pattern to one in which fringe slopes are proportional to derivatives of displacement. This technique enhanced detection and measurement of highly localized shear-strain gradients. Paper was presented at the 1986 SEM Fall Conference on Experimental Mechanics held in Keystone, CO on Nov. 2–5.     

A biaxial stress transducer for fabrics

A transducer^* is developed for measurement of biaxial stresses in fabrics and other flexible sheet materials. The stress transducer responds directly to the far-field stresses by means of strains induced in an elastic inclusion which are monitored by strain gages. Design of the transducer is based on the analytical solution for such inclusions in linear-anisotropic materials, but the concept is applied and tested in nonlinear-orthotropic structural fabrics. Calibration and verification tests are conducted for principal far-field stresses in the yarn directions in vinyl-coated polyester and tefloncoated fiberglass.     

Acoustoelastic determination of stress in slightly orthotropic materials

A method is proposed to determine stresses in acoustoelasticity by making use of orthotropic stress-acoustic relations and the equations of equilibrium. It is derived theoretically that shear stress is determined ny ultrasonic data ofB and ?, which denote a magnitude of acoustic birefringence and its principal direction, respectively. Other stress components are obtained by numerical integration of the equilibrium equation with the shear stress thus determined. Experiments were carried out to show the validity and usefulment of the method. This method was applied to the measurement of stress field on a plate with a circular hole subjected to axial tension. Ultrasonic measurements were made by a Y-cut quartz transducer with 5-MHz fundamental frequency. The specimen was cut out from 1100 aluminum plate of 4-mm thickness, which shows a slight orthotropy due to roll working. The values ofB and ? were measured in both stressed and unstressed state. Then, stress distributions were determined by the method proposed here, and are compared with the known theoretical distributions.     

Determination of the dynamic stress intensity factorsK I d andK II d,for a mixed-mode propagating crack

In this paper, the dynamic propagation problem of a mixed-mode crack was studied by means of the experimental method of caustics. The initial curve and caustic equations were derived under the mixed-mode dynamic condition. A multi-point measurement method for determining the dynamic stress intensity factors,K _I ^d , andK _II ^d , and the position of the crack tip was developed. Several other methods were adopted to check this method, and showed that it has a good precision. Finally, the dynamic propagating process of a mixed-mode crack in the three-point bending beam specimen was investigated with our method.     

Improvement in the frequency response of the electrochemical wall shear stress meter

One of the very few disadvantages of the mass-transfer transducer when compared with the hot-film sensor, is a slightly diminished frequency-response due to the higher Prandtl number encountered. Mass or thermal balance and transfer equations were solved first by Fortuna and Hanratty (1971) for small fluctuations of the wall shear. The solutions allow to make accurate corrections on the frequency spectra and the power of the fluctuations, but in different time. In this paper, the author deduces the frequency response of split rectangular electrodes and shows how a combination of signals improves the response at higher frequencies and makes it comparable to the thermal transducer with the same size, in the same fluid. Two experimental devices are described and compared. With these devices, the measurement of the wall shear fluctuations is improved in real time. Accurate determinations of turbulent power fluctuations and probability density spectra are feasible and illustrate the subject.List of symbols A total area of the electrode - A _j area of the part j of the electrode - a coefficient - C concentration - C bulk concentration - c fluctuation of concentration - D diffusion coefficient - F Faraday's constant - f(n ⁺) transfer function - g gain of the differential electrode - I _j electrolysis current on the part j of the electrode - K transfert coefficient - k fluctuation of K - l electrode length - n frequency - P _r Prandtl number - S wall shear - s fluctuation of the wall-shear - t time - x direction of the flow - y direction normal to the wall - phase delay - v kinematic viscosity A version of this paper was presented at the 11th Symposium on Turbulence, University of Missouri-Rolla, 17–19 October 1988     

Non-Newtonian flow over a wedge with suction

A pseudo-similarity solution is obtained for the flow of an incompressible fluid of second grade past a wedge with suction at the surface. The non-linear differential equation is solved using quasi-linearization and orthonormalization. The numerical method developed for this purpose enables computation of the flow characteristics for any values of the parameters K, a and b, where K is the dimensionless normal stress modulus of the fluid, a is related to the wedge angle and b is the suction parameter. A significant effect of suction on the wall shear stress is observed. The present results match exactly those from an earlier perturbation analysis for Kx^2a ? 0·01 but differ significantly as Kx^2a increases.     

Deformation localization and shear band fracture in strong anisotropy sheet tension

The tensile deformation localization and the shear band fracture behaviors of sheet metals with strong anisotropy are numerically simulated by using Updating Lagrange finite element method, Quasi-flow plastic constitutive theory^[1] and B-L planar anisotropy yield criterion^[2]. Simulated results are compared with experimental ones. Very good consistence is obtained between numerical and experimental results. The relationship between the anisotropy coefficientR and the shear band angle θ is found. The project supported by the National Natural Science Foundation of China and the Excellent Youth Teacher Foundation of the State Education Commission of China     

Tailoring turbulence with an active grid

Using an active grid in a wind tunnel, we generate homogeneous shear turbulence and initiate turbulent boundary layers with adjustable properties. Homogeneous shear turbulence is characterized by a constant gradient of the mean velocity and a constant turbulence intensity. It is the simplest anisotropic turbulent flow thinkable, and it is generated traditionally by equipping a wind tunnel with screens which have a varying transparency and flow straighteners. This is not done easily, and the reachable turbulence levels are modest. We describe a new technique for generating homogeneous shear turbulence using an active grid only. Our active grid consists of a grid of rods with attached vanes which can be rotated by servo motors. We control the grid by prescribing the time-dependent angle of each axis. We tune the vertical transparency profile of the grid by setting appropriate angles of each rod such as to generate a uniform velocity gradient, and set the rods in flapping motion around these angles to tailor the turbulence intensity. The Taylor Reynolds number reached was R _λ = 870, the shear rate S = ∂U/∂y = 9.2 s⁻¹, the nondimensional shear parameter S ^*≡ Sq ²/ε = 12 and u = 1.4 ms⁻¹. As a further application of this idea we demonstrate the generation of a simulated atmospheric boundary layer in a wind tunnel which has tunable properties. This method offers a great advantage over the traditional one, in which vortex-generating structures need to be placed in the wind tunnel to initiate a fat boundary layer.     

New laboratory tests to measure rheological properties of paper coatings in transient and steady-state flows

The paper describes pertinent laboratory tests to characterize the rheological properties of paper coatings with regard to blade coating over a very wide range of shear rates in both transient and steady-state shear flows. Shear rates as high as 10⁶ s^–1 can be reached by means of a gas-driven capillary rheometer. Examples for the evaluation of end effects, wall effects, and coating thixotropy are given. A stiff and fast Couette rheometer is used to determine flow curves and the shear stress overshoot in step shear rate tests. The primary normal stress difference can be measured up to 10⁴ s^–1 by means of a high shear cone-plate rheometer with piezo transducer. A correct evaluation of the measurements has to take into account inertia contributions to the normal force. First results using a sinusoidal modulation of the shear rate are presented.Paper presented at: International Symposium on Pigment Coating Structure and Rheology, Helsinki, Febr. 8–9, 1989     

Study of magnetorheological fluids at high shear rates

The tunable rheological properties of magnetorheological (MR) materials at high shear rates are studied using a piston-driven flow-mode-type rheometer. The proposed method provides measurement of the apparent viscosity and yield stress of MR fluids for a shear rate range of 50 to 40,000 s⁻¹. The rheological properties of a commercial MR fluid, as well as a newly developed MR polymeric gel, and a ferrofluid-based MR fluid are investigated. The results for apparent viscosity and dynamic and static shear stresses under different applied magnetic fields are reported.     

Boundary Regularity of Shear Thickening Flows

This article is concerned with the global regularity of weak solutions to systems describing the flow of shear thickening fluids under the homogeneous Dirichlet boundary condition. The extra stress tensor is given by a power law ansatz with shear exponent p≥ 2. We show that, if the data of the problem are smooth enough, the solution u of the steady generalized Stokes problem belongs to W^{1,(np+2-p)/(n-2)}(W){W^{1,(np+2-p)/(n-2)}(\Omega)} . We use the method of tangential translations and reconstruct the regularity in the normal direction from the system, together with anisotropic embedding theorem. Corresponding results for the steady and unsteady generalized Navier–Stokes problem are also formulated.     

A torsional creep instrument and a torsional pendulum,both with accurate temperature control

Summary This paper describes a torsional pendulum and a torsional creep instrument. With the pendulum shear moduli between 10⁶ and 10¹⁰ N/m² can be measured at frequencies from 0.1 to 20 Hz. The creep instrument is suitable for measurement of shear compliances lower than 10^–7 m²/N in the time range from 1 to 10⁵ seconds. In both instruments, specimens are kept at the right temperature by blowing heated nitrogen gas through a surrounding thermostatic chamber. The signal of a platinum resistance thermometer, provided in each chamber, automatically controls the heating of the gas. Temperatures from –180 to +300 °C can be maintained with an absolute accuracy of ±1 °C and a long term stability of ±0.05 °C. It is shown that one cannot directly compare one and the same shear property, calculated from the shear modulus as measured with the pendulum as well as from the shear compliance as measured with the creep instrument. This is due to differences in the temperature of one thermostatic chamber over against the other. Finally, the paper presents a method to reduce these differences to ±0.1 °C, although the absolute accuracy of temperature control remains ±1 °C.

---

Zusammenfassung Die Arbeit beschreibt ein Torsionspendel und eine Torsions-Kriechapparatur. Mit dem Pendel können Schermoduli zwischen 10⁶ und 10¹⁰ N/m² bei Frequenzen zwischen 0,1 und 20 °C gemessen werden. Die Kriechapparatur ist geeignet für die Messung von Scherkomplianzen kleiner als 10^–7 m²/N in Zeiten zwischen 1 bis 10⁵ sec. In beiden Geräten werden die Proben durch das Einblasen von erhitztem Stickstoff durch eine umgebende thermostatische Kammer bei der richtigen Temperatur gehalten. Die Anzeige eines Platin-Widerstandsthermometers, das in jeder Kammer angebracht ist, kontrolliert automatisch die Erwärmung des Gases. Es können Temperaturen zwischen –180 und +300 °C mit einer absoluten Genauigkeit von ±1 °C und einer Langzeitstabilität von ±0,05 °C eingestellt werden. Es wird gezeigt, daß ein direkter Vergleich der gleichen Schereigenschaft, die zum einen aus dem mit dem Pendel gemessenen Schermodul und zum anderen aus der mit der Kriechapparatur ermittelten Kriechkomplianz errechnet wird, nicht möglich ist. Das beruht auf Temperaturdifferenzen zwischen den thermostatisierten Kammern.Abschließend stellt die Arbeit eine Methode vor, um diese Differenzen auf ±0,1 °C zu senken. Die absolute Genauigkeit der Temperaturkontrolle bleibt bei ±1 °C.

---

---

Paper presented at the Conference on Experimental Rheology, University of Bradford, April 17–19, 1968.     

Flow of viscoelastic fluids through an abrupt contraction

Summary Developing and fully developed velocity profiles were measured for viscoelastic fluids flowing through an abrupt 2 to 1 glass-contraction. An R 16Weissenberg-Rheogoniometer was used to measure the rheological properties of the viscoelastic fluids in the shear rates range of interest in the contraction. The measured entry lengths for the viscoelastic fluids were significantly less than predictions and experimental values for inelastic fluids with the same power-law parameters. Deviations from inelastic entry behaviour ranged from 11.6–100%, were independent ofReynolds number, but were strongly dependent on the ratio of the friction velocity to the shear wave velocity. Increasing the friction velocity relative to the shear wave velocity resulted in an increased deviation from inelastic behaviour. When the friction velocity was of the same order as the shear wave velocity a zero entry length and a fully developed entry velocity profile were observed. Further increase in the friction velocity relative to the shear wave velocity resulted in anomalous entry behaviour accompanied by unusual flow patterns upstream of the contraction.

---

Zusammenfassung Es wurden sich bildende sowie voll ausgebildete Geschwindigkeitsprofile viskoelastischer Flüssigkeiten in einer scharfkantigen Rohrverengung von 2 zu 1 gemessen. Ein Weissenbergsches Rheogoniometer R 16 diente zur Charakterisierung der viskoelastischen Flüssigkeiten im betreffenden Deformationsgeschwindigkeitsbereich.Meßergebnisse für die Einlauflänge viskoelastischer Flüssigkeiten weichen bedeutend von den Voraussagen sowie von Meßergebnissen für unelastische Flüssigkeiten ab, die, nach demOstwald- de Waeleschen Modell berechnet, die gleichen Kenngrößen aufzeigen.Die Abweichung vom viskosen Einlaufverhalten beträgt 11,6 bis 100%. Sie ist unabhängig von der Reynoldschen Zahl, hängt aber sehr stark ab von dem Verhältnis zwischen zwei Geschwindigkeiten u^*=Schubspannungsgeschwindigkeit undu=Scherwellengeschwindigkeit.Eine Erhöhung vonu ^* gegenüberu verursacht eine erhöhte Abweichung vom unelastischen Verhalten. Wenn die zwei Geschwindigkeitenu ^* undu von der gleichen Größenordnung sind, verschwindet die Einlaufsentwicklung und ein vollausgebildetes Geschwindigkeitsprofil tritt schon am Eingang auf. Ein weiteres Erhöhen vonu ^* überu verursacht anomales Einlaufverhalten mit ungewöhnlichem Strömungsbild oberhalb der Verengung.

---

---

On Sabbatical Leave: Dept. of Chemical Engineering, University of Toronto, Toronto 181, Ontario, Canada.     

Tests of a fiber-optic velocity transducer

A novel transducer is developed and tested. The transducer utilizes optical fiber to measure mean and instantaneous flow rates in turbulent flows, and is capable of detecting flow reversal. Calibration of the transducer is conducted in both air and water. The dynamic response of the transducer is tested against hot-wire anemometery in the wake flow of a circular cylinder over a wide range of Reynolds number.List of symbols C _D drag coefficient - D diameter of cylinder - d diameter of fiber - E modulus of elasticity of the fiber - e output voltage - F drag force per unit length of a cylinder - f frequency (Hz) - L length of the fiber cantilever - M magnification factor - m mass per unit length of the fiber - Re Reynolds number - q dynamic pressure (= 1/2 U ²) - U free stream velocity - density - v kinematic viscosity     

Measurement method of complex viscoelastic material properties

A measurement technique of viscoelastic properties of polymers is proposed to investigate complex Poisson’s ratio as a function of frequency. The forced vibration responses for the samples under normal and shear deformation are measured with varying load masses. To obtain modulus of elasticity and shear modulus, the present method requires only knowledge of the load mass, geometrical characteristics of a sample, as well as both the amplitude ratio and phase lag of the forcing and response oscillations. The measured data were used to obtain the viscoelastic properties of the material based on a 2D numerical deformation model of the sample. The 2D model enabled us to exclude data correction by the empirical form factor used in 1D model. Standard composition (90% PDMS polymer + 10% catalyst) of silicone RTV rubber (Silastic^® S2) were used for preparing three samples for axial stress deformation and three samples for shear deformation. Comprehensive measurements of modulus of elasticity, shear modulus, loss factor, and both real and imaginary parts of Poisson’s ratio were determined for frequencies from 50 to 320 Hz in the linear deformation regime (at relative deformations 10^?6 to 10^?4) at temperature 25 °C. In order to improve measurement accuracy, an extrapolation of the obtained results to zero load mass was suggested. For this purpose measurements with several masses need to be done. An empirical requirement for the sample height-to-radius ratio to be more than 4 was found for stress measurements. Different combinations of the samples with different sizes for the shear and stress measurements exhibited similar results. The proposed method allows one to measure imaginary part of the Poisson’s ratio, which appeared to be about 0.04–0.06 for the material of the present study.