Narasimha Marakala 1, Appu Kuttan K K 2 and Ravikiran Kadoli 2

Similar documents
RADIATION THERMOMETRY OF METAL IN HIGH TEMPERATURE FURNACE

MANUFACTURING TOLERANCES AS A CAUSE FOR AUDIBLE NOISE OF INDUCTION MOTORS

Modeling of the Fluid Solid Interaction during Seismic Event

CHAPTER 4 COMPARISON OF PUSH-OUT TEST RESULTS WITH EXISTING STRENGTH PREDICTION METHODS

Analysis the Transient Process of Wind Power Resources when there are Voltage Sags in Distribution Grid

DESIGN AND ANALYSIS OF FLEXIBLE STRUCTURES WITH PARAMETRIC UNCERTAINTIES FOR ACTIVE VIBRATION CONTROL USING ANSYS

Mixed convection in a vertical porous channel: heat transfer enhancement

Figure 1 Siemens PSSE Web Site

Fin shape optimization in tube heat exchangers by means of CFD program

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

Analysis of Active Suspension Systems with Hydraulic Actuators

The Multilayer Impedance Pump Model

Research on sound insulation of multiple-layer structure with porous material and air-layer

Investigation on thermal effects in high power Yb-doped fiber laser. Usha Chakravarty,* B. N. Upadhyaya, A. Kuruvilla, S. M. Oak

Natural Convection of Water-Based CuO Nanofluid Between Concentric Cylinders

Fluid-structure coupling analysis and simulation of viscosity effect. on Coriolis mass flowmeter

Finite Element Analysis of Ferrofluid Cooling of Heat Generating Devices

ME 3560 Fluid Mechanics

Stress Intensity Factors In Two Bonded Elastic Layers Containing Crack Perpendicular on the Interface with Different Elastic Properties

THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER

Volume IV, Issue IX, September 2015 IJLTEMAS ISSN

Chapter 9: Controller design. Controller design. Controller design

A Numerical Study on Mixed Convection of Water Based Cuo Nanofluids in A Lid-Driven Square Enclosure: Effects of Viscosity Models

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization

Design of Two-Channel Low-Delay FIR Filter Banks Using Constrained Optimization

POWER FREQUENCY CONTROL

The Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor

STRAIN LIMITS FOR PLASTIC HINGE REGIONS OF CONCRETE REINFORCED COLUMNS

Heat Transfer Modeling using ANSYS FLUENT

CHAPTER 5. The Operational Amplifier 1

Mixed Convection Characteristics of Ethylene Glycol and Water Mixture Based Al2O3 Nanofluids

Simulation of the Macroscopic Heat Transfer and Flow Behaviours in Microchannel Heat Sinks using Porous Media Approximation

Numerical Investigation of Heat Transfer and Fluid Flow around the Rectangular Flat Plane Confined by a Cylinder under Pulsating Flow

MINIATURIZATION OF INERTIAL MAGNETOELECTRIC VIBRATION VELOCITY SENSOR WITH HIGH FREQUENCY RESPONSE AT LOW FREQUENCY BAND

Fluid Detection in Tight Gas Sand from the Seismic Data

Solved problems 4 th exercise

Heat and mass transfer effects on nanofluid past a horizontally inclined plate

A NEW YIELD CRITERION FOR ORTHOTROPIC SHEET METALS UNDER PLANE-STRESS CONDITIONS

EFFECTS OF NAVIER SLIP AND WALL PERMEABILITY ON ENTROPY GENERATION IN UNSTEADY GENERALIZED COUETTE FLOW OF NANOFLUIDS WITH CONVECTIVE COOLING

SHEAR MECHANISM AND CAPACITY CALCULATION OF STEEL REINFORCED CONCRETE SPECIAL-SHAPED COLUMNS

Critical problem in EB-FRP technique

Operational transconductance amplifier based voltage-mode universal filter

Changes in Fresh and Saltwater Movement in a Coastal Aquifer by Land Surface Alteration

Resonance investigation of pump-turbine during startup process

Theoretical Research on the Calculation of Bending Moment at Arc Vault in Steel Corrugated Plate

ORIGINAL ARTICLE Electron Mobility in InP at Low Electric Field Application

Recent updates in fatigue analysis with LS-DYNA

