Studies on crack growth rate under high temperature creep,fatigue and creep-fatigue interaction—II: On the role of fracture mechanics parameter aeffσg

For high temperature creep, fatigue and creep-fatigue interaction, several authors have recently attempted to express crack growth rate in terms of stress intensity factor $\text{K}{}_{\mathrm{I}}\text{= α}\text{aσ}{}_{\mathrm{g}}$, where a is the equivalent crack length as the sum of the initial notch length a₀ and the actual crack length $\text{a}{}^1$, that is, $\text{a = a}{}_0\text{+ a}{}^1$. On the other hand, it has been shown by Yokobori and Konosu that under the large scale yielding condition, the local stress distribution near the notch tip is given by the fracture mechanics parameter of $\text{aσ}{}_{\mathrm{g}}\text{?(σ}{}_{\mathrm{g}}$), where a is the cycloidal notch length, σ_g is the gross section stress and $\text{?(σ}{}_{\mathrm{g}}$) is a function of σ_g. Furthermore, when the crack growth from the initial notch is concerned, it is more reasonable to use the effective crack length a_eff taking into account of the effect of the initial notch instead of the equivalent crack length a. Thus we believe mathematical formula for the crack growth rate under high temperature creep, fatigue and creep-fatigue interaction conditions may be expressed at least in principle as function of $\text{a}{}_{\mathrm{eff}}\text{σ}{}_{\mathrm{g}}\text{, σ}{}_{\mathrm{g}}$ and temperature.In the present paper, the geometrical change of notch shape from the instant of load application was continuously observed during the tests without interruption under high temperature creep, fatigue and creep-fatigue interaction conditions. Also, the effective crack length a_eff was calculated by the finite element method for the accurate estimation of local stress distribution near the tip of the crack initiated from the initial notch root. Furthermore, experimental data on crack growth rates previously obtained are analysed in terms of the parameter of $\text{a}{}_{\mathrm{eff}}\text{σ}{}_{\mathrm{g}}$ with gross section stresses and temperatures as parameters, respectively.     

Crack tip closure and environmental crack propagation

Crack propagation rate, da/dN, and crack tip closure stress, σ_cc, in part-through crack fatigue specimens of aluminum alloys are drastically affected by gaseous environments. The present studies indicate that the crack closure reflects the influence of the environment on the plastic deformation at the crack tip, and, therefore, on the crack propagation rates. Postulating that da/dN is mainly determined by $\text{ΔK}{}_{\mathrm{eff}}\text{∝ (σ}{}_{\max }\text{?σ}{}_{\mathrm{cc}}\text{)}$ (instead of $\text{ΔK ∝ (σ}{}_{\max }\text{?σ}{}_{\min }\text{)}$, as is done traditionally) leads to the relationship $\text{da/dN = A(ΔK}{}_{\mathrm{eff}}\text{)}{}^{\mathrm{n}}$ in which $\text{A}$ and $\text{n}$ are virtually independent of the gaseous environment. The exponents are $\text{n ≈ 3.3}$ for Al 7075 T651 and $\text{n ≈ 3.1}$ for Al 2024 T351, respectively.     

Strain rate effects on viscoelastic crack bifurcation

The effects of applied strain rate on the viscoelastic crack bifurcation phenomenon in Polymethyl Methacrylate (PMMA) were investigated. It was still verified that the product $\text{σ}{}_{\mathrm{f}}\text{C}{}_{\mathrm{b}}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ was constant, as was already observed by Congleton and Petch, and Anthony, Chubb and Congleton, for brittle elastic materials, for any strain rate, where σ_f = the gross fracture stress and C_b= the main crack length until the bifurcation starts. However, it was found that the higher strain rate increases the main crack length C_b resulting in the decrease in the gross fracture stress σ_f and vice versa. This might be interpreted that the higher stress concentration at the initiation crack tip, which is realized by becoming more brittle due to the higher strain rate owing to the predominance of the elastic element in the viscoelastic material, decreases the gross fracture stress leading to the longer main crack length.     

A model for fatigue crack growth in a threshold region

The prediction of fatigue crack growth at very low ΔK values, and in particular for the threshold region, is important in design and in many engineering applications. A simple model for cyclic crack propagation in ductile materials is discussed and the expression

