Development and experimental verification of a numerical model for optimizing the cleaning performance of the doctor blade press roll tribosystem

Similar documents
3. Using your answers to the two previous questions, evaluate the Mratio

Stability of Smart Beams with Varying Properties Based on the First Order Shear Deformation Theory Located on a Continuous Elastic Foundation

6. Non-uniform bending

Theoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump

A = h w (1) Error Analysis Physics 141

INTRODUCTION AND MATHEMATICAL CONCEPTS

Section A 01. (12 M) (s 2 s 3 ) = 313 s 2 = s 1, h 3 = h 4 (s 1 s 3 ) = kj/kgk. = kj/kgk. 313 (s 3 s 4f ) = ln

Analysis of Static and Dynamic Load on Hydrostatic Bearing with Variable Viscosity and Pressure

Model development for the beveling of quartz crystal blanks

Higher Derivatives. Differentiable Functions

Sample Problems for Exam II

Non-linear Dynamics of Inlet Film Thickness during Unsteady Rolling Process

f a h f a h h lim lim

Finding and Using Derivative The shortcuts

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

AN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES

The effects of shear stress on the lubrication performances of oil film of large-scale mill bearing

Nonlinear correction to the bending stiffness of a damaged composite beam

A general articulation angle stability model for non-slewing articulated mobile cranes on slopes *

Physics 121, April 1, Equilibrium. Physics 121. April 1, Physics 121. April 1, Course Information. Discussion of Exam # 2

Static Response Analysis of a FGM Timoshenko s Beam Subjected to Uniformly Distributed Loading Condition

Tribology in Industry

Derivatives of Exponentials

Prediction of Oil-Water Two Phase Flow Pressure Drop by Using Homogeneous Model

Polynomial Interpolation

Section 15.6 Directional Derivatives and the Gradient Vector

Low Decibel Linear Way E

Derivative as Instantaneous Rate of Change

UNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics

INTRODUCTION AND MATHEMATICAL CONCEPTS

Non-linear Analysis Method of Ground Response Using Equivalent Single-degree-of-freedom Model

Derivatives and Rates of Change

INTRODUCTION DEFINITION OF FLUID. U p F FLUID IS A SUBSTANCE THAT CAN NOT SUPPORT SHEAR FORCES OF ANY MAGNITUDE WITHOUT CONTINUOUS DEFORMATION

Large eddy simulation of turbulent flow downstream of a backward-facing step

Optimization of the thin-walled rod with an open profile

Friction Coefficient s Numerical Determination for Hot Flat Steel Rolling at Low Strain Rate

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator

The Derivative The rate of change

Problem Solving. Problem Solving Process

Chapter 9. τ all = min(0.30s ut,0.40s y ) = min[0.30(58), 0.40(32)] = min(17.4, 12.8) = 12.8 kpsi 2(32) (5/16)(4)(2) 2F hl. = 18.1 kpsi Ans. 1.

= h. Geometrically this quantity represents the slope of the secant line connecting the points

HOMEWORK HELP 2 FOR MATH 151

Study of Convective Heat Transfer through Micro Channels with Different Configurations

June : 2016 (CBCS) Body. Load

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)

Homework Assignment on Fluid Statics

Grade: 11 International Physics Olympiad Qualifier Set: 2

Some Review Problems for First Midterm Mathematics 1300, Calculus 1

UNIT-1 MODERN PHYSICS

pancakes. A typical pancake also appears in the sketch above. The pancake at height x (which is the fraction x of the total height of the cone) has

Continuous Stochastic Processes

Math 34A Practice Final Solutions Fall 2007

Continuity and Differentiability of the Trigonometric Functions

MVT and Rolle s Theorem

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001

Outline. MS121: IT Mathematics. Limits & Continuity Rates of Change & Tangents. Is there a limit to how fast a man can run?

CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY

Work and Energy. Introduction. Work. PHY energy - J. Hedberg

Empirical models for estimating liquefaction-induced lateral spread displacement

Comment on Experimental observations of saltwater up-coning

CFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E.

Krazy Katt, the mechanical cat

Finite Element Analysis of J-Integral for Surface Cracks in Round Bars under Combined Mode I Loading

, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim

Part 2: Introduction to Open-Channel Flow SPRING 2005

Squeeze film effect at elastohydrodynamic lubrication of plain journal bearings

2.1 THE DEFINITION OF DERIVATIVE