EP225 Note No. 5 Mechanical Waves

1. Intensity of Periodic Sound Waves 2. The Doppler Effect

USE OF INTERNET TO DO EXPERIMENTS IN DYNAMICS AND CONTROL FROM ZACATECAS MEXICO IN THE LABORATORY OF THE UNIVERSITY OF TENNESSEE AT CHATANOOGAA.

Numerical Simulation on Unsteady Heat Transfer of a Sphere

Physics 11 HW #9 Solutions

MAE320-HW7A. 1b). The entropy of an isolated system increases during a process. A). sometimes B). always C). never D).

ROOT LOCUS. Poles and Zeros

Bandwidth expansion of a pressure control system for pneumatic anti-vibration apparatuses in presence of dead time

Microwave magnetoelectric fields: An analytical study of topological characteristics

Application of Extended Scaling Law to the Surface Tension of Fluids of Wide Range of Molecular Shapes

Mater. Res. Soc. Symp. Proc. Vol Materials Research Society. Pattern Evolution of Self-Assembled Quantum Dots Under Biaxial Stresses

Simulation of Wound Rotor Synchronous Machine under Voltage Sags

A Mathematical Model of Two Layered Pulsatile Blood Flow through Stenosed Arteries

GEOMETRIC OPTIMISATION OF CONJUGATE COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES

LOAD FREQUENCY CONTROL OF MULTI AREA INTERCONNECTED SYSTEM WITH TCPS AND DIVERSE SOURCES OF POWER GENERATION

Systems Analysis. Prof. Cesar de Prada ISA-UVA

Estimating floor acceleration in nonlinear multi-story moment-resisting frames

Massachusetts Institute of Technology Dynamics and Control II

ECSE 4440 Control System Engineering. Project 1. Controller Design of a Second Order System TA

Modes of Flow in Pneumatic Conveying Systems

Research on DC resistivity for an arbitrarily anisotropic earth using circular scanning measurement

ID-1204 STRUCTURAL-DYNAMIC AND ACOUSTIC DESIGN OF ANISOTROPIC MULTILAYERED COMPOSITE STRUCTURES

Load-deformation Analysis of a Pile in Expansive Soil upon Infiltration

Uncertainty Analysis in Potential Transformer Calibration Using a High Voltage Capacitance Bridge

OVERSHOOT FREE PI CONTROLLER TUNING BASED ON POLE ASSIGNMENT

Active Disturbance Rejection Control for an Electro-statically Actuated MEMS Device

Some Approaches to the Analysis of a Group of Repeated Measurements Experiment on Mahogany Tree with Heteroscedustic Model

55:041 Electronic Circuits

Interaction Diagram - Tied Reinforced Concrete Column (Using CSA A )

Analysis of Friction-Induced Vibration Leading to Eek Noise in a Dry Friction Clutch. Abstract

Steel Fiber-Reinforced Concrete Panels in Shear: Analysis and Modeling

Lecture #5: Introduction to Continuum Mechanics Three-dimensional Rate-independent Plasticity. by Dirk Mohr

Sound waves. Content. Chapter 21. objectives. objectives. When we use Sound Waves. What are sound waves? How they work.

Representation of a Group of Three-phase Induction Motors Using Per Unit Aggregation Model A.Kunakorn and T.Banyatnopparat

Internal Model Control

Efficiency Optimal of Inductive Power Transfer System Using the Genetic Algorithms Jikun Zhou *, Rong Zhang, Yi Zhang

Simple Observer Based Synchronization of Lorenz System with Parametric Uncertainty

UNITS FOR THERMOMECHANICS

OBSERVER-BASED REDUCED ORDER CONTROLLER DESIGN FOR THE STABILIZATION OF LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

Laser Doppler Velocimetry (LDV) Part - 01

One Dimensional Modeling of the Shape Memory Effect

CHAPTER 14 SIGNAL GENERATORS AND WAVEFORM-SHAPING CIRCUITS

Analysis of PI controller by model based tuning method in Real Time Level Control of Conical Tank System using block box modeling

Time Response of Nitinol Ribbons

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS

Uncertainty in long-term behavior and buckling of concrete-filled steel tubular. columns

Reliability of Cu Pillar Bumps for Flip-Chip Packages with Ultra low-k Dielectrics

