Uncertainty is unavoidable in engineering system. Probabilistic approach is the traditional approach
|
|
- Baldwin Summers
- 6 years ago
- Views:
Transcription
1 Interval Finite Elements as a Basis for Generalized Models of Uncertainty in Engineering Analysis Rafi L. Muhanna Georgia Institute of Technology USA Collaborative Research Center 528 German National Science Foundation (DFG) TU Dresden, Germany, June 19, 2007
2 Outline Introduction Interval Arithmetic Interval Finite Elements Element-By By-Element Examples Conclusions
3 Outline Introduction Interval Arithmetic Interval Finite Elements Element-By By-Element Examples Conclusions
4 Introduction- Uncertainty Uncertainty is unavoidable in engineering system structural mechanics entails uncertainties in material, geometry and load parameters (aleatory-epistemic) Probabilistic approach is the traditional approach requires sufficient information to validate the probabilistic model criticism of the credibility of probabilistic approach when data is insufficient (Elishakoff, 1995; Ferson and Ginzburg, 1996; Möller and Beer, 2007)
5 Introduction- Interval Approach Nonprobabilistic approach for uncertainty modeling when only range information (tolerance) is available t = t0 ± δ Represents an uncertain quantity by giving a range of possible values t = [ t δ, t + δ ] 0 0 How to define bounds on the possible ranges of uncertainty? experimental data, measurements, statistical analysis, expert knowledge
6 Introduction- Why Interval? Simple and elegant Conforms to practical tolerance concept Describes the uncertainty that can not be appropriately modeled by probabilistic approach Computational basis for other uncertainty approaches (e.g., fuzzy set, random set, imprecise probability) Provides guaranteed enclosures
7 Outline Introduction Interval Arithmetic Interval Finite Elements Element-By By-Element Examples Conclusions
8 Interval arithmetic Background Archimedes ( B.C.) Α circle of radius one has an area equal to π
9 Interval arithmetic Background Archimedes ( B.C.) Α circle of radius one has an area equal to π 2 < π < 4 r = < π < r=1 π = [ , ]
10 Interval arithmetic Background Modern interval arithmetic Physical constants or measurements g [9.8045, ] Representation of numbers 1/ π / 3 [0.3333, ] 2 [1.4142, ] π [3.1415, ] Rounding errors 1 / / [8.1004, ] R. E. Moore, E. Hansen, A. Neumaier, G. Alefed, J. Herzberger
11 Interval arithmetic Interval number represents a range of possible values within a closed set x [ x, x]: = { x R x x x}
12 Interval Operations Let x = [a, b] and y = [c, d] be two interval numbers 1. Addition x + y = [a, b] + [c, d] = [a + c, b + d] 2. Subtraction x y = [a, b] [c, d] = [a d, b c] 3. Multiplication xy = [min(ac,ad,bc,bd), max(ac,ad,bc,bd)] 4. Division 1 /x= [1/b, 1/a]
13 Properties of Interval Arithmetic Let x, y and z be interval numbers 1. Commutative Law x + y = y + x xy = yx 2. Associative Law x + (y + z) = (x + y) + z x(yz) = (xy)z 3. Distributive Law does not always hold, but x(y + z) ) xy + xz
14 Sharp Results Overestimation The DEPENDENCY problem arises when one or several variables occur more than once in an interval expression f (x) = x x, x= [1, 2] f (x) = [1 2, 2 1] = [ 1, 1] 0 f (x, y) = { f (x, y) = x y x x, y y} f (x) = x (1 1) f (x) = 0 f (x) = { f (x) = x x x x} }
15 Sharp Results Overestimation If a, b and c are interval numbers, then: a (b ± c) ab ± ac If we set a = [ 2,[ 2]; b = [1, 2]; c = [ 2,[ 1], we get a (b + c) = [ 2, 2]([1, 2] + [ 2, 1]) = [ 2, 2] [ 1, 3] = [ 6, 6] However, ab + ac = [ 2, 2 2][1, 2]+[ + [ 2, 2 2][ 2, 21] = [ 4,[ 44]+[ 4] + [ 4, 44] 4] = [ 8, 88] 8]
16 Sharp Results Overestimation Interval Vectors and Matrices An interval matrix is such matrix that contains all real matrices whose elements are obtained from all possible values between the lower and upper bounds of its interval components m n A = { A R A A for i = 1,..., m; j = 1,..., n} ij ij
17 Sharp Results Sharp Results Overestimation Overestimation Let a, b, c and d be independent variables, each with interval [1, 3] b ] 2 2 [ 2] 2 [ 1 1 B, d c b a B = = = 2] 2 [ 2] 2 [ 2] 2 [ 2] 2 [, ,,,, A A b b b b b b b b B b b b b B = = = ] [ ] [ ] [ ] [,, phys phys A A = = = , , * * phys B phys A A B b
18 Outline Introduction Interval Arithmetic Interval Finite Elements Element-By By-Element Examples Conclusions
19 Finite Elements Finite Element Method (FEM) is a numerical method that t provides approximate solutions to differential equations (ODE and PDE)
20 Finite Elements- Uncertainty& Errors Mathematical model (validation) Discretization of the mathematical model into a computational framework Parameter uncertainty (loading, material properties) Rounding errors
21 Interval Finite Elements Uncertain Data Geometry Materials Loads Interval Stiffness Matrix T K = B CB dv Interval Load Vector Element Level K U = F
22 Interval Finite Elements K U = F K = B T C B dv = Interval element stiffness matrix B = Interval strain-displacement matrix C = Interval elasticity matrix F = [F 1,... F i,... F n ] = Interval element tload vector (traction) ti F = N t da i i N i = Shape function corresponding to the i-th DOF t = Surface traction
23 Interval Finite Elements (IFEM) Follows conventional FEM Loads, geometry and material property are expressed as interval quantities i System response is a function of the interval variables and therefore varies in an interval Computing the exact response range is proven NP-hard The problem is to estimate the bounds on the unknown exact response range based on the bounds of the parameters
24 IFEM- Inner-Bound Methods Combinatorial i method (Muhanna and Mullen 1995, Rao and Berke 1997) Sensitivity analysis method (Pownuk 2004) Perturbation (Mc William 2000) Monte Carlo sampling method Need for alternative methods that achieve Rigorousness guaranteed enclosure Accuracy sharp enclosure Scalability large scale problem Efficiency
25 IFEM- Enclosure Linear static finite i element Muhanna, Mullen, 1995, 1999, 2001,and Zhang 2004 Popova 2003, and Kramer 2004 Neumaier and Pownuk 2004 Corliss, Foley, and Kearfott 2004 Dynamic Dessombz, 2000 Free vibration-buckling Modares, Mullen 2004, and Billini and Muhanna 2005
26 Interval Finite Elements Interval Linear System of Equations A x = b x 1 2 [ 1,0] 1.2 = [ 1,0] x Then A A iff A : 2 α = with, [0,1] 2 αβ β
27 Interval Finite Elements x 1 H [0.3, 0.6] A b= S ( A, b ) = [ 0.6, 0.3] S (A, b) P 4 P 3 S sym (A, b) [0.4, 0.6] S sym ( Ab, ) = [ 0.6, 0.4] S sym (A, b) P 1 = (0.3, 0.6) 0.5 P 2 = (0.6, 0.6) P 1 P 2 P H 3 = (0.6, 0.3) but A b S (A,, b) ) [ 0.8, 0.2] P 4 = (0.4, 0.4) 1 [0.2, 0.8] A b= A 1 b 1.0 H A b x 2
28 Outline Introduction Interval Arithmetic Interval Finite Elements Element-By By-Element Examples Conclusions
29 Naïve interval FEA E1, A1, L1 2 p 1 E2, A2, L2 3 p 2 EA / L = k = [0.95, 1.05], EA p / L = k = [1.9, 2.1], = 0.5, p = k1+ k2 k2 u1 p1 [2.85, 3.15] [ 2.1, 1.9] u1 0.5 = = k k u p [ 2.1, 1.9] [1.9, 2.1] u exact solution: u 2 = [1.429, 1.579], u 3 = [1.905, 2.105] naïve solution: u 2 = [-0.052, 3.052], u 3 = [0.098, 3.902] interval arithmetic assumes that all coefficients are independent uncertainty in the response is severely overestimated (1900%)
30 Element-By-Element Element-By-Element (EBE) technique elements are detached no element coupling structure stiffness matrix is block-diagonal l (k 1,, k Ne ) the size of the system is increased u =(u 1,, u Ne ) T need to impose necessary constraints for compatibility and equilibrium u 1 u 3 u 4 u 2 E 1, A 1, L 1 E 2, A 2, L 2 Element-By-Element model
31 Element-By-Element Suppose the modulus of elasticity is interval: E = Eˆ(1 + δ ) δ: zero-midpoint interval The element stiffness matrix can be split into two parts, k = kˆ( I + d) = kˆ+ kˆd kˆ : deterministic part, element stiffness matrix evalued using Eˆ, kd ˆ : interval part d: interval diagonal matrix, diag( δ,..., δ).
32 Element-By-Element : Element stiffness matrix: k = ki ˆ( + d) : Structure stiffness matrix: K = Kˆ( I + D) = Kˆ + KˆD or K k 1 ˆk 1 d1 = = I + ˆ Ne k k N dne e
33 Constraints Impose necessary constraints for compatibility and equilibrium Penalty method Lagrange multiplier method u 1 u 2 u 3 u 4 E 1, A 1, L 1 E 2, A 2, L 2 Element-By-Element model
34 Constraints penalty method Constraint t conditions: cu = 0 Using the penalty method: ( K + Q ) u = p T Q: penalty matrix, Q= c ηc η : diagonal matrix of penalty number η i Requires a careful choice of the penatly number A spring of large stiffness is added to force node 2 and node 3 to have the same displacement. E 1, A 1, L 1 E 2, A 2, L 2 p
35 Constraints Lagrange multiplier Constraint conditions: c u = 0 Using the Lagrange multiplier method: T K c u p = c 0 λ 0 λ : Lagrange multiplier vector, introdued as new unknowns
36 Load in EBE Nodal lload p b = 1 p ( p,..., p ) b T where p i = N φ( ( x ) dx Suppose the surface traction φ( x) is described by j an interval function: φ( x) = a x. p b can be rewritten as p b m j= 0 N j = W F W: deterministic matrix F: interval vector containing the interval coefficients of the surface tractiton e T
37 Fixed point iteration For the interval equation Ax = b, preconditioning: RAx = Rb, R is the preconditioning matrix transform it into g (x * ) = x * : R b RA x 0 + (I RA) x * = x *, x = x * +x 0 Theorem (Rump, 1990): for some interval vector x *, if g (x * ) int (x * ) then A H b x * +x 0 Iteration algorithm: *( l+ 1) *( l) iterate: x = z + G( ε x ) where z = Rb RAx, G = I RA, R = Aˆ, Ax ˆ = bˆ No dependency handling 1 0 0
38 Fixed point iteration, Interval FEA calls for a modified d method which h exploits the special form of the structure equations, ( K + Q ) u= p with K= Kˆ + KˆD R = Kˆ + Q 1 Choose ( ), construct iterations: *( l+ 1) u p K 0 K ˆ *( l ) = Rp u0 RKD ( u0+ε u ) if *( l ) = R R( + Q) u + ( I R( + Q) )( ε u ) ( l ) = Rp u0 RKM Δ ˆ *( l+ 1) *( l) *( l+ 1) ( l) int( ), then + u0 = R RK u u u= u p M Δ: interval vector, Δ=( δ,..., δ ) 1 N e T The interval variables δ,..., δ appear only once in each iteration. i 1 N e Δ
39 Convergence of fixed point The algorithm converges if and only if ρ ( G ) < 1 smaller ρ( G ) less iterations required, and less overestimation in results To minimize ρ( G ): choose R= Aˆ so that G = I RA has a small spectral radius 1 : 1 reduce the overestimation in ˆ 1 ˆ ˆ ˆ G= I R A=I ( K + Q ) ( K + Q + K D ) = RK D G
40 Stress calculation Conventional method: σ = CBu e,(severe overestimation) C: elasticity matrix, B: strain-displacement matrix Present method: E = (1 + δ) E ˆ, C = (1 + δ) C ˆ σ = CBLu ˆ ( l ) =CBL( Rp RCM Δ) = + δ ˆ p ˆ ˆM ( l ) (1 )( CBLR CBLRK Δ) L: Boolean matrix, L u= ue
41 Element nodal force calculation Conventional method: f = ( ku p ), (severe overestimation) T e e e Present method: ( Te) 1( k1( ue) 1 ( pe) 1) in the EBE model, T ( Ku pb ) = ( e) N ( ( ) ( ) ) T k e N u e e N p e e Ne from ( K + Q) u= p + p T( Ku p ) = T( p Qu) Calculate T( p for all elements. c c b b c Qu) to obtain the element nodal forces
42 Outline Introduction Interval Arithmetic Interval Finite Elements Element-By By-Element Examples Conclusions
43 Numerical example Examine the rigorousness, accuracy, scalability, and efficiency of the present method Comparison with the alternative methods the combinatorial method, sensitivity analysis method, and Monte Carlo sampling method these alternative methods give inner estimation x i x x i xo x: exact solution, x: inner bound, x : outer bound o
44 Examples Load Uncertainty Four-bay forty-story frame
45 Examples Load Uncertainty Four-bay forty-story frame Loading A Loading B Loading C Loading D
46 Examples Load Uncertainty Four-bay forty-story frame Total number of floor load patterns = If one were able to calculate kn/m (1.2 kip/ft) 10,000 patterns / s there has not been sufficient time since the creation of the universe (4-8 ) billion years? to solve all load patterns for this simple structure Mt Material ila36, Beams W24 x 55, Columns W14 x m (48 ft)
47 Examples Load Uncertainty Four-bay forty-story frame Four bay forty floor frame - Interval solutions for shear force and bending moment of first floor columns Elements Nodes Combination solution Total number of required combinations = Interval Axial force (kn) [ , 185.7] [ , 0.0] [ , 0.0] solution Shear force (kn) [-5.1, 0.9] [-5.8, 5.0] [-5.0, 5.0] Moment (kn m) [-10.3, 4.5] [-15.3, 5.4] [-10.6, 9.3] [-17, 15.2] [-8.9, 8.9] [-16, 16]
48 Examples Load Uncertainty Ten-bay truss 1 A = m 2 E = kpa kN 20 kn 20kN 20kN 20kN 20kN 20 kn 20kN 20kN 4 = 40 m F = [-4.28, 28.3] kn m F min = -( )( ) 20= kn F max = ( ) 20 =28.3 kn
49 Examples Load Uncertainty Three-Span Beam kn /m m 4 m 4 m 12 m 11.2kNm 4.8 knm 11.2 knm knm knm knm m m 2 m 2 m m m
50 Truss structure P m P 2 P 1 P 3 4.5m 4.5m 4.5m 4.5m A A A,A A A cross-sectional area 2 1, 2, 3 13, 14, 15: [9.95, 10.05] cm (1% uncertainty) 2 of all other elements: [5.97, 6.03] cm (1% uncertainty) modulus of elasticity of all elements: 200,000 MPa p = [190, 210] kn, p = [95,105] kn p 1 2 = [95,105] kn, p = [85.5,94.5] kn (10% uncertainty) 3 4
51 Truss structure - results Table: results of selected responses Method u 5 (LB) u 5 (UB) N 7 (LB) N 7 (UB) Combinatorial Naïve IFEA δ % % 362% 327% Present IFEA δ 0.19% 011% 0.11% 0.19% 0.15% unit: u 5 (m), N 7 (kn). LB: lower bound; UB: upper bound.
52 Truss structure results u 5 (m) Comb. LB Comb. UB IFEA LB IFEA UB % 1% 2% 3% 4% 5% 6% 7% 8% 9% 10% Uncertainty in cross-sectional area for moderate uncertainty ( 5%), very sharp bounds are obtained for relatively large uncertainty, reasonable bounds are obtained in the case of 10% uncertainty: Comb.: u5 = [ , ], IFEM: u5 = [ , ] (relative difference: 2.59%, 1.80% for LB, UB, respectively)
53 Truss with a large number of interval variables p C p p p p D story bay N e N v p p n@l A B A i E i m@l = [0.995,1.005] A, = [0.995,1.005] E for i = 1,..., N 0 0 e
54 Scalability study vertical displacement at right upper corner (node D): v D Table: displacement at node D = a PL EA 0 0 Sensitivity Analysis Present IFEA Story bay LB * UB * LB UB δ LB δ UB wid/d % 0.10% 2.64% % 0.16% 2.84% % 0.22% 3.02% % 0.29% 3.19% % 0.37% 3.37% 37% % 0.45% 3.56% δ = LB- * * = UB- * * = (LB- * * LB LB / LB, δ LB UB / UB, δ LB LB )/ LB
55 Efficiency study Table: CPU time for the analyses with the present method (unit: seconds) Story bay N v Iteratio n t i t r t t i /t t r /t % 78.4% % 80.5% % 89.1% % 89.7% % 90.1% % 90.0% t i : iteration time, t r : CPU time for matrix inversion, t : total comp. CPU time majority of ftime is spent on matrix inversion i
56 Efficiency study CPU time (sec) Computational time: a comparison of the senstitivity analysis method and the present method Sensitivity Analysis method Present interval FEA Number of interval variables Computational time (seconds) N v Sens. Present hr 9 min
57 Plate with quarter-circle cutout thickness: 0.005m Possion ratio: 0.3 load: 100kN/m modulus of elasticity: E = [199, 201]GPa number of element: 352 element type: six-node isoparametric quadratic triangle results presented: u, v, σ and σ at node F A E xx yy
58 Plate with quarter-circle cutout Case 1: the modulus of elasticity for each element varies independently in the interval [199, 201] GPa. Table: results of selected responses Monte Carlo sampling * Present IFEA Response LB UB LB UB u -5 A (10 5 m) v E (10-5 m) σ xx (MPa) σ yy (MPa) * 10 6 samples are made.
59 Outline Introduction Interval Arithmetic Interval Finite Elements Element-By By-Element Examples Conclusions
60 Conclusions Development and dimplementation ti of IFEM uncertain material, geometry and load parameters are described by interval variables interval arithmetic is used to guarantee an enclosure of response Enhanced dependence problem control use Element-By-Element t technique use the penalty method or Lagrange multiplier method to impose constraints modify and enhance fixed point iteration to take into account the dependence problem develop special algorithms to calculate stress and element nodal force
61 Conclusions The method is generally applicable to linear static FEM, regardless of element type Evaluation of the present method Rigorousness: in all the examples, the results obtained by the present method enclose those from the alternative methods Accuracy: sharp results are obtained for moderate parameter uncertainty (no more than 5%); reasonable results are obtained for relatively large parameter uncertainty t (5%~10%)
62 Conclusions Scalability: the accuracy of the method remains at the same level with increase of the problem scale Efficiency: the present method is significantly superior to the conventional methods such as the combinatorial, Monte Carlo sampling, and sensitivity analysis method The present IFEM represents an efficient method to handle uncertainty in engineering applications
63 Center for Reliable Engineering Computing (REC) We handle computations with care
64 Frame structure 7 4 w 3 w 4 Member Shape B 8 3 B 4 C 4 C 5 C 6 w 1 w 2 B 1 5 B 2 C 1 C 2 C m 14.63m m 6.1m C 1 W12 19 C 2 W C 3 W C 4 W10 12 C 5 W C 6 W B 1 W27 84 B 2 W B 3 W18 40 results listed: nodal forces at the left node of member B B 4 W
65 Frame structure case 1 Case 1: load uncertainty w w = [105.8,113.1] kn/m, w = [105.8,113.1] kn/m, 1 2 = [49.255,52.905] ] kn/m, w = [49.255,52.905] ] kn/m, 3 4 Table: Nodal forces at the left node of member B 2 Combinatorial Present IFEA Nodal force LB UB LB UB Axial (kn) Shear (kn) Moment (kn m) exact solution is obtained in the case of load uncertainty
66 Frame structure case 2 Case 2: stiffness uncertainty and load uncertainty 1% uncertainty introduced to AI,, and Eof each element. Number of finterval variables: 34. Table: Nodal forces at the left node of member B 2 Monte Carlo sampling * Present IFEA Nodal force LB UB LB UB Axial (kn) Shear (kn) Moment (kn.m) * 10 6 samples are made.
Interval Finite Element Methods for Uncertainty Treatment in Structural Engineering Mechanics Rafi L. Muhanna Georgia Institute of Technology USA
Interval Finite Element Methods for Uncertainty Treatment in Structural Engineering Mechanics Rafi L. Muhanna Georgia Institute of Technology USA Second Scandinavian Workshop on INTERVAL METHODS AND THEIR
More informationSharp Interval Estimates for Finite Element Solutions with Fuzzy Material Properties
Sharp Interval Estimates for Finite Element Solutions with Fuzzy Material Properties R. L. Mullen, F.ASCE Department of Civil Engineering, Case Western Reserve University, Cleveland, Ohio 44106-7201 rlm@po.cwru.edu
More informationGeometric Misfitting in Structures An Interval-Based Approach
Geometric Misfitting in Structures An Interval-Based Approach M. V. Rama Rao Vasavi College of Engineering, Hyderabad - 500 031 INDIA Rafi Muhanna School of Civil and Environmental Engineering Georgia
More informationNondeterministic Linear Static Finite Element Analysis: An Interval Approach. Hao Zhang
Nondeterministic Linear Static Finite Element Analysis: An Interval Approach A Thesis Presented to The Academic Faculty by Hao Zhang In Partial Fulfillment of the Requirements for the Degree Doctor of
More informationFinite Element Structural Analysis using Imprecise Probabilities Based on P-Box Representation
Finite Element Structural Analysis using Imprecise Probabilities Based on P-Box Representation Hao Zhang 1, R. L. Mullen 2, and R. L. Muhanna 3 1 University of Sydney, Sydney 2006, Australia, haozhang@usyd.edu.au
More informationComparison study of the computational methods for eigenvalues IFE analysis
Applied and Computational Mechanics 2 (2008) 157 166 Comparison study of the computational methods for eigenvalues IFE analysis M. Vaško a,,m.sága a,m.handrik a a Department of Applied Mechanics, Faculty
More informationWorst case analysis of mechanical structures by interval methods
Worst case analysis of mechanical structures by interval methods Arnold Neumaier (University of Vienna, Austria) (joint work with Andrzej Pownuk) Safety Safety studies in structural engineering are supposed
More informationIV B.Tech. I Semester Supplementary Examinations, February/March FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80
www..com www..com Code No: M0322/R07 Set No. 1 IV B.Tech. I Semester Supplementary Examinations, February/March - 2011 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer
More informationFinite Element Structural Analysis using Imprecise Probabilities Based on P-Box Representation
Finite Element Structural Analysis using Imprecise Probabilities Based on P-Box Representation Hao Zhang 1, R. L. Mullen 2, and R. L. Muhanna 3 1 University of Sydney, Sydney 2006, Australia, haozhang@usyd.edu.au
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 informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationFINAL EXAMINATION. (CE130-2 Mechanics of Materials)
UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130- Mechanics of Materials) Problem 1: (15 points) pinned -bar structure is shown in Figure 1. There is an external force, W = 5000N,
More informationDue Tuesday, September 21 st, 12:00 midnight
Due Tuesday, September 21 st, 12:00 midnight The first problem discusses a plane truss with inclined supports. You will need to modify the MatLab software from homework 1. The next 4 problems consider
More informationLecture 7: The Beam Element Equations.
4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite
More informationInterval-Based Parameter Identification for Structural Static Problems
Interval-Based Parameter Identification for Structural Static Problems Naijia Xiao School of Civil and Environmental Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA naijia.xiao@gatech.edu
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationChapter 11. Displacement Method of Analysis Slope Deflection Method
Chapter 11 Displacement ethod of Analysis Slope Deflection ethod Displacement ethod of Analysis Two main methods of analyzing indeterminate structure Force method The method of consistent deformations
More informationDesign of 2D Elastic Structures with the Interval Parameters
Design of 2D Elastic Structures with the Interval Parameters ANDRZEJ POWNUK The University of Texas at El Paso Department of Mathematical Sciences El Paso, Texas USA andrzej@pownuk.com NAVEEN KUMAR GOUD
More informationP.E. Civil Exam Review:
P.E. Civil Exam Review: Structural Analysis J.P. Mohsen Email: jpm@louisville.edu Structures Determinate Indeterminate STATICALLY DETERMINATE STATICALLY INDETERMINATE Stability and Determinacy of Trusses
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 informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationTypes of Structures & Loads
Structure Analysis I Chapter 4 1 Types of Structures & Loads 1Chapter Chapter 4 Internal lloading Developed in Structural Members Internal loading at a specified Point In General The loading for coplanar
More informationMechanical Models with Interval Parameters
In K. Gürlebeck L. Hempel C. Könke (Eds.) IKM2003: Digital Proceedings of 16th International Conference on the Applications of Computer Science and Mathematics in Architecture and Civil Engineering, ISSN
More informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
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 informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached
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 informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationMulti Linear Elastic and Plastic Link in SAP2000
26/01/2016 Marco Donà Multi Linear Elastic and Plastic Link in SAP2000 1 General principles Link object connects two joints, i and j, separated by length L, such that specialized structural behaviour may
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationHowever, reliability analysis is not limited to calculation of the probability of failure.
Probabilistic Analysis probabilistic analysis methods, including the first and second-order reliability methods, Monte Carlo simulation, Importance sampling, Latin Hypercube sampling, and stochastic expansions
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.
D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having
More informationInterval Finite Element Analysis of Thin Plates
Interval Finite Element Analsis of hin Plates M. V. Rama Rao ), R. L. Muhanna ) and R. L. Mullen ) ) Vasavi College of Engineering, Hderabad - 5 India, dr.mvrr@gmail.com ) School of Civil and Environmental
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 informationarxiv: v1 [math.na] 28 Jun 2013
A New Operator and Method for Solving Interval Linear Equations Milan Hladík July 1, 2013 arxiv:13066739v1 [mathna] 28 Jun 2013 Abstract We deal with interval linear systems of equations We present a new
More informationCIVL4332 L1 Introduction to Finite Element Method
CIVL L Introduction to Finite Element Method CIVL L Introduction to Finite Element Method by Joe Gattas, Faris Albermani Introduction The FEM is a numerical technique for solving physical problems such
More informationUNIT IV FLEXIBILTY AND STIFFNESS METHOD
SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : SA-II (13A01505) Year & Sem: III-B.Tech & I-Sem Course & Branch: B.Tech
More informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationStress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy
Stress Analysis Lecture 3 ME 276 Spring 2017-2018 Dr./ Ahmed Mohamed Nagib Elmekawy Axial Stress 2 Beam under the action of two tensile forces 3 Beam under the action of two tensile forces 4 Shear Stress
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationSensitivity and Reliability Analysis of Nonlinear Frame Structures
Sensitivity and Reliability Analysis of Nonlinear Frame Structures Michael H. Scott Associate Professor School of Civil and Construction Engineering Applied Mathematics and Computation Seminar April 8,
More informationME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM?
SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.
More informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationDevelopment of Truss Equations
CIVL 7/87 Chapter 3 - Truss Equations - Part /53 Chapter 3a Development of Truss Equations Learning Objectives To derive the stiffness matri for a bar element. To illustrate how to solve a bar assemblage
More informationBHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I
BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pin-jointed plane frame, shown in Figure Q1, is fixed to rigid supports at nodes and 4 to prevent their nodal displacements. The frame is loaded at nodes 1 and by a horizontal and a vertical
More informationThe University of Melbourne Engineering Mechanics
The University of Melbourne 436-291 Engineering Mechanics Tutorial Four Poisson s Ratio and Axial Loading Part A (Introductory) 1. (Problem 9-22 from Hibbeler - Statics and Mechanics of Materials) A short
More informationStochastic optimization - how to improve computational efficiency?
Stochastic optimization - how to improve computational efficiency? Christian Bucher Center of Mechanics and Structural Dynamics Vienna University of Technology & DYNARDO GmbH, Vienna Presentation at Czech
More informationComposites Design and Analysis. Stress Strain Relationship
Composites Design and Analysis Stress Strain Relationship Composite design and analysis Laminate Theory Manufacturing Methods Materials Composite Materials Design / Analysis Engineer Design Guidelines
More informationMulti-Point Constraints
Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Single point constraint examples Multi-Point constraint examples linear, homogeneous linear, non-homogeneous linear, homogeneous
More information3.4 Analysis for lateral loads
3.4 Analysis for lateral loads 3.4.1 Braced frames In this section, simple hand methods for the analysis of statically determinate or certain low-redundant braced structures is reviewed. Member Force Analysis
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationSabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in
Sabah Shawkat Cabinet of Structural Engineering 17 3.6 Shear walls Walls carrying vertical loads should be designed as columns. Basically walls are designed in the same manner as columns, but there are
More informationMEC-E8001 FINITE ELEMENT ANALYSIS
MEC-E800 FINIE EEMEN ANAYSIS 07 - WHY FINIE EEMENS AND IS HEORY? Design of machines and structures: Solution to stress or displacement by analytical method is often impossible due to complex geometry,
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationME325 EXAM I (Sample)
ME35 EXAM I (Sample) NAME: NOTE: COSED BOOK, COSED NOTES. ONY A SINGE 8.5x" ORMUA SHEET IS AOWED. ADDITIONA INORMATION IS AVAIABE ON THE AST PAGE O THIS EXAM. DO YOUR WORK ON THE EXAM ONY (NO SCRATCH PAPER
More informationFinite Element Method
Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional
More informationProblem d d d B C E D. 0.8d. Additional lecturebook examples 29 ME 323
Problem 9.1 Two beam segments, AC and CD, are connected together at C by a frictionless pin. Segment CD is cantilevered from a rigid support at D, and segment AC has a roller support at A. a) Determine
More informationon the figure. Someone has suggested that, in terms of the degrees of freedom x1 and M. Note that if you think the given 1.2
1) A two-story building frame is shown below. The mass of the frame is assumed to be lumped at the floor levels and the floor slabs are considered rigid. The floor masses and the story stiffnesses are
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 3 Torsion
EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 3 Torsion Introduction Stress and strain in components subjected to torque T Circular Cross-section shape Material Shaft design Non-circular
More informationPrepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING
CHETTINAD COLLEGE OF ENGINEERING AND TECHNOLOGY-KARUR FINITE ELEMENT ANALYSIS 2 MARKS QUESTIONS WITH ANSWER Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING FINITE ELEMENT ANALYSIS
More informationI certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.
NAME: ME 270 Fall 2012 Examination No. 3 - Makeup Please review the following statement: Group No.: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationRODS: THERMAL STRESS AND STRESS CONCENTRATION
RODS: HERML SRESS ND SRESS CONCENRION Example 5 rod of length L, cross-sectional area, and modulus of elasticity E, has been placed inside a tube of the same length L, but of cross-sectional area and modulus
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 informationMethod of Consistent Deformation
Method of onsistent eformation Structural nalysis y R.. Hibbeler Theory of Structures-II M Shahid Mehmood epartment of ivil Engineering Swedish ollege of Engineering and Technology, Wah antt FRMES Method
More informationNonconservative Loading: Overview
35 Nonconservative Loading: Overview 35 Chapter 35: NONCONSERVATIVE LOADING: OVERVIEW TABLE OF CONTENTS Page 35. Introduction..................... 35 3 35.2 Sources...................... 35 3 35.3 Three
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationPh.D. Preliminary Examination Analysis
UNIVERSITY OF CALIFORNIA, BERKELEY Spring Semester 2014 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials Name:......................................... Ph.D.
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
More informationInternal Internal Forces Forces
Internal Forces ENGR 221 March 19, 2003 Lecture Goals Internal Force in Structures Shear Forces Bending Moment Shear and Bending moment Diagrams Internal Forces and Bending The bending moment, M. Moment
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationIf the solution does not follow a logical thought process, it will be assumed in error.
Please indicate your group number (If applicable) Circle Your Instructor s Name and Section: MWF 8:30-9:20 AM Prof. Kai Ming Li MWF 2:30-3:20 PM Prof. Fabio Semperlotti MWF 9:30-10:20 AM Prof. Jim Jones
More informationFINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES
th International Conference on Earthquake Geotechnical Engineering June 5-8, 7 Paper No. 11 FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES A.Halim KARAŞĐN 1, Polat GÜLKAN ABSTRACT Plates on elastic foundations
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More information7. Hierarchical modeling examples
7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More information2 marks Questions and Answers
1. Define the term strain energy. A: Strain Energy of the elastic body is defined as the internal work done by the external load in deforming or straining the body. 2. Define the terms: Resilience and
More informationChapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship
Chapter 5 Elastic Strain, Deflection, and Stability Elastic Stress-Strain Relationship A stress in the x-direction causes a strain in the x-direction by σ x also causes a strain in the y-direction & z-direction
More informationMechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection
Mechanics of Materials II Chapter III A review of the fundamental formulation of stress, strain, and deflection Outline Introduction Assumtions and limitations Axial loading Torsion of circular shafts
More informationMath 175 Common Exam 2A Spring 2018
Math 175 Common Exam 2A Spring 2018 Part I: Short Form The first seven (7) pages are short answer. You don t need to show work. Partial credit will be rare and small. 1. (8 points) Suppose f(x) is a function
More informationSupplement: Statically Indeterminate Trusses and Frames
: Statically Indeterminate Trusses and Frames Approximate Analysis - In this supplement, we consider an approximate method of solving statically indeterminate trusses and frames subjected to lateral loads
More informationCOMPUTING UNCERTAINTIES IN A DELTA WING COMPOSITE STRUCTURE USING INTERVAL-BASED METHODS
COMPUTING UNCERTAINTIES IN A DELTA WING COMPOSITE STRUCTURE USING INTERVAL-BASED METHODS Franck Delcroix 1, Chris Boyer 2, Vincent Braibant 2 1 Laboratoire de Modélisation et Mécanique des Structures,
More informationInterval Finite Element Methods: New Directions
University of Texas at El Paso DigitalCommons@UTEP Departmental Technical Reports (CS) Department of Computer Science 1-1-2006 Interval Finite Element Methods: New Directions Rafi Muhanna Vladik Kreinovich
More information29. Define Stiffness matrix method. 30. What is the compatibility condition used in the flexibility method?
CLASS: III YEAR / VI SEMESTER CIVIL SUBJECTCODE AND NAME: CE 2351 - STRUCTURAL ANALYSIS-II UNIT1 FLEXIBILITY MATRIX METHOD. PART A 1. What is meant by indeterminate structures? 2. What are the conditions
More informationIntroduction to Finite Element Analysis Using Pro/MECHANICA Wildfire 5.0
Introduction to Finite Element Analysis Using Pro/MECHANICA Wildfire 5.0 Randy H. Shih Oregon Institute of Technology SDC PUBLICATIONS Schroff Development Corporation www.schroff.com Better Textbooks.
More informationHydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition
Fluid Structure Interaction and Moving Boundary Problems IV 63 Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition K.-H. Jeong, G.-M. Lee, T.-W. Kim & J.-I.
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationApplications of Eigenvalues & Eigenvectors
Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal
More informationTwo Tier projects for students in ME 160 class
ME 160 Introduction to Finite Element Method Spring 2016 Topics for Term Projects by Teams of 2 Students Instructor: Tai Ran Hsu, Professor, Dept. of Mechanical engineering, San Jose State University,
More information[7] Torsion. [7.1] Torsion. [7.2] Statically Indeterminate Torsion. [7] Torsion Page 1 of 21
[7] Torsion Page 1 of 21 [7] Torsion [7.1] Torsion [7.2] Statically Indeterminate Torsion [7] Torsion Page 2 of 21 [7.1] Torsion SHEAR STRAIN DUE TO TORSION 1) A shaft with a circular cross section is
More informationMTE 119 STATICS FINAL HELP SESSION REVIEW PROBLEMS PAGE 1 9 NAME & ID DATE. Example Problem P.1
MTE STATICS Example Problem P. Beer & Johnston, 004 by Mc Graw-Hill Companies, Inc. The structure shown consists of a beam of rectangular cross section (4in width, 8in height. (a Draw the shear and bending
More informationResearch Article Influence of the Parameterization in the Interval Solution of Elastic Beams
Structures Volume 04, Article ID 3953, 5 pages http://dx.doi.org/0.55/04/3953 Research Article Influence of the Parameterization in the Interval Solution of Elastic Beams Stefano Gabriele and Valerio Varano
More information