MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava

Similar documents
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Application to Plane (rigid) frame structure

One Dimensional Axial Deformations

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

Indeterminate pin-jointed frames (trusses)

DUE: WEDS FEB 21ST 2018

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

FUZZY FINITE ELEMENT METHOD

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS

Increase Decrease Remain the Same (Circle one) (2 pts)

PLATE BENDING ELEMENTS

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS

Finite Element Modelling of truss/cable structures

Transactions of the VŠB Technical University of Ostrava, Mechanical Series. article No. 1907

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

LAB 4: Modulus of elasticity

Second Order Analysis

Module 3: Element Properties Lecture 1: Natural Coordinates

3D beam finite element including nonuniform torsion

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

MECHANICS OF MATERIALS

Numerical Heat and Mass Transfer

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

Grid Generation around a Cylinder by Complex Potential Functions

FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

I have not received unauthorized aid in the completion of this exam.

Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P.

Plate Theories for Classical and Laminated plates Weak Formulation and Element Calculations

In this section is given an overview of the common elasticity models.

Constitutive Modelling of Superplastic AA-5083

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

New Method for Solving Poisson Equation. on Irregular Domains

Lecture 8 Modal Analysis

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

MEMBRANE ELEMENT WITH NORMAL ROTATIONS

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Interconnect Modeling

CHAPTER 9 CONCLUSIONS

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES

Numerical Nonlinear Analysis with the Boundary Element Method

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION

A Timoshenko Piezoelectric Beam Finite Element with Consistent Performance Irrespective of Geometric and Material Configurations

Geometrically exact multi-layer beams with a rigid interconnection

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Description of the Force Method Procedure. Indeterminate Analysis Force Method 1. Force Method con t. Force Method con t

NUMERICAL ANALYSIS AND EXPERIMENTAL MEASUREMENT OF EIGENFREQUENCIES OF ACSR POWER LINES

Computer Based Porosity Design by Multi Phase Topology Optimization

Nordic Insulation Symposium - Nord-IS 13 - Trondheim, Norway, June 9-12, 2013

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Rigid body simulation

Instituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245

A Mechanics-Based Approach for Determining Deflections of Stacked Multi-Storey Wood-Based Shear Walls

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

Chapter - 2. Distribution System Power Flow Analysis

HEAT TRANSFER THROUGH ANNULAR COMPOSITE FINS

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

One-sided finite-difference approximations suitable for use with Richardson extrapolation

Module 11 Design of Joints for Special Loading. Version 2 ME, IIT Kharagpur

Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model

Plan: Fuselages can. multideck

Buckling analysis of single-layered FG nanoplates on elastic substrate with uneven porosities and various boundary conditions

Finite Difference Method

CIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56

Effect of loading frequency on the settlement of granular layer

Tensor Smooth Length for SPH Modelling of High Speed Impact

Army Ants Tunneling for Classical Simulations

Experimental Study on Ultimate Strength of Flexural-Failure-Type RC Beams under Impact Loading

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

Please initial the statement below to show that you have read it

829. An adaptive method for inertia force identification in cantilever under moving mass

Measurement of the shear modulus of wood by asymmetric four-point bending tests

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

Physics 181. Particle Systems

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

THE DEFORMATION OF CYLINDRICAL SHELLS SUBJECTED TO RADIAL LOADS

GEOMETRIC SOLUTION IN PROGRESSIVE COLLAPSE ANALYSIS OF HULL GIRDER

4DVAR, according to the name, is a four-dimensional variational method.

PART 8. Partial Differential Equations PDEs

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp.

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

ANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP

SIMPLE LINEAR REGRESSION

General displacement arch-cantilever element method for stress analysis of arch dam

Torsion Stiffness of Thin-walled Steel Beams with Web Holes

A Hybrid Variational Iteration Method for Blasius Equation

Effects of internal=external pressure on the global buckling of pipelines

6 Supplementary Materials

Modeling of Dynamic Systems

THE EFFECT OF BEAM TO COLUMN CONNECTION IN ARC PORTAL FRAME

Inductance Calculation for Conductors of Arbitrary Shape

Transcription:

