Spectral methods for fuzzy structural dynamics: modal vs direct approach
|
|
- Job Cooper
- 5 years ago
- Views:
Transcription
1 Spectral methods for fuzzy structural dynamics: modal vs direct approach S Adhikari Zienkiewicz Centre for Computational Engineering, College of Engineering, Swansea University, Wales, UK IUTAM Symposium - Dynamical Analysis of Multibody Systems with Design Uncertainties, Stuttgart, Germany
2 Outline of the talk 1 Introduction 2 Brief overview of fuzzy variables 3 Legendre polynomials for fuzzy uncertainty propagation Problem description Legendre polynomials Projection in the basis of Legendre polynomials 4 Fuzzy finite element analysis Formulation of the problem Galerkin approach 5 Fuzzy eigen analysis using spectral projection 6 Numerical example: axial vibration of a rod 7 Summary and conclusion
3 Fuzzy variables The consideration of uncertainties in numerical models to obtain the variability of response is becoming more common for finite element models arising in practical problems. When substantial statistical information exists, the theory of probability and stochastic processes offer a rich mathematical framework to represent such uncertainties. In a probabilistic setting, uncertainty associated with the system parameters, such as the geometric properties and constitutive relations (i.e. Young s modulus, mass density, Poisson s ratio, damping coefficients), can be modeled as random variables or stochastic processes. These uncertainties can be quantified and propagated, for example, using the stochastic finite element method
4 Fuzzy variables The reliable application of probabilistic approaches requires information to construct the probability density functions of uncertain parameters. This information may not be easily available for many complex practical problems. In such situations, non-probabilistic approaches such as interval algebra, convex modelsand fuzzy set based methods can be used. Fuzzy finite element analysis aims to combine the power of finite element method and uncertainty modelling capability of fuzzy variables. In the context of computational mechanics, the aim of a fuzzy finite element analysis is to obtain the range of certain output quantities (such as displacement, acceleration and stress) given the membership of data in the set of input variables. This problem, known as the uncertainty propagation problem, has taken the centre stage in recent research activities in the field.
5 Membership functions of a fuzzy variable Membership functions of a fuzzy variable; (a) symmetric triangular membership function; (b) asymmetric triangular membership function; (c) general membership function. The value corresponding to α = 1 is the crisp (or deterministic) value. The range corresponding to α = 0 is the widest interval. Any intermediate value of 0 < α < 1, yields an interval variable with a finite lower and upper bound.
6 Uncertainty propagation We are interested in the propagation of a m-dimensional vector of Fuzzy variables, which can be expressed as y = f(x) R n (1) where f( ) : R m R n is a smooth nonlinear function of the input fuzzy vector x. The fuzzy uncertainty propagation problem can be formally defined as the solution of the following set of optimisation problems ) y jkmin = min( fj (x αk ) ) j = 1, 2,...,n; k = 1, 2,...,r (2) y jkmax = max( fj (x αk ) Two ways the solution can be made efficient: (1) evaluate f(x) in Eq (1) fast, and/or (2) do a better job with the optimisation in Eq (2). The main idea proposed here is based on efficient and accurate construction of a (hopefully cheaper!) response surface for every α-cut followed by an optimization approach.
7 Transformation of a fuzzy variable Transformation of a Fuzzy variable x to[ ζ [ 1,] 1] for different α-cuts. The transformation ϕ( ) maps the interval, x (α) [ 1, 1]. x (α) l h
8 Transformation of a fuzzy variable Suppose the transformation ϕ( ) maps the interval [ x (α) i l, x (α) i h [ 1, 1]. Here x (α) i l and x (α) i h denotes the lower and upper bound of the Fuzzy variable x i for a given α-cut. Denote ζ i as the variable bounded between [ 1, 1]. [ Considering the variable x (α) i lies between x (α) i l, x (α) i h ], the linear transformation ϕ( ) can be identified as ( ) x (α) ζ i = ϕ(x (α) i x (α) i l 1 i ) = 2 ( ) x (α) i h x (α) 2 i l (3) The inverse transformation of (3) can be obtained as x (α) i = ( x (α) i h x (α) ) i l ζ i + 2 ( ) x (α) i h + x (α) i l 2 (4)
9 Transformation of a fuzzy variable Substituting x (α) i in Eq. (1) for all α and i, we can can formally express the uncertainty propagation problem in terms of the vector valued variable ζ = {ζ i } i=1,,m [ 1, 1] m as y (α) = f (α) (ζ) R n (5) Now we utilise the orthogonal property of Legendre polynomials in [ 1, 1] for this uncertainty propagation problem 1. 1 Adhikari, S. and Khodaparast, H. H., A spectral approach for fuzzy uncertainty propagation in finite element analysis, Fuzzy Sets and Systems, 243[1] (2014), pp
10 Univariate Legendre polynomials In [ 1, 1], Legendre polynomials are orthogonal with respect to the L 2 inner product norm 1 1 L j (ζ)l k (ζ)dζ = 2 2k + 1 δ jk (6) δ jk denotes the Kronecker delta (equal to 1 if j = k and to 0 otherwise) and L j (ζ) is the jth order Legendre polynomial. Each Legendre polynomial L j (ζ) is an kth-degree polynomial and can be expressed using Rodrigues formula as L k (ζ) = 1 d k { 2 k k! dζ k (1 ζ 2 ) k} ; k = 0, 1, 2, (7)
11 Univariate Legendre polynomials Table : One dimensional Legendre Polynomials upto order 10 k Legendre Polynomials of order k: L k (ζ); ζ [ 1, 1] ζ ( ζ 2 1 ) ( ζ 3 3ζ ) ( ζ 4 15ζ ) ( ζ 5 70ζ ζ ) ( ζ 6 315ζ ζ 2 5 ) ( ζ 7 693ζ ζ 3 35ζ ) ( ζ ζ ζ ζ ) ( 12155ζ ζ ζ ζ ζ ) ( 46189ζ ζ ζ ζ ζ 2 63 )
12 Multivariate Legendre polynomials Multivariate Legendre polynomials can be defined as products of univariate Legendre polynomials, similar to that of Hermite polynomials. This can be constructed by considering products of different orders. We consider following sets of integers i = {i 1, i 2, i p }, i k 1 (8) β = {β 1,β 2, β p }, β k 0 (9) The multivariate Legendre polynomial associated with sequence (i,β) as the products of univariate Legendre polynomials p L i,β (ζ) = L βk (ζ k ) (10) The multivariate Legendre polynomials also satisfy the orthogonality condition with respect to the L 2 inner product norm in the respective dimension. k=1
13 Projection approach We follow the idea of homogeneous chaos proposed by Wienner involving Hermite polynomials. If a function f(ζ) is square integrable, it can be expanded in Homogeneous Chaos as f(ζ) = ŷ i0 h 0 + ŷ i1 Γ 1 (ζ i1 ) + i 1 i 1 =1 ŷ i1,i 2 Γ 2 (ζ i1,ζ i2 )+ i 1 i 2 i 1 =1 i 2 =1 i 1 =1 i 2 =1 i 3 =1 + i 1 i 2 i 3 i 1 =1 i 2 =1 i 3 =1 i 4 =1 ŷ i1 i 2 i 3 Γ 3 (ζ i1,ζ i2,ζ i3 ) ŷ i1 i 2 i 3 i 4 Γ 4 (ζ i1,ζ i2,ζ i3,ζ i4 )+..., (11) Γ p (ζ i1,ζ i2, ζ im ) is m-dimensional homogeneous chaos of order p (obtained by products of one-dimensional Legendre polynomials).
14 Projection approach We can concisely write P 1 f(ζ) = y j Ψ j (ζ) (12) j=0 where the constant y j and functions Ψ j ( ) are effectively constants ŷ k and functions Γ k ( ) for corresponding indices. Equation (12) can be viewed as the projection in the basis functions Ψ j (ζ) with corresponding coordinates y j. The number of terms P in Eq. (12) depends on the number of variables m and maximum order of polynomials p as P = p j=0 (m+j 1)! j!(m 1)! ( ) m+p = p (13)
15 Univariate Legendre polynomials Table : Two dimensional Legendre polynomial based homogeneous chaos basis up to 4th order j p Construction of Ψ j Ψ j 0 p = 0 L p = 1 L 1 (ζ 1 ) ζ 1 2 L 1 (ζ 2 ) ζ 2 3 L 2 (ζ 1 ) 3/2ζ 2 1 1/2 4 p = 2 L 1 (ζ 1 )L 1 (ζ 2 ) ζ 1 ζ 2 5 L 2 (ζ 2 ) 3/2ζ 2 2 1/2 6 L 3 (ζ 1 ) 5/2ζ 3 1 3/2ζ1 ) 7 p = 3 L 2 (ζ 1 )L 1 (ζ 2 ) (3/2ζ 2 1 1/2 ζ ( ) 2 8 L 1 (ζ 1 )L 2 (ζ 2 ) ζ 1 3/2ζ 2 2 1/2 9 L 3 (ζ 2 ) 5/2ζ 3 2 3/2ζ2 10 L 4 (ζ 1 ) 35 ζ ζ /8 ) 11 L 3 (ζ 1 )L 1 (ζ 2 ) (5/2ζ 3 1 3/2ζ1 ζ 2 ) 12 p = 4 L 2 (ζ 1 )L 2 (ζ 2 ) (3/2ζ 2 1 1/2 )(3/2ζ 2 2 1/2 ( ) 13 L 1 (ζ 1 )L 3 (ζ 2 ) ζ 1 5/2ζ 3 2 3/2ζ2 14 L 4 (ζ 2 ) 35 8 ζ ζ /8
16 Least square method The unknown coefficients are obtains using a least-square method with respect to the inner product norm in [ 1, 1] m as Here the volume Using this we have <, >= ( )( )dζ 1 dζ 2 dζ m (14) V m }{{} m fold V m = 2 m (15) [ P 1 Ψj (ζ), f(ζ) f(ζ) = ]Ψ Ψj (ζ),ψ j (ζ) j (ζ) (16) j=0 Here f(ζ) is an approximation to the original function f(ζ) for polynomial order upto p.
17 An example with two variables We consider the Camelback function ( 3 x 1 3; 2 x 2 2) f 1 (x) = (4 2.1x x 4 1 /3)x x 1x 2 +(4x 2 2 4)x 2 2 (17) Transform the variables in [ 1, 1]: Omitting the notation α for convenience, we have x 1 = 3ζ 1 and x 2 = 2ζ 2 (18) The function in the transformed variables ( f 1 (ζ) = ) ( ) 10 ζ ζ 1 ζ ζ 1 ζ ζ 2 ζ 2 (19) The values of y j can be obtained as y 1 = , y 4 = , y 5 = 6, y 6 = , y 11 = , y 15 = (20) with all other values 0.
18 Original and fitted function f(ζ 1,ζ 2 ) Approximate f(ζ 1,ζ 2 ) ζ 2 axis ζ 1 axis ζ 2 axis ζ 1 axis (a) The exact function (b) Fitted function using Legendre polynomials
19 Fuzzy finite element (c) SubdomainsD i with Fuzzy parameter a i for i = 1, 2,, M (d) Finite element mesh We define the discretised fuzzy variables a i for all the subdomains as a i = a(r); r D i, i = 1, 2,, M (21)
20 Fuzzy FE formulation A particular subdomain D i is expected to contain several finite elements. The element stiffness matrix can be obtained following the finite element approach as K e = a(r)b (e) T (r)b (e) (r)dr = a i B D e De (e) T (r)b (e) (r)dr; r D i, (22) where B (e) (r) is a deterministic matrix related to the shape function used to interpolate the solution within the element e. Suppose a (α) i l and a (α) i h denotes the lower and upper bound of the Fuzzy variable a i for a given α-cut. We transform a i into the standard variable ζ i [ 1, 1] for every α-cut as ( (α) a a (α) i i = h x (α) ) ( ) a (α) i l i h + a (α) i l ζ i + (23) 2 2
21 Fuzzy FE formulation Substituting this into Eq. (22) we have K (α) e = K (α) e 0 +ζ i K (α) e i (24) where the crisp and fuzzy parts of the element stiffness matrix is given by ( ) a (α) K (α) i h + a (α) i l e 0 = B (e) T (r)b (e) (r)dr (25) 2 D ( ) e a (α) and K (α) i h a (α) i l T e i = (r)b (e) (r)dr; r D i, (26) 2 D e B (e) The global stiffness matrix can be obtained by the usual finite element assembly procedure and taking account of the domains of the discretised fuzzy variables. Similar approach as be used for the mass and damping matrices also.
22 Fuzzy FE formulation The finite element equation for a given α-cut can be expressed as M (α) Ü (α) (t)+c (α) U(α) (t)+k (α) U (α) (t) = F(t) (27) In the preceding equation F R n is the global forcing vector and U (α) (t) is the dynamic response for a given α-cut. The global stiffness, mass and damping matrices can be expressed as K (α) = K (α) m K 0 + i=1 m M +ζ i K (α) i, M (α) = M (α) 0 + +ζ i M (α) i, C (α) = C (α) i=1 m C 0 + i=1 +ζ i C (α) i (28) Considering the stiffness matrix for example, the crisp part can be obtained as K (α) 0 = e K (α) e 0 (29) where e denotes the summation over all the elements.
23 Fuzzy FE formulation The fuzzy parts of the stiffness matrix related to the variable ζ j can be given by = e i ; i = 1, 2,, m K (30) K (α) i e:r D i K (α) where e:r D i denotes the summation over those elements for which the domain belongs to D i. An expression of the stiffness matrix similar to Eq. (28) can alternatively obtained directly for complex systems where different subsystems may contain different fuzzy variables. In that case K (α) i would be block matrices, influenced by only the transformed fuzzy variable ζ i within a particular subsystem.
24 Frequency domain analysis Considering that the initial conditions are zero for all the coordinates, taking the Fourier transformation of Eq. (27) we have [ ω 2 M (α) +iωc (α) + K (α)] u (α) (ω) = f(ω) (31) Here ω is the frequency parameter and f(ω) and u (α) (ω) are the Fourier transforms of F(t) and U (α) (t) respectively. Equation (31) can be conveniently rewritten as where the dynamic stiffness matrix D(ω) (α) u (α) (ω,ζ) = f(ω) (32) M D (α) = D (α) 0 (ω)+ +ζ i D (α) i (ω) (33) Note that u (α) (ω,ζ) is effectivly a function of the frequency ω and well as the variables ζ. i=1
25 Frequency domain analysis The crisp part of the dynamic stiffness matrix is given by and M is the total number of variables so that D (α) 0 = ω 2 M 0 +iωc 0 + K 0 (34) M = m K + m M + m C (35) The matrices D (α) i becomes the deviatoric part of the system matrices for the corresponding values of the index i. The elements of the solution vector u (α) in the fuzzy finite element equation (32) can be viewed as nonlinear functions of the variables ζ j [ 1, 1] for each value of the frequency ω.
26 Galerkin error minimisation We project the solution vector u (α) (ω,ζ) in the basis of Legendre polynomials as P 1 u (α) (ω,ζ) = u (α) j (ω)ψ j (ζ) (36) j=0 The aim to obtain the coefficient vectors u j C N N using a Galerkin type of error minimisation approach. Substituting expansion of u (α) (ω,ζ) in the governing equation (27), the error vector can be obtained as ( M ) P 1 ε (α) (ω) = D (α) i (ω)ζ i (ω)ψ j (ζ) f(ω) (37) i=0 j=1 u (α) j where ζ 0 = 1 is used to simplify the first summation expression.
27 Galerkin error minimisation The coefficients can be obtained by solving D (α) 0,0 D (α) 1,0 (ω) D(α) 0,1 (ω) D(α) 1,1. D (α) P 1,0 (ω) D(α) (ω) D(α).. 0,P 1 (ω) 1,P 1 (ω) (ω) D(α) P 1,1 (ω) D(α) P 1,P 1 (ω) or in a compact notation 0 (ω) f 0 (ω) 1 (ω) f 1 (ω) =.. f P 1 (ω) u (α) u (α) u (α) P 1 (ω) (38) where D (α) (ω) C np np, U (α) (ω),f(ω) C np. D (α) (ω)u (α) (ω) = F(ω) (39) This equation needs to be solved for every α-cut. Once all u (α) j (ω) for j = 0, 1,..., P 1 are obtained, the solution vector can be obtained from (36) for every α-cut.
28 Problem description The generalised eigenvalue problem for a given α-cut can be expressed as K (α) y j = λ (α) j M (α) y (α) j, j = 1, 2, 3, (40) where λ (α) j system. and y (α) j denote the eigenvalues and eigenvectors of the The eigenvalues and eigenvectors can be expressed by spectral decomposition as 2 λ (j) (ξ 1,...ξ M ) = y (j) (ξ 1,...ξ M ) = P λ jk Ψ k (ξ 1,...ξ M ) (41) k=1 P k=1 y (j) k Ψ k(ξ 1,...ξ M ) (42) where λ jk and y (j) k are unknowns and the basis functions considered here are multivariate Legendre s polynomials. 2 Pascual, B., Adhikari, S., Hybrid perturbation-polynomial chaos approaches to the random algebraic eigenvalue problem, Computer Methods in Applied Mechanics and Engineering, [1] (2012), pp
29 Problem description Figure : A cantilever rod subjected to an axial force modelled using two fuzzy variables. L 1 = 1m, L 2 = 1m, crisp values EA 1 = N, EA 2 = N, F = 1kN. The linear (triangular) membership functions for the two variables are shown in the figure.
30 Fuzzy variables At α = 0 we have the maximum variability EA and EA (43) This implies a variability between 10% and 30% for EA 1 and ±20% for EA 2. The maximum and minimum values of the interval for a given α-cut can be described for the two fuzzy variables as EA 1min = EA 1 ( α)and EA 1max = EA 1 ( α) (44) EA 2min = EA 2 ( α)and EA 2max = EA 2 ( α) (45) The stiffness matrix: K = K 0 +ζ 1 K (α) 1 +ζ 2 K (α) 2 (46) Here K 0 R is the stiffness matrix corresponding to the crisp values, K (α) 1 and K (α) 2 are the coefficient matrices corresponding to a given α-cut.
31 Static Results Direct simulation Orthogonal polynomial (1st order) Orthogonal polynomial (2nd order) Orthogonal polynomial (4th order) Direct simulation Orthogonal polynomial (1st order) Orthogonal polynomial (2nd order) Orthogonal polynomial (4th order) α level 0.5 α level Deflection at the mid point (m) x Deflection at the end point (m) x 10 5 (a) Response at the mid point (b) Response at the end point Fuzzy description of the response of the rod at two points computed using different order of Legendre polynomial based homogeneous chaos and direct simulation.
32 Percentage error in the mid point: Static analysis 1st order 2nd order 4th order α-cut min max min max min max
33 Percentage error at the end point: Static analysis 1st order 2nd order 4th order α-cut min max min max min max
34 Frequency response function: mid point (c) α = 0.75 (d) α = 0.5 Variability of the frequency response of the rod at the mid-point computed using Legendre polynomial based homogeneous chaos (4th order) and direct simulation.
35 Frequency response function: mid point (e) α = 0.25 (f) α = 0 Variability of the frequency response of the rod at the mid-point computed using Legendre polynomial based homogeneous chaos (4th order) and direct simulation.
36 Frequency response function: end point (g) α = 0.75 (h) α = 0.5 Variability of the frequency response of the rod at the end-point computed using Legendre polynomial based homogeneous chaos (4th order) and direct simulation.
37 Frequency response function: end point (i) α = 0.25 (j) α = 0 Variability of the frequency response of the rod at the end-point computed using Legendre polynomial based homogeneous chaos (4th order) and direct simulation.
38 Summary and conclusions Uncertainty propagation in complex systems with Fuzzy variables is considered. An orthogonal function expansion approach in conjunction with Galerkin type error minimisation is proposed. The method proposed has three major steps: (a) transformation of a fuzzy variable into a set of interval variables for different α-cuts via the membership function, (b) transformation of interval variables into a standard interval variable between [ 1, 1] for each α-cut, and (c) projection of the response function in the basis of multivariate orthogonal Legendre polynomials in terms of the transformed standard interval variables. The coefficients associated with the basis functions are obtained using a Galerkin type of error minimisation. Depending on the number of basis retained in the series expansion, it is shown that various order of approximation can be obtained.
39 Summary and conclusions A computational method is proposed for fuzzy finite element problems in structural dynamics, where the technique is generalised to vector valued functions with multiple fuzzy variables in the frequency domain. For a given α-cut, the complex dynamic stiffness matrix is expressed as a series involving standard interval variables and coefficient matrices. This representation significantly simplify the problem of response prediction via the proposed multivariate orthogonal Legendre polynomials expansion technique. The coefficient vectors in the polynomial expansion are calculated from the solution of an extended set of linear algebraic equations. A numerical example of axial deformation of a rod with fuzzy axial stiffness is considered to illustrate the proposed methods. The results are compared with direct numerical simulation results.
40 Summary and conclusions The spectral method proposed in this paper enables to propagate fuzzy uncertainty in a mathematically rigorous and general manner similar to what is available for propagation of probabilistic uncertainty. The main computational cost involves the solution of an extended set of linear algebraic equations necessary to obtain the coefficients associated with the polynomial basis functions. Future research is necessary to develop computationally efficient methods to deal with this problem arising in systems with large number of fuzzy variables.
Dynamic response of structures with uncertain properties
Dynamic response of structures with uncertain properties S. Adhikari 1 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Bay Campus, Fabian Way, Swansea, SA1 8EN, UK International
More informationA Vector-Space Approach for Stochastic Finite Element Analysis
A Vector-Space Approach for Stochastic Finite Element Analysis S Adhikari 1 1 Swansea University, UK CST2010: Valencia, Spain Adhikari (Swansea) Vector-Space Approach for SFEM 14-17 September, 2010 1 /
More informationCollocation based high dimensional model representation for stochastic partial differential equations
Collocation based high dimensional model representation for stochastic partial differential equations S Adhikari 1 1 Swansea University, UK ECCM 2010: IV European Conference on Computational Mechanics,
More informationOn the Nature of Random System Matrices in Structural Dynamics
On the Nature of Random System Matrices in Structural Dynamics S. ADHIKARI AND R. S. LANGLEY Cambridge University Engineering Department Cambridge, U.K. Nature of Random System Matrices p.1/20 Outline
More informationIdentification of Damping Using Proper Orthogonal Decomposition
Identification of Damping Using Proper Orthogonal Decomposition M Khalil, S Adhikari and A Sarkar Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk
More informationStochastic structural dynamic analysis with random damping parameters
Stochastic structural dynamic analysis with random damping parameters K. Sepahvand 1, F. Saati Khosroshahi, C. A. Geweth and S. Marburg Chair of Vibroacoustics of Vehicles and Machines Department of Mechanical
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
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 informationProbabilistic Structural Dynamics: Parametric vs. Nonparametric Approach
Probabilistic Structural Dynamics: Parametric vs. Nonparametric Approach S Adhikari School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk URL: http://engweb.swan.ac.uk/
More informationDr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum
STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationDesign of Structures for Earthquake Resistance
NATIONAL TECHNICAL UNIVERSITY OF ATHENS Design of Structures for Earthquake Resistance Basic principles Ioannis N. Psycharis Lecture 3 MDOF systems Equation of motion M u + C u + K u = M r x g(t) where:
More informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationRandom Eigenvalue Problems in Structural Dynamics: An Experimental Investigation
Random Eigenvalue Problems in Structural Dynamics: An Experimental Investigation S. Adhikari, A. Srikantha Phani and D. A. Pape School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk
More informationKinematic Analysis and Inverse Dynamics-based Control of Nondeterministic Multibody Systems
Kinematic Analysis and Inverse Dynamics-based Control of Nondeterministic Multibody Systems Item Type text; Electronic Thesis Authors Sabet, Sahand Publisher The University of Arizona. Rights Copyright
More informationUncertainty Quantification in MEMS
Uncertainty Quantification in MEMS N. Agarwal and N. R. Aluru Department of Mechanical Science and Engineering for Advanced Science and Technology Introduction Capacitive RF MEMS switch Comb drive Various
More informationChapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14
Table of Contents Chapter 1: Research Objectives and Literature Review..1 1.1 Introduction...1 1.2 Literature Review......3 1.2.1 Describing Vibration......3 1.2.2 Vibration Isolation.....6 1.2.2.1 Overview.
More informationDynamic Response of Structures With Frequency Dependent Damping
Dynamic Response of Structures With Frequency Dependent Damping Blanca Pascual & S Adhikari School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk URL: http://engweb.swan.ac.uk/
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 informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 13-14 December, 2017 1 / 30 Forewords
More informationAn Efficient Computational Solution Scheme of the Random Eigenvalue Problems
50th AIAA SDM Conference, 4-7 May 2009 An Efficient Computational Solution Scheme of the Random Eigenvalue Problems Rajib Chowdhury & Sondipon Adhikari School of Engineering Swansea University Swansea,
More informationEstimating functional uncertainty using polynomial chaos and adjoint equations
0. Estimating functional uncertainty using polynomial chaos and adjoint equations February 24, 2011 1 Florida State University, Tallahassee, Florida, Usa 2 Moscow Institute of Physics and Technology, Moscow,
More informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationEfficient Reduced Order Modeling of Low- to Mid-Frequency Vibration and Power Flow in Complex Structures
Efficient Reduced Order Modeling of Low- to Mid-Frequency Vibration and Power Flow in Complex Structures Yung-Chang Tan Graduate Student Research Assistant Matthew P. Castanier Assistant Research Scientist
More informationReliability-Based Microstructural Topology Design with Respect to Vibro-Acoustic Criteria
11 th World Congress on Structural and Multidisciplinary Optimisation 07 th -12 th, June 2015, Sydney Australia Reliability-Based Microstructural Topology Design with Respect to Vibro-Acoustic Criteria
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 informationAn Empirical Chaos Expansion Method for Uncertainty Quantification
An Empirical Chaos Expansion Method for Uncertainty Quantification Melvin Leok and Gautam Wilkins Abstract. Uncertainty quantification seeks to provide a quantitative means to understand complex systems
More informationAlgorithms for Uncertainty Quantification
Algorithms for Uncertainty Quantification Lecture 9: Sensitivity Analysis ST 2018 Tobias Neckel Scientific Computing in Computer Science TUM Repetition of Previous Lecture Sparse grids in Uncertainty Quantification
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
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 informationPolynomial chaos expansions for sensitivity analysis
c DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for sensitivity analysis B. Sudret Chair of Risk, Safety & Uncertainty
More informationResearch Article Multiresolution Analysis for Stochastic Finite Element Problems with Wavelet-Based Karhunen-Loève Expansion
Mathematical Problems in Engineering Volume 2012, Article ID 215109, 15 pages doi:10.1155/2012/215109 Research Article Multiresolution Analysis for Stochastic Finite Element Problems with Wavelet-Based
More informationOn the Application of Intervening Variables for Stochastic Finite Element Analysis
On the Application of Intervening Variables for Stochastic Finite Element Analysis M.A. Valdebenito a,, A.A. Labarca a, H.A. Jensen a a Universidad Tecnica Federico Santa Maria, Dept. de Obras Civiles,
More informationUncertainty quantification for sparse solutions of random PDEs
Uncertainty quantification for sparse solutions of random PDEs L. Mathelin 1 K.A. Gallivan 2 1 LIMSI - CNRS Orsay, France 2 Mathematics Dpt., Florida State University Tallahassee, FL, USA SIAM 10 July
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationUNCERTAINTY ASSESSMENT USING STOCHASTIC REDUCED BASIS METHOD FOR FLOW IN POROUS MEDIA
UNCERTAINTY ASSESSMENT USING STOCHASTIC REDUCED BASIS METHOD FOR FLOW IN POROUS MEDIA A REPORT SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT
More informationPolynomial Chaos and Karhunen-Loeve Expansion
Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the
More informationRandom Matrix Eigenvalue Problems in Probabilistic Structural Mechanics
Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationMatrix Iteration. Giacomo Boffi.
http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 12, 2016 Outline Second -Ritz Method Dynamic analysis of MDOF
More informationDamping Modelling and Identification Using Generalized Proportional Damping
Damping Modelling and Identification Using Generalized Proportional Damping S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk Generalized
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationRandom Eigenvalue Problems Revisited
Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationOPTIMUM PRE-STRESS DESIGN FOR FREQUENCY REQUIREMENT OF TENSEGRITY STRUCTURES
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm OPTIMUM PRE-STRESS DESIGN FOR FREQUENCY REQUIREMENT OF TENSEGRITY STRUCTURES Seif Dalil
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationVirtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing example
AMAS Workshop on Smart Materials and Structures SMART 03 (pp.313 324) Jadwisin, September 2-5, 2003 Virtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing
More informationUncertainty Propagation and Global Sensitivity Analysis in Hybrid Simulation using Polynomial Chaos Expansion
Uncertainty Propagation and Global Sensitivity Analysis in Hybrid Simulation using Polynomial Chaos Expansion EU-US-Asia workshop on hybrid testing Ispra, 5-6 October 2015 G. Abbiati, S. Marelli, O.S.
More informationOutline. Structural Matrices. Giacomo Boffi. Introductory Remarks. Structural Matrices. Evaluation of Structural Matrices
Outline in MDOF Systems Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano May 8, 014 Additional Today we will study the properties of structural matrices, that is the operators that
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationParadigms of Probabilistic Modelling
Paradigms of Probabilistic Modelling Hermann G. Matthies Brunswick, Germany wire@tu-bs.de http://www.wire.tu-bs.de abstract RV-measure.tex,v 4.5 2017/07/06 01:56:46 hgm Exp Overview 2 1. Motivation challenges
More informationTWO-STAGE ISOLATION FOR HARMONIC BASE EXCITATION Revision A. By Tom Irvine February 25, 2008
TWO-STAGE ISOLATION FOR HARMONIC BASE EXCITATION Revision A By Tom Irvine Email: tomirvine@aol.com February 5, 008 Introduction Consider a base plate mass m and an avionics mass m modeled as two-degree-of-freedom.
More informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
More informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More informationResearch Article A Pseudospectral Approach for Kirchhoff Plate Bending Problems with Uncertainties
Mathematical Problems in Engineering Volume 22, Article ID 7565, 4 pages doi:.55/22/7565 Research Article A Pseudospectral Approach for Kirchhoff Plate Bending Problems with Uncertainties Ling Guo, 2 and
More informationDAMPING MODELLING AND IDENTIFICATION USING GENERALIZED PROPORTIONAL DAMPING
DAMPING MODELLING AND IDENTIFICATION USING GENERALIZED PROPORTIONAL DAMPING S. Adhikari Department of Aerospace Engineering, University of Bristol, Queens Building, University Walk, Bristol BS8 1TR (U.K.)
More informationChapter 4 Analysis of a cantilever
Chapter 4 Analysis of a cantilever Before a complex structure is studied performing a seismic analysis, the behaviour of simpler ones should be fully understood. To achieve this knowledge we will start
More informationTime domain analysis of viscoelastic models
Time domain analysis of viscoelastic models L. Rouleau 1, J.-F. Deü 1 1 LMSSC, Conservatoire national des arts et métiers, 2D6R10, 292 rue Saint-Martin, 75141 Paris cedex 03, France e-mail: lucie.rouleau@lecnam.net
More informationStructural Matrices in MDOF Systems
in MDOF Systems http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 9, 2016 Outline Additional Static Condensation
More informationTime integration and the Trefftz Method, Part II Second-order and hyperbolic problems
ime integration and the refftz Method, Part II Second-order and hyperbolic problems JA eixeira de Freitas Departamento de Engenharia Civil e Arquitectura, Instituto Superior écnico echnical University
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationRandom Vibrations & Failure Analysis Sayan Gupta Indian Institute of Technology Madras
Random Vibrations & Failure Analysis Sayan Gupta Indian Institute of Technology Madras Lecture 1: Introduction Course Objectives: The focus of this course is on gaining understanding on how to make an
More informationLecture 02 Linear Algebra Basics
Introduction to Computational Data Analysis CX4240, 2019 Spring Lecture 02 Linear Algebra Basics Chao Zhang College of Computing Georgia Tech These slides are based on slides from Le Song and Andres Mendez-Vazquez.
More informationVector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.
Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +
More informationCOMPLEX MODES IN LINEAR STOCHASTIC SYSTEMS
Proceedings of VETOMAC-I October 25-27, 2, Bangalore, INDIA COMPLEX MODES IN LINEAR STOCHASTIC SYSTEMS S. Adhikari Cambridge University Engineering Department Trumpington Street Cambridge CB2 1PZ (U.K.)
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationSHOCK RESPONSE OF MULTI-DEGREE-OF-FREEDOM SYSTEMS Revision F By Tom Irvine May 24, 2010
SHOCK RESPONSE OF MULTI-DEGREE-OF-FREEDOM SYSTEMS Revision F By Tom Irvine Email: tomirvine@aol.com May 4, 010 Introduction The primary purpose of this tutorial is to present the Modal Transient method
More informationReduction of Random Variables in Structural Reliability Analysis
Reduction of Random Variables in Structural Reliability Analysis S. Adhikari and R. S. Langley Department of Engineering University of Cambridge Trumpington Street Cambridge CB2 1PZ (U.K.) February 21,
More informationPARALLEL COMPUTATION OF 3D WAVE PROPAGATION BY SPECTRAL STOCHASTIC FINITE ELEMENT METHOD
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 569 PARALLEL COMPUTATION OF 3D WAVE PROPAGATION BY SPECTRAL STOCHASTIC FINITE ELEMENT METHOD Riki Honda
More informationCOMPARATIVE STUDIES ON SEISMIC INCOHERENT SSI ANALYSIS METHODOLOGIES
Transactions, SMiRT-22 COMPARATIVE STUDIES ON SEISMIC INCOHERENT SSI ANALYSIS METHODOLOGIES Dan M. Ghiocel 1 1 Ghiocel Predictive Technologies, Inc., Rochester, New Yor, USA (dan.ghiocel@ghiocel-tech.com)
More informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationStochastic Dimension Reduction
Stochastic Dimension Reduction Roger Ghanem University of Southern California Los Angeles, CA, USA Computational and Theoretical Challenges in Interdisciplinary Predictive Modeling Over Random Fields 12th
More informationBenjamin L. Pence 1, Hosam K. Fathy 2, and Jeffrey L. Stein 3
2010 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 30-July 02, 2010 WeC17.1 Benjamin L. Pence 1, Hosam K. Fathy 2, and Jeffrey L. Stein 3 (1) Graduate Student, (2) Assistant
More informationClassical and Quantum Complexity of the Sturm-Liouville Eigenvalue Problem
Classical and Quantum Complexity of the Sturm-Liouville Eigenvalue Problem A. Papageorgiou 1 and H. Woźniakowski 1, Department of Computer Science, Columbia University, New York, USA Institute of Applied
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationPolynomial chaos expansions for structural reliability analysis
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for structural reliability analysis B. Sudret & S. Marelli Incl.
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationNotation. General. Notation Description See. Sets, Functions, and Spaces. a b & a b The minimum and the maximum of a and b
Notation General Notation Description See a b & a b The minimum and the maximum of a and b a + & a f S u The non-negative part, a 0, and non-positive part, (a 0) of a R The restriction of the function
More informationMechanics of Irregular Honeycomb Structures
Mechanics of Irregular Honeycomb Structures S. Adhikari 1, T. Mukhopadhyay 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Singleton Park, Swansea SA2 8PP, UK Sixth International
More informationAdvanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One
Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to
More informationCHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES
CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may
More informationReduction in number of dofs
Reduction in number of dofs Reduction in the number of dof to represent a structure reduces the size of matrices and, hence, computational cost. Because a subset of the original dof represent the whole
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationInterval model updating: method and application
Interval model updating: method and application H. Haddad Khodaparast, J.E. Mottershead, K.J. Badcock University of Liverpool, School of Engineering, Harrison-Hughes Building, The Quadrangle, L69 3GH,
More informationStructural Design under Fuzzy Randomness
NSF workshop on Reliable Engineering Computing, September 15-17, 2004 Center for Reliable Engineering Computing, Georgia Tech Savannah, GA Structural Design under Fuzzy Randomness Michael Beer 1) Martin
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationDinesh Kumar, Mehrdad Raisee and Chris Lacor
Dinesh Kumar, Mehrdad Raisee and Chris Lacor Fluid Mechanics and Thermodynamics Research Group Vrije Universiteit Brussel, BELGIUM dkumar@vub.ac.be; m_raisee@yahoo.com; chris.lacor@vub.ac.be October, 2014
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationSparse polynomial chaos expansions in engineering applications
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Sparse polynomial chaos expansions in engineering applications B. Sudret G. Blatman (EDF R&D,
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More information