1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).

The entransy dissipation minimization principle under given heat duty and heat transfer area conditions

Click here to see an animation of the derivative

A Multiaxial Variable Amplitude Fatigue Life Prediction Method Based on a Plane Per Plane Damage Assessment

Notes: Most of the material in this chapter is taken from Young and Freedman, Chap. 12.

Department of Mechanical Engineering, Azarbaijan Shahid Madani University, Tabriz, Iran b

Desalination by vacuum membrane distillation: sensitivity analysis

An optimum design of robotic food handling by using Burger model

Analysis of Stress and Deflection about Steel-Concrete Composite Girders Considering Slippage and Shrink & Creep Under Bending

Shear Thinning and Hydrodynamic Friction of Viscosity Modifier- Containing Oils. Part II: Impact of Shear Thinning on Journal Bearing Friction

Bulk Deformation Processes

Development of new and validation of existing convection correlations for rooms with displacement ventilation systems

Film thickness Hydrodynamic pressure Liquid saturation pressure or dissolved gases saturation pressure. dy. Mass flow rates per unit length

lecture 26: Richardson extrapolation

Taylor Series and the Mean Value Theorem of Derivatives

arxiv:math/ v1 [math.ca] 1 Oct 2003

Consider the element shown in Figure 2.1. The statement of energy conservation applied to this element in a time period t is that:

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

Week #15 - Word Problems & Differential Equations Section 8.2

Derivatives. By: OpenStaxCollege

Lecture 10: Carnot theorem

Carnot Factor of a Vapour Power Cycle with Regenerative Extraction

HT TURBULENT NATURAL CONVECTION IN A DIFFERENTIALLY HEATED VERTICAL CHANNEL. Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008

Extracting Atomic and Molecular Parameters From the de Broglie Bohr Model of the Atom

University of Alabama Department of Physics and Astronomy PH 101 LeClair Summer Exam 1 Solutions

The Basics of Vacuum Technology

The derivative function

Journal of Engineering Science and Technology Review 7 (4) (2014) 40-45

Polynomial Interpolation

Bending analysis of a functionally graded piezoelectric cantilever beam

THE MODEL OF DRYING SESSILE DROP OF COLLOIDAL SOLUTION

Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using the differential transform method

1 2 x Solution. The function f x is only defined when x 0, so we will assume that x 0 for the remainder of the solution. f x. f x h f x.

Transcription:

Development and experimental verification of a numerical model for optimizing te cleaning performance of te doctor blade press roll tribosystem M. Rodríguez Ripoll, B. Sceicl,, D. Bianci, B. Jakab, F. Franek, A²T Researc Gmb, Austrian Excellence entre for Tribology, Viktor Kaplan Strasse, 700 Wiener Neustadt, Austria Vienna University of Tecnology, Institute of Fluid Mecanics and eat Transfer, Resselgasse /E, 040 Vienna, Austria Vienna University of Tecnology, Institute of Sensor and Actuator Systems, Floragasse 7/, 040 Vienna, Austria Email: ripoll@act.at Abstract: In te paper production, doctor (scraping) blades are placed in contact wit press rolls during wet pressing so as to purge te surface of te rolls from processing water, contamination, and stickies. Te contact is acieved by mounting te blade on a older, wic is tilted around a rotation axis until te blade tip contacts wit te roll. Te contact force is determined by te supply pressure of te air forced troug te tube tat is placed at te bottom of te older. Due to contact, te blade wears off and needs to be replaced periodically. Our aim is to optimize te cleaning performance of te system by modelling te tribological contact between doctor blade and press roll, in order to acieve an optimum cleaning performance, tus increasing te blade lifetime and reducing energy consumption. Te model is susceptible to an inextensive numerical evaluation. A first comparison wit experimental findings is encouraging. Key words: Paper industry, doctor blade, ydrodynamic lubrication, wear. INTRODUTION Modern papermaking factories are divided in four main sections, namely te forming section, te press section, te drying section and te calender section. Troug tese sections, te paper pulp is processed until a omogeneous paper seet is produced []. In te present paper, we focus on te wet pressing section, were te paper pulp gradually reduces its water content by passing between series of rolls and by being pressed against tem. As te paper pulp passes over te press roll and moves to te next one (Figure ), a certain amount of water containing contamination particles remains on te roll s surface. In order to keep its surface clean from particles and avoid contamination of te continuously flowing pulp, a scraping blade known as doctor blade is pressed against te roll s surface. Te blade contact force against te roll is determined by te air pressure applied to a polymer ose. A large pressure flow causes a