ODING OF ASTO-STATICS OF POW INS BY NW COPOSIT BA FINIT NT urín Justín 1 rabovský Jura 1 Gogola oman 1 utš Vladmír 1 Paulech Jura 1 1 Insttute of Automotve echatroncs FI STU n Bratslava Ilkovčova 3 812 19 Bratslava -mal: ustn.murn@stuba.sk eceved 7 a 215; accepted 15 a 215 1. Introducton The power lne s from mechancal pont of vew a 3D sstem. It can be loaded n longtudnal horontal and vertcal drecton. so the torsonal loadng s possble as well. But n the techncal calculatons t s mostl smplfed to the one dmensonal sstem. In the lterature e. g [1] an analtcal method s presented for the power lne elasto-statc analss n ts vertcal plane. The analtcal methods are not much effectve for the general spatal analss. The more effectve are the numercal methods over all the fnte element method. For the smple elasto-statc analss the geometrcall nonlnear lnk fnte element can be used that s able to anale the tensonal forces and stresses and the elongaton of the lne e.g. [2]. For dnamc analss the beam fnte element s preferable [3]. The heterogeneous cross-secton are of several constructon e.g. as shown n Fg.1 [4]. Because the materal of the power lne s nhomogeneous the homogenaton of materal propertes s needed. Therefore the aal fleural and torsonal stffness must be stated. For the mechancal analss also the sold fnte elements are avalable but modelng of the complcated geometr and heterogenet s a ver demandng procedure. In ths paper the results of elasto-statc analss of the sngle and bundle power lnes are presented. Fg. 1. Constructon of power lne cross-secton [4]. For modelng and smulaton of the problem a new 3D composte beam fnte element s used whch was developed at our nsttute [3]. The second order beam theor has been used for the fnte element stffness matr formulaton. The effectve tensonal fleural and torsonal stffness of chosen power lnes are consdered [3]. The results are calculated evaluated and compared wth the ones obtaned b the standard fnte element software [2]. 2. Composte fnte beam element equatons et us consder a 3D straght fnte beam element (Tmoshenko beam theor and Sant-Venant torson theor) of doubl smmetrc cross-secton Fg. 2. The nodal degrees 375