Salina Aktar 1, Md.Abdul. Alim 2, Ali J. Chamkha 3

Reliability Analysis of Embedded System with Different Modes of Failure Emphasizing Reboot Delay

3D Metamaterial Based on a Regular Array of Resonant Dielectric Inclusions

Thermal Resistance Measurements and Thermal Transient Analysis of Power Chip Slug-Up and Slug-Down Mounted on HDI Substrate

Lecture 10. Erbium-doped fiber amplifier (EDFA) Raman amplifiers Have replaced semiconductor optical amplifiers in the course

NONLINEAR CONTROLLER DESIGN FOR A SHELL AND TUBE HEAT EXCHANGER AN EXPERIMENTATION APPROACH

Transcription:

IOSR Journal o Mechanical and Civil Engineering (IOSR-JMCE) ISSN: 78-1684, PP: 63-68 www.iorjournal.org Exerimental and Theoretical Invetigation o Combined Eect o Fluid and Thermal Induced Vibration on Vertical Thin Slender Tube Naraimha Marakala 1, Au uttan and Ravikiran adoli 1( Deartment o Mechanical Engineering, NMAM Intitute o Technology, Nitte-574110.) ( Deartment o Mechanical Engineering, National Intitute o Technology arnataka Srinivanagar-57505, Mangalore, India). Abtract: The aer reent an exerimental and theoretical analyi on the tudy o dynamic reone o a thin lender cantilever ie conveying air at dierent reure and. A dynamic characteritic o a ie line conveying internal luid exerience mechanical load due to inertia eect o luid, corioli orce, luid rotation kinetic orce due to luid low velocity, dynamic load due to inertia eect and thermal load due to hot air on the ie. One dimenional beam inite element i ued or invetigating the dynamic behavior o the thin lender ie with aroriate numerical comutational method or analyi. The undamental natural requency o vibration i correlated with exerimental and numerical method. eyword- Fluid and thermal induced vibration, inite element method, natural requency, lender tube I Introduction Thin walled tube are widely ound in variou tructural engineering alication. Static and dynamic reone o ie tructure to luid load and thermal load i very imortant or ae deign and oeration. Fluid lowing through a ie can imoe reure on the ie wall which delect the ie, lead to luid induced vibration. Non uniorm heat traner due to dierence in convective heat traner over the urace o the tube lead to thermal induced vibration. Many reearcher have tudied thermal and luid induced vibration. Studie have been conducted to ae the relationhi between luid low velocity, reure o the luid and it inluence on the ree vibration. Among the aer available in the literature, Paidoui et al. [1] and Rouelet and Herrmann [] reented a aer on analytical tudie to invetigate dynamic tability o a lender and long tube conveying luid. Numerical invetigation on the vibration behavior o the ie conveying luid wa reented by Lin and Tai [3] and Lee and Oh [4]. Bar-Avi [5] ha conducted an exerimental and analytical tudy on the dynamic reone o a vertical cantilever ie conveying luid which i ubjected to ea wave orce. uier and Metrikine [6] invetigated dynamic tability o a ubmerged vertical cantilever ie. Paidoui and Semler [7] examined the lanar dynamic o a luid conveying cantilever ie with a mall ma attached at the ree end by theoretically and exerimentally. Sinha et al. [8] tudied the reult o a modal exeriment on an oen ended cantilever ie conveying luid and their numerical imulation uing inite element method. Zou et al. [9] tudied analytically luid induced vibration o comoite natural ga ieline and invetigated the inluence o luid velocity, internal reure and initial tenion on ieline vibration requency. Tran et al. [10] analyed the tructural intenity attern o thin ientroic and laminated comoite late which are ubjected to thermally induced vibration. idawa-ukla [11] ha roved the requency o the thermally induced beam vibration i a multile o harmonic motion requency o the heat ource and then reonance can occur in the ytem. Hong [1] analyed the thermally induced vibration o a thermal leeve uing comutational Generalized Dierential Quadrature(GDQ) method or calculating the natural requency, dilacement and thermal tree. ukla [13] analyed the ditribution and tranvere vibration o beam induced by uing the roerty o the Green unction. Blandino and Thornton [14] carried out the detailed tudy o a thermally induced vibration caued by internal heating. From the literature review it i oberved that not much work ha been reorted on the reone o tructure due to combined loading o thermal, luid, luid inertia and tructure inertia. Hence the reent work ocuing on thi toic et out the ollowing objective. To viualize the roblem, an exerimental and theoretical analyi i contemlated. (a) An exerimental etu i abricated in the laboratory to imulate luid a well thermal load or the thin tube. (b) Exeriment were conducted or variou luid and reure (c)obervation on the tranient and teady tate reone o thin ie due to thee load will be analyzed. Second International Conerence on Emerging Trend in Engineering (SICETE) Dr.J.J.Magdum College o Engineering, Jayingur 63 Page