iger contact pressure at te blade tip, wic improves te cleaning performance at expenses of iger energy consumption and an increased wear rate of te blade. Figure. Press section of a paper macine. Te doctor blade keeps te press roll surface clean from contaminations and splas water. Due to te interest of te industry to save energy costs and an increasing pressure of lawmakers to reduce te immense energy consumption at paper mills, tere is great potential for energy saving if a more efficient cleaning strategy could be acieved. Te large number of subsequent rolls used in te press section and teir large dimensions provide room for improvement in order to reduce energy costs, avoid pollution and increase profit. Nowadays te optimum process parameters set at te doctor blade press roll tribosystem are based on te accumulated experience during years. Te aim of tis work is to develop a compreensive model of te tribosystem in order to optimise te cleaning performance of doctor blades troug a systematic and scientific approac.. NOTATION In wat follows, a list of te symbols and variables used in te paper is provided (see also Figure ). R F : Reaction force acting on te blade tip F : Force perpendicular to te blade tip F : Force parallel to te blade tip F µ : Friction force F N : Normal force v M : Roll speed vector at te contact point, wose components are v and v R FA : Reaction force acting on te rotating axis M : Momentum of te blade older r : Vector tat links te end of te older wit te blade tip r 0 : Initial position of r c0 before rotation F : Force per unit lengt of te pressurized ose r : Vector tat links te rotation axis wit te center of te ose p : Pressure of te ose E : Young s Modulus of te ose L : Lengt (perimeter) of te relaxed ose d : eigt of te compressed ose r: Vector tat links te rotation axis wit te blade tip x: Displacement of te blade due to bending R x : Displacement of te ose arm due to pressure canges γ : Rotation angle of te older due to bending of te blade β: Angle between te tangent to te mill and te blade µ: oefficient of friction E: Young s Modulus of te doctor blade I: Second moment of inertia of te doctor blade d: Tickness of te doctor blade l: Free lengt of te doctor blade l r U: Sliding speed at te surface of te press roll δ: Opening angle at te blade tip due to bending η: Viscosity of te lubricating fluid L: ontact lengt of te doctor blade 0 : Gap eigt at fluid intake : Gap eigt at fluid outtake So: Sommerfeld number K: onvergence ratio W: Blade widt in te out of plane direction α: Angle of te blade tip : ardness of te doctor blade k: Arcard wear coefficient k D : Arcard dimension wear coefficient V: Removed wear volume

guarantied as long as te film tickness between te doctor blade and te press roll is kept small. (a) (b) In a recently model proposed by te autors [], te doctor blade is considered as a pad bearing sliding over te roll surface. Assuming equilibrium conditions between ydrodynamic and contact forces, a non-dimensional group involving te key parameters was obtained. Its optimum value can be calculated by imposing te onset of ydrodynamic sliding conditions. In te present work, a mecanical model governing te deflection of te blade yields a relationsip between te linear contact force at its tip and te imposed air pressure. Te cange of te blade tip due to wear is taken into account by applying a geometrical wear model. Under tese assumptions, te influence of te process parameters on te cleaning performance can be systematically analysed... ydrodynamic Model A toroug description of te ydrodynamic model is found in []. In wat follows, only te most salient features are described for te sake of completeness. Te doctor blade is considered as a one dimensional beam, wose deflection is governed by te Euler-Bernoulli elastic beam teory Figure. (a) Scematic representation of te doctor blade and its older. (b) Detail of te doctor blade tip. (c) Detail of te pressurized ose.. MODELLING.. Modelling Approac Te modelling approac is based on te assumption tat an optimum cleaning performance can be acieved by forcing te doctor blade to operate under weak ydrodynamic conditions. In tis case, friction and wear are minimised, wile an optimum cleaning performance is (c) F Ed δ, () 6l were E is te Young s Modulus of te blade, l and d are te blade lengt and tickness and δ is te opening angle, 0 L sinδ ~ δ. () Imposing equilibrium conditions at te blade tip between ydrodynamic and contact forces, te following relationsip is obtained F cos β, () F N were β is te positioning angle of te blade (blade angle).

