Robust Reference Tension Optimization in Winding Systems Using Wound Internal Stress Calculation

Size: px
Start display at page:

Download "Robust Reference Tension Optimization in Winding Systems Using Wound Internal Stress Calculation"

Transcription

1 01 AIMCAL Europe and USA Web Coating Conference & Web Handling Conference Robust Reference Tension Optimization in Winding Systems Using Wound Internal Stress Calculation By 1 M hamed Boutaous & Patrick Bourgin CETHIL, INSA Lyon - France Mhamed.boutaous@insa-lyon.fr Dominique Knittel Université de Strasbourg, France

2 OUTLINE Introduction Optimization of the winding tension - Offline optimization - Online optimization Influence of the dynamic gauge Effect of the reference tension variations rate Optimization of the admissible tension dispersion Conclusion

3 Introduction : Winding flexible media Main defects - Wrinkles, - Interlayer slipping, - Waviness Why? - Technology or Process (geometry, design, tension, winding velocity, nip force,...) - Product (mechanical properties, surface topography ) Challenge in winding : good quality rolls. Mathematical models of the stress within the roll. The role of the TENSION REFERENCE in stress state. 3 New proposal: control the tension in terms of computed internal stress

4 Role of the tension reference Influence of the nominal tension on the stress components. Stress Winding tension (Kg) Stress 8 x Tangentiel stress Radius Radius Effect of tension fluctuation.

5 Optimization of the winding tension J Optimum reference tension = industrial know-how. Model of stress computation Criterion J The reference tension which minimizes criterion J is optimized using an algorithm based on the principle of the simplex (Nelder and Mead, 1965). R = max R roll (σ θ (T w ) σ θmean (r)) g(σ θ,r)) dr 10 x 10 Tangential stresses and stress gauge g (σ θ, r) weight function defined by : 9 8 g (σ θ, r) = 1 if is in the gauge g (σ θ, r) >>1 if else. Stress (Pa) (a) (b) Normalized Radius 5 Optimization gauge for tangential stress

6 Offline optimization of the reference tension Distributed control: Sensorless measurement of the tension with control based on fuzzy logics and neuronal networks.(ex: Wolfermann.) Offline optimization of the reference tension - No perfect follow-up of the instruction - Fluctuations due to the imperfections of the winding chain. - Induced variations of the stress within the roll, as compared with the one which would result from an ideal control. NEED TO MODIFY THE REFERENCE TENSION AFTER EACH LAYER: A NEW CONCEPT 6

7 Online optimization of the reference tension Winding tension reference Winding Tension of the global wound J R = max ( σ ( T ) σ ( r)) g( σ, r)) dr R roll w mean Pack 1... Pack i measured tension Corrector Winding Plant Pack 1... Pack i measures g (σ, r) weight function : g (σ, r) = 1 if σ is in the gauge g (σ, r) >>1 if else. J= min (J T, J R ) Pack i+1... Pack n Radius Prediction algorithm Stress gauge Predictions Pack i+1... Pack n Computation of stress state Computation and minimization of J Principle: The reference tension is modified (adjusted) layer after layer according to the criterion J At time ti, which corresponds to the winding of pack i. Packs 1 to i-1 are already wound: we know their tension reference and their actual values (measured). The instruction for i is also clearly defined and it is now, possible to calculate the instruction from i+1 to N, then to compare the values of the tension measured with the ones predicted by the model for a fictive roll. By the gauge of the tangential stress, we will calculate the weight function. The optimization of this function by minimizing it according to the principle of the simplex makes it possible to find the optimal online tension reference for the remainder of the roll: i+1 to N. The next instructions will be thus updated as the roll is being wound. 7

8 Examples: Comparison Winding web tension Online measurements Winding web tension Prediction of the reference tension Offline On / Off Initial tension Final tension : Time t i Tf on Tf off Offline optimized reference tension Pack 1 Radius Pack Pack i Pack i+1 Pack n Offline and online reference tension 8 Dividing the roll radius into several segments, the tension reference is computed and corrected for each range of roll radius values, by using the predictive model for the stress within the roll. The adjusted tension is updated step by step, following the optimization principle as described above and it will be considered as the new tension reference value for the coming layers

9 Examples: Control with simulated noise in tension measurement (Constant gauge) 8 x 10 Curve Offline, real and optimized control tension with random noise Radial stresses 3 1 Curve The tangential stresses are more affected than the radial stresses - The perturbation affects only the wound part, but not the remainder part Tangential stress at several time during the winding process 9

10 local perturbation and constant gauge - Before the perturbation, the offline and the online tensions are the same. - We observe the different steps of the online optimization 10

11 Effect of a dynamic gauge The gauge is fixed during all the process of winding for the offline optimization. Periodic disturbance in the reference tension Corresponding radial and tangential stresses The gauge is assumed 10 times higher from the pack N 0 to 55 for the online one, for both the radial and the tangential stresses. 11

12 Effect of a dynamic gauge The gauge is fixed during all the process of winding for the offline optimization, Periodic disturbances in the reference tension The gauge is assumed 10 times higher from the pack N 0 to 55 in the online one but just for the radial stresses. Corresponding radial and tangential stresses 1 The result is that the tangential stresses in the central zone are lower compared to the precedent case.

13 Effect of a dynamic gauge Corresponding radial and tangential stresses Periodic disturbance in the reference tension The gauge is assumed 10 times lower from the pack of layers N 0 to 55 in the online one but just for the radial stresses. 13 For example, in the case of risk to have negative tangential stresses, this optimization is a way to prevent it.

