Microsoft Word - numero_50_art_30_2603 G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 354 Focused on the research activities of the Greek Society of Experimental Mechanics of Materials Probabilistic geotechnical engineering analysis based on first order reliability method Georgios Belokas University of West Attica, Greece gbelokas@uniwa.gr, https://orcid.org/0000-0003-2674-7877 ABSTRACT. Limit equilibrium engineering analysis needs the incorporation of probabilistic approaches for the determination of soil strength statistical measures for deterministic and probabilistic analyses. For the commonly used Mohr–Coulomb strength model, the uncertainty and the characteristic values of cohesion (c) and angle of shearing resistance (φ) determination is not straight- forward. This paper applies the first order reliability method (FORM) to estimate these values from the direct shear and the typical triaxial tests. The method is verified with test data. Furthermore, the FORM is applied to the planar failure limit equilibrium problem to determine the statistical measures of the safety margin (SM) and safety factor (SF). It is observed that the critical slip surface for the best estimate of the mean (SMm, SFm), for a 5% probability of exceedance (SMp=5%, SFp=5%) and for the characteristic value (SMk, SFk) do not coincide. It is interesting that the maximum probability of having a SM<0 or SF<1 does not correspond to the minimum best estimate of the SM or SF. FORM can be a very useful tool for complete probabilistic analyses. Further- more, probabilistic approaches applied to soil properties estimation can set a framework for the selection of their characteristic values for deterministic analyses. KEYWORDS. Soil strength; Geotechnical engineering; Probabilistic approaches; Reliability; Characteristic values. Citation: Belokas, G., Probabilistic geotech- nical engineering analysis based on first order reliability method, Frattura ed Integrità Strut- turale, 50 (2019) 354-369. Received: 16.02.2019 Accepted: 22.05.2019 Published: 01.10.2019 Copyright: © 2019 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION he limit equilibrium analysis for the design of geotechnical works requires the application of reliability theory principles for: a) the determination of soil properties in deterministic analyses and b) the uncertainty estimation of the safety margin (SM) or safety factor (SF) in probabilistic analyses. Modern codes of practice (e.g. Eurocode 7 – EC7) provide the alternative of both types of analysis, which may include one of the following methods: a) Monte Carlo, b) Point Estimate, and c) First Order Reliability. The application of First Order Reliability Method (FORM) for geo- technical engineering limit states analyses (e.g. [1,2]), a method for the estimation of error propagation, has been commonly T http://www.gruppofrattura.it/VA/50/2603.mp4 G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 355 applied to the uncertainty estimation of the experimental results of laboratory tests (e.g. GUM: 1995, see [3]) and can be applied directly wherever a closed form of analytical solution exists. Under the framework of EC7, the reliability analysis (uncertainty calculation) for a limit equilibrium problem can be per- formed with respect to the safety margin (SM), which is expected to have a value of SM≥0.0 for a certain level of con- fidence. That requires the knowledge of the uncertainty of the parameters that affect the value of SM. These parameters normally include the value and the uncertainty of the external and internal loads (permanent and mobile) and of the strength constants, as well as the spatial variability and uncertainty of the model (see Fig.1 based on Kulhawy [4]). One of the most critical components on the overall SM uncertainty is the strength uncertainty both due to spatial variability and model uncertainty. On the other hand, a deterministic analysis requires the best estimate of the loading conditions and the materials properties, which correspond to a specific level of confidence. Inherent soil  variability SOIL IN SITU or LAB  MEASUREMENT TRANSFORMATION  MODEL ESTIMATED SOIL  PROPERTY Data scatter Statistical  uncertainty Measurement  error Model  uncertainty Sampling  disturbance Figure 1: Factors affecting the property uncertainty (based on Kulhawy [4]). Therefore, for the Mohr – Coulomb strength parameters, cohesion (c) and angle of shearing resistance (φ), which are used in limit state analysis, an estimation of their mean value and their corresponding variation (or uncertainty) is required. These measures may be calculated either by direct application of statistical methods (e.g. for the direct shear test) or by an error propagation method (e.g. FORM for the typical triaxial test). The present work explores on the application of the FORM for the statistical evaluation of the strength parameters and for the slope stability analytical solution of a wedge failure mechanism. Issues with respect to the design and characteristic strength are also discussed, as well as the capability to apply the FORM into a general limit equilibrium slope stability problem. STATISTICAL MEASURES OF SOIL PROPERTIES deterministic analysis requires the knowledge of the best estimates of its individual components, i.e. loading conditions and material properties, for a specific confidence level. These correspond to the characteristic values of actions (Fk) and of soil parameters (Xk) defined in EC7 (ΕΝ-1997-1). Focusing on soil parameters, for any specific parameter, X, that affects the development of the limit state condition, its characteristic value Χk is defined as a cautious estimate of the mean value, (i.e. of the best estimate) of the mean, Χm (see ΕΝ-1997-1). The selection of this characteristic value has to be representative of the volume involved in the considered failure mechanism and it can depend on the type of the failure mechanism (e.g. local vs generalized failure). When the sample size, n, is large enough to apply statistical methods, the characteristic value corresponds to a worse value governing the occurrence of the soil parameter with a calculated probability not greater than 5% (it is 90% confidence interval, see ΕΝ-1997-1), which for a single variable model is given by Eq.(1).  , ,1 /       k m d X k d X mX X k S X k V V S X (1) where k is the confidence level coefficient for a given probability distribution, Sd,X is the sample standard deviation and V the variation coefficient. For a specific sample with unknown standard deviation, the Sd,X is the corrected standard deviation, which is related to the corrected – unbiased sample variance (s2) according to Eq.(2).    2 22 , 1 1 var( ) 1        n X d X i m i X s S X X n (2) A G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 356 The use of the corrected variance (sX2=var(X), Eq.(2)) instead the uncorrected sample variance (σX2), implies that there is little confidence that σX2 is a close estimate of the population variance, σ2. Had it been σX2σ2, the sample would closely follow a normal distribution, which is not the case in geotechnical engineering investigations as sample size is usually small, despite that population may follow a normal distribution. Therefore, the direct use of the statistical methods may not be applicable, because they may lead to non-representative values for the soil mass (see [2]). In order to account for this difficulty on the statistical error, direct values (e.g. [2,5]) or semi-empirical methods (e.g. [5, 6]) have been proposed for the estimation of the variation coefficient V. The V coefficient can also be used to include other types of errors, such as the error uncertainty due to ground spatial variability, the measurement error and the trans- formation uncertainty of the empirical equations application ([6]). However, a good knowledge of the statistical background can provide a better understanding on the selection of the cautious estimate (see [5,7]), or even to apply Bayesian statistical methods (e.g. [2,8]), not only in cases of a small sample but also of complex uncertainties. Concerning the characteristic values used in deterministic analysis, due to the different types of uncertainties involved (e.g. inherent soil variability, sampling disturbance, Fig.1), engineering judgement is also recommended for their selection (e.g. [9]). Such judgement should evaluate the relative importance of the following uncertainties (see also [6]):  Sampling quality (sampling type and sample condition).  The extent of the in-situ and laboratory investigation (the number and spatial distribution of samples and in-situ and laboratory tests).  The quality of the laboratory tests (accreditation and uncertainty of the laboratory measurements)  The spatial variability of the parameters and samples distribution with respect to the extent of the considered mechanism of geotechnical soil model. In addition, the following should also be taken into account: a) existing experience and data on similar soil units (including their uncertainty) and b) the failure mechanism with respect to the geotechnical profile (e.g. generalized vs local failure, short term vs long term conditions, small – large strains). For a single variable model (e.g. undrained strength, Su) and a small sample size (n) of laboratory data with unknown standard deviation of the population (the usual case for geotechnical engineering) and assuming a normally distributed population, the resulting estimated distribution follows the Student t-distribution. The estimated characteristic value (Xk), which corresponds to a probability P(Xk<μ)=1-p=1-α/2 (i.e. in a certain percentage of the cases the expected value of the true mean, μ, is greater than Xk), is then given by Eq.(3), in which SEX is the sample standard error given by Eq.(4). , 1 , 1 ,    k m p n X m p n d XX X t SE X t S n (3) ,X d XSE S n (4) where tp,n-1 is the n-1 degrees of freedom student distribution confidence parameter for the one-sided 1-p lower confidence limit of the true mean (μ). From Eqs.(1,3) we get the coefficient k=tp,n-1/n0.5. Schneider [10] proposed the approximate re- lationship Χk=Xm–0.5Sd,X, which corresponds to p=5% and n=14. As already mentioned, Eqs.(3,4) correspond to the single variable model (e.g. strength determined from the unconfined compression test, in which qu=2Su, where Su undrained shear strength). The Mohr – Coulomb failure criterion used in limit state analyses, is a two variables linear model that includes two constants (cohesion, c, and angle of shearing resistance, φ) and two variables. The for the cases of: a) the direct shear test is given by Eq.(5), in which σn is the imposed normal stress (the nonrandom or independent variable) and τ the measured – observed shear stress (the random or dependent variable) and b) the typical triaxial test is given by Eq.(6), in which σ3 is the imposed cell pressure – radial stress (the non- random or independent variable) and σ1 the measured – observed vertical stress (the random or dependent variable). τ=c+σntanφ= a+b·σn (5) σ1=2c·cosφ/(1–sinφ)+σ3(1+sinφ)/(1–sinφ)=a+b·σ3 (6) With regards to the experiments, if xi is the imposed value and yi is the corresponding observed value of a sample size n, the linear model is given by Eq.(7), where εi is the error and a and b the model constants. Eq.(7) gives a best estimate of the mean of y for a given x. The measurement error εi is the deviation of the yi from its best (deterministic) estimate y=a+xib (Fig.2). In simple linear regression, the best estimate of the constants may be derived from minimizing the error square sum. yi = a + xib + εi (7) G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 357 Figure 2: Graphical representation of the regression model. Assuming that: a) xi is an accurate observation, b) each xi is an independent observation, c) the error εi has a constant variation for each xi and d) the uncertainties of the yi observations are equivalent (otherwise weight coefficients are required), the best mean estimators for the linear regression coefficients are given by Eqs.(8,9). Observation xi can be a predetermined imposed loading (e.g. the σn in direct shear and the σ3 in the typical triaxial), while observation yi a measured reaction (e.g. the τ in direct shear and the σ1 in the typical triaxial).          ,1 1 1 1 22 ,2 1 1 1 1 , ,1 n n n n i i i i i i d yi i i i m xyn n n d x i i i i i i x x y y x y y x SCov x yn b r Var x y S x x x x n                               (8) m ma y b x    (9) In Eqs.(8,9) (xi, yi) are the data measurements of the two dimensional sample, ,y x are their mean values, Sd,x and Sd,y are the sample standard deviation of x and y measurements (Eqs.(10,11)) and rxy the Pearson sample correlation coefficient (given by Eq.(12)). The rxy is sensitive only to a linear relationship between two variables (|rxy|≤1, when rxy=1 the correlation is a perfect direct, i.e. increasing). Moreover, an unbiased estimate of the variance of y(x) with n-2 degrees of freedom is given by Eq.(13). The standard error estimators SEb and SEa of the b and a regression coefficients are given by Εqs.(14,15), re- spectively. Some applications in civil and geotechnical engineering of the two variables linear model have been presented by Baecher & Christian [1], Pohl [8] and Kottegoda & Rosso [11], as for instance the case of a variation with depth. A classic example is the increasing undrained shear strength with depth. The application of this model in the Mohr – Coulomb strength failure criterion has some individual characteristics that will be presented later.  2, 1 1 1      n d x i i S x x n             (10)  2, 1 1 1      n d y i i S y y n             (11)    2 22 2                     i i i i xy i i i i n x y x y r n x x n y y         (12)     2 1 , 1 | 2 2        n i i dvar y x S n S n n   , ˆˆi i iy a bx          (13) G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 358     2 1 2 2 1 1 var( | ) 1 2            n i i b n n i i i i x y SE n x x x x        (14)     2 2 2 1 2 1 1 1 1 1 2                  n in n i a b i i n i i i i x SE SE x n n n x x        (15) STATISTICAL MEASURES FOR THE DIRECT SHEAR TEST – TWO VARIABLES MODEL he direct shear test gives a direct determination of c and tanφ of the linear Mohr – Coulomb criterion (Eq.(5)), with σn=x and τ=y. For this occasion: a) the best estimates of the mean are cm=am and (tanφ)m=bm given by Eqs.(9,8), respectively (with xi=σni and yi=τi being the data measurements) and b) the standards error estimators SE(tanφ)=SEb and SEc=SEa are given by Εqs.(14,15). Concerning the confidence intervals of the linear regression coefficients estimators, the standard method relies on the normality assumption, which is justified if either: a) the errors in the regression are normally distributed (this leads to a t-statistic) or b) the number of observations n is sufficiently large (in this case the estimator is approximately normally distributed). Applying a statistical t-test, the linear regression random variables follow a student’s t-distribution with n-2 degrees of freedom (dof), i.e. Tc=(cm–c)/SEc ~ tn-2, T(tanφ)=[(tanφ)m–tanφ)/SE(tanφ) ~ tn-2, where c and tanφ represent the true mean values (or population mean values). The t-test includes the assumptions that the sample is representative of a specific soil unit, the observations are independent, while the variation of x=σn depends only on the uncertainty of the laboratory measurement. The σn variation has a negligible influence on the total uncertainty. Ignoring the influence of sampling disturbance and spatial variability, the characteristic values of cohesion, ck, and angle of shearing resistance, (tanφ)k are given by Eqs.(16,17) and their standard errors by Eqs.(14,15), respectively. The best estimates and their standard errors can be used for probabilistic analyses, by incorporating the standard errors as quantitative measures of the corresponding uncertainties (uc=SEtanφ=Sd,tanφ/n0.5 and utanφ=SEtanφ=Sd,tanφ/n0.5). , 2 k m p n cc c t SE (16)     , 2 (tan )tan tan   p nk m t SE   (17) Eqs.(16,17) may be used to give the characteristic Mohr – Coulomb failure envelope: τ=ck+σn(tanφ)k. Alternatively, an estimate of the characteristic failure envelope (see also [1,11]) can be obtained by incorporating the shear stress estimates for a specific probability (given by Eq.(18)). The resulting from Eq.(18) failure envelope is non-linear and needs to be ap- proximated by a linear regression to get the characteristic values of the Mohr – Coulomb failure criterion parameters.         2 2 2 2 1 1 1 1 tan 2                         n n m n n j nm j i i x τ c t n n x x     (18) Table 1 shows a fictional example for the application of the abovementioned relationships for a direct shear test. Typically, three specimens are obtained from each sample, which sample corresponds to a specific depth and location. Applying Eqs. (8,9) we obtain the best estimates of strength parameters and their corresponding uncertainties respectively: cm=23.95 kPa, (tanφ)m=0.37950 (i.e. φm=20.8ο) SEc=uc=5.3 kPa and SEtanφ=utanφ=0.01561. Applying Eqs.(16,17) for a probability p=5% we get the following characteristic values: ck1=14.5 kPa and (tanφ)k1=0.35186 (see Fig.3, characteristic 1). Applying a linear regression on the characteristic envelope derived from Eq.(18) we get the following characteristic values: ck2=16.50 kPa and (tanφ)k2= 0.38500 (see Fig.3, characteristic 2), which is less conservative than the characteristic 1. T G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 359 Sample 1 2 3 4 5 Specimen 1.1 1.2 1.3 2.1 2.2 2.3 3.1 3.2 3.3 4.1 4.2 4.3 5.1 5.2 5.1 σn (kPa) 100 300 500 100 300 500 100 300 500 100 300 500 100 300 500 τ (kPa) 56 140 218 56 130 208 69 151 222 48 143 241 76 134 196 Table 1: Example for direct shear test. 0 50 100 150 200 250 0 50 100 150 200 250 300 350 400 450 500 550 S h e a r  st re ss  τ  ( k P a ) Normal stress σn' (kPa) Exp Data Best estimate Characteristic 1 Characteristic 2 Characteristic 3 Figure 3: Failure envelopes, best estimate and characteristic envelopes. The typical test procedure of applying on the three different specimens of each soil sample the σn = 100, 300 and 500 kPa normal stresses respectively, may raise a question concerning the independence of the observations. An alternative is to consider a single variable model for τ on each one of the three σn values (in our case n=5 samples) and then apply a t-test on the observed τ values, which in our example leads to Table 2. Applying a linear regression on σn, τp=5% pairs we get the following characteristic values: ck3=14.02 kPa and (tanφ)k3= 0.38056 (characteristic 3). This approach, overcomes the issues concerning the independence of observation and gives similar results to characteristic 1 (Fig.3), however, it does not give a standard error of the parameters. For this specific case, the maximum standard error on τ lies close to the previously determined standard error of the cohesion. Specimen 1 2 3 σn (kPa) 100 300 500 τm (kPa) 61 139.6 212.8 SD (kPa) 11.27 8.14 10.83 n 5 5 5 SE (kPa) 5.0398413 3.641428 4.841487 a/2 0.05 0.05 0.05 tn-1 -2.131847 -2.13185 -2.13185 τp=5%,min 50.255831 131.837 202.4787 Table 2: Application of t – test on the values of the observed τ for each σn. Another approach often used in engineering practice is to apply a t-test on the determined c and tanφ pairs from each sample (again n=5 for the examined case). This approach gives cm=23.95 kPa, (tanφ)m=0.37950 (i.e. same as before), SEc=uc=6.6 kPa and SEtanφ=utanφ=0.02183 (i.e. greater than before). Applying Eqs.(15,16) for a probability p=5% we get the following cautious estimate of strength parameters: ck4=9.8 kPa and (tanφ)k4=0.33296 (characteristic 4), which are unrealistically conservative, as can been in Fig.4, due to the higher uncertainties. G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 360 0 50 100 150 200 250 0 50 100 150 200 250 300 350 400 450 500 550 S h e a r  st re ss τ  (k P a ) Normal stress σn' (kPa) Experimental Data Best estimate Characteristic 1 Best estimate (alternative) Characteristic 4 (alternative) Figure 4: Best estimates and corresponding characteristic envelopes comparison. STATISTICAL MEASURES FOR THE TYPICAL TRIAXIAL COMPRESSION TEST BY FORM n the typical triaxial compression test the Mohr – Coulomb failure criterion parameters are determined indirectly. For this test, a constant horizontal radial stress (the cell pressure) σr=σc=σ3 is applied to a cylindrical specimen, while the reaction of the axial stress is measured (Δσa, σ1=σ3+Δσa). We can consider that the cell pressure is an accurate ob- servation (i.e. the non – random). The Mohr – Coulomb failure criterion in terms of principal stresses is given by Eq.(6), in which the statistical measures of constants a and b can be determined directly from the two variables models, while the best estimates of the mean for c and tanφ constants are calculated indirectly by the transformation of Eqs.(19,20), respectively. sinφ=(b–1)/(1+b)  φ=asin[(b–1)/(1+b)] and tanφ=tan{asin[(b–1)/(1+b)}] (19) c=a(1–sinφ)/(2cosφ) (20) A way to calculate the uncertainties of c and tanφ is to apply the FORM, which method makes use of the second moment statistics (the mean and the standard deviation) of the random variables and assumes a linearized form of their perform- ance function (e.g. z=g(X1,…, Xn)) at the mean values of the random variables and independency between all variables. Truncating at the linear terms the Taylor expansion of the performance function about the mean, it is possible to obtain the first order approximation of the variance (σz2) of the true mean (μz) of z. Assuming uncorrelated non – random variables X1, …, Xn, the approximation of the variance is given by Eq.(21), an equation commonly used to estimate the uncertainties by error propagation for laboratory tests results.   2 2 2 2 , 1 1 var( )                     i n n z i d z X i ii i g g X SE SE X X (21) The sample variance var(Xi)=sXi2=(Sd,Xi)2 of a Xi variable (where Xi is tanφ or c) relates to the standard error SEXi by Eq.(4), while SEXi is a quantitative measure of the corresponding uncertainty uXi (i.e. uXi=SEXi). The variances stanφ2 and sc2 of the Mohr – Coulomb constants may be now calculated by applying Eq.(21) and considering Eqs.(19,20) as the perform- ance functions (i.e. c=g(a, φ) and tanφ=g(b) or φ=g(b)):   2 2      bSE SE b   ,    2 2 tan tan        bSE SE b   (22)     22 22              c a c c SE SE SE a       (23) I G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 361 where: 2 2 2 2 ( 1) 1 ( 1) / ( 1) )         b bb b ,    3/ 22 2 3 2 ( 1) 1 1 ( a 1 t n )            b b b b (24) 1 – /( ) ( )2     sin cos c a           (25) 21 )2(/( )    c a sin cos              (26) Next, the above equations are applied to laboratory results from the Herakleion marl (see Table 3, data from [12]), with 27 specimens and 9 samples (3 specimens per sample). Treating each specimen separately (i.e. sample size n=27), the best estimates of a and b constants of Eq.(6) are determined from the linear regression (Eqs.( 8,9)), which gives am=229.6 kPa and bm=3.02553. Applying these values into Eqs.(20,19) we get the best estimates cm=66.0 kPa and (tanφ)m=0.58225. The FORM (by application of Eqs.(23,22,15,14)) gives the corresponding uncertainties SEc=uc=15.5 kPa and SEtanφ=utanφ=0.0347. A t- student distribution for n-2 dof (Eqs.(16,17)) to account for the error propagation leads to the characteristic values ck1= 39.60 kPa and (tanφ)k1=0.52297. The best estimate of the mean failure envelope, as well as the estimated characteristic en- velopes are represented in Fig.5, in which figure the characteristic 2 and 3 envelopes are explained later in the same section. It is evident from Fig.5 that characteristic 1 can be indeed considered as a conservative estimate of the mean failure envelope. Sample ΓΣ2 (5m) Γ1 (2m) Γ1 (3.5m) ΓΣ13 (17m) Γ1 (18m) σ3 (kPa) 117.72 155.98 219.74 62 110 304 95 158 383 110.85 200.12 245.25 104 210 427 σ1 (kPa) 448.32 591.54 736.73 375 606 1488 475 506 1383 593.51 884.86 1044.77 514 785 1339 s (kPa) 283 374 478 219 358 896 285 332 883 352 542 645 309 498 883 t (kPa) 165 218 258 157 248 592 190 174 500 241 342 400 205 288 456 Table 3.1: Typical triaxial compression data from Herakleion marl ([12]). Sample Γ1 (6m) ΓΣ2 (9m) ΓΣ13 (24m) Γ3 (12.5m) σ3 (kPa) 79 146 265 119.68 258.98 346.29 73.57 141.26 162.85 52 139 222 σ1 (kPa) 365 618 997 640.59 1056.54 1339.07 429.68 587.62 750.47 672 819 1004 s (kPa) 222 382 631 380 658 843 252 364 457 362 479 613 t (kPa) 143 236 366 260 399 496 178 223 294 310 340 391 Table 3.2: Typical triaxial compression data from Herakleion marl ([12]). 0 200 400 600 0 200 400 600 800 1000 1200 1400 1600 σ 3 (k P a ) σ1 (kPa) Experimental results Best estimate Characteristic 1 (FORM) Characteristic 2 (envelope for p>5%) Characteristic 2 (linear regression) Figure 5: Statistical evaluation of the experimental results of Table 3. G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 362 Treating each sample separately we get the Mohr – Coulomb failure envelope constants of Table 4 (sample size n=9), which are represented in Fig.6 by the thin coloured lines, in the classic τ – σn Mohr – Coulomb diagram. Again the remaining char- acteristic envelopes of Fig.6 are explained next. Alternatively, for the characteristic envelope we can apply a t-student distribution into the y(x) estimate (i.e. the predicted σ1 for given σ3, given by Eq.(27)) for a p=5% probability (i.e. tp,n-2=1.70814) and n-2 dof. In Fig. 5 the characteristic 2 line shows the σ1, σ3 graph of the characteristic envelope from Eq.(27) for a p=5%. This characteristic envelope is non – linear and a linear regression leads to characteristic values ak=190.36 kPa and bk=2.92550 for Eq.(5) and then ck2=55.6 kPa and (tanφ)k2=0.56287 for the Mohr – Coulomb failure criterion (Figs.(5,6), characteristic 2). Taking into account the results for the direct shear test, it is obvious that Eq.(18) or Eq.(27) systematically gives a more optimistic characteristic envelope than the one computed based on mean values and uncertainty calculation of c and tanφ. However, this approach does not take into account the error propagation to calculate the uncertainties. Sample Mean ΓΣ2 (5m) Γ1 (2m) Γ1 (3.5m) ΓΣ13 (17m) Γ1 (18m) Γ1 (6m) ΓΣ2 (9m) ΓΣ13 (24m) Γ3 (12.5m) c (kPa) 64.34 40.99 22.34 21.24 60.41 77.75 29.69 76.83 48.20 201.58 Tanφ 0.59415 0.53191 0.83698 0.64201 0.64099 0.4862 0.64606 0.59134 0.63149 0.34038 φ (ο) 30.42 28.01 39.93 32.70 32.66 25.93 32.86 30.60 32.27 18.78 Table 4: Application of t – test on the values of the observed τ for each σn. Figure 6: Characteristic and best estimate Mohr – Coulomb failure envelopes superimposed on the Mohr – Coulomb failure envelopes of each sample.       2 3 32 1 3 2 2 1 3 3 1 1 1 2                         n m m n j n j i i a b t n n        (27) An alternative for the characteristic envelope is to compute first the characteristic values of a and b constants by applying a t-student distribution (i.e. ak=am-tp,n-2SEa and bk=bm-tp,n-2SEb) for p=5% and n-2 dof and then apply Eqs.(20,19) to compute the Mohr – Coulomb characteristic constants ck and (tanφ)k. This approach, which does not take into account the error propagation, leads to ak=139.86 kPa and bk=2.58981 constants of Eq.(6) (Fig.5, characteristic 3) and to ck3=43.4 kPa and (tanφ)k3=0.493949 constants of the Mohr – Coulomb failure criterion (Fig.6, characteristic 3). G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 363 Finally, as it was applied for the direct shear test, a common approach in engineering practice for the estimation of all stat- istical measures of c and tanφ is first to derive the Mohr – Coulomb constants (c, tanφ) for each sample separately (see Table 4) and then apply a statistical t-test on each constant independently for n/3-1 dof (where n the complete number of specimens and n/3 the number of specific locations of soil sampling). For the Herakleion marl case the n/3-1=8 dof lead to: a) mean values: cm=64.34kPa, (tanφ)m=0.59415, b) standard deviation: Sd,c=55.69kPa, Sd,tanφ=0.13577, c) standard error: Figure 7: Characteristic and best estimate Mohr – Coulomb failure envelopes superimposed on the Mohr – Coulomb failure envelopes of each sample.   Best estimate cm (kPa) (tanφ)m SEc (kPa) SEtan(φ) am, bm from linear regression of all σ1, σ3 data points cm, (tanφ)m from Eqs.(20,19) SEc and SEtanφ from FORM error propagation (this gives characteristic 1) 66.00 0.58225 15.50 0.0347 cm, (tanφ)m from mean values of Table 4 SEc and SEtanφ from single variable model (this gives characteristic 4) 64.34 0.59415 18.56 0.0423 Table 5: Mohr – Coulomb failure criterion constants (mean values and uncertainties).   Characteristic values ck (kPa) (tanφ)k Characteristic 1 (SE from FORM error propagation): ck=cm-tp,n-2SEc, tan(φ)k= tan(φ)m-tp,n-2SEtan(φ) 39.60 0.52297 Characteristic 2: ak, bk from linear regression on σ1pred, σ3 data points (σ1pred=am+σ3bm±tn-2SEσ1) ck, (tanφ)k from Eqs.(20,19) 55.65 0.56288 Characteristic 3: ak=am-tp,n-2SEa, bk=bm-tp,n-2SEb ck, (tanφ)k from Eqs.(20,19) 43.43 0.49395 Characteristic 4 (SE from one variable model, Eq.(4), for c and tanφ): cm and tanφ from mean values of Table 4 ck=cm-tp,n-2SEc, tan(φ)k= tan(φ)m-tp,n-2SEtan(φ) 29.82 0.50990 Table 6: Mohr – Coulomb failure criterion constants (characteristic values). G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 364 SEc=18.56kPa, SEtanφ=0.04526 and d) characteristic values: ck4=29.82kPa, (tanφ)k4=0.5099. The mean value from this alternate approach (characteristic 4, Fig.7) is almost the same with the previously calculated. Fig.7 compares the character istic value from this approach (characteristic 4) with the characteristic 1 from the FORM. The alternate approach is very conservative because it gives a higher SE than the FORM. All the above results are summarized into Tables 5 for the mean values and the corresponding uncertainties and Table 6 for the characteristic values.     APPLICATION OF THE FORM TO A SIMPLE PLANAR FAILURE PROBLEM urther extending the FORM application to engineering problems calculations, a limit state analysis may be performed in terms of either the safety margin SM=R-E (R is resistance and Ε action as defined in EC7) in terms of the safety factor SF=R/E. More specifically, in FORM applications the SM application is preferable compared to the FS because the actions Ε in the denominator of the SF enhances the non-linearity effects in the error propagation. For the safety margin the reliability index is β=SMm/Sd,SM, where SMm is the best estimate of the mean and Sd,SM the standard deviation, while by definition it is βSM=1/VSM, where VSM the variation coefficient. Since resistances and actions describe different types of random variables, they are expected to be uncorrelated, and their covariance, ρR,D, can be considered zero. Moreover, the variation of permanent loads is generally small compared and should not greatly affect the SM or SF value. The FORM is applied herein to a planar failure problem, which can be adapted to any type of failure surface (e.g. [13]) of limit equilibrium problems (e.g. method of slices). For the planar wedge failure problem considered herein Fig.8 shows the geometry. W is the weight of the wedge, H is the height of the slope, β is the angle of the slope to the horizontal, θ is the angle of the plane of failure with respect to horizontal and N and T are the normal and shear reaction forces on the plane of failure. The safety margin is then determined from the equilibrium equations that lead to Eq.(28). On the other hand the SF is given by Eq.(29).   212= c sin cos tan sin sin sinSM H γH β θ θ φ θ β θ      (28)     = 2c / γ sin cos tan sin sin sin sinSF H β θ θ φ β H β θ θ β          (29) H β θ W Wcosθ Wsinθ N T Α Β Γ ΑΒ=Η/sinθ Figure 8: Geometry of two dimensional planar failure. Applying the FORM on Eq.(28) the standard deviation of the SM is given by Eqs 30 to 33, in which uc, utanφ and uγ are the uncertainties of c, tan(φ) and γ respectively. Likewise, applying the FORM on Eq.(29) the standard deviation of the SF is given by Eqs.(34-37).   2 22 2 2 2 tan tan                      c SM SM SM u SM u u u c    (30) = sin   SM H c θ           (31) F G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 365  2 sin1 1= cos tan sin 2 sin         β θSM γH θ θ β           (32)    2 sin1 1 = cos tan sin sin 2 sin         β θSM H θ φ θ γ θ β (33)   2 22 2 2 2 tan tan c SF SF SF u SF u u u c                         (34)   2 1 = sin sin sin SF c H β θ θ β              (35) 1 = tan tan SF                      (36)  2 2 1 = - sin sin sin SF c γ γ H β θ θ β        = 2c / γ sin cos tan sin sin sin sinSF H β θ θ φ β H β θ θ β            (37) A slope geometry example of β=60ο and Η=25m and the material properties from Tables 5 and 6 are considered to exhibit the safety probability calculation. A deterministic calculation of the best estimate of the mean of SMm and SFm from Eqs.(28,29) respectively is computed first, by the best estimates of the mean of the soil properties (i.e. cm, tan(φ)m and γm). Then a normal distribution for SM and SF is applied, in order to estimate their value for a probability not greater than 5% (i.e., Eq.(38) with k=1.64485). The uncertainties u(SM) and u(SF) are calculated by Eqs.(30,34) applying the best estimates of the mean and the corresponding uncertainties of c, tan(φ) and γ (error propagation considers the standard error of the mean, i.e. u=SE). The values of SM and SF for probability p=5% and the corresponding probability p for SM<0 and SF<0, all calculated for the critical plane θcr that gives the minimum SMm, are presented in Tables 7 and 8. The two sets of statistical measures from Table 5 for soil strength properties have been used. ( )  mSM SM k u SM ,  ( )mSF SF k u SF   (38) cm (kPa) (tanφ)m uc (kPa) utanφ SMm= min(SMm) (kPa) θcr (o) uSM (kPa) SMp=5%= SMm - kuSM (kPa) VSM p(SM<0) (%) 66.00 0.58225 15.46 0.03470 1435.08 48 527.49 567.43 2.721 0.33 64.34 0.59415 18.56 0.04526 1396.92 48 631.85 357.62 2.211 1.35 Table 7: Calculations of SM for p=5% and p for SM<0 based on the minimum best estimate of the SMm and the corresponding θcr. cm (kPa) (tanφ)m uc (kPa) utanφ SFm= min(SFm) (kPa) θcr (o) uSF (kPa) SFp=5%= SFm - kuSF (kPa) VSM p(SF<1) (%) 66.00 0.58225 15.46 0.03470 1.639 40 0.241 1.243 6.800 0.40 64.34 0.59415 18.56 0.04526 1.6300 40 0.284 1.163 5.740 1.33 Table 8: Calculations of SF for p=5% and p for SF<1 based on the minimum best estimate of the SFm and the corresponding θcr. G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 366 It is obvious that the minimum values of the estimated means SMm=1435.08kPa and SFm=1.639 develop for different critical angle θcr of the plane of failure. However, the corresponding values of SMp=5% and SFp=5% (calculated in Tables 7 and 8 by applying Eq.(38)) are not the minimum ones, nor are pSM<0 and pSF<1 the maximum ones, as their minimum and maximum values develop for different critical slopes θcr (presented in Tables 9 and 10). The differences for the problem examined are of the order of 4% for the SM with respect to the min(SMp=5%) value and of the order of 1 to 2% for the SF with respect to the min(SFp=5%). Calculations based on min SM for p=5% cm (kPa) (tanφ)m uc (kPa) utanφ min(SMp=5%) (kPa) θcr (o) uSM (kPa) SMm (kPa) VSM pSM<0 66.00 0.58225 15.46 0.03470 554.18 46 546.51 1453.11 2.659 0.39 64.34 0.59415 18.56 0.04526 340.75 46 654.42 1417.49 2.165 1.54 Table 9.1: Calculated minimum SM for p=5% and the corresponding θcr, uSM and SMm values. Calculations based on max p for SM<0 cm (kPa) (tanφ)m uc (kPa) utanφ max(pSM<0) (%) θcr (o) uSM (kPa) SMm (kPa) VSM SMp=5% (kPa) 66.00 0.58225 15.46 0.03470 0.40 45 556.72 1475.12 2.650 559.40 64.34 0.59415 18.56 0.04526 1.54 45 667.07 1440.88 2.160 343.64 Table 9.2: Calculated maximum p for SM<0 and the corresponding θcr, uSM and SMm values. Calculations based on min SF for p=5% cm (kPa) (tanφ)m uc (kPa) utanφ min(SFp=5%) (kPa) θcr (o) uSF (kPa) SFm (kPa) VSF pSF<0 66.00 0.58225 15.46 0.03470 1.231 42.5 0.260 1.659 6.383 0.56 64.34 0.59415 18.56 0.04526 1.141 43.5 0.317 1.662 5.243 1.83 Table 10.1: Calculated minimum SF for p=5% and the corresponding θcr, uSF and SFm values. Calculations based on max p for SF<0 cm (kPa) (tanφ)m uc (kPa) utanφ max(pSF<0) (%) θcr (o) uSF (kPa) SFm (kPa) VSF SFp=5% (kPa) 66.00 0.58225 15.46 0.03470 0.62 45 0.287 1.718 5.976 1.245 64.34 0.59415 18.56 0.04526 1.90 45 0.338 1.701 5.035 1.145 Table 10.2: Calculated maximum p for SF<0 and the corresponding θcr, uSF and SFm values.. The second set of soil strength parameters lead to more conservative results. The comparison of the results of Tables 9.1 and 9.2 with Table 7 and of Tables 10.1 and 10.2 with Table 8 shows that for the probabilistic analyses it is preferable to determine the critical surface that corresponds to the minimum calculated SM or SF for a probability of exceedance 5% (or 90% confidence level, i.e. minSMp=5% or minSFp=5%), instead of calculating the critical surface of the minimum mean value SMm or SFm first and then the corresponding SMp=5% or SFp=5%. This happens because there is no linear relationship between SMm (or SFm) and uSM (or uSF) for a monotonically increasing or reducing failure plane angle. Concerning the influence of the individual uncertainties Fig.9 shows their influence on the minimum SM and Fig.10 on the probability of having an SM<0. Fig.11 shows their influence on the minimum SF and Fig.12 on the probability of having an SF<1. It is apparent that for this specific problem, the most influential factor on the probabilistic SM or SF is the uncertainty of cohesion. This is important because cohesion generally has a greater uncertainty from the angle of shearing resistance. The dependence SM and SF with the various uncertainties present the same trends. The above results will now be compared with the results from the deterministic analysis, for which the application of Eurocode 7 has been considered and more specifically the design analysis 3 (DA-3). DA-3 has become the national choice G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 367 of many countries for the case of overall stability of natural and cut slopes. According to this, for the drained conditions, the following components are used for deterministic analysis: a) the design strength cd=ck/(γΜγm), (tanφ)d=(tanφ)k/(γΜγm) with γm=1.0 (irrespectively of the type of analysis) and γΜ=1.25, b) the permanent actions Gd=γGGk, with γG=1.0 and c) the mobile design loads are Qd=γQQk, with γQ=0 or 1.3 for favourable or unfavourable loads respectively. Table 11 shows the results for the deterministic analyses, in which the four characteristic strength values of Table 6 have been used, factored by γΜ=1.25 for the design strength values. The minimum safety margin (SM) and safety factor (SF)   Figure 9: Influence of each uncertainty coefficient on the value of SM=SMm-kuSM.   Figure 10: Influence of each uncertainty coefficient on the probability of SM<0. Figure 11: Influence of each uncertainty coefficient on the value of SF=SFm-kuSF. G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 368 Figure 12: Influence of each uncertainty coefficient on the probability of SF<1. and their corresponding critical slip surface (θcr) are presented. The material unit weight has been taken γsoil,m=22kN/m3 with uncertainty uγ=2kN/m3. According to the calculations, only characteristic 2 gives acceptable safety margin (SM>0) and safety factor (SF>1.00) and. The critical failure plane determined by SM and SF do not coincide. Characteristic values 1 and 4 were derived from the statistical measures used in the above analyses. It is obvious that probabilistic and deterministic analyses can give the opposite result. In our case the deterministic analysis was more conservative, which has to do with the selection of the partial safety factors. i cki (kPa) tanφki cdi (kPa) tanφdi SMi θcr,SM,i (o) FSi θcr,FS,i (o) 1 39.60 0.52297 31.38 0.41837 -133.35 41 0.947 42 2 55.65 0.56288 44.52 0.45030 431.74 43 1.174 40.5 3 43.43 0.49395 34.74 0.39516 -86.60 40.5 0.967 41 4 29.82 0.50999 23.86 0.40799 -466.99 40 0.814 43.5 Table 11: Deterministic stability calculation with design parameters.. CONCLUSIONS Reliability methods can be incorporated for the determination of soil strength statistical measures and for probabilistic limit equilibrium analyses. Regarding soil strength, the determination of the uncertainties of the Mohr – Coulomb failure criterion constants c and tanφ from results of the typical triaxial test can be calculated through an error propagation technique, such as the FORM, applied to an appropriate performance function. This function relates the measured (σ1) and applied (σ3) quantities to the material constants (c, φ) of the M-C failure criterion by a nonlinear relation- ship. The FORM was applied to triaxial data and the results from different approaches for calculating the best estimate of the mean, the uncertainties and the characteristic values of the M-C constants were compared. For the approaches con- sidered, the direct application of the FORM gave a lower uncertainty, since the consideration of each single specimen as independent increased the sample size. The uncertainties together with the corresponding best estimates were also used to estimate the characteristic failure envelope. Regarding the characteristic failure, envelope four different approaches were presented with different results. The FORM gave the more reliable results. In any case, engineering judgement on the results is necessary. The results were applied to a simple planar failure problem, which allows for a direct comparison of the various statistical measures of the SM and SF, since a single failure surface is considered. It was observed that the maximum probability of having a SM<0 or SF<1 does not correspond to the minimum best estimate of the SM or SF. This is influenced by the non – linear relationship of SM and uSM (or SF and uSF) function and needs further exploring. The sensitivity analysis showed that the only influential uncertainty was that of cohesion, which generally happens to have the greater variability with regards to other soil constants. A sensitivity analysis is recommended in probabilistic analyses in order to determine the influence of each uncertainty separately. R G. Belokas, Frattura ed Integrità Strutturale, 50 (2019) 354-369; DOI: 10.3221/IGF-ESIS.50.30 369 The deterministic analyses gave in general more conservative results and only one gave acceptable SM and SF values. These analyses greatly depend on the selection of the partial safety factor selection that determines the characteristic value. The partial factors are indirectly related to an empirical probability of exceedance and the confidence of material properties and method of analysis. On the contrary, all probabilistic analyses gave acceptable SM and SF values, which means that probabilistic analyses could be used for an appropriate selection of partial safety factors. Therefore, the application of statistical methods can also set a framework for the selection of the characteristic mechanical properties, for deterministic analyses. REFERENCES [1] Baecher G. and Christian J. (2003). Reliability and Statistics in Geotechnical Engineering, Wiley, p. 618. [2] Orr T. and Breysse D. (2008). Eurocode 7 and reliability-based design. Reliability-Based Design in Geotechnical Engineering. Computations and Applications, Kok-Kwang Phoon(edr), Taylor & Francis, pp. 298–343. [3] ISO/IEC Guide 98-3:2008. Uncertainty of measurement – Part 3: Guide to the expression of uncertainty in measure- ment (GUM: 1995). [4] Kulhawy H. (1992). On the evaluation of static soil properties, Stability and Performance of Slopes and Embankments – II, Proceedings of a Specialty Conference, Seed R.B. and Boulanger R.W. (eds), ASCE, pp. 95–115. [5] Duncan J. M. (2000). Factors of safety and reliability in geotechnical engineering, Journal of Geotechnical and Geo- environmental Engineering, 126(4), pp. 307–316. [6] Schneider H.R. and Fitze P. (2013). Characteristic shear strength values for EC7: Guidelines based on a statistical framework. Proc. 15th European Conference on Soil Mechanics and Geotechnical Engineering, 4, pp. 318–324. [7] Fellin W. (2005). Assessment of characteristic shear strength parameters of soil and its implication in geotechnical design, Analyzing Uncertainty in Civil Engineering, Fellin, W., Lessmann, H., Oberguggenberger, M., Vieider, R. (eds.), Springer, pp. 1–15. [8] Pohl C. (2011). Determination of characteristic soil values by statistical methods. ISGSR 2011 - Vogt, Schuppener, Straub & Bräu (eds), pp. 427–434. [9] Frank R., Bauduin C., Driscoll R., Kavvadas M., Krebs Ovensen N., Orr T. and Schuppener B. (2004). Designers' Guide to EN 1997-1 Eurocode 7: Geotechnical Design - General Rules, Thomas Telford, p. 216. [10] Kottegoda N. T. and Rosso R. (2008). Applied Statistics for Civil and Environmental Engineers, 2nd Edition, Wiley- Blackwell, p. 736. [11] Tsiambaos G. (1988). Technicogeological features of the Herakleion Crete Marls. PhD, University of Patras, p. 358 (in Greek). [12] Schneider H.R. (1997). Definition and determination of characteristic soil properties. Proc. 14th International Conf- erence on Soil Mechanics and Foundation Engineering, 4, pp. 2271–2274. [13] Wu T. (2008). Reliability analysis of slopes, Reliability-Based Design in Geotechnical Engineering, Computations and Applications, Kok-Kwang Phoon (ed.), Taylor & Francis, pp. 413–447. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice