Statistical evaluation of low cycle stress-strain curves parameters for alloyed structural steels weld metals at room and elevated temperature/Legiruotuju plienu suvirinimo siuliu medziagu deformavimo diagramu parametru statistinis ivertinimas kambario ir aukstesneje temperaturoje.
Daunys, M. ; Stulpinaite, A. ; Sniuolis, R. 等
1. Introduction
To improve the durability of the critical importance constructions,
it is necessary to analyze how working conditions and material
properties are influencing the strength and reliability of these
constructions. In the construction the fatigue damage under low cycling
loading is one of the most frequent failure mode. Because the
accumulation of plastic strain generally occurs in small volumes of the
material, particularly in parts and structures near the stress
concentrators, such as geometrical discontinuities, shoulders, keyways,
oil holes, welded joints, technological, welded and foundry defects,
termed notches [1-4]. Those components are surrounded by the elastically
deformed material. That is why the conditions of loading with limited
strain in these areas are similar to real constructions and machines.
For low cycle durability calculation the parameters [A.sub.1],
[alpha], [[bar.S].sub.T], characterizing the low cycle stress-strain
curves, are used. The cyclic proportionality limit stress
[[bar.S].sub.T] and parameters [alpha] and [A.sub.1] are determined at
symmetric tension-compression low cycle loading. This loading cycle is
the most universal and most correctly shows the characteristics of the
materials, because at elastic plastic cyclic deformation almost evenly
develop tension and compression deformations. Because it is complicated
and rather expensive to determine precise values of these
characteristics, particularly at elevated temperature, considering
temperature control and the cyclic stress-strain curves recording, the
dependences of the parameters [A.sub.1], [alpha], [[bar.S].sub.T] on
main mechanical characteristics of the material were investigated. In
our previous work detailed statistical analysis showed that those low
cycle stress-strain curves characteristics correlate the best with
modified plasticity criterion ([[sigma].sub.u]/[[sigma].sub.0.2])Z, i.e.
the parameter depending on ultimate tensile and yield strengths and
reduction of the area at fracture [5, 6].
In this work the materials are additionally specified in to three
groups by chemical composition of the main metal, such as 22K, 19MN5,
10ChN1M, 10GN2MFA, 12Ch2NMFA, 15Ch2MF, 15Ch2MFA, 15Ch2MFAA, 15Ch3NMFAA
and others. Analytical dependences of the cyclic stress-strain
parameters [alpha], [A.sub.1] and the proportionality limit stress
[[bar.S].sub.T] on the modified plasticity criterion
([[sigma].sub.u]/[[sigma].sub.0.2])Z for Cr-Ni, Cr-Ni-Mo, Cr-Ni-Mo-V
steels and theirs weld metals at room and elevated (250[degrees]C -
350[degrees]C) temperature are considered. The data of investigation of
75 welded joint materials at room temperature and 40 at elevated
temperature are obtained by the methods and testing equipment with
computer control and recording of stress-strain curves in Low Cycle
Fatigue laboratory of Kaunas University of Technology during 30 years.
2. Determination of low cycle stress-strain curves characteristics
At the best usable relationship it is assumed the relationship
between standard mechanical characteristic [[sigma].sub.pl],
[[sigma].sub.u], [[sigma].sub.0.2], [[sigma].sub.f], Z (where
[[sigma].sub.pl] is the proportional limit, [[sigma].sub.f] is the
fracture stress and Z is the reduction of area at fracture) and cyclic
characteristics. The low cycle loading experiments were carried out
under symmetrical tension compression and strain controlled conditions.
Low cycle loading with k semicycle curve is described by the equation
[5]
[[bar.[epsilon]].sub.k] = [[bar.S].sub.k] + [A.sub.1]
([[bar.e].sub.0] - [[bar.S].sub.T]/2)[k.sup.[alpha]] (1)
where [[bar.[epsilon]].sub.k] is cyclic strain; [[bar.S].sub.k] is
stress amplitude for k semi cycle; [[bar.e].sub.0] is initial stain;
[[bar.S].sub.T] is the cyclic proportionality limit stress and
[A.sub.1], [alpha] are constants.
Under strain limited conditions [[bar.[epsilon]].sub.k] =
2[[bar.e].sub.0] = const. To calculate the low cycle strain and stresses
the following units are used: [[bar.S].sub.k] =
[S.sub.k]/[[sigma].sub.pl]; [[bar.S].sub.1] =
[S.sub.1]/[[sigma].sub.pl]; [[bar.S].sub.T] =
[S.sub.T]/[[sigma].sub.pl]; [[bar.e].sub.0] = [e.sub.0]/[e.sub.pl];
[bar.[epsilon]] = [epsilon]/[e.sub.pl]; [bar.[delta]] =
[delta]/[e.sub.pl].
The parameters [A.sub.1] and [alpha] of low cycle loading under
limited strain are calculated by the k semicycle diagram using the
equations
[alpha] = 1/logk log [bar.[epsilon]] -
[[bar.S].sub.k]/[A.sub.1]([[bar.e].sub.0] - [[bar.S].sub.T]/2) (2)
or graphically from the equation
[alpha] = log[[bar.[delta]].sub.k] - log[[bar.[delta]].sub.1]/logk
(3)
and when k = 1
[A.sub.1] = [bar.[epsilon]] - [[bar.S].sub.1]/[[bar.e].sub.0] -
[[bar.S].sub.T]/2 (4)
or graphically from the equation
[A.sub.1] = [[bar.[delta]].sub.1]/[[bar.e].sub.0] -
[[bar.S].sub.T]/2 (5)
Due to unsettled of cyclic stress-strain curves for 1-9 semicycles
the values of [[bar.[delta]].sub.1] - [[bar.[delta]].sub.9] were
rejected. The parameters [A.sub.1] and [alpha] were determined from Eqs.
(1) - (5) and experimental results of all materials tested under low
cycle straining. When [alpha] < 0 - the materials cyclic harden, when
[alpha] > 0 - the materials cyclic soften and when [alpha] = 0 - the
materials are stable.
3. Statistic evaluation of cyclic stress-strain curves parameters
The machines used in nuclear power energetic, metallurgy and other
industries are operating at different temperatures. Therefore the
results of the materials investigated in work [5] were divided according
in-service temperature into groups: 1) at room temperature; 2) at
elevated temperature. In this work the materials are specified in to
three groups by chemical composition: 1) Cr-Ni; 2) Cr-Ni-Mo; 3)
Cr-Ni-Mo-V. The numbers of percentage of the investigated materials are
given in Table 1. As we can see 65.2% of the analyzed of alloyed steels
were tested at room temperature. In this work the largest part of the
investigated materials fall in to chemical classified Cr-Ni-Mo-V group
at room and elevated temperature.
Kurtosis coefficient and skewness coefficient are rating numerical
characteristics, which describe empirical distribution asymmetry and
flatness of the parameter's data comparing with normal
distribution. Histogram of stress-strain curves parameter [A.sub.1] for
Cr-Ni-Mo steels weld metals at room temperature is represented at Fig.
1. It has the left skewness compared with normal distribution ([A.sub.s]
= 0.25). The kurtosis coefficient ([E.sub.K] = -0.33) evidence that
analyzed parameters of [A.sub.1] for Cr-Ni-Mo of alloyed structural
steels weld metals at room temperature are spread in a wider interval
comparing with normal distribution.
To reduce the variance of the results, to define the more precise
mean value and obtain the better correlation relationship between the
analyzed data, clearly distinct results were eliminated, considering the
graphic view, using N. Smirnov criterion and the quartile width method.
Rectangular diagram in Fig. 2 shows that the scatter interval of the
results of the parameter [[bar.S].sub.T] for Cr-Ni-Mo-V steels weld
metals at elevated temperature is not wide (within limits [x.sub.min] -
[x.sub.max]) and there are no clearly distinct results. The represented
median value [x.sub.me] of investigated n number of materials divides
the scatter of the results into two equal parts. Defined area (within
quartiles limits [x.sub.0.25] - [x.sub.0.75]) describes the 50% scatter
of the middle values. Statistical characteristics of the proportionality
limit stress [[bar.S].sub.T] and the cyclic stress-strain curves
parameters [alpha] and [A.sub.1] at room (20[degrees]C) and elevated
(250[degrees]C-350[degrees]C) temperatures are given in Table 2. There
is no strongly outstanding result of the materials because mean values
of the parameters are similar to median values.
[FIGURE 1 OMITTED]
[FIGURE 2 OMITTED]
In previous works [5, 6] the accomplished statistical analysis
conformed that the results of the parameters [[bar.S].sub.T], [alpha]
and [A.sub.1] of cyclic stress-strain curves correlate the best with the
modified plasticity criterion ([[sigma].sub.u]/[[sigma].sub.0.2])Z at
room and elevated temperatures. Accordantly the dependences of the
cyclic stress-strain parameters a on the modified plasticity criterion
([[sigma].sub.u]/[[sigma].sub.0.2])Z for integrated group of Cr-Ni-Mo-V
steels weld metals and 95% confidence interval boundaries (dotted line)
for the theoretical regression line at elevated temperature are depicted
in Fig. 3. The same dependences of the proportionality limit stress
[[bar.S].sub.T] on the mechanical characteristic for integrated group of
Cr-Ni-Mo-V steels weld metals at elevated temperature and 95% confidence
interval ranges (dotted line) to theoretic line are represented in Fig.
4. As we can see from Figs. 3 and 4 the scatter of the results,
dependences for the characteristics [alpha] and [[bar.S].sub.T] of
cyclic stress-strain curves on the modified plasticity criterion, that
fall in to 95% confidence interval boundaries to theoretic line are
better for separated group than for integrated group at elevated
temperature. The number of percentage for the results inside 95%
confidence interval for cyclic stress-strain curves parameters a on the
modified plasticity criterion of Cr-Ni-Mo-V steels weld metals is 53.3%,
while for the integrated group it includes only 21.2%. Analogical
results are for the dependence of the characteristics [[bar.S].sub.T] on
the modified plasticity criterion ([[sigma].sub.u]/[[sigma].sub.0.2])Z.
It is 60% for Cr-Ni-Mo-V welded joints and 40% for integrated group.
Therefore we can say that with 95% guaranty the scatter of the results
of proportionality limit stress [[bar.S].sub.T] and low cycle
stress-strain curve's parameters [alpha], [A.sub.1] on the modified
plasticity criterion is better for separated group than the scatter of
results for integrated group.
Correlation analysis is statistical relation strength between
analyzed variables, which is expressed by correlation coefficient.
Pearson's correlation coefficient describes the strength of linear
dependence between the random and normally distributed results.
Correlation analysis is used to determine linear correlations [7]. If
linear model is not adequate, it is necessary to adapt nonlinear model.
It other works it was determined, that the best approximation choice for
the analyzed parameters is linear regression, accordingly the
accomplished statistical analysis conformed that the low cycle
stress-strain curve's parameters correlate the best with the
modified plasticity criterion ([[sigma].sub.u]/[[sigma].sub.0.2])Z at
room and elevated temperatures. The result given in Table 3 confirmed,
that the parameters [A.sub.1], [alpha] and [[bar.S].sub.T] of low cycle
stress-strain curves for the grouped weld metals are linearly dependent
on the modified plasticity criterion
([[sigma].sub.u]/[[sigma].sub.0.2])Z at room and elevated temperatures.
Pearson correlation coefficient vary from the minimum value [absolute
value of -0.24] for the parameter a of integrated group of weld metals
at room temperature to the maximum value [absolute value of 0.99] for
Cr-Ni steelsweld metals of alloyed structural steels at elevated
temperature. The correlation coefficient for the classified by the
chemical composition groups of the weld metals is better than for
integrated group at room and elevated temperature. Therefore we can say,
that the grouping of the random and normally distributed results by the
chemical composition influences the stronger linear relationship between
the parameters [A.sub.1], [alpha], [[bar.S].sub.T] of low cycle
stress-strain curves and the modified plasticity criterion
([[sigma].sub.u]/[[sigma].sub.0.2])Z for weld metals at room and
elevated temperatures.
The correlation analysis of weld metals, grouping them by chemical
composition, showed, that cyclic characteristics and the modified
plasticity criterion are significant by correlated by linear regression.
Accordingly analytical linear dependences of these parameters for weld
metals of Cr-Ni, Cr-Ni-Mo, Cr-Ni-Mo-V alloyed structural steels at room
and elevated temperature are given in Table 4.
[FIGURE 3 OMITTED]
[FIGURE 4 OMITTED]
The experimental and calculated results were distributed according
normal law in the intervals: [bar.x] [+ or -] 1.96 s with the
propobility P [approximately equal to] 0.95 (95% of the normal curve
area); [bar.x] [+ or -] s with the propobility P [approximately equal
to] 0.68; [bar.x] [+ or -] 1.96s with the propobility P [approximately
equal to] 0.95. Here [bar.x] is the mean value of the cyclic
stress-stain curves characteristics [[bar.S].sub.T], [alpha] and
[A.sub.1]; s is standard deviation [8]. The comparison of experimental
[[bar.S].sub.T.sup.exp] and calculated [[bar.S].sub.T.sup.cal]
parameters of integrated and grouped Cr-Ni-Mo-V steelsweld metals at
room temperature are shown in Fig. 5. The intervals are narrower for
Cr-Ni-Mo-V group than for the integrated group (because the correlation
cofficients are better) therefore the standart deviation is smaller.
However, in the intervals [bar.x] [+ or -] 0.675 x s, [bar.x] [+ or
-] s, [bar.x] [+ or -]1.96 x s the scatter of results is better for the
group of the materials, classified by chemical composition, than for the
integrated group at room temperature. The comparison of the experimental
[[bar.S].sup.exp.sub.T] and calculated [[bar.S].sup.cal.sub.T]
parameters at room temperature showed, that the scatter of the results,
within the interval [bar.x] [+ or -] 1.96 s is 4% greater for the
integrated group of weld metals of alloyed structural steels comparing
with the results of grouped Cr-Ni-Mo-V steels weld metals. The
analogical results are obtained for parameters [alpha] and [A.sub.1] of
weld metals at room temperature. The scatter of the results at elevated
temperature for all grouped weld metals in interval [bar.x] [+ or -]
1.96 s with the probilistic P [approximately equal to] 0.95, is 100% .
However, in the narrower intervals [bar.x] [+ or -] 0.675s and [bar.x]
[+ or -] s with the probabilities P [approximately equal to] 0.50 and P
[approximately equal to] 0.68 accordingly, the scatter of results is
better for all the groups of materials, classified by chemical
composition, than for the integrated group. Therefore, it is possible to
conclude that grouping of the random, normally distributed data by the
chemical composition gives more valuable information, because the
stronger linear relationship between the low cycle stress-strain
curves' parameters and the modified plasticity criterion is
obtained for different groups at room and elevated temperature.
[FIGURE 5 OMITTED]
4. Conclusions
1. Correlation coefficient of the parameters [alpha], [A.sub.1] and
[[bar.S].sub.T] for integrated group is smaller [absolute value of 0.02]
- [absolute value of 0.49] than for Cr-Ni, Cr-Ni-Mo, Cr-Ni-Mo- V steels
weld metals at room and elevated temperature.
2. With 95% guaranty the scatter of the results of proportionality
limit stress [[bar.S].sub.T] and cyclic stress-strain parameters
[alpha], [A.sub.1] on the modified plasticity criterion for separated
groups is smaller than the scatter of results for integrated group.
3. The scatter of the results, within the interval [bar.x] [+ or -]
1.96 s is 4% greater for integrated group of weld metals comparing with
the results of grouped Cr-Ni, Cr-Ni-Mo, Cr-Ni-Mo-V steels weld metals at
room temperature. The number of the results for all weld materials'
groups at elevated temperature between the experimental
[[alpha].sup.exp], [A.sup.exp.sub.1], [[bar.S].sup.exp.sub.T] and
calculated [[alpha].sup.cal], [A.sup.cal.sub.1], [[bar.S].sup.cal.sub.T]
parameters that fall in to the interval [bar.x] [+ or -]1.96 x s
boundaries is 100%. But in the narrower intervals [bar.x] [+ or -]
0.675- s and [bar.x] [+ or -] s the scatter of results is from 10% to
52% better for all the groups of materials, classified by chemical
composition, than for the integrated group.
4. Proposed analytical dependencies are recommended for preliminary
evaluation of the low cycle stress-strain parameters and durability.
Accomplished analysis showed that additional grouping of the materials
by the chemical composition has influence while defining the parameters
by the modified plasticity criterion at room and elevated temperatures.
Received July 01, 2010 Accepted October 22, 2010
References
[1.] Zeng, Z., Fatemi, A. Elasto-plastic stress and stain behavior
at notch roots under monotonic and cyclic loadings. -Journal of Strain
Analysis, 2001, vol.36, no.3, p.287-300.
[2.] Daunys, M., Bazaras, Z., Timofeev, B.T. Low cycle fatigue of
materials in nuclear industry. -Mechanika. -Kaunas: Technologija, 2008,
Nr.5(73), p.12-17.
[3.] Serdjuks, D., Rocens K., Pakrastins L. Hybrid composite cable
with an increased specific strength for tensioned structures. -The
Baltic Journal of Road and Bridge Engineering. -Vilnius: Technika, 2008,
vol.III, No3, p.129-136.
[4.] Kuranovas, A., Goode, D., Kvedaras, A. K., Zhong, S.
Load-bearing capacity of concrete-filled steel columns. -Journal of
Civil Engineering and Management International research and
Achievements. -Vilnius: Technika, 2009, vol.15, No.1, p.21-33.
[5.] Daunys, M., Sniuolis, R. Statistical evaluation of low cycle
loading curves parameters for structural materials by mechanical
characteristics. -Nuclear Engineering and Design, 2006, vol.236, no.13,
p.1352-136.
[6.] Daunys, M. Strength and fatigue life under low cycle
non-stationary loading. -Vilnius: Mokslas. 1989, p.46-64 (in Russian).
[7.] Ctepnov, M.N. Statistical Evaluation Methods of the Mechanical
Research Results. Catalog. -Moscow: Mashinostroenie, 1985.-232p. (in
Russian).
[8.] Zaks, L. Statistical Estimation. Statistics. -Moscow, 1976,
p.138-153 (in Russian).
M. Daunys *, A. Stulpinaite **, R. Sniuolis ***
* Kaunas University of Technology, Kestucio str. 27, 44312 Kaunas,
Lithuania, E-mail: Mykolas.Daunys@ktu.lt
** Siauliai University, Vilniaus str. 141, 76353 Siauliai,
Lithuania, E-mail: agette@gmail.com
*** Siauliai University, Vilniaus str. 141, 76353 Siauliai,
Lithuania, E-mail: rsrs@tf.su.lt
Table 1
Number and percentage of the material's groups
Group of
the material Room temperature Elevated temperature
Number % Number %
Cr-Ni 20 26.7 10 25.0
Cr-Ni-Mo 27 36.0 12 30.0
Cr-Ni-Mo-V 28 37.3 18 45.0
Total 75 100.0 40 100.0
Table 2
Statistical characteristics of stress-strain curves parameters
[A.sub.1], [alpha], [[bar.s].sub.T] of weld metals at room and
elevated temperatures
Parameters Cr-Ni
Room temperature
[[bar.s].sub.T] [alpha] [A.sub.1]
Number of materials 15 17 17
Mean value 1.66 0.013 1.67
Median value 1.62 0.010 1.71
Minimum value 1.55 -0.002 1.42
Maximum value 1.87 0.040 1.84
Kurtosis coefficient -0.14 0.73 0.50
Skewness coefficient 0.89 1.12 -1.05
Elevated temperature
Number of materials 5 3 6
Mean value 1.74 -0.003 1.57
Median value 1.72 -0.005 1.58
Minimum value 1.64 -0.009 1.43
Maximum value 1.85 0.004 1.69
Kurtosis coefficient 0.53 -- -1.42
Skewness coefficient 0.41 1.06 -0.25
Parameters Cr-Ni-Mo
Room temperature
[[bar.s].sub.T] [alpha] [A.sub.1]
Number of materials 22 22 21
Mean value 1.74 0.016 1.57
Median value 1.77 0.017 1.57
Minimum value 1.38 -0.002 1.39
Maximum value 2.11 0.036 1.78
Kurtosis coefficient -0.78 -1.32 -0.33
Skewness coefficient -0.41 -0.07 0.25
Elevated temperature
Number of materials 9 10 10
Mean value 1.71 -0.003 1.38
Median value 1.70 0.001 1.37
Minimum value 1.64 -0.034 1.10
Maximum value 1.78 0.018 1.71
Kurtosis coefficient -0.17 -0.25 -1.33
Skewness coefficient 0.16 -0.60 0.20
Parameters Cr-Ni-Mo-V
Room temperature
[[bar.s].sub.T] [alpha] [A.sub.1]
Number of materials 23 22 21
Mean value 1.84 0.022 1.63
Median value 1.83 0.024 1.63
Minimum value 1.6 -0.003 1.49
Maximum value 2.25 0.055 1.85
Kurtosis coefficient -0.32 0.10 -0.38
Skewness coefficient 0.56 0.14 0.45
Elevated temperature
Number of materials 15 15 15
Mean value 1.70 0.018 1.40
Median value 1.70 0.018 1.36
Minimum value 1.61 0.001 1.16
Maximum value 1.78 0.042 1.73
Kurtosis coefficient -0.63 1.12 -0.76
Skewness coefficient 0.025 0.52 0.52
Table 3
Correlation analysis of the parameters [[bar.s].sub.T], [alpha],
[A.sub.1] and the modified plasticity criterion ([[sigma].sub.u]/
[[sigma].sub.0.2]) Z at room and elevated temperature
Chemical
composition
(groups) Pearson correlation coefficient r
Room temperature
[[bar.s].sub.T] [alpha] [A.sub.1]
Cr-Ni 0.73 -0.63 0.63
Cr-Ni-Mo 0.43 -0.47 0.48
Cr-Ni-Mo-V 0.44 0.51 0.36
Integrated 0.37 -0.24 0.34
group
Chemical
composition
(groups) Pearson correlation coefficient r
Elevated temperature
[[bar.s].sub.T] [alpha] [A.sub.1]
Cr-Ni 0.93 0.99 -0.72
Cr-Ni-Mo 0.58 0.76 -0.59
Cr-Ni-Mo-V 0.65 0.58 -0.69
Integrated 0.53 0.50 -0.37
group
Table 4
Analytical dependences of the parameters [alpha], [A.sub.1],
[[bar.s].sub.T] on the modified plasticity criterion
([[sigma].sub.u]/[[sigma].sub.0.2])Z at room and elevated
temperature
Weld metals of Cr-Ni alloyed
structural steels
Room temperature
[alpha] = 0.053- 0.040 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[A.sub.1] = 1.26 + 0.426 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[[bar.s].sub.T] = 1.27 + 0.395 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
Elevated temperature
[alpha] = -0.073 + 0.053 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[A.sub.1] = 2.08 - 0.369 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[[bar.s].sub.T] = 1.20 + 0.381 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
Weld metals of Cr-Ni-Mo alloyed
structural steels
Room temperature
[alpha] = 0.051 - 0.036 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[A.sub.1] = 1.30 + 0.284 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[[bar.s].sub.T] = 1.18 + 0.575 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
Elevated temperature
[alpha] = -0.071 + 0.070 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[A.sub.1] = 2.53 - 1.156 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[[bar.s].sub.T] = 1.60 + 0. 121 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
Weld metals of Cr-Ni-Mo-V
alloyed structural steels
Room temperature
[alpha] = -0.023 + 0.044 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[A.sub.1] = 1.45 + 0.179 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[[bar.s].sub.T] = 1.43 + 0.410 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
Elevated temperature
[alpha] = -0.013 + 0.027 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[A.sub.1] = 2.04 - 0.588 ([[sigma].sub.u]/[[sigma].sub.0.2])Z
[[bar.s].sub.T] = 1.55 + 0.136 ([[sigma].sub.u]/[[sigma].sub.0.2])Z