14 Maximum Reference Tension Gauge? T max? and T min? such as T min < T ref < T max σ (T ref ) located between fixed stress gauge Problem to solve: T ref +/- C: C maximum =? T ref T max C can also be C(r) T min C=? 1

15 Effect of the reference tension variations rate during winding The control of winding systems is generally based on practical experience and the tension reference does not change or decreases according to a more or less complicated function of the radius. The optimum tension is that which guarantees that the stresses within the roll still confined in a desired gauge The practical question is to know the limits of variations of the reference tension, i. e. witch slope of decreasing to impose to the reference tension without exceeding prescribed stresses values. 10 x T decreases from 18 to N: Tmin T decreases from 0 to 6 N: T max T decreases from 18 to 6 N: T T decreases from 0 to N: T1 Tangential stresses ( Pa) Layers number

16 Effect of the reference tension variations rate during winding In the precedent figures, surprisingly, the obtained results are different from that expected. Indeed, the stresses σ(t 1 ) calculated with T 1 as reference, is not confined between σ( T max ) and σ( T min ). The same observation has been made for tangential stresses σ(t ) computed with T as reference tension. Conclusion: The slope of the reference tension influences the stresses as well as the tension value itself. The magnitude of the tension is a key parameter, but not sufficient. The rate of the tension decrease plays an important role too. 16

17 Effect of the fluctuations with maximum slope around the reference tension Consider now the extreme case: T varies alternatively, with a maximum slope, between T max and T min : for each new wound layer, T passes alternatively from T max to T min. Hence, the applied reference tension, named T alt is given by T alt = (a.r +b) +/- c (curves in the figure below are given with c = 1 or ) 1.8 x 105 T varies alternatively between Tmax and Tmin Tension (N) Layers Number Alternative variations of the tension between maximal and minimal values

18 Effect of the fluctuations with maximum slope around the reference tension Tangential Stresses (Pa) 10 x T decreases from 18 to N T decreases from 0 to 6N T = (a* r+b ) +/- 1 T = (a* r+b ) +/- 1- Any reference tension which lies between T min (r) and T max (r) gives stresses with values confined between those calculated for a tension which varies alternatively between T min (r) and T max (r) Layers number - If the dispersion (T max (r) - T min (r) ) is too large, the stresses can overcross their gauge, and one risks to have defects in the roll. It is thus important to calculate this authorized maximum dispersion. 18

19 Effect of the fluctuations with maximum slope around the reference tension If the fluctuations become more and more (or less and less) important with the radius. The parameter c becomes function of the radius, and the reference tension presents an attenuation (or an exaggeration) of its fluctuations during winding. Reference Tension (N) x 105 Tw = (a*r +b)+-c c increases with the radius r c decreases with the radius r Layers Number Radial stresses (Pa) Tangential stresses (Pa) 0 x 10 Radial Stresses x 10 Tangential Stresses 5 c increasing with the radius c decreasing with the radius c increasing with the radius c decreasing with the radius Layers Number The stresses remain globally comparable. But the risk of web wrinkling due to negative tangential stresses, when the fluctuation decreases (triangle symbols), is more important than in the inverse case (circle symbols). The same phenomena is confirmed when the roll radius is more important. 19

20 Effect of the fluctuations with maximum slope around the reference tension Radial stresses (Pa) x 105 Radial stresses c decreases with the radius c increases with the radius r Tangential stresses (Pa) 8 x 105 Tangential stresses Layers Number c decreases with the radius r lower limit of the gauge c increases with the radius r 0 The same phenomena is confirmed when the roll radius is more important.

21 Optimization of the ADMISSIBLE TENSION DISPERSION Tref = ar + b ± C The maximum value for C (admissible dispersion) is obtained by minimization of the criterion J using the simplex optimization algorithm (Nelder and Mean 1967) J = ( σ σ ( T )) g( σ ) g (σ θ ) weight function defined by : g (σ θ ) = 1 if θ critical θ ref θ σ σ σ θmin θ( T ref ) θmax g (σ θ ) >>1 else 1

22 Optimization of the ADMISSIBLE TENSION DISPERSION C is choosen constant 9 Optimized winding tension 8 7 Tension (N) x 10 Optimized tangential stresses Number of layers C = Tangential stresses (Pa) Number of layers

23 Optimization of the ADMISSIBLE TENSION DISPERSION C can be C (r) C (r) = α r + β α =?, β =? 9 Optimized winding tension 8 Tension (N) Number of layers Tangential stresses (Pa) x 10 Optimized tangential stresses 5 3 C(r) max = Number of layers

24 CONCLUSION -Well-known observation: the winding tension strongly affects the roll internal stress. -BASIC IDEA: control the tension with respect to the stress generated within a roll. -NEW CONCEPT: ONLINE ADAPTATION OF THE REFERENCE TENSION -The new optimization strategy called online optimization associated with a dynamic gauge is used and compared to the classical offline one. The comparison shows interesting improvements, in providing an internal stress state compatible with elastic deformations of the web within the roll. -A second approach consists to optimize the tension reference gauge in order to make the choice of the winding reference confined in a constraining gauge, with respect to the stress state within the roll.