Exerimental and Theoretical Invetigation o Combined Eect o Fluid and Thermal Induced (d)exerimental and numerically obtained natural requency i comared. II Exerimental etu An exerimental etu i abricated in the laboratory to demontrate the luid thermal interaction in a thin lender tube. The etu conit o two narrow tube o diameter 3 mm. The tube are connected through an oriice aembly. The oriice aembly roduce an additional load due to luid exanion. The oriice ize can be changed to obtain dierent hole diameter. Exeriment were conducted or 1 mm and 1.5 mm oriice hole ize. The one end o the irt tube i or the entry o the comreed air. The other end o thi tube i connected to another tube through oriice aembly. The other end o the econd tube i exoed to atmohere. Thee two tube aemblie are in uch away that one end i ixed and other end i ree with oriice aembly. The etu i mounted inide a tet tand. The tet tand act a a mean to hold the two end o the tube a hown in Figure 1. The uly reure i controlled by the reure regulator. The air i aed through tainle teel heating box o cylindrical hae o diameter 00 mm and the arrangement i uitable or generating reurized air at required. The tainle teel heating box contain 00 watt heating coil connected through 0 V, 5 am AC ource with a voltage regulator. The voltage regulator i required to vary the ower uly to the heating coil in order to vary the uly air. A the air low, the heat rom the comreed air i lot by conduction through the tube and by convection rom the tube urace to atmohere. A thermocoule wa mounted inide the tet box to meaure the ambient which wa recorded by mean o digital thermometer. The data acquiition rogram in LabVIEW through a cloed loo controlled the acquiition o dilacement data rom the tet ecimen. Due to conveying o comreed air reurized air, luid load on the tructure and an additional load due to udden exanion o the air through an oriice i reent. The air enter into one chamber o the udden exanion aembly and low through an oriice late having mall hole. Thu, low reult in udden exanion o air and hence additional load at the tube ti. From the econd chamber, air exit through the econd tube to atmohere. To obtain the tranient and teady tate reone o the tube, a iezoelectric accelerometer AD11 N605 with enitivity 10 mv/g i ued. In order to acquire the ignal NIPXI 1050 chai containing NIPXI 447 DAQ with 4 bit having 8 analog inut channel and analog outut channel i ued a hown in Figure 1. Then, intrument i connected to CPU and monitor to dilay o the outut. The rogramming wa done to acquire reone ignal through LabVIEW 8.0. III Finite element ormulation Conider a hot luid element and ie element with variou orce and moment acting on it due to low o hot luid through ie a hown in igure 3 and 4. The luid element i ubjected to reure orce, reaction orce o the ie on the luid normal to the luid element F, and tangential to it i q S, Second International Conerence on Emerging Trend in Engineering (SICETE) Dr.J.J.Magdum College o Engineering, Jayingur 64 Page

