THE EQUATIONS OF THE SECOND ORDER LINEAR MODEL OF SURFACE GRIDS

Similar documents
4WCE * 5 * : GEO. Air Products and Chemicals, Inc., 2009

JOINT-STOCK COMPANY GIDROPRIVOD. RADIAL PISTON PUMPS OF VARIABLE DISPLACEMENT type 50 НРР

KbdKaz 500 layout tables

uninsta un in sta 9 weights & italics 5 numeral variations Full Cyrillic alphabet

«, 68, 55, 23. (, -, ).,,.,,. (workcamps).,. :.. 2

ON CERTAIN INEQUALITIES IN THE LINEAR SHELL PROBLEMS

2x x l. Module 3: Element Properties Lecture 4: Lagrange and Serendipity Elements

Infusion Pump CODAN ARGUS 717 / 718 V - Release Notes. Firmware V

Myriad Pro Light. Lining proportional. Latin capitals. Alphabetic. Oldstyle tabular. Oldstyle proportional. Superscript ⁰ ¹ ² ³ ⁴ ⁵ ⁶ ⁷ ⁸ ⁹,.

OFFER VALID FROM R. 15 COLORS TEXT DISPLAYS SERIES RGB12-K SERIES RGB16-K SERIES RGB20-K SERIES RGB25-K SERIES RGB30-K

Analysis of Continuous Beams in General

Operating Manual version 1.2

OFFER VALID FROM R. TEXT DISPLAYS SERIES A SERIES D SERIES K SERIES M

Module I. Unit 5. д е far ед е - near д м - next to. - е у ед е? railway station

, «Ruby»..,

Module 6: FEM for Plates and Shells Lecture 6: Finite Element Analysis of Shell

Town Development Fund

А. Љто нќша кђшка That s our cat

THE MATHEMATICAL MODEL OF AN OPERATOR IN A HUMAN MACHINE SYSTEMS. PROBLEMS AND SOLUTIONS

R s s f. m y s. SPH3UW Unit 7.3 Spherical Concave Mirrors Page 1 of 12. Notes

SHALLOW SPHERICAL LATTICE TYPE SHELL SUBJECTED TO NORMAL POINT LOAD FUNDAMENTAL SOLUTION

Hermite Splines in Lie Groups as Products of Geodesics

UCA Chart Help. Primary difference. Secondary Difference. Tertiary difference. Quarternary difference or no difference

NON-PROFIT ORGANIZATION CHARITY FUND

Title of your Paper AUTHOR NAME. 1 Introduction. 2 Main Settings

Optimization Methods: Integer Programming Integer Linear Programming 1. Module 7 Lecture Notes 1. Integer Linear Programming

National Women Commission

Accounting for the Use of Different Length Scale Factors in x, y and z Directions

NUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS

3D vector computer graphics

Electrical analysis of light-weight, triangular weave reflector antennas

Lecture #15 Lecture Notes

TEX Gyre: The New Font Project. Marrakech, November 9th 11th, Bogusław Jackowski, Janusz M. Nowacki, Jerzy B. Ludwichowski

ц э ц эр е э вс свэ эч эр э эвэ эч цэ е рээ рэмц э ч чс э э е е е э е ц е р э л в э э у эр це э в эр э е р э э э в э ес у ч р Б эш сэ э в р э ешшв р э

INDEPENDENT AUDITORIS REPORT. in the procedures

S.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION?

THE COCHINEAL FONT PACKAGE

Romanization rules for the Lemko (Ruthenian or Rusyn) language in Poland

R E N E W A B L E E N E R G Y D E V E L O P M E N T I N K A Z A K H S T A N MINISTRY OF ENERGY OF THE REPUBLIC OF KAZAKHSTAN

Design of Structure Optimization with APDL

Mathematics 256 a course in differential equations for engineering students

Finite Element Analysis of Rubber Sealing Ring Resilience Behavior Qu Jia 1,a, Chen Geng 1,b and Yang Yuwei 2,c

September 13, 1963 Letter from the worker of Donetsk metallurgy plant Nikolai Bychkov to Ukrainian Republican Committee of Peace Protection, Donetsk

Line geometry, according to the principles of Grassmann s theory of extensions. By E. Müller in Vienna.