25 CONCLUSION Two applications can be aimed by this approach: - the estimation of the dispersions of the measured tension with regard to the reference one, let to deduce if the wound roll stress state respects the stress gauge, without recalculating the stresses. - for a defined stress gauge, one can calculate the reference tension gauge, with respect to the stress state within the roll. By defining au Stress Gauge We estimate the T ref max (T ref optimal or not) 5

Robust Reference Tension Optimization in Winding Systems Using Wound Internal Stress Calculation

Robust Reference Tension Optimization in Winding Systems Using Wound Internal Stress Calculation obust eference Tension Optimization in Winding Systems Using Wound Internal Stress Calculation M hamed BOUTAOUS (), Patrick BOUGIN (), Dominique KNITTEL () () Centre de Thermique de Lyon (CETHIL), INSA,

More information

Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases

Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of

More information

Disturbance Rejection in Parameter-varying Web-winding Systems

Disturbance Rejection in Parameter-varying Web-winding Systems Proceedings of the 17th World Congress The International Federation of Automatic Control Disturbance Rejection in Parameter-varying Web-winding Systems Hua Zhong Lucy Y. Pao Electrical and Computer Engineering

More information

H unwinding web tension control of a strip processing plant using a pendulum dancer

H unwinding web tension control of a strip processing plant using a pendulum dancer 2009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 10-12, 2009 WeB07.3 H unwinding web tension control of a strip processing plant using a pendulum dancer Vincent Gassmann,

More information

Supplementary Material

Supplementary Material 1 2 3 Topological defects in confined populations of spindle-shaped cells by G. Duclos et al. Supplementary Material 4 5 6 7 8 9 10 11 12 13 Supplementary Note 1: Characteristic time associated with the

More information

MOTION OF A CHARGED PARTICLE IN A UNIFORM MAGNETIC FIELD

MOTION OF A CHARGED PARTICLE IN A UNIFORM MAGNETIC FIELD MOTION OF A CHARGED PARTICLE IN A UNIFORM MAGNETIC FIELD When the velocity of a charged particle is perpendicular to a uniform magnetic field, the particle moves in a circular path in a plane perpendicular

More information

Real-Time Feasibility of Nonlinear Predictive Control for Semi-batch Reactors

Real-Time Feasibility of Nonlinear Predictive Control for Semi-batch Reactors European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. Real-Time Feasibility of Nonlinear Predictive Control

More information

FCP Short Course. Ductile and Brittle Fracture. Stephen D. Downing. Mechanical Science and Engineering

FCP Short Course. Ductile and Brittle Fracture. Stephen D. Downing. Mechanical Science and Engineering FCP Short Course Ductile and Brittle Fracture Stephen D. Downing Mechanical Science and Engineering 001-015 University of Illinois Board of Trustees, All Rights Reserved Agenda Limit theorems Plane Stress

More information

2010 F=ma Solutions. that is

2010 F=ma Solutions. that is 2010 F=ma Solutions 1. The slope of a position vs time graph gives the velocity of the object So you can see that the position from B to D gives the steepest slope, so the speed is the greatest in that

More information

Resolution of scheduling problems by using a generic optimisation software

Resolution of scheduling problems by using a generic optimisation software Resolution of scheduling problems by using a generic optimisation software Christian HOULBERT Centre National d Etudes Spatiales (CNES) 18 avenue Edouard Belin 31401 Toulouse France Tél. 05 61 27 33 38

More information

Robust discrete time control

Robust discrete time control Robust discrete time control I.D. Landau Emeritus Research Director at C.N.R. Laboratoire d utomatique de Grenoble, INPG/CNR, France pril 2004, Valencia I.D. Landau course on robust discrete time control,

More information

Physics 53. Dynamics 2. For every complex problem there is one solution that is simple, neat and wrong. H.L. Mencken

Physics 53. Dynamics 2. For every complex problem there is one solution that is simple, neat and wrong. H.L. Mencken Physics 53 Dynamics 2 For every complex problem there is one solution that is simple, neat and wrong. H.L. Mencken Force laws for macroscopic objects Newton s program mandates studying nature in order

More information

7.6 Journal Bearings

7.6 Journal Bearings 7.6 Journal Bearings 7.6 Journal Bearings Procedures and Strategies, page 1 of 2 Procedures and Strategies for Solving Problems Involving Frictional Forces on Journal Bearings For problems involving a

More information

Robust Control 9 Design Summary & Examples

Robust Control 9 Design Summary & Examples Robust Control 9 Design Summary & Examples Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/6/003 Outline he H Problem Solution via γ-iteration Robust stability via

More information

7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment

7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment 7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment à It is more difficult to obtain an exact solution to this problem since the presence of the shear force means that

More information

Aircraft Stress Analysis and Structural Design Summary

Aircraft Stress Analysis and Structural Design Summary Aircraft Stress Analysis and Structural Design Summary 1. Trusses 1.1 Determinacy in Truss Structures 1.1.1 Introduction to determinacy A truss structure is a structure consisting of members, connected

More information

Shape Optimization of Oldham Coupling in Scroll Compressor

Shape Optimization of Oldham Coupling in Scroll Compressor Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 24 Shape Optimization of Oldham Coupling in Scroll Compressor In Hwe Koo LG Electronics

More information

Physics 8 Monday, October 28, 2013

Physics 8 Monday, October 28, 2013 Physics 8 Monday, October 28, 2013 Turn in HW8 today. I ll make them less difficult in the future! Rotation is a hard topic. And these were hard problems. HW9 (due Friday) is 7 conceptual + 8 calculation

