Interval Finite Element Methods for Uncertainty Treatment in Structural Engineering Mechanics Rafi L. Muhanna Georgia Institute of Technology USA
|
|
- Winfred Shelton
- 6 years ago
- Views:
Transcription
1 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 APPLICATIONS August 25-27, 27, 2003,Technical University of Denmark, Copenhagen, Denmark
2 Outline Introduction Interval Finite Elements Element-By By-Element Examples Conclusions
3 Acknowledgement Professor Bob Mullen Dr. Hao Zhang
4 Center for Reliable Engineering Computing (REC) We handle computations with care
5 Outline Introduction Interval Finite Elements Element-by by-element Examples Conclusions
6 Introduction- Uncertainty Uncertainty is unavoidable in engineering system structural mechanics entails uncertainties in material, geometry and load parameters 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)
7 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
8 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) Provides guaranteed enclosures
9 Introduction- Finite Element Method Finite Element Method (FEM) is a numerical method that provides approximate solutions to partial differential equations
10 Introduction- Uncertainty & Errors Mathematical model (validation) Discretization of the mathematical model into a computational framework Parameter uncertainty (loading, material properties) Rounding errors
11 Outline Introduction Interval Finite Elements Element-by by-element Examples Conclusions
12 Interval Finite Elements Uncertain Data Geometry Materials Loads Interval Stiffness Matrix T K B CB dv Interval Load Vector = F = N t da i i Element Level K U = F
13 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 load vector (traction) F = N t da i i N i = Shape function corresponding to the i-th DOF t = Surface traction
14 Interval Finite Elements (IFEM) Follows conventional FEM Loads, geometry and material property are expressed as interval quantities 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
15 IFEM- Inner-Bound Methods Combinatorial 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
16 IFEM- Enclosure Linear static finite 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
17 Interval arithmetic Interval number: Interval vector and interval matrix, e.g., T x = ( x, x ) = ([0,1], [ 2,1]) 1 2 midpoint, width, absolute value: defined componentwise Notations x = [ x, x] ( midpoint: x = ( x+ x) / 2, width: wid( x) = x x, absolute value: x = max{ x, x }. intervals: boldface, e.g., x, b, A real: non-boldface, x x, A A T
18 Linear interval equation Linear interval equation Ax = b ( A A, b b) Solution set Σ(A, b) = {x R A A b b: Ax = b} Hull of the solution set Σ(A, b) A H b := Σ(A, b)
19 Linear interval equation Example 2 [ 1, 0] x1 1.2 = [ 1, 0] 2 x Hull x Solution set 1.0 Enclosure x 2
20 Outline Introduction Interval Finite Elements Element-by by-element Examples Conclusions
21 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 1 = [1.429, 1.579], u 2 = [1.905, 2.105] naïve solution: u 1 = [-0.052, 3.052], u 2 = [0.098, 3.902] interval arithmetic assumes that all coefficients are independent uncertainty in the response is severely overestimated
22 Element-By-Element Element-By-Element (EBE) technique elements are detached no element coupling structure stiffness matrix is block-diagonal (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
23 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 + kd ( ( k : deterministic part, element stiffness matrix evalued using E, ( kd : interval part d: interval diagonal matrix, diag( δ,..., δ).
24 Element-By-Element : : ( Element stiffness matrix: k = k( I + d) Structure stiffness matrix: ( ( ( K = K( I + D) = K + KD or K ( k 1 k 1 d 1 = O = O I + O ( N k k e N N e d e
25 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
26 Constraints penalty method Constraint conditions: cu = 0 Using the penalty method: ( K + Q) u = p T Q: penalty matrix, Q= c ηc η: diagonal matrix of penalty number η Requires a careful choice of the penatly number i E 1, A 1, L 1 E 2, A 2, L 2 A spring of large stiffness is added to force node 2 and node 3 to have the same displacement.
27 Constraints Lagrange multiplier Constraint conditions: cu = 0 Using the Lagrange multiplier method: T K c u p = c 0 λ 0 λ : Lagrange multiplier vector, introdued as new unknowns
28 Load in EBE Nodal load applied by elements p b T where p = N φ( x) dx i p = ( p,..., p ) Suppose the surface traction φ( x) is described by j an interval function: φ( x) = a x. p b can be rewritten as b p 1 b W: deterministic matrix F: interval vector containing the interval coefficients of the surface tractiton m j= 0 N j = WF e T
29 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 ) ( ( ( 1 where z = Rb RAx0, G = I RA, R= A, Ax0 = b No dependency handling
30 Fixed point iteration, Interval FEA calls for a modified method which exploits the special form of the structure equations ( ( ( K + Q) u= p with K = K + KD ( 1, Choose R= ( K + Q), construct iterations: *( l+ 1) *( l) u = Rp R( K + Q) u0 + ( I R( K + Q) )( ε u ) ( *( l) = Rp u0 RKD ( u0+ε u ) ( () l = Rp u0 RKM ( + + if u int( u ), then u= u + u = Rp RKM *( l 1) *( l) *( l 1) ( l) 0 : interval vector, =( δ,..., δ ) 1 1 N e The interval variables δ,..., δ appear only once in each iteration. Most sources of dependence are eliminated. N e T
31 Convergence of fixed point The algorithm converges if and only if To minimize ρ( G ): ( 1 choose R = A so that G = I RA : 1 ρ ( G ) < 1 smaller ρ( G ) less iterations required, and less overestimation in results has a small spectral radius reduce the overestimation in G ( ( ( ( 1 G= I RA=I ( K + Q) ( K + Q+ KD) = RKD
32 Stress calculation Conventional method: σ = CBu Present method: e, (severe overestimation) C: elasticity matrix, B: strain-displacement matrix ( ( E = (1 + δ) E, C = (1 + δ) C σ = CBLu ( () l =CBL( Rp RKM ) ( ( ( = + δ p M L: Boolean matrix, L u= ue () l (1 )( CBLR CBLRK )
33 Element nodal force calculation Conventional method: f = ( ku p ), (severe overestimation) T e e e Present method: ( Te) 1( k1( ue) 1 ( p ) 1) e in the EBE model, T ( Ku pb ) = M ( 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
34 Outline Introduction Interval Finite Elements Element-by by-element Examples Conclusions
35 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
36 Truss structure P m P 2 P 1 P 3 4.5m 4.5m 4.5m 4.5m 2 1, 2, 3 13, 14, 15: [9.95, 10.05] cm (1% uncertainty) A A A,A A A cross-sectional area 2 of all other elements: [5.97, 6.03] cm (1% uncertainty) modulus of elasticity of all elements: 200,000 MPa p p = [190, 210] kn, p = [95,105] kn 1 2 = [95,105] kn, p = [85.5,94.5] kn (10% uncertainty) 3 4
37 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 δ Present IFEA δ % % % % % % % % unit: u 5 (m), N 7 (kn). LB: lower bound; UB: upper bound.
38 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 = [ , ], I u5 FEM: = [ , ] (relative difference: 2.59%, 1.80% for LB, UB, respectively)
39 Frame structure w 3 w B 8 3 B 4 C 4 C 5 C 6 w 1 w B 1 5 B 2 C 1 C 2 C m 14.63m 4.57m 6.1m Member Shape C 1 W12 19 C 2 C 3 W C 4 W C 5 W10 12 C 6 W B 1 W W27 84 B 2 W B 3 W18 40 W27 94 results listed: nodal forces at the left node of member B 2 B 4
40 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 Nodal force Axial (kn) Shear (kn) Moment (kn m) Combinatorial LB UB Present IFEA LB UB exact solution is obtained in the case of load uncertainty
41 Frame structure case 2 Case 2: stiffness uncertainty and load uncertainty 1% uncertainty introduced to AI,, and Eof each element. Number of interval variables: 34. Table: Nodal forces at the left node of member B 2 Nodal force Axial (kn) Shear (kn) Moment (kn.m) Monte Carlo sampling * LB UB Present IFEA LB UB * 10 6 samples are made.
42 Truss with a large number of interval variables p p p p C D p p n@l p B A A i E i m@l = [0.995,1.005] A, = [0.995,1.005] E for i = 1,..., N 0 0 e story bay N e N v
43 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 ba y 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% % 0.45% 3.56% δ LB = LB- LB * / LB *, δ LB = UB- UB * / UB *, δ LB = (LB- LB * )/ LB *
44 Efficiency study Table: CPU time for the analyses with the present method (unit: seconds) Story bay N v Iteratio n t i t t i /t t r /t t r % 12.4% 10.2% 9.7% 9.1% 9.2% 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 time is spent on matrix inversion
45 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
46 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: ua, ve, σxx and σyy at node F
47 Plate case 1 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 * Response LB UB u A (10-5 m) v E (10-5 m) σ xx (MPa) σ yy (MPa) * 10 6 samples are made. Present IFEA LB UB
48 Plate case 2 Case 2: each subdomain has an independent modulus of elasticity. E i = [199, 201] GPa, for i = 1,..., m 0.025m 0.025m 0.025m E 1 E2 E 5 E m E 3 E 4 E 7 E m
49 Plate case 2 Table: results of selected responses Response u A (10-5 m) v E (10-5 m) σ xx (MPa) σ yy (MPa) Combinatorial LB UB Present IFEA LB UB
50 Outline Introduction Interval Finite Elements Element-by by-element Examples Conclusions
51 Conclusions Development and implementation 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 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
52 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 (5%~10%)
53 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
Uncertainty is unavoidable in engineering system. Probabilistic approach is the traditional approach
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
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationAlvaro F. M. Azevedo A. Adão da Fonseca
SECOND-ORDER SHAPE OPTIMIZATION OF A STEEL BRIDGE Alvaro F. M. Azevedo A. Adão da Fonseca Faculty of Engineering University of Porto Portugal 16-18 March 1999 OPTI 99 Orlando - Florida - USA 1 PROBLEM
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 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 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 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 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 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 information4 Finite Element Method for Trusses
4 Finite Element Method for Trusses To solve the system of linear equations that arises in IPM, it is necessary to assemble the geometric matrix B a. For the sake of simplicity, the applied force vector
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 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 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 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 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 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 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 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 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 informationEsben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer
Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics
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 informationPractice Final Examination. Please initial the statement below to show that you have read it
EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use
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 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 15: Revisiting bars and beams.
3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
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 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 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 informationStiffness Matrices, Spring and Bar Elements
CHAPTER Stiffness Matrices, Spring and Bar Elements. INTRODUCTION The primary characteristics of a finite element are embodied in the element stiffness matrix. For a structural finite element, the stiffness
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 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 information6. Bending CHAPTER OBJECTIVES
CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where
More informationINTRODUCTION TO STRAIN
SIMPLE STRAIN INTRODUCTION TO STRAIN In general terms, Strain is a geometric quantity that measures the deformation of a body. There are two types of strain: normal strain: characterizes dimensional changes,
More informationCHAPTER OBJECTIVES CHAPTER OUTLINE. 4. Axial Load
CHAPTER OBJECTIVES Determine deformation of axially loaded members Develop a method to find support reactions when it cannot be determined from euilibrium euations Analyze the effects of thermal stress
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 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 informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
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 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 informationEFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM
Proceedings of the International Conference on Mechanical Engineering 2007 (ICME2007) 29-31 December 2007, Dhaka, Bangladesh ICME2007-AM-76 EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE
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 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 informationDHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)
DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types
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 informationOptimal Preconditioning for the Interval Parametric Gauss Seidel Method
Optimal Preconditioning for the Interval Parametric Gauss Seidel Method Milan Hladík Faculty of Mathematics and Physics, Charles University in Prague, Czech Republic http://kam.mff.cuni.cz/~hladik/ SCAN,
More informationChapter 4-b Axially Loaded Members
CIVL 222 STRENGTH OF MATERIALS Chapter 4-b Axially Loaded Members AXIAL LOADED MEMBERS Today s Objectives: Students will be able to: a) Determine the elastic deformation of axially loaded member b) Apply
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 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 informationBasic Energy Principles in Stiffness Analysis
Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting
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 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 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 informationStructural Analysis I Chapter 4 - Torsion TORSION
ORSION orsional stress results from the action of torsional or twisting moments acting about the longitudinal axis of a shaft. he effect of the application of a torsional moment, combined with appropriate
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 informationJ. Harish* and R.P. Rokade + *Graduate scholar Valliammai Engineering College, Kattankulathur, India + Principal Scientist
International Journal of Scientific & Engineering Research Volume 8, Issue 6, June-2017 109 Stochastic finite element analysis of simple hip truss Abstract- In this paper, nodal displacements and member
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 informationFREE VIBRATIONS OF FRAMED STRUCTURES WITH INCLINED MEMBERS
FREE VIBRATIONS OF FRAMED STRUCTURES WITH INCLINED MEMBERS A Thesis submitted in partial fulfillment of the requirements for the degree of Bachelor of Technology in Civil Engineering By JYOTI PRAKASH SAMAL
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 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 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 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 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 informationUnit 15 Shearing and Torsion (and Bending) of Shell Beams
Unit 15 Shearing and Torsion (and Bending) of Shell Beams Readings: Rivello Ch. 9, section 8.7 (again), section 7.6 T & G 126, 127 Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics and Engineering
More informationAluminum shell. Brass core. 40 in
PROBLEM #1 (22 points) A solid brass core is connected to a hollow rod made of aluminum. Both are attached at each end to a rigid plate as shown in Fig. 1. The moduli of aluminum and brass are EA=11,000
More informationSolution: T, A1, A2, A3, L1, L2, L3, E1, E2, E3, P are known Five equations in five unknowns, F1, F2, F3, ua and va
ME 323 Examination # 1 February 18, 2016 Name (Print) (Last) (First) Instructor PROBLEM #1 (20 points) A structure is constructed from members 1, 2 and 3, with these members made up of the same material
More informationTHE USE OF DYNAMIC RELAXATION TO SOLVE THE DIFFERENTIAL EQUATION DESCRIBING THE SHAPE OF THE TALLEST POSSIBLE BUILDING
VII International Conference on Textile Composites and Inflatable Structures STRUCTURAL MEMBRANES 2015 E. Oñate, K.-U.Bletzinger and B. Kröplin (Eds) THE USE OF DYNAMIC RELAXATION TO SOLVE THE DIFFERENTIAL
More informationCIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives
CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a
More informationFolded-plate structures with openings analysis and optimization
Folded-plate structures with openings analysis and optimization Bogdan Szybiński, Andrzej P. Zieliński Institute of Mechanics & Machine Design, Cracow University of Technology Marek Karaś Institute of
More information