of freedom at node are: the dsplacements u v w n the local as drecton and the cross-sectonal area rotatons -. The degrees of freedom at the node are denoted n a smlar manner. The nternal forces at node are: the aal force N the transversal forces and the bendng moments and and the torson moment. The frst dervatve wth respect to of the relevant varable s denoted wth an apostrophe. Fg. 2: The local nternal varables and loads. Furthermore n n s the aal force dstrbuton q q and q q transversal and lateral force dstrbutons m m m m and m m are the are the dstrbuted moments A s the cross-sectonal area I and I are the second area moments I p I I s the area polar moment. The effectve homogened and longtudnall varng N N N stffness reads: A A s the aal stffness ( s the effectve elastct modulus for aal loadng) I I s the effectve elastct modulus for bendng about as ) I I s the fleural stffness about the as ( sm bendng about as ) GA Gk A G G s the effectve shear modulus and s the fleural stffness about the -as ( s the effectve elastct modulus for s the reduced shear stffness n drecton ( sm k s the average shear correcton factor n sm drecton [5]) GA G k A s the reduced shear stffness n drecton ( s the effectve shear modulus and sm k G G s the average shear correcton factor n drecton [5]) G I T s the effectve torsonal stffness ( G G s the torsonal elastct modulus and I T s the torson constant). Detaled dervaton of the effectve materal propertes s presented n [3] and [4]. stablshng of the local composte beam fnte element equatons s presented n [6]: 376

377 w v u w v u Y T Y S N N 1212 1111 11 911 99 812 88 77 612 68 66 511 59 55 41 44 311 39 35 33 212 28 26 22 17 11 (1) In (1) the terms contan the lnear and the lneared geometrc non-lnear stffness terms contanng the aal force effect on the fleural beam stffness. Shear correcton s accounted as well. The global stffness matr of the beam structures can be done b classcal methods. stablshng of the local and global stffness matrces as well as the whole soluton procedure were coded b the software ATATICA [7]. 3. Numercal eperments The smmetrc power lne marked as 456 (3+9 steel and 11+17 alumnum wres) whch s loaded n ts ntal state b the self-weght n -drecton has been consdered. Span of the power lne s = 3 m the mamal deflecton s m ma 3966 and the average nternal aal force s N II = 49.566 kn. The power lne s subsequentl loaded at ts mdpont b forces 1 F N and 1 F N. The dameter of the alumnum wres s d = 45 mm and the dameter of the steel wres s d = 28 mm. The effectve crosssectons of the power lne parts are: A = 7389 mm 2 A = 44532 mm 2 and the effectve cross-sectonal area of the power lne s A = 51921 mm 2. Fg. 3: power lne n ntal state consequentl loaded b concentrated forces. ateral propertes of the components are constant and ther values are: alumnum the elastct modulus GPa 7 the Posson s rato 32 the mass denst 3 27 kgm ; steel the elastct modulus GPa 21 the Posson s rato 28 the mass denst 3 785 kgm. For the elasto-statc analss of the sngle power lnes the followng effectve materal propertes have been used [3]: the effectve elastct modulus for aal loadng - Pa 8992376 N the fleural stffness - Pa 9714655

the effectve shear modulus - G G 345866 Pa the torsonal elastct modulus - G 28928 Pa The followng calculatons were done wth our 3D FG beam fnte element (NF). The same problem has been solved usng a commercal F program ANSYS wth 3 of beam fnte elements (BA188). The global dsplacements v [m] n ( drecton) and w [m] n ( drecton) at dstances of 5 1 and 15 m from the left lne end calculated wth ANSYS as well as wth the NF are presented n Table 1. The average relatve dfference [%] between dsplacements calculated b our NF-method and the ANSYS soluton has been evaluated as well. Tab. 1. Spatal dsplacements at the selected ponts (DSP) of the power lne. 5 1 15 DSP NF ANSYS [%] v 5-22229 -22255 12 w 5 499 497 4 v 1-35776 -3582 7 w 1 999 994 5 v 15-4.619-4627 2 v 15 1492 1491 7 The comparson of total dsplacement of the power lne calculated b our new approach and commercal F program ANSYS s shown n Fgure 4. Fg. 4: Total dsplacement of the power lne. The local nternal forces and moments n the power lne can be calculated from the local dsplacements at the fnte elements nodal ponts. Accordng to the gven loads and the low power lne stffness for bendng the relevant nternal force s the aal tensonal force N. Its mamal value s at the left and rght lne end: N ma 5324 kn. The average normal N stress s N/ A and the normal stran / at the -poston. Because of dfferent elastct modulus of the alumnum and the steel a dfferent normal stress n the materal components wll be arse. Poston of the crtcal cross secton (were the mamal normal force arses) s at the left and rght end of the power lne. In our case the mamal 75 Pa and for steel stress for the alumna s ma 51 226 53 Pa. ma 378

4. Conclusons New composte beam fnte element was used for the elasto-statc analss of the power lne. omogenaton of the heterogeneous materal propertes was made b the reference volume method (V). For comparson of the effectveness and accurac of the new fnte element the same problem was solved b the commercal F software ANSYS. A ver good agreement of both results has been obtaned. The most advantage of the new beam fnte element s that the stffness matr contans the effect of aal force and all the stffness (tensonal fleural and torsonal) are calculated b a consstent manner. Acknowledgement Ths paper has been supported b the Slovak Grant Agences VGA No. 1/453/15 and APVV-246-12 and VGA No. 1/228/14. Ths contrbuton s the result of the proect mplementaton: Industral esearch Centre for operatng lfetme of selected components of power plants (ITS: 26242281) supported b the esearch & Development Operatonal Programme funded b the DF. eferences: [1] Š. cko et. al.: Vonkaše elektrcké vedena enesans s.r.o. Bratslava ISBN 8-8942-35-9 (21). [2] ANSYS Swanson Analss Sstem Inc. 21 Johnson oad ouston PA 15342/13 USA. [3] J. urín J. rabovský. Gogola G. Gálk: odal analss of the power lnes b fnte element methods. Časops Vol. 2. NO 5/S 22-26 (214). [4] Ocelovo-hlníkové vodče ( lana) pro venkovní elektrcké vedení. Avalable on web: http://www.acword.c/portal/hromosvodov-materal/produkt/alfe-lana/alfelana-1644.htm [5] J. urín J. rabovský V. utš J. Paulech: Shear correcton functon dervaton for the FG beams. In: 2nd Internatonal Conference on ult-scale Computatonal ethods for sold and Fluds. 1. 6 12. 6.215 Saraevo Bosna and ercegovna (215). [6] J. urín. Amnbagha J. rabovský V. utš J. Pulech: 3D beam fnte element for elasto-statc analss of the FG structures. In: 8ICC 29.9. 2.1.215 Opata Croata (215). [7] S. Wolfram ATATICA 5 Wolfram research Inc. (23. 379