More information

Regulating Web Tension in Tape Systems with Time-varying Radii

Regulating Web Tension in Tape Systems with Time-varying Radii Regulating Web Tension in Tape Systems with Time-varying Radii Hua Zhong and Lucy Y. Pao Abstract A tape system is time-varying as tape winds from one reel to the other. The variations in reel radii consist

More information

Mechanics Answers to Examples B (Momentum) - 1 David Apsley

Mechanics Answers to Examples B (Momentum) - 1 David Apsley TOPIC B: MOMENTUM ANSWERS SPRING 2019 (Full worked answers follow on later pages) Q1. (a) 2.26 m s 2 (b) 5.89 m s 2 Q2. 8.41 m s 2 and 4.20 m s 2 ; 841 N Q3. (a) 1.70 m s 1 (b) 1.86 s Q4. (a) 1 s (b) 1.5

More information

Introduction to Control (034040) lecture no. 2

Introduction to Control (034040) lecture no. 2 Introduction to Control (034040) lecture no. 2 Leonid Mirkin Faculty of Mechanical Engineering Technion IIT Setup: Abstract control problem to begin with y P(s) u where P is a plant u is a control signal

More information

Supplementary Figure 1 Extracting process of wetting ridge profiles. a1-4, An extraction example of a ridge profile for E 16 kpa.

Supplementary Figure 1 Extracting process of wetting ridge profiles. a1-4, An extraction example of a ridge profile for E 16 kpa. Supplementary Figure 1 Extracting process of wetting ridge profiles. a1-4, An extraction example of a ridge profile for E 16 kpa. An original image (a1) was binarized, as shown in a2, by Canny edge detector

More information

Practice Final Examination. Please initial the statement below to show that you have read it

Practice 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 information

PLATE GIRDERS II. Load. Web plate Welds A Longitudinal elevation. Fig. 1 A typical Plate Girder

PLATE GIRDERS II. Load. Web plate Welds A Longitudinal elevation. Fig. 1 A typical Plate Girder 16 PLATE GIRDERS II 1.0 INTRODUCTION This chapter describes the current practice for the design of plate girders adopting meaningful simplifications of the equations derived in the chapter on Plate Girders

More information

For an imposed stress history consisting of a rapidly applied step-function jump in

For an imposed stress history consisting of a rapidly applied step-function jump in Problem 2 (20 points) MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0239 2.002 MECHANICS AND MATERIALS II SOLUTION for QUIZ NO. October 5, 2003 For

More information

Physics 121, Final Exam Do not turn the pages of the exam until you are instructed to do so.

Physics 121, Final Exam Do not turn the pages of the exam until you are instructed to do so. , Final Exam Do not turn the pages of the exam until you are instructed to do so. You are responsible for reading the following rules carefully before beginning. Exam rules: You may use only a writing

More information

Rheological and Engineering Properties of Orange Pulp. Elyse Payne Juan Fernando Muñoz José I. Reyes De Corcuera

Rheological and Engineering Properties of Orange Pulp. Elyse Payne Juan Fernando Muñoz José I. Reyes De Corcuera Rheological and Engineering Properties of Orange Pulp Elyse Payne Juan Fernando Muñoz José I. Reyes De Corcuera September 20, 2012 Acknowledgements Industry Dr. Paul Winniczuk Mr. Marcelo Bellarde Mr.

More information

Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS

Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS of the 18 th International Conference on Process Control Hotel Titris, Tatranská

More information

TOPIC D: ROTATION EXAMPLES SPRING 2018

TOPIC D: ROTATION EXAMPLES SPRING 2018 TOPIC D: ROTATION EXAMPLES SPRING 018 Q1. A car accelerates uniformly from rest to 80 km hr 1 in 6 s. The wheels have a radius of 30 cm. What is the angular acceleration of the wheels? Q. The University

More information

General Physics I. Lecture 8: Rotation of a Rigid Object About a Fixed Axis. Prof. WAN, Xin ( 万歆 )

