Constitutive Relations of Stress and Strain in Stochastic Finite Element Method
|
|
- Erik Barber
- 6 years ago
- Views:
Transcription
1 American Journal of Computational and Applied Mathematics 15, 5(6): DOI: 1.593/j.ajcam Constitutive Relations of Stress and Strain in Stochastic Finite Element Method Drakos Stefanos International Centre for Computational Engineering, Rhodes, Greece Abstract The analysis and design in structural and geotechnical engineering problems requires the calculation of stress and strain which is generally a difficult task because of the uncertainty and spatial variability of the properties of soil materials. This paper presents a procedure of conducting Stochastic Finite Element Analysis using Polynomial Chaos in order to propagate the uncertainties of input to constitutive relation of stress and strain. The problem is dominated by highly non linearity. Among other methods the procedure leads to an efficient computational cost for real practical problems. This is achieved by polynomial chaos expansion displacement field, stress and strain also. An example of a plane-strain strip load on a semi-infinite elastic foundation is presented and the results of settlement are compared to those obtained from the closed form solution method. A close matching of the two is observed. The constitutive relation of stress and strain is presented as result of the Polynomial Chaos expansion and Monte Carlo method. A close matching of the two method is observed also. Keywords Stochastic Finite Element Method, Constitutive Relations, Polynomial Chaos, Quantification of Uncertainty 1. Introduction The analysis and design in structural and geotechnical engineering problems requires the calculation of stress and strain which is generally a difficult task because of the uncertainty and spatial variability of the materials s properties. Various forms of uncertainties arise which depend on the nature of geological formation or construction method, the site investigation, the type and the accuracy of design calculations etc. In recent years there has been considerable interest amongst engineers and researchers in the issues related to quantification of uncertainty as it affects safety, design as well as the cost of projects. A number of approaches using statistical concepts have been proposed in engineering in the past 5 years or so. These include the Stochastic Finite Element Method (SFEM) [1-3], and the Random Finite Element Method (RFEM) [4-8]. The RFEM involves generating a random field of soil or structure properties with controlled mean, standard deviation and spatial correlation length, which is then mapped onto a finite element mesh. However the number of works on the stochastic stress and strain calculation and their statistical moments are limited. An essential paper on the field is presented by Ghosh & Farhat [9] where the constitutive relation of stress and strain calculated by different approaches. * Corresponding author: stefanos.drakos@gmail.com (Drakos Stefanos) Published online at Copyright 15 Scientific & Academic Publishing. All Rights Reserved In the past SFEM has been developed using different expansions of stochastic variables. In this paper we present SFEM [11-13] using the method of Generalized Polynomial Chaos (GPC) [14]. To descretise the stochastic process of material the Karhunen-Loeve Expansion was used and it is presented. The constitutive relation of stress and strain calculated using the Generalized Polynomial Chaos and verified against Monte Carlo simulation which is treated as the exact solution based on a series of computational experiment. A numerical example of foundation settlement given in the last part of the paper and the results of settlement compared with those arises from closed form solution. The two methods of stress and strains constitutive relation compared also and the results are presented.. Problem Description and Model Formulation Let us consider a general boundary value problem of computation of probable deformation of a body of arbitrary shape having randomly varying material properties caused by the application of a randomly varying load as shown in Fig. 1. According to the elasticity theory a boundary value problem can be described as follow: σσ iiii,jj (xx, ωω) = ff(xx, ωω) iiii DD ΩΩ σσ iiii (xx, ωω) = CC iiiiiiii (xx, ωω)εε kkkk (xx, ωω) iiii DD ΩΩ (1) uu(xx, ωω) = gg DD iiii BB DD σσ iiii (xx, ωω)nn jj = gg NN iiii BB DD
2 American Journal of Computational and Applied Mathematics 15, 5(6): And in the weak form as: aa(uu, vv) = ll(vv) () aa(uu, vv) = εε ΤΤ (vv)cc(xx, ωω)εε(uu) dddd (3) DD ll(vv) = ff(xx, ωω) vvvvvv + gg NN vvvvvv DD BB NN εε ΤΤ (vv)cc(xx, ωω)εε(uu) dddd gg BB DD DD (4) Figure 1. Body of arbitrary shape In order to model the problem assuming the sample space (Ω, F, P) where F is the σ-algebra and is considered to contain all the information that is available, P is the probability measure and the spatial domain of the soil or the structure is DD R. The Elasticity modulus {EE( xx, ωω) : DD ΩΩ} and the external load {ff( xx, ωω) : DD ΩΩ} considered as second order random fields and their functions are determined EE, ff: DD ΩΩ R VV = LL ΩΩ, LL (DD) and characterized by specific distribution where in our case as lognormal. Considering as μμ kk, σσ kk and vv kk = σσ kk the μμ kk mean value the standard deviation and the coefficient of variation of Elasticity modulus, the lognormal distribution is given [8]: EE = exp(μμ llllll + σσ llllll ZZ(ωω)) (5) Where the mean values and the variance of the distribution are equal to: σσ llllll = ln(1 + vv kk ) μμ llllll = ln(μμ kk ) 1 σσ (6) llllll And ωω ΩΩ, ΖΖ ΝΝ(,1) To separate the deterministic part from the stochastic part of the formulation the Karhunen-Loeve expansion has been used. It is considered as the most efficient method for the discretization of a random field, requiring the smallest number of random variables to represent the field within a given level of accuracy. Based on that the stochastic process of Young modulus over the spatial domain with a known mean value E (x) and covariance matrix CCCCCC(x 1, x ) assuming lognormal distribution is given by: EE xx, ξξ(ωω) = exp(εε (xx) + κκ λλ κκ ww κκ (xx)ξξ κκ (ωω)) (7) In practice, calculations were carried out over a finite number of summations (for example 1-5) so the approximate stochastic representation is given by the trancuated part of expansion: EE xx, ξξ(ωω) = exp(εε (xx) + ΚΚ κκ=1 λλ κκ ww κκ ξξ κκ (ωω)) (8) λλ κκ : are the eingenvalues of the covariance function ww κκ (xx): are the eingenfunctions of the covariance function CCCCCC(x 1, x ) x D and ω Ω and ξξ = [ξξ 1, ξξ 1,, ξξ ΜΜ ]: ΩΩ R ΜΜ = 1 1 ΜΜ The pairs of eingenvalues and eingenfunctions arised by the equation: CC(xx DD 1, xx )φφ κκ (xx ) = λλ κκ ww κκ (xx 1 ) (9) Using the Karhunen-Loeve expansion the stochastic elasticity tensor is given by: CC iiiiiiii (xx, yy) = EE(xx)CC iiiiiiii (xx), ii, jj, kk, ll = 1,,3 (1) CC iiiiiiii (xx): is expressed in terms of (deterministic) Poisson's ratio as CC iiiiiiii (xx) = vv δδ (1+vv) iiii δδ kkkk + 1 (δδ (1+vv) iiiiδδ jjjj + δδ llll δδ jjjj ) (11) To compute the statistical moments of the calculations output we perform a change of variable yy kk : = ξξ kk (ωω) and yy = [yy 1, yy, yy MM ] [1]. If the random variables are independent and ρρ ii denote the density of ξξ ii then the joint density is given by: ρρ(yy) = ρρ 1 (yy 1 )ρρ (yy ) ρρ MM (yy MM ) (1) Drakos & Pande [1, 13] developed a new algorithm of Stochastic Finite method using the Generalized Polynomial Chaos (appendix A) and proved that the problem formulation has the final form: mm KK mm = qq ff + tt gggg mm KK BBBB gg dd (13) mm = ρρ(yy)ψψ κκ (yy)ψψ pp (yy)exp( ΚΚ λλ κκ φφ κκ yy κκ KK mm = BB TT exp(εε (xx))cc (xx)bbbbbbbbbb DD tt gggg = κκ=1 ) qq = ρρ(yy)ψψ pp (yy)ψψ 1 (yy) dddd KK BBBB = BBBB ff = φφ ΤΤ ff(xx)dddd DD φφ ΤΤ gg NN dddd BBBB B is strain displacement matrix. φφ is the hat function. ψψ is the Polynomial Chaos dddd BB TT exp(εε (xx))cc (xx)bbbbbbbbbb (14) (15)
3 166 Drakos Stefanos: Constitutive Relations of Stress and Strain in Stochastic Finite Element Method 3. Constitutive Relations of Stress and Strain The calculation of the constitutive relation of stress and strain in the case of stochastic problems is quite complicated especially when the invariant of them are needed where the equations become highly nonlinear. In [9] several numerical integration schemes to evaluate the statistical moments of strains and stresses in a random system is presented. In the current work the propagation of the input uncertainty to the stress and strain relation is modelled by the polynomial Chaos expansion and verified against Monte Carlo simulation which is treated as the exact solution of the problems. The computational implementation of the Monte Carlo Method leads to the random field generation and the requested function u k (x) gets a new value for each realization. At the end of all simulations the statistical εε iiii (xx, yy) = 1 uu (kk) ii (xx)ψψ κκ (yy) moment are calculated. The expected value and the variance are given by: EE uu(xx) = 1 KK uu KK kk(xx) VVVVVV(uu(xx)) = 1 (uu (16) KK KK 1 kk(xx) EE(uu(xx)) In an elastostatic problem of homogeneous isotropic body one of the field equations that must be satisfied at all interior points of the body is the Strain-Displacement relations:,jj εε iiii (xx, yy) = 1 uu(xx, yy) ii,jj + uu(xx, yy) jj,ii ii, jj = 1,,3 (17) Using the displacement polynomial chaos expansion uu(xx, yy) = uu kk (xx)ψψ κκ (yy) (18) Q and ψ are given in appendix A The equation 17 leads to: + uu (kk) jj (xx)ψψ κκ (yy) εε iiii (xx, yy) = 1 uu (kk) ii,jj (xx)ψψ κκ (yy) + uu (kk) jj,ii (xx)ψψ κκ (yy) = εε (kk) iiii ψψ κκ (yy) 3.1. Expected Value of Strains According to the polynomial chaos expansion of strains the expected value can be evaluated by the following: EE[εε iiii (xx, yy)] = EE 1 uu (kk) ii,jj (xx)ψψ κκ (yy) + uu (kk) jj,ii (xx)ψψ κκ (yy) = 1 () uuii,jj (xx) EE[ψψ (yy)] + uu (kk) ii,jj (xx)ee[ψψ κκ (yy)] + uu () jj,ii (xx) EE[ψψ (yy)] + uu (kk) jj,ii (xx)ee[ψψ κκ (yy)] 1 1 EE εε iiii (xx, yy) = 1 uu () (xx) ii,jj + uu () () jj,ii (xx) = εε iiii () 3.. Variance of Strains Respected to the expected value evaluation and the Chaos Polynomial characteristics the variance of the strain tensor can be calculated as: σσ = EE εε iiii (xx, yy) (EE εε iiii (xx, yy ) PP = [εε (κκ) iiii (xx)] ρρ(yy)ψψ kk (yy)ddyy [εε () iiii (xx)] = [εε () iiii (xx)] ρρ(yy)ψψ (yy)ddyy + [εε (κκ) iiii (xx)] ρρ(yy)ψψ kk (yy)ddyy 1,ii [εε () iiii (xx)] σσ = (κκ) (xx)] < ψψ kk (yy) > (1) [εε iiii (19)
4 American Journal of Computational and Applied Mathematics 15, 5(6): Expected Value of Stress Tensor The constitutive relation of stress and strain given by the well known equation of Hooke s law equation: Using the polynomial chaos expansion of strains we get: σσ iiii (xx, yy) = EE(xx, yy)cc iiiiiiii εε mmmm () σσ iiii (xx, yy) = EE(xx, yy)cc iiiiiiii εε (kk) mmmm ψψ κκ (yy) (3) According to the elasticity modulus distribution and the strain tensor Chaos Polynomial expansion the expected value of stress tensor takes the following form: < σσ iiii (xx, yy) >=< EE(xx, yy)cc iiiiiiii εε mmmm > =< eeee pp(mmmm llllll + σσ llllll yy) CC iiiiiiii = ee mmmm llllll CC iiiiiiii εε (kk) mmmm ψψ kk (yy) > εε (kk) mmmm ρρ(yy)ee σσ llllll yy ψψ kk (yy)ddyy ) < σσ iiii (xx, yy) >= ee mmmm llllll CC iiiiiiii εε (kk) mmmm < ee σσ llllll yy ψψ kk (yy) > (4) 3.4. Variance of Stress Tensor Having calculated the expected value of stress tensor and knowing its stochastic equation of by the Hook s law constitutive relation, the variance of stress tensor can be calculated as: σσ = EE σσ iiii (xx, yy) (EE σσ iiii (xx, yy ) =< EE(xx, yy)cc iiiiiiii PP εε (kk) mmmm ψψ κκ (yy) > [ee mmmm llllll CC iiiiiiii Given the elasticity modulus distribution the variance of stress tensor is equal to: PP εε (kk) mmmm < ee σσ llllll yy ψψ kk (yy) > σσ =< exp(mmmm llllll + σσ llllll yy) CC iiiiiiii εε (kk) mmmm ψψ κκ (yy) > ee mmmm llllll (kk) CC iiiiiiii εε mmmm ρρ(yy)ee σσ llllll yy ψψ kk (yy)ddyy (5) ] ) 3.5. Pore Pressure Calculation A major issue in a wide range of geotechnical and geomechanics problem is the estimation of the excess pore pressure in the ground. Using the Chaos Polynomial expansion the statistical moments of pore pressure can be evaluated as following. The pore pressure is given by the following equation: pp = KK aa εε vv = KK aa (εε 11 + εε + εε 33 ) (5) εε vv : is the volumetric strain KK aa : is shown [15] as Assuming a minimum value of KK aa = EE (1 vv) KK aa EE (1 vv) Expected Value of Pore Pressure The expected value of pore pressure is the result of the summation of the expected value of the strain s Chaos expansion multiplied by the stochastic fluid modulus: EE[pp] = EE[KK aa (εε 11 + εε + εε 33 )] (7) Replacing the longnormal distribution o elasticity modulus we get: EE[pp] = EE eeeeee(mmmm llllll + σσ llllll yy) (1 vv) εε (kk) 11 ψψ κκ (yy) + εε (kk) ψψ κκ (yy) + εε (kk) 33 ψψ κκ (yy) PP (6)
5 168 Drakos Stefanos: Constitutive Relations of Stress and Strain in Stochastic Finite Element Method And finally = eeeeee(mmmm llllll ) (1 vv) EE ee σσ llllll yy εε (kk) 11 ψψ κκ (yy) + εε (kk) ψψ κκ (yy) + εε (kk) 33 ψψ κκ (yy) EE[pp] = eeeeee (mmmm llllll ) (1 vv) [εε (kk) ] vv ee σσ llllll yy ψψ κκ (yy)dddd (8) Variance of Pore Pressure Using the outcome of the pore pressure expected value its variance is given by: VVVVVV[pp] = EE[pp ] (EE[pp]) EE[pp ] = EE eeeeee (mmmm llllll ) ee σσyy (1 vv) (EE[pp]) = eeeeee (mmmm llllll ) (1 vv) [εε (kk) vv ] εε (kk) vv ψψ κκ (yy) ee σσ llllll yy ψψ κκ (yy)dddd 3.6. Invariants of Stress Tensor It is well known that to solve engineering problem the model itself should be defined independent of the coordinate system attached to the material. Thus, it is necessary to define the model in terms of stress invariants which are, by definition, independent of the coordinate system selected. The physical content of a stress tensor is reflected exclusively in the stress invariants (II 1, II, II 3 ). (9) II 1 = 1 3 σσ kkkk (3) II = 1 (σσ kkkk ) σσ iiii σσ iiii (31) ΙΙ 3 = εε iiiiii σσ ii1 σσ jj σσ kk3 (3) Solving a stochastic problem the statistical moments of stress invariants are needed. In the following paragraphs the calculation of the expected and variance value of each of invariants are presented Expected value of II 11 According to the linearity of the excepted value of II 1 can be calculated: ΕΕ[II 1 ] = ΕΕ[ 1 3 σσ qqqq ] and EE σσ qqqq = ΕΕ[II 1 ] = 1 3 EE σσ qqqq (kk) CC qqqqqqqq εε mmmm ρρ(yy)ee σσ llllll yy ψψ κκ (yy)dddd (33) CC = exp(mmmm llllll ) CC (34) Variance of II 11 Knowing the mean value of II 1 the variance is: VVVVVV(II 1 ) = EE[II 1 ] (EE[II 1 ]) But due to linearity: VVVVVV(II 1 ) = 1 9 (VVVVVV σσ qqqq ) = EE ee σσ llllll yy CC qqqqqqqq εε (kk) mmmm ψψ κκ (yy) EE[II 1 ] = 1 EE σσ 9 qqqq (EE[II 1 ]) = 1 EE[σσ 9 qqqq ] = (kk) CC qqqqqqqq εε mmmm (35) ρρ(yy)ee σσσσ ψψ κκ (yy)dddd
6 American Journal of Computational and Applied Mathematics 15, 5(6): Expected Value of II As shown before: EE[II ] = EE[ 1 σσ qqqq σσ iiii σσ iiii ] This gives EE[II ] = 1 EE eeσσ llllll yy CC qqqqqqqq εε (kk) mmmm ψψ κκ (yy) EE ee σσ llllll yy CC iiiiiiii εε (kk) mmmm ψψ κκ (yy) ee σσ llllll yy CC iiiiiiii εε (kk) mmmm ψψ κκ (yy) (36) Variance of II The variance of II due to the stress product on its equation become highly nonlinear. Thus VVVVVV(II ) = EE[II ] (EE[II ]) EE[II ] = 1 EE ee σσ llllll yy CC qqqqqqqq εε (kk) mmmm ψψ κκ (yy) 4 EE ee σσ llllll yy CC iiiiiiii εε (kk) mmmm ψψ κκ (yy) ee σσ llllll yy CC iiiiiiii εε (kk) mmmm ψψ κκ (yy) (EE[II ]) = 1 EE ee σσ llllll yy CC qqqqqqqq εε (kk) mmmm ψψ κκ (yy) EE ee σσ llllll yy CC iiiiiiii εε (kk) mmmm ψψ κκ (yy) ee σσ llllll yy CC iijjjjjj εε (kk) mmmm ψψ κκ (yy) (37) Expected Value of ΙΙ 33 The high non linearity presented also in the statistical moments of II 3. EE[ΙΙ 3 ] = EE[ee iiiiii σσ ii1 σσ jj σσ kk3 ] This leads to: EE[ΙΙ 3 ] = ee iiiiii EE ee 3σσ llllll yy CC ii1mmmm εε (kk) mmmm ψψ κκ (yy) CC jj mmmm εε (kk) mmmm ψψ κκ (yy) CC kk3mmmm εε (kk) mmmm ψψ κκ (yy) (38) Variance of ΙΙ 33 Similarly as documented above the variance of II 3 is equal to: VVVVVV(II 3 ) = EE[II 3 ] (EE[II 3 ]) EE ee iiiiii ee 3σσ llllll yy CC ii1mmmm εε (kk) mmmm ψψ κκ (yy) CC jj mmmm εε (kk) mmmm ψψ κκ (yy) CC kk3mmmm εε (kk) mmmm ψψ κκ (yy) EE ee iiiiii ee 3σσ llllll yy CC ii1mmmm εε (kk) mmmm ψψ κκ (yy) CC jj mmmm εε (kk) mmmm ψψ κκ (yy) CC kk3mmmm εε (kk) mmmm ψψ κκ (yy) (39) 4. Numerical Example A shallow foundation problem for various values of variation s coefficient vv ee is solved taken to account the randomness of the ground. To estimate the statistical moments of the soil deformation the numerical algorithm of SFEM using the Generalized Polynomial Chaos as described in the previous paragraphs is applied and the results are compared to those obtained by the closed form solution. To avoid the negative values of the elastic modulus assumed to have lognormal. It is known that the settlement beneath a foundation with uniform but random elastic modulus is given by the equation [8]: ss = ss dddddd μμ EE (4) EE ss dddddd is the deterministic value of settlement with EE = μμ ΕΕ everywhere. Assuming lognormal distribution for the settlements the mean values is equal to μμ llll (ss) = llll(ss dddddd ) + 1 σσ llll (EE) (41) The geometry of the finite elements used for the simulation of the problem presented in Fig.. The input data of the problem is the random field modulus with a constant
7 17 Drakos Stefanos: Constitutive Relations of Stress and Strain in Stochastic Finite Element Method average value equal to 1 Mpa and a fixed Poisson ratio equal to.5. Calculations have been made for ten different coefficients vv ee = σσ ΕΕ of the elastic modulus with a minimum μμ ΕΕ value of.1 and then with step.1 to a maximum value equal to 1. The randomness of Elasticity modulus in Fig. 3 is shown. For SFEM one dimensional Hermite GPC with order 5 [14] were used. In the Fig. B1 the results of SFEM method comparatively with the closed form solution are shown and they present great accuracy. It is observed that for of v e =.5 the error is equal to.8%. In the figures B-B13 the strains and stress components, the pore pressure and the stress tensor invariants are presented as resulted by the Chaos Polynomial expansion and compared with those raised by the Monte Carlo Method. Simulations of 1-5 samples were carried and the convergence of the outcomes decreases as the number of Monte Carlo simulations increases. presented. An algorithm of Stochastic Finite Element using Polynomial Chaos has been developed. An analysis of settlement of a plane strain strip load on an elastic foundation has been given as an example of the proposed approach. It is shown that the results of SFEM using polynomial chaos compare well with those obtained from closed form solution. The stress and strain constitutive relation the pore pressure and the stress invariants are modeled by the polynomial Chaos expansion and verified against Monte Carlo simulation which is treated as the exact solution of the problems. The main advantage in using the proposed methodology is that a large number of realizations which have to be made for (Random Finite Element Method) avoided, thus making the procedure viable for realistic practical problems. Appendix A Figure. Finite element mesh Figure 3. Modulus of Random Elasticity of two different realizations 5. Conclusions To propagate the uncertainties of input parameters to constitutive relations of strain and stress where arises due to spatial variability of mechanical parameters of soil/rock in geotechnical and geomechanics problems, a procedure of conducting a Stochastic Finite Element Analysis has been Galerkin approximation and Generalized Polynomial of chaos In order to solve the problem 1 we have to create the new space LL pp (, ΗΗ 1 (DD)). For that reason the subspace SS kk LL pp () is considered as [1]. SS kk = ssssssss{ψψ 1, ψψ,, ψψ κ } (A.1) Using the dyadic product of the space VV h, SS kk the space LL pp (, ΗΗ 1 (DD)) created. Thus VV hkk = VV h VV kk = ssssssss{φφ ii ψψ jj, i = 1 NN, j = 1, } (A.) The space VV hkk has dimension QN and regards the test function v. In the case where exists NN BB finite element supported by boundaries condition then the subspace of solution belongs is: WW hkk = VV hkk ssssssss{φφ NN+1, φφ NN+,, φφ ΝΝ+NNNN } (A.3) kk Assuming that the SS ii represents a space of univariate orthonormal polynomial of variable yy ii ιι R with order k or lower and: SS kk ii = ssssssss PP ii aaii (yy ii ), aa ii =,1,, kk, ii = 1, MM (A.4) kk The tensor product of the M SS ii subspace results the space of the Generalized Polynomial Chaos: SS kk = SS 1 SS SS MM (A.5) And using (A4) SS kk MM = ssssssss PP aa ii ii (yy ii ): aa ii =,1, kk, ii = 1 MM, ii=1 (A.6) aa kk Where aa = And MM ii=1 aa ii = dim(ss kk ) = (MM+kk)! (A.7) MM!kk! Xiu & Karniadakis [14] show the application of the method for different kind of orthonormal polynomials and in the current paper the Hermite polynomial was used with the following characteristics:
8 American Journal of Computational and Applied Mathematics 15, 5(6): PP = 1, < PP ii >=, ii > < PP mm PP nn >= PP mm (yy)pp nn (yy)ρρ(yy)ddyy = γγ nn δδ mmmm (A.8) γγ nn =< PP nn > : are the normalization factors, δδ mmmm is the Kronecker delta ρρ(yy) = 1 yy : is the density function and ππ ee PP nn = ( 1) nn ee yy ddnn ddyy nn ee yy (A.9) For a 3 rd order of one dimension of uncertainty the Hermite Polynomial Chaos is given by: ψψ οο (yy) = PP (yy) = 1, ψψ 1 (yy) = PP 1 (yy) = yy, ψψ (yy) = PP (yy) = yy 1 Appendix B Results of Numerical Example 1 1 Std[I1] Figure B3. Standard deviation of stress tensor complements. (MC 1 E[strain] E[ε],Chaos E[ε],Random 1 E[ε11],Chaos E[ε11],Random std[σ], Chaos std[σ],random std[σ11], Chaos std[σ11],random E[Stress] Figure B1. Closed form solution and SFEM results E[σ],GPC E[σ],Random E[σ11],GPC E[σ11],Random Figure B4. Expected value of strain tensor complements. (MC 1 Std[Srtain] Std[ε],Chaos Std[ε],Random 1 Std[ε11], Chaos Std[ε11],Random Figure B. Expected value of stress tensor complements (MC 1 Figure B5. Standard deviation of strain tensor complements. (MC 1
9 17 Drakos Stefanos: Constitutive Relations of Stress and Strain in Stochastic Finite Element Method E[p] Figure B6. Expected value of pore pressure. (MC 1 std[p] Figure B7. Standard deviation of pore pressure. (MC 1 E[I1] Figure B8. Expected value of stress tensor invariant II 1. (MC 1, 3 E[p],chaos E[p],random-1 samples std[p],chaos std[p],random 1 samples E[I1],chaos E[I1],random 1 samples E[I1],random 3 samples std[i1] Figure B9. Standard deviation of stress tensor invariant II 1. (MC 3 E[I] std[i1], Chaos std[i1], Random 3 samples Figure B1. Expected value of stress tensor invariant II. (MC 1, 5 Std[I] Figure B11. Standard deviation of stress tensor invariant II. (MC 1, 5 E[I],chaos E[I],random, 1 samples E[I],random, 5 samples std[i],chaos std[i],random, 1 samples std[i],random,5 samples
10 American Journal of Computational and Applied Mathematics 15, 5(6): E[I3] Figure B1. Expected value of stress tensor invariant II 3. (MC 1, 5 Std[I3] Figure B13. Standard deviation of stress tensor invariant II 3. (MC 1, 5 REFERENCES E[I3],Chaos E[I3],Random 1 E[I3],Random [1] Phoon, K., Quek, S., Chow, Y. and Lee, S (199), Reliability analysis of pile settlements. Journal of Geotechnical Engineering, ASCE, 116(11), std[i3],chaos std[i3],random 1 std[i3],random 5 [] Mellah, R., Auvinet, G. and Masrouri, F. (), Stochastic finite element method applied to non-linear analysis of embankments. Probabilistic Engineering Mechanics, 15 (3), [3] Eloseily, K., Ayyub, B. and Patev, R. (), Reliability assessment of pile groups in sands. Journal of Structural Engineering, ASCE, 18(1), [4] Fenton, G. A., and Vanmarcke, E. H. (199), Simulation of random fields via local average subdivision. J. Eng. Mech., 116(8), [5] Paice, G. M., Griffiths, D. V., and Fenton, G. A. (1996), Finite element modeling of settlements on spatially random soil. J. Geotech. Eng., 1(9), [6] Fenton, G. A., and Griffiths, D. V. (), Probabilistic foundation settlement on spatially random soil. J. Geotech. Geoenviron. Eng.,18(5), [7] Fenton, G. A., and Griffiths, D. V. (5), Three-dimensional probabilistic foundation settlement. J. Geotech. Geoenviron. Eng., 131(), [8] Fenton, G. A., and Griffiths, D. V. (8), Risk assessment in geotechnical engineering, Wiley, Hoboken, NJ. [9] Ghosh Debraj and Farhat Charbel (8), Strain and stress computations in stochastic finite element methods, International Journal for Numerical Method In Engineering, Vol. 74, Iss. 8, [1] Lord G., Powel C., Shardlow T. (14), An Introduction to Computational Stochastic PDEs, Cambridge Texts in Applied Mathematics. [11] Ghanem, R. G., and Spanos, P. D. (1991), Stochastic finite elements: A spectral approach, Springer-rlag, New York. [1] Drakos I. S. & Pande G. N. Quantitative of uncertainties in Earth Structures, World Congress on Advances in Structural Engineering and Mechanics, Incheon Korea 5-9 August, 15. [13] Drakos I. S. & Pande G. N. Stochastic Finite Element Analysis using Polynomial Chaos Studia Geotechnica et Mechanica Journal.(Under Review). [14] D Xiu, G Em Karniadakis. (3), Modeling uncertainty in steady state diffusion problems via generalized polynomial chaos Computer Methods in Applied Mechanics and Engineering 191 (43), [15] Griffiths, D.V (1985). The effect of pore-fluid compressibility on failure loads in elasto-plastic soil. Int J Numer Anal Methods Geomech, vol.9, no.3, pp
Lecture No. 5. For all weighted residual methods. For all (Bubnov) Galerkin methods. Summary of Conventional Galerkin Method
Lecture No. 5 LL(uu) pp(xx) = 0 in ΩΩ SS EE (uu) = gg EE on ΓΓ EE SS NN (uu) = gg NN on ΓΓ NN For all weighted residual methods NN uu aaaaaa = uu BB + αα ii φφ ii For all (Bubnov) Galerkin methods ii=1
More information14- Hardening Soil Model with Small Strain Stiffness - PLAXIS
14- Hardening Soil Model with Small Strain Stiffness - PLAXIS This model is the Hardening Soil Model with Small Strain Stiffness as presented in PLAXIS. The model is developed using the user-defined material
More informationA Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions
Lin Lin A Posteriori DG using Non-Polynomial Basis 1 A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions Lin Lin Department of Mathematics, UC Berkeley;
More informationIndependent Component Analysis and FastICA. Copyright Changwei Xiong June last update: July 7, 2016
Independent Component Analysis and FastICA Copyright Changwei Xiong 016 June 016 last update: July 7, 016 TABLE OF CONTENTS Table of Contents...1 1. Introduction.... Independence by Non-gaussianity....1.
More informationAngular Momentum, Electromagnetic Waves
Angular Momentum, Electromagnetic Waves Lecture33: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay As before, we keep in view the four Maxwell s equations for all our discussions.
More information(1) Introduction: a new basis set
() Introduction: a new basis set In scattering, we are solving the S eq. for arbitrary VV in integral form We look for solutions to unbound states: certain boundary conditions (EE > 0, plane and spherical
More informationModule 4 (Lecture 16) SHALLOW FOUNDATIONS: ALLOWABLE BEARING CAPACITY AND SETTLEMENT
Topics Module 4 (Lecture 16) SHALLOW FOUNDATIONS: ALLOWABLE BEARING CAPACITY AND SETTLEMENT 1.1 STRIP FOUNDATION ON GRANULAR SOIL REINFORCED BY METALLIC STRIPS Mode of Failure Location of Failure Surface
More informationTECHNICAL NOTE AUTOMATIC GENERATION OF POINT SPRING SUPPORTS BASED ON DEFINED SOIL PROFILES AND COLUMN-FOOTING PROPERTIES
COMPUTERS AND STRUCTURES, INC., FEBRUARY 2016 TECHNICAL NOTE AUTOMATIC GENERATION OF POINT SPRING SUPPORTS BASED ON DEFINED SOIL PROFILES AND COLUMN-FOOTING PROPERTIES Introduction This technical note
More informationVariations. ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra
Variations ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra Last Time Probability Density Functions Normal Distribution Expectation / Expectation of a function Independence Uncorrelated
More informationModule 7 (Lecture 25) RETAINING WALLS
Module 7 (Lecture 25) RETAINING WALLS Topics Check for Bearing Capacity Failure Example Factor of Safety Against Overturning Factor of Safety Against Sliding Factor of Safety Against Bearing Capacity Failure
More informationWorksheets for GCSE Mathematics. Algebraic Expressions. Mr Black 's Maths Resources for Teachers GCSE 1-9. Algebra
Worksheets for GCSE Mathematics Algebraic Expressions Mr Black 's Maths Resources for Teachers GCSE 1-9 Algebra Algebraic Expressions Worksheets Contents Differentiated Independent Learning Worksheets
More informationThick Shell Element Form 5 in LS-DYNA
Thick Shell Element Form 5 in LS-DYNA Lee P. Bindeman Livermore Software Technology Corporation Thick shell form 5 in LS-DYNA is a layered node brick element, with nodes defining the boom surface and defining
More informationThe Elasticity of Quantum Spacetime Fabric. Viorel Laurentiu Cartas
The Elasticity of Quantum Spacetime Fabric Viorel Laurentiu Cartas The notion of gravitational force is replaced by the curvature of the spacetime fabric. The mass tells to spacetime how to curve and the
More informationON CALCULATION OF EFFECTIVE ELASTIC PROPERTIES OF MATERIALS WITH CRACKS
Materials Physics and Mechanics 32 (2017) 213-221 Received: November 7, 2017 ON CALCULATION OF EFFECTIVE ELASTIC PROPERTIES OF MATERIALS WITH CRACKS Ruslan L. Lapin 1, Vitaly A. Kuzkin 1,2 1 Peter the
More information(2) Orbital angular momentum
(2) Orbital angular momentum Consider SS = 0 and LL = rr pp, where pp is the canonical momentum Note: SS and LL are generators for different parts of the wave function. Note: from AA BB ii = εε iiiiii
More informationModule 7 (Lecture 27) RETAINING WALLS
Module 7 (Lecture 27) RETAINING WALLS Topics 1.1 RETAINING WALLS WITH METALLIC STRIP REINFORCEMENT Calculation of Active Horizontal and vertical Pressure Tie Force Factor of Safety Against Tie Failure
More informationSupport Vector Machines. CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington
Support Vector Machines CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington 1 A Linearly Separable Problem Consider the binary classification
More informationPHY103A: Lecture # 4
Semester II, 2017-18 Department of Physics, IIT Kanpur PHY103A: Lecture # 4 (Text Book: Intro to Electrodynamics by Griffiths, 3 rd Ed.) Anand Kumar Jha 10-Jan-2018 Notes The Solutions to HW # 1 have been
More informationCHAPTER 5 Wave Properties of Matter and Quantum Mechanics I
CHAPTER 5 Wave Properties of Matter and Quantum Mechanics I 1 5.1 X-Ray Scattering 5.2 De Broglie Waves 5.3 Electron Scattering 5.4 Wave Motion 5.5 Waves or Particles 5.6 Uncertainty Principle Topics 5.7
More informationCreeping Movement across a Long Strike-Slip Fault in a Half Space of Linear Viscoelastic Material Representing the Lithosphere-Asthenosphere System
Frontiers in Science 214, 4(2): 21-28 DOI: 1.5923/j.fs.21442.1 Creeping Movement across a Long Strike-Slip Fault in a Half Space of Linear Viscoelastic Material Representing the Lithosphere-Asthenosphere
More informationProf. Dr.-Ing. Armin Dekorsy Department of Communications Engineering. Stochastic Processes and Linear Algebra Recap Slides
Prof. Dr.-Ing. Armin Dekorsy Department of Communications Engineering Stochastic Processes and Linear Algebra Recap Slides Stochastic processes and variables XX tt 0 = XX xx nn (tt) xx 2 (tt) XX tt XX
More informationGrover s algorithm. We want to find aa. Search in an unordered database. QC oracle (as usual) Usual trick
Grover s algorithm Search in an unordered database Example: phonebook, need to find a person from a phone number Actually, something else, like hard (e.g., NP-complete) problem 0, xx aa Black box ff xx
More informationHaar Basis Wavelets and Morlet Wavelets
Haar Basis Wavelets and Morlet Wavelets September 9 th, 05 Professor Davi Geiger. The Haar transform, which is one of the earliest transform functions proposed, was proposed in 90 by a Hungarian mathematician
More informationChapter 22 : Electric potential
Chapter 22 : Electric potential What is electric potential? How does it relate to potential energy? How does it relate to electric field? Some simple applications What does it mean when it says 1.5 Volts
More informationSoil Uncertainty and Seismic Ground Motion
Boris Jeremić and Kallol Sett Department of Civil and Environmental Engineering University of California, Davis GEESD IV May 2008 Support by NSF CMMI 0600766 Outline Motivation Soils and Uncertainty Stochastic
More informationTopics. Module 3 Lecture 10 SHALLOW FOUNDATIONS: ULTIMATE BEARING CAPACITY NPTEL ADVANCED FOUNDATION ENGINEERING-I
Topics Module 3 Lecture 10 SHALLOW FOUNDATIONS: ULTIMATE BEARING CAPACITY 1.1 THE GENERAL BEARING CAPACITY EQUATION Bearing Capacity Factors General Comments 1.2 EFFECT OF SOIL COMPRESSIBILITY 1.3 ECCENTRICALLY
More information(1) Correspondence of the density matrix to traditional method
(1) Correspondence of the density matrix to traditional method New method (with the density matrix) Traditional method (from thermal physics courses) ZZ = TTTT ρρ = EE ρρ EE = dddd xx ρρ xx ii FF = UU
More informationReview for Exam Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa
57:020 Fluids Mechanics Fall2013 1 Review for Exam3 12. 11. 2013 Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa 57:020 Fluids Mechanics Fall2013 2 Chapter
More informationOn The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time Delays
Journal of Mathematics and System Science 6 (216) 194-199 doi: 1.17265/2159-5291/216.5.3 D DAVID PUBLISHING On The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time
More informationSECTION 5: CAPACITANCE & INDUCTANCE. ENGR 201 Electrical Fundamentals I
SECTION 5: CAPACITANCE & INDUCTANCE ENGR 201 Electrical Fundamentals I 2 Fluid Capacitor Fluid Capacitor 3 Consider the following device: Two rigid hemispherical shells Separated by an impermeable elastic
More informationJasmin Smajic1, Christian Hafner2, Jürg Leuthold2, March 23, 2015
Jasmin Smajic, Christian Hafner 2, Jürg Leuthold 2, March 23, 205 Time Domain Finite Element Method (TD FEM): Continuous and Discontinuous Galerkin (DG-FEM) HSR - University of Applied Sciences of Eastern
More informationMacroelement Modeling for Single Vertical Piles
6 th International Conference on Earthquake Geotechnical Engineering 1-4 November 2015 Christchurch, New Zealand Macroelement Modeling for Single Vertical Piles Z. Li 1,2, P. Kotronis 1, S. Escoffier 2,
More informationControl of Mobile Robots
Control of Mobile Robots Regulation and trajectory tracking Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Organization and
More information7.3 The Jacobi and Gauss-Seidel Iterative Methods
7.3 The Jacobi and Gauss-Seidel Iterative Methods 1 The Jacobi Method Two assumptions made on Jacobi Method: 1.The system given by aa 11 xx 1 + aa 12 xx 2 + aa 1nn xx nn = bb 1 aa 21 xx 1 + aa 22 xx 2
More informationON INFLUENCE OF SHEAR TRACTION ON HYDRAULIC FRACTURE PROPAGATION
Materials Physics and Mechanics 32 (2017) 272-277 Received: November 2, 2017 ON INFLUENCE OF SHEAR TRACTION ON HYDRAULIC FRACTURE PROPAGATION A. M. Linkov 1, 2 1 Rzeszow University of Technology, Poland
More informationGradient expansion formalism for generic spin torques
Gradient expansion formalism for generic spin torques Atsuo Shitade RIKEN Center for Emergent Matter Science Atsuo Shitade, arxiv:1708.03424. Outline 1. Spintronics a. Magnetoresistance and spin torques
More informationIntroduction to Kinetic Simulation of Magnetized Plasma
Introduction to Kinetic Simulation of Magnetized Plasma Jae-Min Kwon National Fusion Research Institute, Korea 018 EASW8 July 30 Aug 3, 018 1 Outline Introduction to kinetic plasma model Very brief on
More informationStochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations
Stochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations Boris Jeremić 1 Kallol Sett 2 1 University of California, Davis 2 University of Akron, Ohio ICASP Zürich, Switzerland August
More informationMath, Stats, and Mathstats Review ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD
Math, Stats, and Mathstats Review ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD Outline These preliminaries serve to signal to students what tools they need to know to succeed in ECON 360 and refresh their
More informationModelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases
Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of
More informationWorksheets for GCSE Mathematics. Quadratics. mr-mathematics.com Maths Resources for Teachers. Algebra
Worksheets for GCSE Mathematics Quadratics mr-mathematics.com Maths Resources for Teachers Algebra Quadratics Worksheets Contents Differentiated Independent Learning Worksheets Solving x + bx + c by factorisation
More informationMathematics Ext 2. HSC 2014 Solutions. Suite 403, 410 Elizabeth St, Surry Hills NSW 2010 keystoneeducation.com.
Mathematics Ext HSC 4 Solutions Suite 43, 4 Elizabeth St, Surry Hills NSW info@keystoneeducation.com.au keystoneeducation.com.au Mathematics Extension : HSC 4 Solutions Contents Multiple Choice... 3 Question...
More informationLast Name _Piatoles_ Given Name Americo ID Number
Last Name _Piatoles_ Given Name Americo ID Number 20170908 Question n. 1 The "C-V curve" method can be used to test a MEMS in the electromechanical characterization phase. Describe how this procedure is
More informationWave Motion. Chapter 14 of Essential University Physics, Richard Wolfson, 3 rd Edition
Wave Motion Chapter 14 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 Waves: propagation of energy, not particles 2 Longitudinal Waves: disturbance is along the direction of wave propagation
More informationNumerical Solution of Fredholm and Volterra Integral Equations of the First Kind Using Wavelets bases
The Journal of Mathematics and Computer Science Available online at http://www.tjmcs.com The Journal of Mathematics and Computer Science Vol.5 No.4 (22) 337-345 Numerical Solution of Fredholm and Volterra
More informationGeneral Strong Polarization
General Strong Polarization Madhu Sudan Harvard University Joint work with Jaroslaw Blasiok (Harvard), Venkatesan Gurswami (CMU), Preetum Nakkiran (Harvard) and Atri Rudra (Buffalo) May 1, 018 G.Tech:
More informationAdvanced data analysis
Advanced data analysis Akisato Kimura ( 木村昭悟 ) NTT Communication Science Laboratories E-mail: akisato@ieee.org Advanced data analysis 1. Introduction (Aug 20) 2. Dimensionality reduction (Aug 20,21) PCA,
More informationReview for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa
Review for Exam3 12. 9. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University
More informationIntegrating Rational functions by the Method of Partial fraction Decomposition. Antony L. Foster
Integrating Rational functions by the Method of Partial fraction Decomposition By Antony L. Foster At times, especially in calculus, it is necessary, it is necessary to express a fraction as the sum of
More informationWork, Energy, and Power. Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition
Work, Energy, and Power Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 With the knowledge we got so far, we can handle the situation on the left but not the one on the right.
More informationMutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene and carbon tetrachloride
- 1 - PLMMP 2016, Kyiv Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene and carbon tetrachloride Jadran Vrabec Tatjana Janzen, Gabriela Guevara-Carrión,
More informationChapter 10. Simple Harmonic Motion and Elasticity
Chapter 10 Simple Harmonic Motion and Elasticity 10.1 The Ideal Spring and Simple Harmonic Motion F Applied x = k x spring constant Units: N/m FF SSSSSSSSSSSS = kkkk 10.1 The Ideal Spring and Simple Harmonic
More informationConharmonically Flat Vaisman-Gray Manifold
American Journal of Mathematics and Statistics 207, 7(): 38-43 DOI: 0.5923/j.ajms.207070.06 Conharmonically Flat Vaisman-Gray Manifold Habeeb M. Abood *, Yasir A. Abdulameer Department of Mathematics,
More informationSolving Fuzzy Nonlinear Equations by a General Iterative Method
2062062062062060 Journal of Uncertain Systems Vol.4, No.3, pp.206-25, 200 Online at: www.jus.org.uk Solving Fuzzy Nonlinear Equations by a General Iterative Method Anjeli Garg, S.R. Singh * Department
More informationA reduced-order stochastic finite element analysis for structures with uncertainties
A reduced-order stochastic finite element analysis for structures with uncertainties Ji Yang 1, Béatrice Faverjon 1,2, Herwig Peters 1, icole Kessissoglou 1 1 School of Mechanical and Manufacturing Engineering,
More informationApproximate Second Order Algorithms. Seo Taek Kong, Nithin Tangellamudi, Zhikai Guo
Approximate Second Order Algorithms Seo Taek Kong, Nithin Tangellamudi, Zhikai Guo Why Second Order Algorithms? Invariant under affine transformations e.g. stretching a function preserves the convergence
More informationSolar Photovoltaics & Energy Systems
Solar Photovoltaics & Energy Systems Lecture 3. Solar energy conversion with band-gap materials ChE-600 Kevin Sivula, Spring 2014 The Müser Engine with a concentrator T s Q 1 = σσ CffT ss 4 + 1 Cff T pp
More informationAn Enhanced Assumed Strain (EAS) Solid Element for Nonlinear Implicit Analyses
An Enhanced Assumed Strain (EAS) Solid Element for Nonlinear Implicit Analyses Fredrik Bengzon ϯ, Thomas Borrvall ϯ and Ushnish Basu ϯϯ ϯ DYNAmore Nordic AB ϯϯ Livermore Software Technology Corporation
More informationPHY103A: Lecture # 9
Semester II, 2017-18 Department of Physics, IIT Kanpur PHY103A: Lecture # 9 (Text Book: Intro to Electrodynamics by Griffiths, 3 rd Ed.) Anand Kumar Jha 20-Jan-2018 Summary of Lecture # 8: Force per unit
More informationCharge carrier density in metals and semiconductors
Charge carrier density in metals and semiconductors 1. Introduction The Hall Effect Particles must overlap for the permutation symmetry to be relevant. We saw examples of this in the exchange energy in
More informationS. Freitag, B. T. Cao, J. Ninić & G. Meschke
SFB 837 Interaction Modeling Mechanized Tunneling S Freitag, B T Cao, J Ninić & G Meschke Institute for Structural Mechanics Ruhr University Bochum 1 Content 1 Motivation 2 Process-oriented FE model for
More informationPHL424: Nuclear Shell Model. Indian Institute of Technology Ropar
PHL424: Nuclear Shell Model Themes and challenges in modern science Complexity out of simplicity Microscopic How the world, with all its apparent complexity and diversity can be constructed out of a few
More informationSECTION 8: ROOT-LOCUS ANALYSIS. ESE 499 Feedback Control Systems
SECTION 8: ROOT-LOCUS ANALYSIS ESE 499 Feedback Control Systems 2 Introduction Introduction 3 Consider a general feedback system: Closed-loop transfer function is KKKK ss TT ss = 1 + KKKK ss HH ss GG ss
More informationFermi Surfaces and their Geometries
Fermi Surfaces and their Geometries Didier Ndengeyintwali Physics Department, Drexel University, Philadelphia, Pennsylvania 19104, USA (Dated: May 17, 2010) 1. Introduction The Pauli exclusion principle
More informationCALCULATION OF A SHEET PILE WALL RELIABILITY INDEX IN ULTIMATE AND SERVICEABILITY LIMIT STATES
Studia Geotechnica et Mechanica, Vol. XXXII, No. 2, 2010 CALCULATION OF A SHEET PILE WALL RELIABILITY INDEX IN ULTIMATE AND SERVICEABILITY LIMIT STATES JERZY BAUER Institute of Mining, Wrocław University
More informationReview for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa
Review for Exam2 11. 13. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University
More informationMath 171 Spring 2017 Final Exam. Problem Worth
Math 171 Spring 2017 Final Exam Problem 1 2 3 4 5 6 7 8 9 10 11 Worth 9 6 6 5 9 8 5 8 8 8 10 12 13 14 15 16 17 18 19 20 21 22 Total 8 5 5 6 6 8 6 6 6 6 6 150 Last Name: First Name: Student ID: Section:
More informationLecture 3. STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher
Lecture 3 STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher Previous lectures What is machine learning? Objectives of machine learning Supervised and
More informationGeneral Strong Polarization
General Strong Polarization Madhu Sudan Harvard University Joint work with Jaroslaw Blasiok (Harvard), Venkatesan Gurswami (CMU), Preetum Nakkiran (Harvard) and Atri Rudra (Buffalo) December 4, 2017 IAS:
More informationAnalytical procedure for the design of PFRP-RC hybrid beams including shear interaction effects
Analytical procedure for the design of PFRP-RC hybrid beams including shear interaction effects Catalin A. Neagoe *, Lluís Gil Department of Strength of Materials and Structural Engineering, Universitat
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationClassical RSA algorithm
Classical RSA algorithm We need to discuss some mathematics (number theory) first Modulo-NN arithmetic (modular arithmetic, clock arithmetic) 9 (mod 7) 4 3 5 (mod 7) congruent (I will also use = instead
More informationMATHEMATICAL MODEL OF ULTRASOUND TRANSMISSION IN GRADIENT MATERIALS SYNTHESIZED IN ARC SURFACING UNDER CONTACT LOADING
Materials Physics and Mechanics 32 (207) 86-93 Received: May 7, 207 MATHEMATICAL MODEL OF ULTRASOUND TRANSMISSION IN GRADIENT MATERIALS SYNTHESIZED IN ARC SURFACING UNDER CONTACT LOADING V.D. Sarychev,
More informationProblem 3.1 (Verdeyen 5.13) First, I calculate the ABCD matrix for beam traveling through the lens and space.
Problem 3. (Verdeyen 5.3) First, I calculate the ABCD matrix for beam traveling through the lens and space. T = dd 0 0 dd 2 ff 0 = dd 2 dd ff 2 + dd ( dd 2 ff ) dd ff ff Aording to ABCD law, we can have
More information2.4 Error Analysis for Iterative Methods
2.4 Error Analysis for Iterative Methods 1 Definition 2.7. Order of Convergence Suppose {pp nn } nn=0 is a sequence that converges to pp with pp nn pp for all nn. If positive constants λλ and αα exist
More informationUsing the Application Builder for Neutron Transport in Discrete Ordinates
Using the Application Builder for Neutron Transport in Discrete Ordinates C.J. Hurt University of Tennessee Nuclear Engineering Department (This material is based upon work supported under a Department
More informationComparison of Slope Reliability Methods of Analysis
Comparison of Slope Reliability Methods of Analysis D.V. Griffiths,F ASCE, Jinsong Huang,M ASCE and Gordon A. Fenton 3,M ASCE Professor, Colorado School of Mines, 60 Illinois Street, Golden, CO 8040; d.v.griffiths@mines.edu
More informationA Study of Probabilistic FEMs for a Slope Reliability Analysis Using the Stress Fields
Send Orders for Reprints to reprints@benthamscience.ae 96 The Open Civil Engineering Journal, 05, 9, 96-06 Open Access A Study of Probabilistic FEMs for a Slope Reliability Analysis Using the Stress Fields
More informationRadiation. Lecture40: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay
Radiation Zone Approximation We had seen that the expression for the vector potential for a localized cuent distribution is given by AA (xx, tt) = μμ 4ππ ee iiiiii dd xx eeiiii xx xx xx xx JJ (xx ) In
More informationExam 2 Fall 2015
1 95.144 Exam 2 Fall 2015 Section instructor Section number Last/First name Last 3 Digits of Student ID Number: Show all work. Show all formulas used for each problem prior to substitution of numbers.
More informationModeling Damage in Composite Materials Using an Enrichment Based Multiscale Method
AD Technical Report ARWSB-TR-15002 Modeling Damage in Composite Materials Using an Enrichment Based Multiscale Method Michael F. Macri Andrew G. Littlefield March 2015 ARMAMENT RESEARCH, DEVELOPMENT AND
More informationNon-Parametric Weighted Tests for Change in Distribution Function
American Journal of Mathematics Statistics 03, 3(3: 57-65 DOI: 0.593/j.ajms.030303.09 Non-Parametric Weighted Tests for Change in Distribution Function Abd-Elnaser S. Abd-Rabou *, Ahmed M. Gad Statistics
More informationCHAPTER 4 Structure of the Atom
CHAPTER 4 Structure of the Atom Fall 2018 Prof. Sergio B. Mendes 1 Topics 4.1 The Atomic Models of Thomson and Rutherford 4.2 Rutherford Scattering 4.3 The Classic Atomic Model 4.4 The Bohr Model of the
More informationUncertainty analysis of large-scale systems using domain decomposition
Center for Turbulence Research Annual Research Briefs 2007 143 Uncertainty analysis of large-scale systems using domain decomposition By D. Ghosh, C. Farhat AND P. Avery 1. Motivation and objectives A
More informationPLANE WAVE AND FUNDAMENTAL SOLUTION IN THERMOPOROELASTIC MEDIUM
Materials Physics and Mechanics 35 (08 0- Received: June 7, 05 PLANE WAVE AND FUNDAMENTAL SOLUTION IN THERMOPOROELASTIC MEDIUM R. Kumar, S. Kumar *, M.G. Gourla 3 Department of Mathematics, Kurukshetra
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationTime Domain Analysis of Linear Systems Ch2. University of Central Oklahoma Dr. Mohamed Bingabr
Time Domain Analysis of Linear Systems Ch2 University of Central Oklahoma Dr. Mohamed Bingabr Outline Zero-input Response Impulse Response h(t) Convolution Zero-State Response System Stability System Response
More informationMaterials & Advanced Manufacturing (M&AM)
Modeling of Shear Thickening Fluids for Analysis of Energy Absorption Under Impulse Loading Alyssa Bennett (University of Michigan) Nick Vlahopoulos, PhD (University of Michigan) Weiran Jiang, PhD (Research
More informationUnit II. Page 1 of 12
Unit II (1) Basic Terminology: (i) Exhaustive Events: A set of events is said to be exhaustive, if it includes all the possible events. For example, in tossing a coin there are two exhaustive cases either
More informationPhysics 371 Spring 2017 Prof. Anlage Review
Physics 71 Spring 2017 Prof. Anlage Review Special Relativity Inertial vs. non-inertial reference frames Galilean relativity: Galilean transformation for relative motion along the xx xx direction: xx =
More informationLaterally Loaded Piles. Rocscience 2018
Laterally Loaded Piles Rocscience 2018 Contents 1 Soil Response Modelled by p-y Curves... 3 2 Governing Differential Equation... 4 3 Finite Element Method... 5 4 Pile Bending Stiffness... 5 5 Soil Models...
More informationNUMERICAL ANALYSIS OF BEAMS ON RANDOM ELASTIC FOUNDATIONS
NUMERICAL ANALYSIS OF BEAMS ON RANDOM ELASTIC FOUNDATIONS D.. Griffiths Jumpol Paiboon Jinsong Huang d.v.griffiths@mines.edu jpaiboon@mines.edu jhuang@mines.edu Geomechanics Research Center, Colorado School
More informationIntegrated Attitude Determination and Control System via Magnetic Measurements and Actuation
Proceedings of the th WSEAS International Conference on Automatic Control, Modelling and Simulation Integrated Attitude Determination and Control System via Magnetic Measurements and Actuation MOHAMMAD
More informationLecture No. 1 Introduction to Method of Weighted Residuals. Solve the differential equation L (u) = p(x) in V where L is a differential operator
Lecture No. 1 Introduction to Method of Weighted Residuals Solve the differential equation L (u) = p(x) in V where L is a differential operator with boundary conditions S(u) = g(x) on Γ where S is a differential
More informationReliability of Traditional Retaining Wall Design
Reliability of Traditional Retaining Wall Design by Gordon A. Fenton 1, D. V. Griffiths 2, and M. B. Williams 3 in Géotechique, Vol. 55, No. 1, pp. 55-62, 2005 Keywords: retaining walls, earth pressure,
More informationSECTION 5: POWER FLOW. ESE 470 Energy Distribution Systems
SECTION 5: POWER FLOW ESE 470 Energy Distribution Systems 2 Introduction Nodal Analysis 3 Consider the following circuit Three voltage sources VV sss, VV sss, VV sss Generic branch impedances Could be
More informationDiscrete and Three-Parameter Models of Hydraulic Fracture Crack Growth
Discrete and Three-Parameter Models of Hydraulic Fracture Crack Growth SERGEY KANAUN Department of Mechatronics, Monterrey Institute of Technology and Higher Education, State of Mexico Campus, Carretera
More informationRotational Motion. Chapter 10 of Essential University Physics, Richard Wolfson, 3 rd Edition
Rotational Motion Chapter 10 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 We ll look for a way to describe the combined (rotational) motion 2 Angle Measurements θθ ss rr rrrrrrrrrrrrrr
More informationThermoelastic Behavior of a Small Bearing Cage
Clemson University TigerPrints All Theses Theses 12-2009 Thermoelastic Behavior of a Small Bearing Cage Eric Planting Clemson University, eplanti@clemson.edu Follow this and additional works at: https://tigerprints.clemson.edu/all_theses
More informationADAPTIVE NEURAL NETWORK CONTROLLER DESIGN FOR BLENDED-WING UAV WITH COMPLEX DAMAGE
ADAPTIVE NEURAL NETWORK CONTROLLER DESIGN FOR BLENDED-WING UAV WITH COMPLEX DAMAGE Kijoon Kim*, Jongmin Ahn**, Seungkeun Kim*, Jinyoung Suk* *Chungnam National University, **Agency for Defense and Development
More information