A New Approach For the Ranking of Fuzzy Sets With Different Heights

Harmonic Coordinates for Character Articulation PIXAR

OCCASION DISCLAIMER FAIR USE POLICY CONTACT. Please contact for further information concerning UNIDO publications.

MEGAFLASHROM SCC+ SD

MANAGEMENT LETTER. Implemented by. Prepared by: (igdem KAHVECi Treasury Controller

Circuit Analysis I (ENGR 2405) Chapter 3 Method of Analysis Nodal(KCL) and Mesh(KVL)

Type-2 Fuzzy Non-uniform Rational B-spline Model with Type-2 Fuzzy Data

Solitary and Traveling Wave Solutions to a Model. of Long Range Diffusion Involving Flux with. Stability Analysis

-. (Mr R. Bellin),.: MW224r

Support Vector Machines

6.854 Advanced Algorithms Petar Maymounkov Problem Set 11 (November 23, 2005) With: Benjamin Rossman, Oren Weimann, and Pouya Kheradpour

The XCharter Font Package

Lecture 5: Probability Distributions. Random Variables

Lobachevsky State University of Nizhni Novgorod. Polyhedron. Quick Start Guide

Geometric Error Estimation

Classification / Regression Support Vector Machines

S1 Note. Basis functions.

Feature Reduction and Selection

Usage guidelines. About Google Book Search

Multiblock method for database generation in finite element programs

Р у к о в о д и т е л ь О О П д -р б и о л. н а у к Д. С. В о р о б ь е в 2017 г. А Г И С Т Е Р С К А Я Д И С С Е Р Т А Ц И Я

GSLM Operations Research II Fall 13/14

Introduction to Geometrical Optics - a 2D ray tracing Excel model for spherical mirrors - Part 2

12/2/2009. Announcements. Parametric / Non-parametric. Case-Based Reasoning. Nearest-Neighbor on Images. Nearest-Neighbor Classification

Barycentric Coordinates. From: Mean Value Coordinates for Closed Triangular Meshes by Ju et al.

Kinematics of pantograph masts

Structured Grid Generation Via Constraint on Displacement of Internal Nodes

Parallelism for Nested Loops with Non-uniform and Flow Dependences

i v v 6 i 2 i 3 v + (1) (2) (3) (4) (5) Substituting (4) and (5) into (3) (6) = 2 (7) (5) and (6) (8) (4) and (6) ˆ

P07303 Series Customer Display User Manual

Modeling of Fillets in Thin-walled Beams Using Shell/Plate and Beam Finite. Elements

APPLICATION OF MULTIVARIATE LOSS FUNCTION FOR ASSESSMENT OF THE QUALITY OF TECHNOLOGICAL PROCESS MANAGEMENT

LOOP ANALYSIS. The second systematic technique to determine all currents and voltages in a circuit

y and the total sum of

A NOTE ON FUZZY CLOSURE OF A FUZZY SET

Intra-Parametric Analysis of a Fuzzy MOLP

User Manual. June 2008 Revision 1.7. P07303 Series Customer Display

Assignment # 2. Farrukh Jabeen Algorithms 510 Assignment #2 Due Date: June 15, 2009.

Ц пл йъ п К л л я с д д ь с лп с И Ц щ пй л И й л мц Ц л ё чý и ф Ч иь Ерс ч и шь Ц ф с Й Г х Ч й М п д м цв фн пп йл ч йй м

Active Contours/Snakes

ON SOME ENTERTAINING APPLICATIONS OF THE CONCEPT OF SET IN COMPUTER SCIENCE COURSE

AUDIT OF THE FINANCIAL STATEMENT OF THE "EMERGENCY NATIONAL POVERTY TARGETING PROGRAM (PROJECT)" ENPTP

Hyperelastic Materials Characterization by Planar Tension Tests and Full-field Strain Measurement

User Manual. Revision v1.2 November 2009 EVO-RD1-VFD

User Manual September 2009 Revision 1.9. P07303 Series Customer Display

MIVOICE OFFICE 400 MITEL 6863 SIP / MITEL 6865 SIP

Visual Curvature. 1. Introduction. y C. IEEE Conf. on Computer Vision and Pattern Recognition (CVPR), June 2007

OPERATING RULES FOR CLEARING OF INTERNATIONAL PAYMENTS

Cell Count Method on a Network with SANET

SBOX-Max Operating Manual

ROBOT KINEMATICS. ME Robotics ME Robotics

Classifier Selection Based on Data Complexity Measures *

Topology Design using LS-TaSC Version 2 and LS-DYNA

FEATURE EXTRACTION. Dr. K.Vijayarekha. Associate Dean School of Electrical and Electronics Engineering SASTRA University, Thanjavur

On Some Entertaining Applications of the Concept of Set in Computer Science Course

Mixed-mode delamination in 2D layered beam finite elements

Transcription:

MECHANKA TEORETYCZNA STOSOWANA 2/3 2 (983) THE EQUATIONS OF THE SECOND ORDER LINEAR MODEL OF SURFACE GRIDS PIOTR W I Ś N I A K O W S KI Poltechnka Warszawska. Introducton The paper deals wth the method of fndng the governng equatons for a surface structure havng a form of a dense and regular grd made of bars. The lateral deformaton of elements of the system are taken nto account. It s assumed that the materal of the structure s elastc homogeneous and sotropc. The problem of statcs s analysed wthn the lnear theory. The numercal methods employed to solve the problems related to the consdered systems were based on dscret representaton of the structure (see among others [ 2 3]) and lead to a system of algebrac equatons wth a large number of unknowns. The dmensons of nodes ther deformablty and the lateral deformablty of structure's bars were not taken nto account. The applcaton of a contnuum model of a structure conssts n an approxmaton of the mult connected geometry of the system by a certan smply connected and contnuous model (see among others [4 5 6]). The advantage of the dscussed approach over the prevous one les n the fact that the analytcal methods can be employed. The negatve aspects are: a) consderable naccuracy of results for not suffcently dense grds b) the requred geometrcal symmetry of the structure. An nterestng dea of a contnuum model of such structures based on the concept of a contnuum wth nternal mcrostructure and hgher order nternal reactons s presented n [6]. In the present paper Cz. Woźnak's model wll be appled to obtan equatons of the second order theory. An energetc approach dfferent from the prevously consdered one whch wll be employed makes t possble to descrbe n the explct from all propertes of the contnuum model. As a specal case (n whch the hgher order effects are neglected) equatons of the frst approxmaton wll be obtaned. 2. Basc assumptons It s assumed that the structure consst of (homogeneously) deformable cubcod nodes connected by means of the prsmatc lnks of rectangular cross sectons (and subject to homogeneous deformaton n ther plane) and consttute the regular and orthogonal

22 P. WlŚ NIAKOWSKI surface grd made of bars (Fg. ). The lengths of the elements of structure are small as compared wth the lengths of the surface and ts curvature rad. A system of JC x 2 coordnates on the л surface on whch the structure s shaped and a z coordnate n the drecton normal to surface л were chosen n such a way that x' x 2 z axes represent a rght hand system of coordnates. It was assumed that the geometrc centers of the nodes le at ntersectons of parametrc lnes.y' = const x 2 = const z = 0 and that the axes of the lnks concde wth drectons of these parametrc lnes. A typcal segment of such a structure s shown n Fg.. x (x'=const) D JC T~ A t ; г f Щ Щ I lx I к I 4 Д у Д / a v 72 H ' I (xr=const) Fg. We shall ntroduce twelve contnuous suffcently smooth functons defned on the surface л of the structure. These functons represent translatons rotatons deformatons along the coordnate axes and the shape deformatons. The forementoned functons consttute unknown quanttes of the model and have a physcal sense only at the node centres. In every net mesh they can be treated as lnear nature. 3. The analyss of the structure components Node. When a structure s loaded a typcal node s subjected to a homogeneous deformaton havng 2 degrees of freedom. Let (u x u y u : ) be the dsplacements (d x #>. the components of an ndependent vector of rotaton (o) x w y co : ) and (ш х у co X2 oo yz ) the lnear and devatorc components of a homogeneous deformaton respectvely. Denotng by w x w u> dsplacements wthn the node area n drectons X y z respectvely

SECOND ORDER SURFACE GRIDS 23 the followng formula hold w x (xyz) = u x co x x (w xy e z )y (& y (o x: )z (3.) w y (xyz) = u y (&. co xy )x co y y (co yz '& x )z w z (xyz) = u z (oj xz & y )x (& x co yz )y(o z z. After applyng the prncple of deal constrants we can arrve at 2 equatons descrbng the node equlbrum wth 6 generalzed nternal forces and 2 generalzed external forces. From the equatons of the lnear theory of elastcty the general consttutve relatons can be obtaned together wth a formula for the stran energy of a node. Lnk. Let us take nto account a typcal element connectng the th and the j th nodes stuated on the x 2 = const parametrc lne (see Fg. ). Let w x w y w. represent dsplacements of the lnk area n drectons of a local coordnates x y z (see Fg. 2). It s assumed that the lateral cross sectons of a lnk are subjected to homogeneous deformatons n ther planes as well as to the rgd dsplacements (9 degrees of freedom). Hence: w x (xyz) = v x (x) ycp z (x)z(p y (x) (3.2) w y (xyz) = v y (x)yy y (x)z y y (x) 9> x (x)j w : (xyz) = v z (x)y Јy yw(»)c)(*)j2y«(») where v x v y v z are dslocatons <p x q> y q> z rotatons y y y. y yz deformatons of the cross secton of an element along the.v coordnate. 'ytvy f P 4 Fg. 2 The state of lnk area dsplacements s descrbed by 9 functons of the varable x beng the Lagrange's generalzed dsplacements. The assumpton (3.2) can be called the hypothess of a flat homogeneously deformable cross secton wth ndependent rotatons. Ths s a generalzaton of the well known hypothess of Bernoull and Tmoshenko for the classcal model of a bar. The constrants for stresses are assumed n the form (3.3) = o 0. Ths assumpton smplfes consderably the formulae gven below. After applyng the prncple of deal constrants of the statc and knematc types we shall obtan 9 equatons 7 Mech. Teorct Stos. 2 3/83

24 P. WlŚ NIAKOWSKl descrbng the equlbrum of the lnk wth 0 generalzed nternal forces and 9 generalzed external forces. Takng nto account the known equatons of the lnear elastcty the generalzed consttutve equatons and the formulae defnng the stran energy of a lnk can be found. Node lnk node system. A system consstng of the / th node the j th node and the ( j) th lnk connectng these two nodes s presented below (see Fg. 2). From the equlbrum and consttutve equatons descrbng the lnk the dfferental equatons for the generalzed dsplacements can be obtaned. The knematc boundary condtons result from the assumpton that the dsplacements of the approprate boundares of the / th and v' th nodes have to be competble wth the dsplacements of the sutable boundares of the lnk stuated between them. In ths way we obtan functons v x v y ) (p x <p y q> z y y y z y y whch are expressed as the functons of the / th node and the A operator defned as follows: (3.4) zj() =()' ()' These functons can be understood as certan shape functons of the bar treated as a three dmensonal body. The total elastc energy of the (/ j) th lnk s a functon of parameters attrbuted to the / th node and the A operator defned above. The analogous procedure can be appled to lnk stuated on the x = const parametrc lne and connectng the / th and the A: th nodes. Instead of the A operator we deal now wth the A operator defned as follows: (3.5) ~A{ )=()* ( )' 4. Bovcrnng equatons Accordng to the forementoned assumptons the parameters descrbng the dsplacement deformaton strans and stresses as well as the elastc modul are descrbed by the contnuous suffcently regular functons of arguments x x 2. These functons have a physcal sense only n certan ponts of the surface. The dsplacement state of a structure s defned by 2 parameters for each node. A contnuous suffcently regular extenson of these dscrete functons leads to the relatons: U(x} X 2 ) = u x Ux x 2 ) = Uy u(x\ x 2 ) = u z ^(.x x 2 ) = 9 X 2 {x\ x 2 ) = y {x\ x 2 ) = a> lt (x l x 2 ) = co x ш 22(х x 2 ) = co y w(x l x 2 ) = w 2 (o 2 (x\ x 2 ) = co 2 (x x 2 ) = co xy co^x x 2 ) = co xz a> 2 (x l x 2 ) = co yz. The stran energy of a typcal structure segment (.e. the energy of the ( )) th and the (/ /<) th lnks energy of the th node Д уf energy of the 7 th node of energy of the Ar th node) s related to ABCD surface segment wth / l 2 dmensons (see Fg. ). Ths energy s a functon of parameters assgned to the th node and nvolve A and A operators.

SECOND ORDER SURFACE GRIDS 25 Assumng that: (4 2) _ = _ L ( ) = _ L = ( b and that the densty of elastc energy a s equal to the densty of energy a () n the th node the basc relaton of the contnuum model of the consdered structure n ts explct form was found: ao = ^ C^ >»y.kly.mn ~ A^y K y L G^ 'x K L y M F^ y K L r M N D K y K co L H KL»x KL co M R«I M y K r L M l G KLMFL T K Lr MJV I G*4 K r L у F KL o> K to L у Am 2 C KLM co K r LM А к т к(о у A KLMN y KL y MN (4.3) ~ C KL y.k x E KLM K K y LM B K N y KL co MN D««x K r LMN pklm.. JKLMNP _ pklmn. ' r KLMN. у Zl*' LMPJS T KI A T P?s 'y4 KLMP/S T KLM T^jjs C^^^W^LTj^pn f tf"c XL a> D KLM OJ KL т м B K L M T K L M M H K L M N r K L M r N. KL MN P R S = 2. The densty of work of the external forces can be olso defned. The relaton (4.3) was orgnally expressed n the Cartesan coordnate system and then generalzed to a curvlnear orthogonal system of coordnates n terms of whch the surface system s descrbed. The parameters: y KL y K x KL x K r KLM = r KML r KL r K a) KL = co LK co K OJ consttute generalzed components of the state of deformaton wth the geometrc relatons takng the form YKL = L\K b LK 4 e LK & y K = u\ K b I { l L e KL <& L (4.4) T KLM U>LM\K~Ь щ К () b K m M *KL = b>l\c bk(o N L b LK m r K Ы (o\ K 2b%co N where b KL e KL ( ) K represent the components of the second metrc tensor of the surface Rcc's bvector and the symbol of covarant dfferentaton on the surface respectvely. The functons A K L M P R S... C K L M... A K A stand for the tensor of elastc modul of the structure and descrbe ts geometrc and physcal propertes. The components of the stress state of the structure are gven by the formulae

26 P. WlŚ NIAKOWSKl (4.5) [cont.] v J s klm 8a... 8а. 8а = = s = s* = 8r KL 8r KL dv 8a K 8a 8a. dft) KŁ ' 8co K ' dco We see that 2 from 30 ntroduced aqove components of the stress state s of the force type (p KL p K r KL r K r) and the remanng 8 s of the couple type (m KL m K s KLM s KL s K ). Wth the ad of the prncple of vrtual work the equlbrum equatons and the boundary condtons for the contnuum model of the structure can be obtaned n the form p KL \K b L Kp K q L = 0 p K \ K b LK p KL q = 0 m KL \ K bkm K e' I tp K h L = 0 m K \ K e KL p KL b LK m KL h = 0 (4.6) s KLM \ K ~ j(b L Ks KM b^s KL ) r LM f LM = 0 s KL \ K 2b MKS KLM 2b L Ks K r L f L = 0 s*\ K b LK s KL rf = 0. (4.7) p KL n K p L 0 o r "_*=«_; Р к п к р = 0 or u = u * * * m KL n K m L = 0 or # L = L ; m K n K m = 0 or SKL U K_*LM Ш o o r W ł m = 4 A (; sr^n F s Ł = 0 or s K «K * = 0 or w = *J where «K represents the components of a unt normal vector to the boundary 8Q of the structure q K q h K hj' KL f K f are the denstes of the surface type external stress; p pm K tn s KL s K s are the denstes of boundary stresses u L u # Ł co LS OJ LOJ are the gven values of generalzed dsplacements wthn the 8Q. The equlbrum equatons (4.6) and the boundary condtons (4.7) together wth the consttutve (4.5) and geometrc (4.4) relatons form the basc system of equatons descrbng the contnuous model of the structure. Ths system enables us to calculate the dsplacement dstrbuton n the lnk and node areas as well as the stress dstrbuton. It must be stressed that parameters y K y.kl r KL K (ol [see (4.4)] defne the components of the plate lke deformaton state whle y K L co KL x K r KLM the components of the plane lke deformaton state. The components yxl. w *. *к» _ L m a r e defned exactly as n the rst order model (see [6]) however the parameters whch do not appear n that model.e. T K L M T K L T K <X>kl cok ш result from the deformablty of a node and the deformablty of lnk's lateral cross sectons. The analyss of the nfluence of the second order parameters on the nternal forces together wth the sutable numercal calculaton wll be the subject of separate papers.

SECOND ORDER SURFACE GRIDS 27 References. W. NOWACKI Problems related to the theory of flat grds. Arch. Mech. Stos. 6 954. (n Polsh) 2. W. GUTKOWSKI Grd surface structures. Mech. Teor. Stos. 3 3 965. (n Polsh) 3. H. FRĄ CKIEWICZ Deformaton of a dscrete set of ponts. Arch. Mech. Stos. 3 8 966. 4. M.T. HUBER Probleme der Statk technsch wehlger orthotroper Platten Zurch 929. 5. C. WOŹ NIAK Introducton to mechancs of fbrous meda Arch. Mech. Stos. 5 6 964. 6. С. WOŹ NIAK Lattce type surface structures PWN Warsaw 970. (n Polsh). Р е з ю ме У Р А В Н Е НЯ И Л И Н Е Й Н Й О Т Е О Р И В Т О Р О О Г Р Я ДА У П Р У Г Х И П О В Е Р Х Н О С Т Х Н Ы Р О С Т В О РК О В р а б ое т п р е д с т а в о л ел ни н е й не ы у р а в н е я н ис т а т ии к у п р у гх и п о в е р х н о с х т нр ыа с т в о к р ои м е ю щ их п л о т ню у и р е г у л я рю н ус е т ку э л е м с о е д и н е н е н мы е ж ду с о б й о р пи п о м о щи п р и з м а т и ч ех с кс ти е р ж ч е н и. е е н т окво т о р х ы д е ф о р м и р о в а н нп ыр яе м о у г о л е ь нуыз лы нй е и м е ю щ х и п р я м о у г о л е ь сн ое П р и н и мя а и с х о д не ы д а н н ы: еу р а в н е я н ил и н е й нй о т е о р и и у п р у г о с т а и т а к же п о д х о д я е щ и к и н е м а т и ч е сгк ии п о т ы е зп о л у ч о е нв а р и а ц и о нм н мы е т о дм о у р а в н е я н ис п л о ш Р а б оа т с о д е р жт ио б о б щ в о л я ю ще иу ч и т ы е е н ит е о р и и В р о в а не и п о п е р е ч но осге ч е ня и с т е р ж н о о мг о д е я л п р о г о н. а о з н я к ав ы х о д я ще из а п р е д еы л т е о р и и г о р я д а п о з в ь а эт ф ф е кы т в ы с ш х и р я д о " в (р а з м еы р у з л о в и х д е ф о р м и р о в а нд ие еф о р м и нй ес о е д и н я ю х щ уи з л ы. ) Streszczene RÓWNANIA LINIOWEJ TEORII DRUGIEGO RZĘ DU SPRĘ Ż YSTYC H POWIERZCHNIOWYCH RUSZTÓW W pracy wyprowadzono lnowe równana statyk sprę ż ystyc h rusztów powerzchnowych o gę stej regularnej satce elementów których odksztalcalne prostopadloś cenne wę zły połą czone są za pomocą pryzmatycznych prę tów o przekroju prostoką tnym. Przyjmując za punkt wyjś ca równana lnowej teor sprę ż ystoś coraz zakładając odpowedne hpotezy knematyczne otrzymano na drodze waracyjnej równana cą głego modelu dź wgara. Praca zawera uogólnene teor Woź naka wykraczają ce poza teorę I go rzę du zezwalają ce na uwzglę dnene efektów wyż szych rzę dów" (wymary wę złów ch odkształcalnoś ć odksztalcalność przekrojów poprzecznych prę tów łą czą cyc h wę zły). Praca została złoż ona w Redakcj dna 3 marca 983 roku