General Physics I. Lecture 8: Rotation of a Rigid Object About a Fixed Axis. Prof. WAN, Xin ( 万歆 ) General Physics I Lecture 8: Rotation of a Rigid Object About a Fixed Axis Prof. WAN, Xin ( 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ New Territory Object In the past, point particle (no rotation,

More information

EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES

EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES OF WOOD OR PLYWOOD CONSIDERED AS ORTHOTROPIC MATERIALS Information Revied and Reaffirmed March 1956 No. 1510 EFFECT OF ELLIPTIC

More information

Let the plant and controller be described as:-

Let the plant and controller be described as:- Summary of Fundamental Limitations in Feedback Design (LTI SISO Systems) From Chapter 6 of A FIRST GRADUATE COURSE IN FEEDBACK CONTROL By J. S. Freudenberg (Winter 2008) Prepared by: Hammad Munawar (Institute

More information

6. Find the net torque on the wheel in Figure about the axle through O if a = 10.0 cm and b = 25.0 cm.

6. Find the net torque on the wheel in Figure about the axle through O if a = 10.0 cm and b = 25.0 cm. 1. During a certain period of time, the angular position of a swinging door is described by θ = 5.00 + 10.0t + 2.00t 2, where θ is in radians and t is in seconds. Determine the angular position, angular

More information

Investigation Of The Parameters Affecting Crankshaft And Rotor Interference Fit

Investigation Of The Parameters Affecting Crankshaft And Rotor Interference Fit Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2014 Investigation Of The Parameters Affecting Crankshaft And Rotor Interference Fit Naz?m

More information

SELECTED PROBLEMS OF SHORT CIRCUIT WITHSTANDABILITY Section II - POWER TRANSFORMER October 2004, Vigo - Spain

SELECTED PROBLEMS OF SHORT CIRCUIT WITHSTANDABILITY Section II - POWER TRANSFORMER October 2004, Vigo - Spain Dr. Władysław Pewca Institute of Power Engineering, Transformer Division (IenOT( IenOT), Poland SELECTED PROBLEMS OF SHORT CIRCUIT WITHSTANDABILITY Section II - POWER TRANSFORMER 28-30 October 2004, Vigo

More information

M340(921) Solutions Problem Set 6 (c) 2013, Philip D Loewen. g = 35y y y 3.

M340(921) Solutions Problem Set 6 (c) 2013, Philip D Loewen. g = 35y y y 3. M340(92) Solutions Problem Set 6 (c) 203, Philip D Loewen. (a) If each pig is fed y kilograms of corn, y 2 kilos of tankage, and y 3 kilos of alfalfa, the cost per pig is g = 35y + 30y 2 + 25y 3. The nutritional

More information

Iterative Learning Control for Tailor Rolled Blanks Zhengyang Lu DCT

Iterative Learning Control for Tailor Rolled Blanks Zhengyang Lu DCT Iterative Learning Control for Tailor Rolled Blans Zhengyang Lu DCT 2007.066 Master Internship Report Supervisors: Prof. Maarten Steinbuch (TU/e) Camile Hol (Corus) Eindhoven University of Technology Department

More information

Physics 111: Mechanics Lecture 9

Physics 111: Mechanics Lecture 9 Physics 111: Mechanics Lecture 9 Bin Chen NJIT Physics Department Circular Motion q 3.4 Motion in a Circle q 5.4 Dynamics of Circular Motion If it weren t for the spinning, all the galaxies would collapse

More information

Lecture 34 Minimization and maximization of functions

Lecture 34 Minimization and maximization of functions Lecture 34 Minimization and maximization of functions Introduction Golden section search Parabolic interpolation Search with first derivatives Downhill simplex method Introduction In a nutshell, you are

More information

Switched systems: stability

Switched systems: stability Switched systems: stability OUTLINE Switched Systems Stability of Switched Systems OUTLINE Switched Systems Stability of Switched Systems a family of systems SWITCHED SYSTEMS SWITCHED SYSTEMS a family

More information

UNIT 2 KINEMATICS OF LINKAGE MECHANISMS

UNIT 2 KINEMATICS OF LINKAGE MECHANISMS UNIT 2 KINEMATICS OF LINKAGE MECHANISMS ABSOLUTE AND RELATIVE VELOCITY An absolute velocity is the velocity of a point measured from a fixed point (normally the ground or anything rigidly attached to the

More information

N = N A Pb A Pb. = ln N Q v kt. = kt ln v N

N = N A Pb A Pb. = ln N Q v kt. = kt ln v N 5. Calculate the energy for vacancy formation in silver, given that the equilibrium number of vacancies at 800 C (1073 K) is 3.6 10 3 m 3. The atomic weight and density (at 800 C) for silver are, respectively,

More information

BUTT SPLICE HINGING. KEVIN COLE, PhD Senior Web Handling Development Engineer Optimation Technology Incorporated

BUTT SPLICE HINGING. KEVIN COLE, PhD Senior Web Handling Development Engineer Optimation Technology Incorporated BUTT SPLICE HINGING BY KEVIN COLE, PhD Senior Web Handling Development Engineer Optimation Technology Incorporated Introduction Splicing is a process used to join the tail of an expiring roll to the start

More information

Ch 7 Homework. (a) Label physical quantities in this problem using letters you choose.

Ch 7 Homework. (a) Label physical quantities in this problem using letters you choose. Ch 7 Homework Name: Homework problems are from the Serway & Vuille 10 th edition. Follow the instructions and show your work clearly. 1. (Problem 7) A machine part rotates at an angular speed of 0.06 rad/s;

More information

8.8 Applications of Taylor Polynomials

8.8 Applications of Taylor Polynomials 8.8 Applications of Taylor Polynomials Mark Woodard Furman U Spring 2008 Mark Woodard (Furman U) 8.8 Applications of Taylor Polynomials Spring 2008 1 / 14 Outline 1 Point estimation 2 Estimation on an

More information

I. D. Landau, A. Karimi: A Course on Adaptive Control Adaptive Control. Part 9: Adaptive Control with Multiple Models and Switching

I. D. Landau, A. Karimi: A Course on Adaptive Control Adaptive Control. Part 9: Adaptive Control with Multiple Models and Switching I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 1 Adaptive Control Part 9: Adaptive Control with Multiple Models and Switching I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 2 Outline

More information

Excerpt from the Proceedings of the COMSOL Conference 2010 Boston

Excerpt from the Proceedings of the COMSOL Conference 2010 Boston Excerpt from the Proceedings of the COMSOL Conference 21 Boston Uncertainty Analysis, Verification and Validation of a Stress Concentration in a Cantilever Beam S. Kargar *, D.M. Bardot. University of

More information

Pressure corrected SPH for fluid animation

Pressure corrected SPH for fluid animation Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional

More information

Integer Sequences and Circle Chains Inside a Tangential Quadrilateral

Integer Sequences and Circle Chains Inside a Tangential Quadrilateral Sangaku Journal of Mathematics (SJM c SJM ISSN 534-956 Volume (018, pp.31-40 Received August 018. Published on-line 10 August 018 web: http://www.sangaku-journal.eu/ c The Author(s This article is published

More information

The Learning Objectives of the Compulsory Part Notes:

The Learning Objectives of the Compulsory Part Notes: 17 The Learning Objectives of the Compulsory Part Notes: 1. Learning units are grouped under three strands ( Number and Algebra, Measures, Shape and Space and Data Handling ) and a Further Learning Unit.

More information

Beam Design and Deflections

Beam Design and Deflections Beam Design and Deflections tation: a = name for width dimension A = name for area Areq d-adj = area required at allowable stress when shear is adjusted to include self weight Aweb = area of the web of

More information

Multi-Objective Trajectory Planning of Mobile Robots Using Augmented Lagrangian

Multi-Objective Trajectory Planning of Mobile Robots Using Augmented Lagrangian CSC 007 6-8 May Marraech Morocco /0 Multi-Objective rajectory Planning of Mobile Robots Using Augmented Lagrangian Amar Khouhi Luc Baron and Mare Balazinsi Mechanical Engineering Dept. École Polytechnique

More information

Effect of taper tension profile on the telescoping in a winding process of high speed roll to roll printing systems

Effect of taper tension profile on the telescoping in a winding process of high speed roll to roll printing systems Journal of Mechanical Science and Technology 23 (2009) 3036~3048 Journal of Mechanical Science and Technology www.springerlink.com/content/1738-494x DOI 10.1007/s12206-009-0906-2 Effect of taper tension

More information

PHYSICS 221 Fall 2016 EXAM 2: November 02, :15pm 10:15pm. Name (printed): Recitation Instructor: Section #:

PHYSICS 221 Fall 2016 EXAM 2: November 02, :15pm 10:15pm. Name (printed): Recitation Instructor: Section #: PHYSICS 221 Fall 2016 EXAM 2: November 02, 2016 8:15pm 10:15pm Name (printed): Recitation Instructor: Section #: INSTRUCTIONS: This exam contains 25 multiple-choice questions, plus 2 extra-credit questions,

More information

CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER

CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER 114 CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER 5.1 INTRODUCTION Robust control is a branch of control theory that explicitly deals with uncertainty in its approach to controller design. It also refers

More information

Supporting Information

Supporting Information Supporting Information Calculation of the swelling ratios of DN hydrogels To estimate χ parameter for PAAm water system, eqs 1a and 1b were solved using the experimentally determined swelling ratios and

More information

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator. MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY

More information

Basis of Structural Design

Basis of Structural Design Basis of Structural Design Course 2 Structural action: cables and arches Course notes are available for download at http://www.ct.upt.ro/users/aurelstratan/ Structural action Structural action: the way

More information

Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f

Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f Material undergoes plastic deformation when stress exceeds yield stress σ 0 Permanent strain results from

More information

Evolutionary Multiobjective. Optimization Methods for the Shape Design of Industrial Electromagnetic Devices. P. Di Barba, University of Pavia, Italy

Evolutionary Multiobjective. Optimization Methods for the Shape Design of Industrial Electromagnetic Devices. P. Di Barba, University of Pavia, Italy Evolutionary Multiobjective Optimization Methods for the Shape Design of Industrial Electromagnetic Devices P. Di Barba, University of Pavia, Italy INTRODUCTION Evolutionary Multiobjective Optimization

More information

MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD

MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used

More information

Computational Methodology for the Prediction of Functional Requirement Variations Across the Product Life-Cycle

Computational Methodology for the Prediction of Functional Requirement Variations Across the Product Life-Cycle 1 Computational Methodology for the Prediction of Functional Requirement Variations Across the Product Life-Cycle Guillaume Mandil (PhD Student) Pr. Alain Desrochers (Université de Sherbrooke) Pr. Alain

More information

H -Optimal Control and Related Minimax Design Problems

H -Optimal Control and Related Minimax Design Problems Tamer Başar Pierre Bernhard H -Optimal Control and Related Minimax Design Problems A Dynamic Game Approach Second Edition 1995 Birkhäuser Boston Basel Berlin Contents Preface v 1 A General Introduction

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland 1/21 Outline 1 Lecture 4: Modeling Tools for Mechanical Systems

More information

Ch. 10: Fundamental of contact between solids

Ch. 10: Fundamental of contact between solids Ch. 10: Fundamental of contact between solids Actual surface is not smooth. At atomic scale, there are always defects at surface, such as vacancies, ledges, kinks, terraces. In micro or macro scale, roughness

More information

Acceleration due to Gravity

Acceleration due to Gravity Acceleration due to Gravity 1 Object To determine the acceleration due to gravity by different methods. 2 Apparatus Balance, ball bearing, clamps, electric timers, meter stick, paper strips, precision

More information

Physics 12. Unit 5 Circular Motion and Gravitation Part 1

Physics 12. Unit 5 Circular Motion and Gravitation Part 1 Physics 12 Unit 5 Circular Motion and Gravitation Part 1 1. Nonlinear motions According to the Newton s first law, an object remains its tendency of motion as long as there is no external force acting

More information

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 10 Columns

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 10 Columns EMA 370 Mechanics & Materials Science (Mechanics of Materials) Chapter 10 Columns Columns Introduction Columns are vertical prismatic members subjected to compressive forces Goals: 1. Study the stability

More information

Rotation. PHYS 101 Previous Exam Problems CHAPTER

Rotation. PHYS 101 Previous Exam Problems CHAPTER PHYS 101 Previous Exam Problems CHAPTER 10 Rotation Rotational kinematics Rotational inertia (moment of inertia) Kinetic energy Torque Newton s 2 nd law Work, power & energy conservation 1. Assume that

More information

How materials work. Compression Tension Bending Torsion

How materials work. Compression Tension Bending Torsion Materials How materials work Compression Tension Bending Torsion Elemental material atoms: A. Composition a) Nucleus: protons (+), neutrons (0) b) Electrons (-) B. Neutral charge, i.e., # electrons = #

More information

Chapter 3. Load and Stress Analysis

Chapter 3. Load and Stress Analysis Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3

More information

DEVIL PHYSICS BADDEST CLASS ON CAMPUS IB PHYSICS

DEVIL PHYSICS BADDEST CLASS ON CAMPUS IB PHYSICS DEVIL PHYSICS BADDEST CLASS ON CAMPUS IB PHYSICS OPTION B-1A: ROTATIONAL DYNAMICS Essential Idea: The basic laws of mechanics have an extension when equivalent principles are applied to rotation. Actual

More information

Final Exam April 30, 2013

Final Exam April 30, 2013 Final Exam Instructions: You have 120 minutes to complete this exam. This is a closed-book, closed-notes exam. You are allowed to use a calculator during the exam. Usage of mobile phones and other electronic

More information

Application of the immersed boundary method to simulate flows inside and outside the nozzles

Application of the immersed boundary method to simulate flows inside and outside the nozzles Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,

More information

Radial Compressors. Damian Vogt Course MJ2429. Nomenclature

Radial Compressors. Damian Vogt Course MJ2429. Nomenclature Turbomachinery Lecture Notes 1 007-10-04 Radial Compressors Damian Vogt Course MJ49 Nomenclature Subscripts Symbol Denotation Unit c Absolute velocity m/s h Enthalpy J/kg m& Mass flow rate kg/s r Radius

More information

Some improvements of Xfem for cracked domains

Some improvements of Xfem for cracked domains Some improvements of Xfem for cracked domains E. Chahine 1, P. Laborde 2, J. Pommier 1, Y. Renard 3 and M. Salaün 4 (1) INSA Toulouse, laboratoire MIP, CNRS UMR 5640, Complexe scientifique de Rangueil,

More information

Adaptive Dual Control

Adaptive Dual Control Adaptive Dual Control Björn Wittenmark Department of Automatic Control, Lund Institute of Technology Box 118, S-221 00 Lund, Sweden email: bjorn@control.lth.se Keywords: Dual control, stochastic control,

More information

Find the magnitude of F when t = 2. (9 marks)

Find the magnitude of F when t = 2. (9 marks) Condensed M2 Paper These questions are all taken from a Mechanics 2 exam paper, but any intermediate steps and diagrams have been removed, leaving enough information to answer the question, but none of

More information

VIII. Rubber Elasticity [B.Erman, J.E.Mark, Structure and properties of rubberlike networks]

VIII. Rubber Elasticity [B.Erman, J.E.Mark, Structure and properties of rubberlike networks] VIII. Rubber Elasticity [B.Erman, J.E.Mark, Structure and properties of rubberlike networks] Using various chemistry, one can chemically crosslink polymer chains. With sufficient cross-linking, the polymer

More information

Quiz #7. T f k m A gsinθ = m A a N m A gcosθ = 0 f k = µ k N m B g T = m B a

Quiz #7. T f k m A gsinθ = m A a N m A gcosθ = 0 f k = µ k N m B g T = m B a Quiz #7 Vector The method used in 2 Dimensions is exactly the same in 3D; just keep one more components (Pythagorean theorem also still holds in higher dimensions as long as the space is Euclidean). 1)

More information

Generally, there exists an optimum tip-speed-ratio, λ that maximized C p. The exact λ depends on the individual wind turbine design

Generally, there exists an optimum tip-speed-ratio, λ that maximized C p. The exact λ depends on the individual wind turbine design Summary Chapter 6-End 1 Wind Turbine Control The control system on a wind turbine is designed to: 1. seek the highest efficiency of operation that maximizes the coefficient of power, C p, 2. ensure safe

More information

Flin Flon Mining Belt

Flin Flon Mining Belt EOSC433: Geotechnical Engineering Practice & Design Lecture 7: Stress Analysis around Underground Openings 1 of 40 Erik Eberhardt UBC Geological Engineering EOSC 433 (2007) Flin Flon Mining Belt Since

More information

Critical Borehole Orientations Rock Mechanics Aspects

Critical Borehole Orientations Rock Mechanics Aspects Critical Borehole Orientations Rock Mechanics Aspects By R. BRAUN* Abstract This article discusses rock mechanics aspects of the relationship between borehole stability and borehole orientation. Two kinds

More information

General Physics I. Lecture 8: Rotation of a Rigid Object About a Fixed Axis. Prof. WAN, Xin ( 万歆 )

General Physics I. Lecture 8: Rotation of a Rigid Object About a Fixed Axis. Prof. WAN, Xin ( 万歆 ) General Physics I Lecture 8: Rotation of a Rigid Object About a Fixed Axis Prof. WAN, Xin ( 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ New Territory Object In the past, point particle (no rotation,

More information

Optimization with COMSOL Multiphysics COMSOL Tokyo Conference Walter Frei, PhD Applications Engineer

Optimization with COMSOL Multiphysics COMSOL Tokyo Conference Walter Frei, PhD Applications Engineer Optimization with COMSOL Multiphysics COMSOL Tokyo Conference 2014 Walter Frei, PhD Applications Engineer Product Suite COMSOL 5.0 Agenda An introduction to optimization A lot of concepts, and a little

More information

This chapter focuses on the study of the numerical approximation of threedimensional

This chapter focuses on the study of the numerical approximation of threedimensional 6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation

More information

On the Inherent Robustness of Suboptimal Model Predictive Control

On the Inherent Robustness of Suboptimal Model Predictive Control On the Inherent Robustness of Suboptimal Model Predictive Control James B. Rawlings, Gabriele Pannocchia, Stephen J. Wright, and Cuyler N. Bates Department of Chemical & Biological Engineering Computer

More information

Fundamentals of Structural Design Part of Steel Structures

Fundamentals of Structural Design Part of Steel Structures Fundamentals of Structural Design Part of Steel Structures Civil Engineering for Bachelors 133FSTD Teacher: Zdeněk Sokol Office number: B619 1 Syllabus of lectures 1. Introduction, history of steel structures,

More information

Exam 3 PREP Chapters 6, 7, 8

Exam 3 PREP Chapters 6, 7, 8 PHY241 - General Physics I Dr. Carlson, Fall 2013 Prep Exam 3 PREP Chapters 6, 7, 8 Name TRUE/FALSE. Write 'T' if the statement is true and 'F' if the statement is false. 1) Astronauts in orbiting satellites

More information

Chapter 12 Elastic Stability of Columns

Chapter 12 Elastic Stability of Columns Chapter 12 Elastic Stability of Columns Axial compressive loads can cause a sudden lateral deflection (Buckling) For columns made of elastic-perfectly plastic materials, P cr Depends primarily on E and

More information

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA MECHANICAL PRINCIPLES AND APPLICATIONS NQF LEVEL 3 OUTCOME 1 - LOADING SYSTEMS TUTORIAL 3 LOADED COMPONENTS

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA MECHANICAL PRINCIPLES AND APPLICATIONS NQF LEVEL 3 OUTCOME 1 - LOADING SYSTEMS TUTORIAL 3 LOADED COMPONENTS EDEXCEL NATIONAL CERTIICATE/DIPLOMA MECHANICAL PRINCIPLES AND APPLICATIONS NQ LEVEL 3 OUTCOME 1 - LOADING SYSTEMS TUTORIAL 3 LOADED COMPONENTS 1. Be able to determine the effects of loading in static engineering

More information

Analytical and Numerical Solution for Collapse Load in Oil Pipelines under Compressive Load

Analytical and Numerical Solution for Collapse Load in Oil Pipelines under Compressive Load Proceedings of the World Congress on Engineering 01 Vol I WCE 01, July 4-6, 01, ondon, U.K. Analytical and Numerical Solution for Collapse oad in Oil Pipelines under Compressive oad uciano M. Bezerra,

More information

Influence of the filament winding process variables on the mechanical behavior of a composite pressure vessel

Influence of the filament winding process variables on the mechanical behavior of a composite pressure vessel Influence of the filament winding process variables on the mechanical behavior of a composite pressure vessel G. Vargas 1 & A. Miravete 2 1 Universidad Pontificia Bolivariana, Facultad de Ingeniería Mecánica,

More information

A Luminosity Leveling Method for LHC Luminosity Upgrade using an Early Separation Scheme

A Luminosity Leveling Method for LHC Luminosity Upgrade using an Early Separation Scheme LHC Project Note 03 May 007 guido.sterbini@cern.ch A Luminosity Leveling Method for LHC Luminosity Upgrade using an Early Separation Scheme G. Sterbini and J.-P. Koutchouk, CERN Keywords: LHC Luminosity

More information

Simulation in Computer Graphics Elastic Solids. Matthias Teschner

Simulation in Computer Graphics Elastic Solids. Matthias Teschner Simulation in Computer Graphics Elastic Solids Matthias Teschner Outline Introduction Elastic forces Miscellaneous Collision handling Visualization University of Freiburg Computer Science Department 2

More information

Spatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media

Spatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Spatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media Karim G. Sabra, School of Mechanical Engineering,

More information

Imperfect Gases. NC State University

Imperfect Gases. NC State University Chemistry 431 Lecture 3 Imperfect Gases NC State University The Compression Factor One way to represent the relationship between ideal and real gases is to plot the deviation from ideality as the gas is

More information

1 Numerical optimization

1 Numerical optimization Contents Numerical optimization 5. Optimization of single-variable functions.............................. 5.. Golden Section Search..................................... 6.. Fibonacci Search........................................

More information

7.2.1 Seismic waves. Waves in a mass- spring system

7.2.1 Seismic waves. Waves in a mass- spring system 7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-

More information

Inhomogeneous elastic response of amorphous solids

Inhomogeneous elastic response of amorphous solids Inhomogeneous elastic response of amorphous solids Jean-Louis Barrat Université de Lyon Institut Universitaire de France Acknowledgements: Anne Tanguy, Fabien Chay Goldenberg, Léonforte, Michel Tsamados

More information