$\text{da}\text{dN}\text{=}\text{2}{}^{\mathrm{1+n}}\text{(1?2v)(ΔK}{}^2{}_{\mathrm{eff}}\text{?ΔK}{}^2{}_{\mathrm{c,eff}}\text{)}\text{4(1+n)π σ}{}^{\mathrm{1?n}}{}_{\mathrm{yc}}\text{E}{}^{\mathrm{1+n}}\text{?}{}^{\mathrm{1+n}}{}_{\mathrm{f}}$

developed. Here, n is the cyclic strain hardening exponent, σ_yc is cyclic yield, and ε_f is the true fracture strain. The model is successfully used in the analysis of fatigue data BS 4360-50D steel.     

Fatigue crack propagation work coefficient—a material constant giving degree of resistance to fatigue crack growth

Study on fatigue crack growth in steels was carried out from energetic point of view, i.e. taking account of plastic work around the fatigue crack. Based on the examination of the relation between fatigue crack growth rate (da/dN) and the plastic work around the fatigue crack tip (W_0.02 in SUS304, Fe-3Si and HT 60 steels, a material constant-fatigue crack propagation work coefficient-Q_0.02 is proposed. It is the ratio of W_0.02 to da/dN and means the degree of the resistance to fatigue crack growth. Numerical expression of Q_0.02 by mechanical properties was derived, which is given by

$\text{Q}{}_{0.02}\text{=}\text{9.3x10}{}^1\text{(σ}{}_{\mathrm{y}}\text{+σ}{}_{0.2}\text{)}\text{σ}{}_{\mathrm{y}}{}^{1.3}$

Comparison of Q_0.02 of various steels showed that Q_0.02 of high strength steels is very small compared with that of low strength steels. Graphical representation of the relation between Q_0.02 and da/dN at various values of ΔK/σ_y for steels revealed that da/dN at given value of ΔK/σ_y increase with decreasing Q_0.02. It is shown that fatigue crack growth behaviour of a steel (da/dN-ΔK relation) can be obtained from the Q_0.02-da/dN diagram by knowing the mechanical properties. Discussion on design stress level of the steels is also given.     

Fracture behaviour of dynamically loaded ship building plate material

The fracture toughness values of ship building mild steel measured over a temperature range ? 196°C to 28°C and crack tip strain rates ranging from 10^?5/sec to 10^?1/sec are examined in the light of the models recently proposed by Malkin and Tetelman. The effect of a change in inclusion morphology brought about by electroslag refining on the fracture toughness of the steel is also evaluated. It is found that the stress-induced fracture criterion ofthe model applies for the case where the ratio $\text{σ}{}^1{}_{\mathrm{f}}\text{σ}{}_{\mathrm{YS}}\text{? 3.94}$. This ratio is independent of the strain rate. In the strain induced fracture region of the model, the critical strain near the crack tip, ?_f(Rβ) is a function of the yield stress irrespective of temperature and strain rate. Electroslag refining reduces significantly the size and volume fraction of the inclusions and changes their shape from prolate ellipsoid to spherical. Apparently the electroslag refining does not improve fracture toughness significantly if the fracture toughness of the as received material measured with the major axis of the inclusions perpendicular to the crack front, is taken as a basis of comparison.     

Properties of bone at low temperatures and its potential as cryostat support material

The ultimate compressive strength, $\text{σ}{}_{\mathrm{<mtext>f</mtext>}}$, Young's modulus, E, and the integrated thermal conductivity k_T1–T2, of bone have been measured between 20 and 80 K. The two figures of merit indicating the quality for transmission of forces to low temperature apparatus are determined as $\text{σ}{}_{\mathrm{<mtext>f</mtext>}}\text{/}\text{k}{}_{\mathrm{4–77}}\text{= 42.5 ± 4 Ms m}{}^{\mathrm{?2}}$, and E/k_4–77 = 3.5 ± 0.2 Gs m^?2. According to these figures bone is comparable or superior to the best glass-fibre composites. Some observations on creep strength are added.     

Fracture toughness characterization of several aluminum alloys in semi-brittle fracture

A method has recently been developed for determining a nonlinear fracture toughness parameter defined by the relation $\text{G}\text{?}{}_{\mathrm{c}}\text{=}\text{C}\text{?}\text{G}{}_{\mathrm{c}}$ where G_c is the critical elastic strain energy rate as defined by Irwin. The $\text{C}\text{?}$ term is a function of the nonlinearity of the load-displacement test record and has been evaluated using the three parameter Ramberg-Osgood approach, although other curve fitting techniques could be applied as well. The method is quite straightforward and is applicable to plane stress, plane strain and mixed mode testing although only plane stress conditions are considered in this paper. For the case of a linear load-displacement record $\text{C}\text{?}\text{→ 1}$ and $\text{G}\text{?}{}_{\mathrm{c}}$ reduces to the linear elastic result.The toughness parameter $\text{G}\text{?}{}_{\mathrm{c}}$ has been evaluated for a number of high strength aluminum alloys and compared with published G_c values for these materials. The tests were conducted on center-cracked sheets of 2014-T6, 2024-T81, 7075-T6 and 7475-T61 aluminum alloys under conditions of varying specimen geometry and displacement gage length. It was found that the values of $\text{G}\text{?}{}_{\mathrm{c}}$ obtained from displacement readings with a gage length of 2 in. generally agreed with published values of $\text{G}{}_{\mathrm{c}}\text{=}\text{K}{}_{\mathrm{c}}{}^2\text{E}$. The $\text{G}\text{?}{}_{\mathrm{c}}$ values were found to vary inversely with gage length and a/w ratios. The variation in values for $\text{G}\text{?}{}_{\mathrm{c}}$ is of the same order of magnitude as the scatter in published values for G_c. However, $\text{G}\text{?}{}_{\mathrm{c}}$ appears to be less sensitive than G_c to changes in a/w.     

Effect of mechanical properties of material on rate of fatigue crack propagation

Many experimental and analytical equations on a rate of a fatigue crack propagation have been proposed. However, it seems that they can not fully express its complex behavior. There are still many problems remaining to be solved in order to clarify its mechanism. One of them is to clarify the relation between the rate of the crack propagation and the mechanical properties of material. In this paper, the rate of the crack propagation is analysed to clarify this problem. This analysis is based on the observation results of the fatigue crack propagation behavior previously by the authors. The analytical result is compared with the experimental one to make sure that they agree with each other. The conclusion obtained is; the rate of fatigue crack propagation is expressed by using the stress intensity factors as

$\text{dl}\text{dN}\text{= {}\text{c}\text{[Y}{}^2\text{F}{}^{\mathrm{a}}\text{E}{}^{\mathrm{a(1?n)}}\text{]}\text{} (K}{}_{\max }\text{)}{}^2\text{(K}{}_{\mathrm{a}}\text{)}{}^{\mathrm{a(2?n)}}$

. where C is a constant; E, Young's modulus; F, plastic coefficient; Y, yield stress; K_max and K_a, maximum and amplitude of the stress intensity factor, and α and n, exponents of the Manson-Coffin's law and work-hardening.     

Some problems in experimental fracture mechanics

A standard procedure for the determination of fracture toughness K_IC is discussed. The insufficiency of the existing K_ic determination confidence criteria is stressed and the following criteria are proposed instead: φ_max ? 1.5%; $\text{σ}{}_{\mathrm{fr}}{}^{\mathrm{net}}\text{σ}{}_{\mathrm{0.2\; ?\; 0.8}}$, in conjunction with the old criterion $\text{P}{}_{\max }\text{P}{}_{\mathrm{Q}}\text{? 1.1}$. Determination of K_IC from P_max should be used instead of from P_Q.A method for the determination of a point on the “force-displacement” diagram corresponding to crack growth initiation is set forth. The method is based on specimen compliance tests under repeated load-relief cycles. The crack growth initiation point is used to determine both the critical crack opening and plane strain fracture toughness. The indefinite effect of the growing crack (in the ease of crack opening or Cherepanov-Rice integral calculations) is thereby eliminated. Necessity is emphasized to determine the share of the J-integral which contributes to fracture process. A method for plotting the elastic displacement diagram is proposed which allows on the basis of preliminary estimates to determine fracture toughness of small-sized specimens without using special setups. The area ratio between the plastic and elastic strain diagrams is proposed to be adopted as fracture type criterion. Certain experiments to determine crack resistance of material specimens are described and discussed.     

A new parameter for prediction of creep crack growth rate at high temperature

By analysis based on a series of experimental data obtained by continuous observations using high temperature microscope during the creep test without interruption in vacuum of 10^?5 mm Hg for the purpose of the crack length measurements, a new mathematical equation for prediction of high temperature creep crack growth rate has been proposed in terms of disposable parameters, that is $\text{α}\text{a}{}_{\mathrm{<mtext>eff</mtext>}}\text{σ}{}_{\mathrm{g}}\text{,σ}{}_{\mathrm{g}}$ and temperature for 304 stainless steel within the range of $\text{α}{}_{\mathrm{g}}$ and temperature concerned. It can be seen that it is the best one to fit the experimental data among any other formula proposed hitherto.The new parameter proposed herein

$\text{8.48 × 10}{}^3\text{t}\text{log}{}_{10}\text{α}\text{α}{}_{\mathrm{<mtext>eff</mtext>}}\text{α}{}_{\mathrm{<mtext>g</mtext>}}\text{4.66 × 10}{}^2\text{+ 5.46}\text{log}{}_{10}\text{α}{}_{\mathrm{<mtext>g</mtext>}}$

where

$\text{α = 1.98 +0.36}\text{a}\text{w}\text{? 2.12}\text{a}\text{w}{}^2\text{+ 3.42}\text{a}\text{w}{}^3\text{, a≦0.7w}$

may be used for characterizing the creep crack growth rate just similar as Larson-Miller parameter for the creep life.     

A comparison of the geometry dependence of several nonlinear fracture toughness parameters

A number of fracture toughness tests on compact tension specimens have been performed for the purpose of comparing several nonlinear fracture toughness methods; including the nonlinear energy ($\text{G}\text{?}{}_{\mathrm{I}}$), J-integral ($\text{J}{}_{\mathrm{I}}$), COD ($\text{G}{}_{\mathrm{\delta }}$), and linear (–$\text{G}{}_{\mathrm{I}}$) approaches. The effect of variations in specimen thickness ($\text{B}$) and width ($\text{w}$) on the fracture toughness was examined for 7075-T651, 2124-T851, 2048-T35I, and 2048-T851 aluminum alloys, Ti-6Al-4V, and 4340 steel. Fracture toughness values were evaluated at both the initiation of stable crack growth and the onset of unstable fracture (peak load).It was found that the peak load toughness values are quite geometry sensitive at thicknesses below the requirement for plane strain fracture. At the initiation of stable crack growth, the toughness values are constant over a much larger range of specimen thickness. However, the nonlinearity of the load displacement curve is quite limited at this point and the associated fracture toughness is only 30–50% of the peak-load values.     

G̃c and R-curve fracture toughness values for aluminum alloys under plane stress conditions

The effect of specimen geometry and subcritical crack growth on the nonlinear energy fracture toughness, $\text{G}\text{?}{}_{\mathrm{c}}$, has been examined for thin, center-cracked sheets of 2024-T3 and 7075-T6 aluminum alloys. The procedure followed was to independently vary the specimen length, L, width, w, andd crack length-to-specimen width ratio and to determine the toughness both at the onset of subcritical crack growth and at the initiation of unstable fracture. Comparisons were also made with the R-curve toughness, G_R, evaluated at unstable fracture from which it was found that both $\text{G}\text{?}{}_{\mathrm{c}}$ and G_R displayed the same trend of change with geometrical variables, with $\text{G}\text{?}{}_{\mathrm{c}}$ consistently higher than G_R. When the nonlinear energy fracture toughness was evaluated at the onset of subcritical crack growth, it was found that the geometry dependence essentially disappeared.Scanning electron microscopic examination of some typical fracture surfaces showed that stable crack growth was accompanied by a gradual change of fracture mode from plane strain to plane stress. An analysis of possible errors in the experimental procedure showed that the scatter observed in $\text{G}\text{?}{}_{\mathrm{c}}$ values was not due to experimental errors, but apparently due to inhomogeneities in the materials. Several techniques were also introduced for the purpose of more directly incorporating crack growth into the $\text{G}\text{?}{}_{\mathrm{c}}$ determination, but it was found that they did not cause significant variation in the toughness values.     

Effect of loading history on stress corrosion cracking in high strength steel

The effect of preloading on crack nucleation time was examined with compact tension specimens having various notch radius in 0.1N-H²SO⁴ aqueous solution for 200°C tempered AISI 4340 steel. Crack nucleation time t_n increases by preloading for a given apparent stress intensity factor K_p2. The curve K_?2 vs. t_n deviates upward from the curve for the non preloading case. A linear relationship between the crack nucleation time and parameter is seen in semi-log diagram, where is taken as the value at t_n=α due to preloading. The apparent threshold stress intensity factor increases with K_?2 which is the apparent stress intensity factor of preloading. A detached crack is nucleated at some distance from the notch root and extends in a form of circle. This distance increases with increasing K_?2. The effect of load reduction during crack growth was examined. When the K-value was reduced from K₁ to K₂, an incubation time was observed before the crack started growing under the K₂-value. The incubation time t_m tends to increase with increasing ΔK = K₁-K₂. The threshold stress intensity factor was also found to increase for high load reduction.In order to explain these experimental results, a new dislocation model is proposed on the basis of stress induced diffusion of hydrogen in high stress region ahead of the notch root or a crack. This model suggests that the change in the crack nucleation time and the increase of the incubation time due to preloading or load reduction are caused by reducing the hydrostatic pressure and by spreading the hydrogen saturated region which requires more time for the hydrogen accumulation due to preloading or load reduction. The theory predicts the experimentally observed relations between and t_n and between log t_in and ΔK.     

Electronic and spin configurations of Co3+ and Ni3+ ions in oxides of K2NiF4 structure: A magnetic susceptibility study

In La₄LiCoO₈, Li⁺ and Co³⁺ ions are ordered in two dimensions and Co³⁺ ions undergo transitions from the low-spin to the intermediate as well as the high-spin states. Both Sr₄TaCoO₈ and Sr₄NbCoO₈ exhibit low to intermediate-spin state transitions of Co³⁺ ions. In the system LaSr_1?xBa_xNiO₄, the e_g electrons are essentially in extended states forming a band. With increase in x, the band width decreases accompanying an increase in unit cell volume; high-spin Ni³⁺ ions are formed to a small extent with increasing x, but there is no spin-state transition. In LaSrAl_1?xNi_xO₄, at small x, there is a small proportion of high-spin Ni³⁺; when x ≈ 0.6, there is an abrupt decrease in the $\text{c}$/$\text{a}$ ratio, signalling the formation of the band. In LnSrNiO₄, the $\text{c}$/$\text{a}$ ratio decreases sharply between Ln = La and Nd; this is likely to be accompanied by a broadening of the band.     

Fatigue crack propagation from a crack inclined to the cyclic tensile axis

Cyclic stresses with stress ratio R = 0.65 were applied to sheet specimens of aluminium which have an initial crack inclined to the tensile axis at angles of 30°, 45°, 72° or 90°. The threshold condition for the non-propagation of the initial crack was found to be given by a quadratic form of the ranges of the stress intensity factors of modes I and II. The direction of fatigue crack extension from the inclined crack was roughly perpendicular to the tensile axis at stress ranges just above the threshold value for non-propagation. On the other hand, at stress ranges 1.6 times higher than the threshold values the crack grew in the direction of the initial crack. The rate of crack growth in the initial crack direction was found to be expressed by the following function of stress intensity factor ranges of mode I, K₁, and mode II, K₂: $\text{dc}\text{dN}\text{= C(K}{}_{\mathrm{eff}}\text{)sum}$, where $\text{K}{}_{\mathrm{eff}}\text{= [K}{}_1{}^4\text{+ 8K}{}_2{}^4\text{]}{}^{\mathrm{<mtext>1</mtext><mtext>4</mtext>}}$. This law was derived on the basis of the fatigue crack propagation model proposed by Weertman.     

Nucleation mechanism of stress corrosion cracking from notches

Crack nucleation mechanism of hydrogen assisted cracking at notched cracks in aqueous solutions is investigated, using the compact type specimens with various notch radius in low-tempered 4340 steel. A detached crack initiates at some distance ahead of the notch root. The crack nucleation at the notched root is determined by the electrical potential method. When the crack initiates, the voltage difference starts to increase. The crack nucleation site is examined by SEM. The time for crack nucleation increases with the notch root radius, ρ, and decreases with the apparent stress intensity factor K_ρ. A linear relationship between the crack nucleation time, t_n, and the parameter ${}^2\text{K}{}_{\mathrm{\rho }}\text{/(πρ)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{-(2K}{}_{\mathrm{\rho }}\text{/(πρ)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}{}_{\mathrm{th}}\text{}}$ is seen in semi-log diagram, where $\text{(2K}{}_{\mathrm{\rho }}\text{/(πρ)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}{}_{\mathrm{th}}$ is almost equal to the yield shear strength.In order to explain these experimental results, a new model of micromechanics is proposed on the basis of stress induced diffusion of hydrogen in the high stress region ahead of the notch root. This model suggests that the detached crack initiates at the elasto-plastic boundary where the hydrogen concentration is from 2 to 5 times higher than that of the notch root surface. The theory agrees with experiments with respect to $\text{{2K}{}_{\mathrm{\rho }}\text{/(πρ)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{-(2K}{}_{\mathrm{\rho }}\text{/(πρ)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}{}_{\mathrm{th}}\text{}}$ vs t_n and t_n vs ρ.The empirical equation holds under constant t_n, K_ρ = K_o(ρ/ρ_eff)^m where K₀ is the stress intensity factor with ρ ≈ 0 under the present environment, ρeff is the effective notch radius and m is constant. The value of m is 0.25 for the crack nucleation time (t_n)_th corresponding to the threshold stress intensity factor (K_ρ)_th, 0.5 for t_n < (t_n)_th and 0 for ρ ≦ ρ_eff. The above equation agrees with the theoretical equation proposed by Tanaka and Mura for any t_n and ρ_eff.     

Cumulative damage in creep

The effect of two step stress variation and intermittent loading on the creep behaviour of commercially pure aluminium and stainless steel has been investigated. In the two step stress variation, the first stress σ₁ was applied for a given time t₁ and the stress level was switched over to σ₂. The resultant creep rate ?_s2⁺ and the failure time t₂ have been observed. Under the intermittent loading programme, the stress cycle was applied in the order σ₁-zero-σ₁ and the average creep rate on each reloading has been observed. The experimental data appear to give a cumulative damage rule in the form

$\text{t}{}_1\text{t}{}_{\mathrm{r1}}\text{+}\text{t}{}_2\text{t}{}_{\mathrm{r2}}\text{e}\text{?}{}_{\mathrm{s2}}{}^{\mathrm{+}}\text{e}\text{?}{}_{\mathrm{s2}}\text{=1}$

where t_r is the creep rupture time corresponding to a given stress and ?_s2 is the creep rate under the second stress in normal creep.     

Bounds on the effective thermal conductivity of a dispersion of fully penetrable spheres

We evaluate upper and lower bounds on the effective thermal conductivity K_e of a model of two-phase composite materials in which one of the phases consists of spherical inclusions (or voids) of conductivity K₂ and volume fraction φ₂, dispersed randomly throughout a matrix phase of conductivity K₁ and volume fraction φ₁. Our evaluations compare third-order bounds of Beran and of Brown, which utilize the three-point matrix probability function of the model, with bounds of De Vera and Strieder (which apply only to the aforementioned “fully penetrable sphere” model) and of Hashin and Shtrikman. The comparisons are made over extended ranges of values of both φ₂ and $\text{α =}\text{K}{}_2\text{K}{}_1$ and reveal that the best bounds among those we have tested (generally those of Beran) are sharp enough to give quantitatively useful estimates of K_e for $\text{0.1 ≤}\text{K}{}_2\text{K}{}_1\text{≤ 10}$ over a wide range of φ₂ values. They are sharp at high φ₂ values (i.e., φ₂ = 0.9) and very sharp at low φ₂ values (e.g. φ₂ = 0.1) where they remain useful for $\text{K}{}_2\text{K}{}_1\text{≈ 100}$. They are less sharp at intermediate values (e.g. φ₂ = 0.5). As is well known, such results immediately translate into equivalent results for the electrical conductivity, dielectric constant, or magnetic permeability of composites.