Exerimental and Theoretical Invetigation o Combined Eect o Fluid and Thermal Induced aociated with the wall-hear treq and gravity orce m g. The ie element i ubjected to the longitudinal tenion T o, tranvere hear orce Q T, bending moment M and daming due to riction with the w urrounding luid a B. Heat traner rom hot luid to ie i h A t T T, heat traner rom hot ie T T, where h i heat to atmohere at comreed ide i 1 and at elongated ide i traner coeicient, T i o hot luid, T i o ie and T i ambient. The thermal moment due to change in rate o heat traner i M T. The baic equation o motion conveying luid [15] can be modiied to include the thermal load [14] and hence the governing equation or ie conveying hot luid could be derived a: 4 E I w m 4 U T A1 w m w ( ) w w 0, U m m E T x x x t t x (1) where E i the modulu o elaticity o the ie, I the area moment o inertia o the ie, A the luid low area o the ie, the denity o luid, m i the luid ma er unit length, and m i the ie ma er unit length, U i the low velocity o luid i the average reure,t i the alied tenion, i thermal coeicient o exanion o material and T i. The lat term, thermal moment w MT E T conider deormation due to thermal load (Blandino and Thornton, 001) which include x conduction and convection through the tube. Uing the Galerkin aroach or the equation (1), the inite element equation o motion can be obtained a: t th ( e) ( e) ( e) ( e) ( e) ( e) ( e) ( e) e ( e) M w G w w 0, () where, Stine matrix 1 6l 1 6l e EI 6l 4l 6l l 3 l 1 6l 1 6l 6l l 6l 4l Fluid element matrix 3 6 3l 3 6 3l e MU 3l 4l 3l l 3 6 3l 3 6 3l 3l l 3l 4l Pie tenion matrix 3 6 3l 3 6 3l e T 3l 4l 3l l t 3 6 3l 3 6 3l 3l l 3l 4l Fluid reure matrix Second International Conerence on Emerging Trend in Engineering (SICETE) Dr.J.J.Magdum College o Engineering, Jayingur 65 Page

Exerimental and Theoretical Invetigation o Combined Eect o Fluid and Thermal Induced 3 6 3l 3 6 3l e A1 3l 4l 3l l 3 6 3l 3 6 3l 3l l 3l 4l Fluid daming matrix 3 0 6l 3 0 6l e MU 6l 0 6l l G 30 3 0 6l 3 0 6l 6l l 6l 0 Element ma matrix 1 5 6 l 5 4 1 3l e M ml l 4l 1 3l 3l M 40 5 4 1 3l 1 5 6 l 1 3l 3l l 4l Element thermal matrix 3 6 3l 3 6 3l e 3l 4l 3l l E T th 3 6 3l 3 6 3l 3l l 3l 4l Aemble the element matrice o tine, daming, and ma matrice and rereent in a tate ace orm (Zou et.al, 005) a:. Θ Ψ where 0 0 M 0,, 0 G Eigenvalue are comuted uing matlab or the equation (3), e e e e e. t th (3) M T w B t M T Q T T 0 qs m g F o Q QT Fig. Hot luid element 1 T T M M T o M MT T F A T qs A( ) m g Fig.3 Pie element T Second International Conerence on Emerging Trend in Engineering (SICETE) Dr.J.J.Magdum College o Engineering, Jayingur 66 Page

Amlitude in meter Magnitude Exerimental and Theoretical Invetigation o Combined Eect o Fluid and Thermal Induced IV Reult and dicuion Figure 4 and 5 indicate the tranient reone o the tube and undamental natural requency at uly air reure o bar and at ambient 31 0 C. Exeriment were reeated or 1 bar, 3 bar, 4 bar and 5 bar with o air being 58 0 C and 78 0 C. Table 1 how the exerimental dilacement o the tube or variou reure and. Table how the exerimental and theoretical undamental requencie o the ie conveying hot air with variou reure and. From the reone analyi it i oberved that dilacement varie rom. mm to 5.3 mm a the reure increae rom bar to 5 bar at ambient o air lowing internally. When the o the air being conveyed through the tube i at 78 0 C and 1 bar reure the dilacement o the tube variou rom. mm to 4.0 mm. The requency o vibration or a tube o 800mm length increae rom.614 Hz to 6.17 Hz or an increae in reure rom 1 bar to 5 bar. Similarly requency o vibration varie rom.614 Hz to 6.317 Hz or an increae in air rom ambient to 78 0 C. Thi how that a the reure increae the reone become more ocillatory and take more time to attain teady tate. Frequency o tube increae with reure and it decreae with. Thi may be due to the internal reure o air being tranlated into a tenile load along the ie, hence the natural requency increae. With the hot air being conveyed, there would likely be the otening eect that can caue the requency to decreae. 40000 35000 0.0 30000 5000 0.00 0000-0.0 15000 10000-0.04 5000 0-1.0-0.5 0.0 0.5 1.0 1.5.0.5 3.0 Time in econd Fig.4 Time reone 0.0 0.5 1.0 1.5.0.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 Frequency Fig. 5 Frequency reone Table 1 Exerimental Dilacement o tube or variou reure and Preure in bar Ambient (31 o C) 58 0 C 78 0 C 1.3.1 4.0 3.1 4.0 4.7 3 3.8 4.7 5.4 4 4.5 5.6 6. 5 5.3 6.5 7.1 Preure in bar Table Exerimental and theoretical natural requency or variou and reure Fundamental natural requency (Hz) Exerimental reult Theoretical reult Ambient 31 0 C o the air 58 0 C o the air 78 0 C Ambient 31 0 C o the air 58 0 C 1.614.481.134.511.48.156 3.834 3.356 3.4 3.656 3.51 3.36 3 4.91 3.935 3.434 4.41 3.91 3.351 4 5.811 4.365 3.564 5.1 4.41 3.641 5 6.317 5.83 4.458 6.516 5.974 4.59 o the air 78 0 C Second International Conerence on Emerging Trend in Engineering (SICETE) Dr.J.J.Magdum College o Engineering, Jayingur 67 Page