Te normal load F N : can be calculated analytically for a pad bearing using te Reynolds equation [] 6UηL So, (4) F N were U denotes te sliding velocity, η te viscosity, L te contact lengt and te minimum film tickness, i.e. te gap eigt at fluid discarge. So is te Sommerfeld number for te lubricated wedge flow and ence a function of te convergence ratio K 0 (5) solely: So So(K). Upon substituting eq. () and eq. (4) into eq. (), a linear furter relationsip between So and K is establised: So Ed. (6) 6UηL l cos β ( K ) K From tis relationsip one determines a unique value of K (Figure ). Figure. Relationsip between te Sommerfeld Number So and te convergence ratio K, for two different contact lengts L. An optimum cleaning performance is acieved for a small gap and small converging ratios ( 0 ~ ), i.e. wen K 0 in eq. (6). By tis means, te following condition for a non-dimensional group involving te key parameter is obtained: 6UηL Ed l cos β. (7).. Geometrical canges due to wear During te paper production, doctor blades need to be replaced periodically due to wear. During te wear process, te geometry of te blade canges, wic modifies te contact conditions at te tip. Based on purely geometrical relations, te actual blade lengt l and its contact lengt L (i.e. at a certain point of time) can be calculated for a given wear volume V [] V sinα L, (8) W sin β sin ( π α β ) ( cotα cot β ) V + l, (9) W were W is te blade widt in te out of plane direction. Te removed wear volume at every differential increment dv can be simply imposed or derived from te operational parameters by using a wear law, suc as te one proposed by Arcard [4]: FN U dv k dt kdfnudt, (0) were k D is a dimension wear coefficient to be fitted to experimental results, F N is te normal load, U is te sliding speed, te material ardness, and dt is te differential time increment..4. Mecanical Model Te ydrodynamic model and te relationsip between te geometrical key quantities at te blade tip, wic were recalled in te previous sections, rely in te fact tat te contact force at te blade tip is known. owever, in a paper mill, te only parameter readily available to te operator is te air pressure applied to te polymer ose. In tis section, a simple mecanical model aims to link bot pysical quantities. Te doctor blade is assumed to be a noncompressible one dimensional beam connected to a fixed older. Te older is free to rotate around

an axis in order to contact te press roll. Te contact force depends on te air pressure flowing troug a ose placed below te older. Te reaction force at te blade tip can be decomposed in two components, one component perpendicular F and one component parallel F to te blade tip contact line: F ( F F ) R +. () Te friction force F µ between te blade and te roll is given as v M FN Fµ µ FN µfn ˆv, () v F M N were v is a unit vector tangent to te roll at te contact point and F N and F µ are given by F F F, () + N F +. (4) µ FB F Te momentum at te blade older M is given by: M r R. (5) F We assume tat te forces and torques acting on te blade older are equal in magnitude and wit opposite sign to tose acting on te blade. Tereby, te following balance applies R F R 0. (6) FA + F By imposing conservation of momentum at te rotation axis M r F r R M 0 (7) so tat A F r F r R + M. (8) F We aim to relate te contact force at te blade tip wit te force imparted by te pressure ose. Based on te follow relation F + F F + F (9) N µ and substituting and combining eq. () and () we obtain tat te reaction force at te blade tip is R ( F vˆ + v ) ˆ F N µfn (0) and te related momentum ( F vˆ + v ) ˆ N M rc N µf. () Substituting te latter two equations in te balance of momentum (eq. 8) we finally obtain a relation between te force imparted by te pressure ose F and te normal load at te blade tip F N ( r + r ) ( F vˆ + v ) ˆ N r F c N µf. () So far, bending of te doctor blade upon contact wit te roll as been neglected in te mecanical model. In te conventional leading-order approximation, te deflection angle at te blade tip is given by F l α arctan. () EI According to tis angle, te displacement of te blade is r r 0 + x (4) wit x F l Fˆ. (5) EI Given a rotation γ, te relation between F N and F becomes: ( r Rx( γ )) F (( r + r 0 + x) R x( γ )) ( FN v ( γ ) + µfn v ( γ )), (6) wic according to eq. (9) can be updated for any wear rate. In a refined mecanical model, te elasticity of te pressure ose can be included readily. Tus, te pressure can be calculated from te required force component F according to p π 4F d L E W + + L d EW ( ) d π L (7)

4. RESULTS Te proposed model is applied in order to calculate te film tickness, wen typical operational parameters are applied (Table ). Table. Reference case parameters Parameter l d E F N Value 40 mm.5 mm 5 MPa 00 N/m Te results sow te gap eigt at fluid intake and at fluid outtake for tree different contact lengts L, namely 0.,.0 and.0 mm (Figure 4). For a rater new doctor blade (L 0. mm), te value of 0 and is small and a satisfactory cleaning performance can be expected. As te blade tip wears off, te cleaning performance degrades and for a severely worn doctor blade tip, te water film bot at flow intake and outtake becomes ticker. Figure 5. Gap eigt at fluid intake 0 and at fluid outtake for a doctor blade wit a free lengt of 0 mm. Instead of sortening te blade lengt, an even better cleaning performance can be acieved by increasing te blade tickness d up to 4.5 mm. Figure 6. Gap eigt at fluid intake 0 and at fluid outtake for a doctor blade wit a tickness of 4.5 mm. Figure 4. Gap eigt at fluid intake 0 and at fluid outtake for te parameters sown in Table. If te free lengt of te doctor blade l is sortened up to 0 mm, te model predicts an improvement of te cleaning performance for all values of L. 5. EXPERIMENTAL Dry and water-lubricated tribological model tests were performed on a pin-on-disc tribometer. Te pin was replaced by doctor blade samples wit a widt of 8 mm. During te tests, te blade was pressed against a 00r6 steel disc wit a constant normal load F N of 5 N, a blade angle β of 8 and a sliding velocity of 6 m/s. During te tests, wear was measured as te vertical displacement of te blade. After te test, te Arcard wear coefficient k d was calculated from te measured wear rate. Te results sow iger wear rates for glass fiber reinforced plastic blades wen compared to carbon fiber reinforced

polymer blades (Figure 7). In bot cases, wear was more severe under dry contact conditions. to predict te contact lengt as a function of time, provided tat te Arcard wear coefficient is obtained using independent experiments under similar wear conditions. 6. ONLUSIONS A matematical model for describing te cleaning performance and te lifetime of te doctor blade press roll tribosystem was presented. Figure 7. Arcard wear coefficient for different doctor blade materials An additional test was performed using a scale paper mill component test, wic was designed to performed tests on doctor blades wit a widt of 90 mm. Dry contact conditions were selected in order to accelerate wear and keep te duration of te component test witin a reasonable time. Te blades slid against a mill wit 0.7 m radius wit a sliding velocity of 0 m/s and a contact force of 95 N. Te latter value was set in order to ave te same contact force per unit lengt as in te pin-ondisc tests. Te blade angle was set in tis case to 0. Te test was interrupted regularly in order to measure te contact lengt of te doctor blade. Te results obtained are sown in Figure 8. Figure 8. Evolution of te contact lengt L as a function of time. Te dots are measured via a component test and te line is te prediction of te model. Te dots are te measurements obtained in te component test, wereas te line is te prediction of te geometrical wear model using te Arcard wear coefficient measured wit te pin-on-disc test. Tis result sows te feasibility of te model Te approac relied on a recently developed ydrodynamic model, wic links in a nondimensional group te most relevant operational parameters during wet pressing. anges of te blade tip due to wear were taken into account by a simple wear model based on Arcard. A mecanical model linked te contact force at te blade tip wit te force of te pressure ose, wic is te actual parameter controlled by operators. Te model allows a systematic analysis of te cleaning performance of te doctor blade as a function of te process parameters. Ticker and sorter blades reduce te gap eigt between te blade tip and te press roll, consequently improving te cleaning performance. Te building blocks of te model can be used independently in order to focus on particular penomena. For instance, te geometrical wear model was successfully applied to predict te result of a component test, by fitting te Arcard wear coefficient using model tests. 7. AKNOWLEDGEMENTS Te autors express teir tanks to te Austrian Researc Promotion Agency (FFG) for financial support. 8. REFERENES [] Papiermacer Tascenbuc. Auflage, Dr. Kurt aefner-verlag, eidelberg, 999, (in German). [] Rodríguez Ripoll, M., Sceicl, Jakab, B., Franek, F., Modelling te doctor blade-roller

tribosystem for improving te cleaning performance during paper production, Trib. Letters, 0, (under review). [] ameron, A., Basic Lubrication Teory ( nd Ed.). Jon Wiley & Sons, 976. [4] Arcard, J.F., irst, W.: Te wear of metals under unlubricated conditions, Royal Soc. of London Proc. Series A 6, 956, 97 40.