Exerimental and Theoretical Invetigation o Combined Eect o Fluid and Thermal Induced V Concluion Table how the concluding remark obtained exerimentally and theoretically on the undamental requencie o the ie conveying hot air. It i oberved that the requency o vibration increae with the increae in reure and decreae with increae in. The increae in reure increae the velocity o the luid low and reduce the daming eect. Temerature ha an eect on dilacement a well a requency, ince thermal contraction and exanion increae due to high heat traner rate at high velocity. A the o the luid increae, the requency decreae due to otening eect o tube. REFERENCES [1] Paidoui M.P.,Luu T.P. and Laithier B.E., Dynamic o inite length tubular beam conveying luid, Journal o Sound and Vibration, 106, 1986, 311-331. [] Rouelt J. and Herrmann G., Dynamic behavior o continuou cantilever ie conveying luid near the critical velocity, Journal o Alied Mechanic, 48, 1981, 943-947. [3] Yih-Hwang Lin and Yau-un Tai, Nonlinear vibration o thimohenko ie conveying luid, Journal o Solid and Structure, 34, 1997,945-96. [4] Uik Lee and Hyuckjin Oh,, The ectral element model or ieline conveying internal teady low, Journal o Engineering Structure, 5, 003, 1045-1055. [5] Patrick Bar-Avi,, Dynamic reone o rier conveying luid, Journal o Ohore Mechanic and Arctic Engineering, 1, 000,188-193. [6] uier G.L. and Metrikine A.V., Exerimental invetigation o dynamic tability o a cantilever ie airating luid, Journal o Fluid and Structure, 4, 008, 541-558. [7] Paidoui M.P. and Semler C., Non-linear dynamic o a luid conveying cantilever ie with a mall ma attached at the ree end, International Journal o Non-linear Mechanic, 33,1998, 15-3. [8] Sinha J.., Singh S. and Rama Rao A., Finite element imulation o dynamic behavior o an oen ended cantilever ie conveying luid, Journal o Sound and Vibration, 40(1), 001, 189-194. [9] Zou G.P.,Cheraghi N.,and Taheri F., Fluid induced vibration o comoite natural ga ieline, Journal o Solid and Structure, 4, 005, 153-168. [10] Tran T.Q.N, Lee H.P. and Lim S.S., Structural intenity analyi o thin laminated comoite late ubjected thermally induced vibration, Journal o Comoite Structure, 78, 007, 70-83. [11] idawa-ukla, Vibration o a beam induced by harmonic motion o a heat ource, Journal o Sound and Vibration, 05(), 1997, 13-. [1] Hong C. C., Liao H. W., Lee L. T., e J. B. and Jane. C., Thermally induced vibration o a thermal leeve with the GDQ method, International Journal o Mechanical Science, 47, 005, 1789 1806. [13] idawa-ukla J., Alication o green unction to the roblem o thermally induced vibration o a beam, Journal o Sound and Vibration, 6, 003, 865-878. [14 Blandino J. R and Thornton E. A., Thermally induced vibration o an internally heated beam, Journal o vibration and acoutic, 13, 001, 67-74. [15] PaidouiM.P., Fluid-Structure interaction ( Academic re, Caliornia,USA, 1998). Second International Conerence on Emerging Trend in Engineering (SICETE) Dr.J.J.Magdum College o Engineering, Jayingur 68 Page