APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS
|
|
- Jasmine Boyd
- 5 years ago
- Views:
Transcription
1 6th ICPT, Sapporo, Japan, July 008 APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS James MAINA Prncpal Researcher, Transport and Infrastructure Engneerng, CSIR Bult Envronment P.O.Box 395, Pretora 000, South Afrca Yoshak OZAWA Department of Cvl and Envronmental Engneerng, Tokyo Denk Unversty Ishaka, Hatoyama Town, Hk-gun, Satama , Japan Kunhto MATSUI Department of Cvl and Envronmental Engneerng, Tokyo Denk Unversty Ishaka, Hatoyama Town, Hk-gun, Satama , Japan ABSTRACT GAMES software s well known for ts capablty to compute responses for unformly dstrbuted load actng on the surface of a mult-layered lnear elastc system. In ths study a method was developed to approxmate rgd plate loadng to be used n GAMES software. When a rgd plate load s actng on the surface of a multlayered system, the resultng nterface contact stress dstrbuton s not unform. Snce contact stress dstrbuton between rgd plate and sem-nfnte medum s well known, ths dstrbuton was approxmated usng unformly dstrbuted multple loads and analyss performed usng GAMES. Results have shown good agreement wth the theory for the case of a sem-nfnte medum. Furthermore, extenson of ths method to multlayered system has shown good results f the rato of the frst layer thckness over the loadng plate radus s greater than. KEY WORDS rgd plate loadng, elastc multlayered system, GAMES, ax-symmetrc analyss 43
2 w (cm) Approxmate Analyss of Rgd Plate Loadng on Elastc Mult-layered Systems INTRODUCTION Rgd plate loadng test s generally used to evaluate load bearng capacty of pavement layers lke subgrade, subbase and base layers. When structures lke pavements are to be analyed, most of the tme the structural system s consdered to be multlayered lnear elastc. Furthermore, unformly dstrbuted crcular load s assumed to act on the surface of the pavement and t s possble to use software lke ELSA, CHEVRON, BISAR ) as well as GAMES ). However, when a test lke rgd plate loadng s performed, the nterface contact stresses are not unformly dstrbuted and software such as the ones lsted above may not be used for analyss. Authors of ths paper were unable, through lterature search, to fnd any document descrbng the applcaton of multlayer lnear elastc theory to the rgd plate loadng. Because of the facts hghlghted above, ths research evaluates the method of approxmatng theoretcal solutons of rgd plate loadng by usng GAMES software. unformly dstrbuted loads Sneddon derved the solutons for a ( P 49kN, E 00 MPa, 0.35 ) sem-nfnte medum subected to crcular rgd body loadng as shown n Fgure 3). In ths case, surface deflecton under the base of the rgd body s constant and normal stresses at the surface but outsde the base become ero. Fgure was drawn to compare surface deflectons from the applcatons of unformly dstrbuted loadng and rgd plate loadng. Ths fgure clearly shows that n the vcnty of the load applcaton area, deflectons from the analyss of unformly dstrbuted loadng are dfferent from results of rgd body loadng. Ths phenomenon s supported by Sant - Venant prncple, whch states that there s hgh nfluence of load dstrbuton near the pont of applcaton and ths nfluence dmnshes as you move away from the pont of load applcaton. It s clear from Fgure that there s a very good match of deflecton results for r a, there s a good match of deflecton results. In ths research, dstrbuton of contact stresses derved by Sneddon was approxmated by multple unformly dstrbuted loads, and responses for sem-nfnte medum were analyed usng GAMES software. Accuracy of ths approach was evaluated by comparng results obtaned by Sneddon for a sem-nfnte layer to results from usng approxmate stress dstrbuton. Extenson of ths approach to a multlayered system was evaluated by checkng whether deflectons under the base of a rgd plate remaned constant when layer modul were vared for a two-layer system. w φ = a r a Fgure Deformaton due to rgd body loadng Rgd plate loadng Unform loads ρ = r /a Fgure Surface deflectons due to rgd plate and 44
3 Mana, Oawa and Matsu Contact stress dstrbuton under rgd plate Theoretcal stress dstrbuton Assume a rgd plate load actng on the surface of a sem-nfnte body to be consdered an ax-symmetrc problem, and then the boundary condtons can be represented as shown below: w (r,0) 0 ra (a) ( r,0) 0 a r (b) s the amount of compresson (surface deflecton) and a s radus of the loadng plate. Accordng to Sneddon, shear modulus G, Posson s rato, and dstrbuton of contact stresses on a sem-nfnte medum p (r) are related as shown n the followng equaton: G p( r) ( ) a r () From equaton (), contact stress at the center of the load s G ( ) a, whle at the load edge, contact stresses becomes. Equaton () can be ntegrated across the crcular loadng secton to obtan the total load P as shown below: P a 4aG p( r) rddr 0 0 (3) From equatons () and (3), the relaton between p (r) and P can be obtaned as follows: P p( r) (4) a a r Furthermore, Sneddon determned surface deflectons usng the followng equaton: 0 r a ( r,0) a sn r a r w (5) Approxmate contact stresses Contact stress p (r) as expressed by equaton (4), becomes P ( a ) at the load center, where r 0 and at the load edge where r a. Such knd of a load dstrbuton can not be handled n any multlayer lnear elastc software such as GAMES. In order for such a contract stress dstrbuton to be analyed, t s mportant to approxmate ts dstrbuton by multple unformly dstrbuted loads. The method that was used to approxmate ths dstrbuton s explaned below. The load was subdvded nto N dfferent crcular loads of dfferent rad wth same centers. The smallest radus s r and the largest radus s r a N. Assumng every radus of each rgd plate crcle to represent the total loadng of P, we can obtan N domans comprsed of the nner most crcle of radus r and rngs formed by dfference between adacent crcles 45
4 Approxmate Analyss of Rgd Plate Loadng on Elastc Mult-layered Systems r r (, ).The value of r s determned by equalty wth the total loadng. Ths means, total load n each doman s equal to P P N. Ths total load n each doman s assumed to be unformly dstrbuted and we can ntegrate equaton (4) from r to r to obtan total th load for the doman. Consderng that total load n each doman can be expressed as P P / N, the followng equaton s obtaned: P P r a r a (6) N Makng use of the relatonshp between r and r shown n Equaton (5), yelds the followng equaton: r a r a N (7) Replacng the contact stresses between r and r by unformly dstrbuted load obtan: p to p N ( r P r ) (8) Assume a load of magntude P s actng on a crcular rgd plate. In order for ths load to be analyed usng GAMES software, t wll be approxmated as multple concentrc crcular th loads whose centers are on the loadng plate. Furthermore, for the load of radus r, responses (.e. stresses and dsplacements) at pont can be represented by ( r ), and for th ( ) load of radus r, responses can be represented by ( r ). When a unformly dstrbuted load of magntude p s actng on rng of wdth ( r r r ), responses at pont can be obtaned from the followng equaton: r ) ( r ) (9) ( th represents responses at pont when a unformly dstrbuted unt load s actng wthn a rng. Supermposng responses on pont due to the effect of the loads actng on all the rngs, the total response S may be obtaned as follows: S P N N p ( r r ) (0) p s the unformly dstrbuted load on rng, whle P s ts total load. Results from multple unformly dstrbuted loadng When multple unformly dstrbuted loads were analyed, there was mprovement n accuracy for bgger number of multple loads but the down sde of t was the ncrease n computatonal tme. Three methods of load subdvsons were used where N 5,0, 5, and 46
5 Mana, Oawa and Matsu results from numercal analyses were compared wth those obtaned by Sneddon. -σ (MPa) Fgure 3 s a comparson of contact stress dstrbutons between theoretcal analyss and approxmated load dstrbuton usng the three dfferent numbers of subdvsons. There s a sgnfcant dfference between contact stresses near the load edge for the theoretcal and approxmate loads but as the number of multple loads ncreases the contact stresses approach theoretcal values. At the load edge, theoretcal contact stress s nfnte whle approxmate contact stress stll mantans a fnte value. Fgure 4 shows comparson of surface deflectons for theoretcal load as determned by equaton (5) and surface deflectons as computed by GAMES software for the three dfferent N load subdvsons. Ths fgure shows that as the number of load subdvsons N ncreases, results of surface deflectons approach the theoretcal values. However, t s not the same tendency for deflectons near the load edge, where there s very lttle mprovement as the number of load subdvsons ncreases. The reason may be attrbuted to the fact t s dffcult for approxmate loadng to smulate theoretcal loadng near the edge, whch approaches nfnte. Deflecton results by fewer load subdvsons, n ths example 5 subdvsons, show poor accuracy near the load center whle results from bgger number of subdvsons, that s N 0, 5 show very good agreement wth theoretcal results. It s -aσ/(gδ ) An approxmate value by 5 ponts An approxmate value by 0 ponts An approxmate value by 5 ponts Theoretcal responses r /a Fgure 3 Theoretcal and approxmate Fgure 4 Subdvsons and surface deflectons contact stresses ( P 49kN, E 00 MPa, 0.35 ) ( P 49kN, 00 MPa, 0.35 ) Theory 0.0 Theory 0. Theory 0.4 Theory 0.6 Theory 0.8 Theory.0 GAMES 0.0 GAMES 0. GAMES 0.4 GAMES 0.6 GAMES 0.8 GAMES r / a * Legend value s a Fgure 5 Comparson of responses between approxmate and theoretcal loads w (cm) An approxmate value by 5 ponts An approxmate value by 0 ponts An approxmate value by 5 ponts Theoretcal responses r /a E therefore recommended that at least 0 subdvsons are mportant to obtan reasonably good results. Fgure 5 shows comparson of normal stresses ( ) at dfferent depths from the surface ( a ) n order to confrm ablty of 0 load subdvsons to approxmate theoretcal analyss. In ths fgure, dmensonless parameters of normal stresses ( ) dvded by G a 47
6 Approxmate Analyss of Rgd Plate Loadng on Elastc Mult-layered Systems was plotted aganst r a. Numbers n the legend of Fgure 5 ndcate a ratos. In ths fgure, sold lnes represent theoretcal results whle approxmate results are represented by the marks as defned n the legend. The fgure shows for surface results a dfference of about 5% between theoretcal and approxmate loads and a very good agreement for ponts deeper than a 0.. EXTENSION TO A TWO-LAYER STRUCTURE Method for applcaton evaluaton Equaton (4) shows dstrbuton of contact stresses resultng from a rgd plate loadng actng on the surface of a sem-nfnte medum and ths dstrbuton s not the same as for a two-layer structure. However, Sant-Venant s prncple states that nfluence of load dstrbuton s n the vcnty of the load applcaton pont and t s assumed that f the thckness of the surface layer s bg enough, contact stress dstrbuton for sem-nfnte medum may be used n two-layer system wthout affectng computatonal accuracy. Fgure 6 shows a two-layer analytcal model. Thckness of the frst layer and the rato of the elastc modul between the frst and second layer were vared and analyses performed n order to evaluate the lmtaton of the approxmate method explaned above. The accuracy of the analyss was evaluated n terms of deflectons under the base of the rgd plate n order to see whether they remaned constant. The root mean square error shown n equaton () between mean deflecton under the rgd plate and each deflecton wthn the doman was used as evaluaton crtera. RMSE () w 5 5 w w where w s the mean surface deflecton and w s the deflecton wthn r 0~4cm. ( 0 r a ) at cm each ( r a ). Good agreement of dsplacements for all the ponts of nterest wll make RMSE equal to ero. Ths wll be an ndcaton that approxmate soluton s equal to the theoretcal soluton. It s clear from Fgure 4 that rrespectve of the number of load subdvsons, there s poor agreement between approxmate soluton and theoretcal soluton at the load edge, and hence results at the load edge poston were not consdered n equaton (). The scope of applcaton of the approxmate loadng was evaluated by usng a two-layer system shown n Fgure 6 by followng the procedure explaned below: ) Thckness of the frst layer h was vared from h a to 3 n 0. steps resultng n thckness varatons. ) For each value of h a, E E was vared between 0.9~0 resultng n varatons. 3) h a and E E varables resulted n combnatons, for whch analyses were performed and example of the results s shown n Fgure 7. 4) Usng equaton (), each combnaton was evaluated and analytcal errors RMSE are gven n Table. 48
7 Mana, Oawa and Matsu Vertcal Load 0 P a 5cm 49 kn 0.04 r / a h varable r Layer modul Posson's rato E E varable 0.35 E 00MPa 0.35 w (cm) Fgure 6 Two layer structure ( h a. 8, E E 3 ) Fgure 7 Surface deflectons E /E Table Least mean square error (LMSE) from Eq.() h / a Table shows errors ( RMSE) for each combnaton consdered. Wth regards to the ratos of layer modul and rrespectve of the h a value, rato 7 gves the largest RMSE value. Smaller ratos of layer modul, where the two-layer system s close to sem-nfnte medum, result n smaller errors n the surface deflectons. Practcal applcaton of the approxmate method should be for all combnatons wth h a and RMSE less than %. All the errors fallng outsde ths crteron have been grayed out n Table. Fgure 7 shows surface deflecton w for the acceptable regon of applcaton wth the hghest error ( h a. 8, E E 3). Surface deflecton at the center of the load s cm, whle at the load edge where ( r a ) the surface deflecton s 0.08 cm, showng consderably smaller error. Rgd plate versus unformly dstrbuted loads An example for the above case wth largest RMSE was assumed. Thckness of the frst layer was h 30cm and layer modulus was E 700 MPa. Elastc modulus of the second layer was E 00 MPa. Radus of the load was a 5cm and analyses for a 49 kn was performed consderng a rgd plate load and unformly dstrbuted load. Results were compled and presented n Cartesan coordnates where Fgure 8(a) shows contour plot of x due to rgd plate load and Fgure 8(b) s for unformly dstrbuted load. At the boundary between the frst and second layer where 30cm and ust drectly under the load, both 49
8 Approxmate Analyss of Rgd Plate Loadng on Elastc Mult-layered Systems (a) Rgd plate load (b) Unformly dstrbuted load Fgure 8 Contour plots for n X-Z plane loads result n tensle strans but results from the rgd plate loadng cover a relatvely smaller area. For the case of unformly dstrbuted load, there are hgh strans, whch are transmtted downwards whle for the rgd plate loadng strans are transmtted from the load edge downwards. CONCLUSIONS An approxmate method was developed whereby GAMES software, whch s capable of analyng effect of unformly dstrbuted load on a multlayered lnear elastc system, was utled to analyse a multlayered lnear elastc system under the acton of rgd plate loadng. Results obtaned are summared as follows: ) Wth an excepton of rgd plate edge, very accurate results were obtaned when load subdvson s greater or equal to 0. ) It was possble to accurately evaluate a multlayered lnear elastc system usng rgd approxmate contact stresses for a sem-nfnte medum. 3) There were very hgh strans developed at the layer nterface between the frst and second layers for the case of unformly dstrbuted load. However, because of better load dsperson from the rgd plate edge downward, relatvely low strans were developed at the layer nterface rght below the load center. 4) For h a, contact stresses for rgd plate loadng on a sem-nfnte medum may be used to analye a multlayered lnear elastc system rrespectve of the magntude of layer modul ratos. REFERENCES ) De Jong, D.E. Peut, M.G.F. and Korswagen, A.R.:COMPUTER PROGRAM BISAR, Layered Systems under Normal and Tangental Surface Loads, Konnklke/Shell -Laboratorum, Amsterdam, 979. ) Mana, J.W. and Matsu, K.: Development of Software for Elastc Analyss of Pavement Structures due to Vertcal and Horontal Surface Loadngs, Journal of Trasportaton Research Record No.896, TRB, pp.07-8, ) Sneddon, I. FOURIER TRANSFORMS, McGraw-Hll Book Company, pp , 95. x 50
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationFUZZY FINITE ELEMENT METHOD
FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationGEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE
GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE Prof. J. N. Mandal Department of cvl engneerng, IIT Bombay, Powa, Mumba 400076, Inda. Tel.022-25767328 emal: cejnm@cvl.tb.ac.n Module - 9 LECTURE - 48
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationA Robust Method for Calculating the Correlation Coefficient
A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal
More informationOne Dimensional Axial Deformations
One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the
More informationFINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION
VOL. 6, NO. 3, MARCH 0 ISSN 89-6608 006-0 Asan Research Publshng Network (ARPN). All rghts reserved. FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION Adel A. Al-Azzaw and Al S. Shaker
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationThermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)
hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00
ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationChapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods
Chapter Eght Energy Method 8. Introducton 8. Stran energy expressons 8.3 Prncpal of statonary potental energy; several degrees of freedom ------ Castglano s frst theorem ---- Examples 8.4 Prncpal of statonary
More informationTorsion Stiffness of Thin-walled Steel Beams with Web Holes
Torson Stffness of Thn-walled Steel Beams wth Web Holes MARTN HORÁČEK, JNDŘCH MELCHER Department of Metal and Tmber Structures Brno Unversty of Technology, Faculty of Cvl Engneerng Veveří 331/95, 62 Brno
More informationSTATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION
STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationNON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE
October 1-17, 008, Bejng, Chna NON LINEAR ANALYSIS OF SRUCURES ACCORDING O NEW EUROPEAN DESIGN CODE D. Mestrovc 1, D. Czmar and M. Pende 3 1 Professor, Dept. of Structural Engneerng, Faculty of Cvl Engneerng,
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationOFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES
ICAMS 204 5 th Internatonal Conference on Advanced Materals and Systems OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES VLAD LUPĂŞTEANU, NICOLAE ŢĂRANU, RALUCA HOHAN, PAUL CIOBANU Gh. Asach Techncal Unversty
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationSCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.
SCALARS AND ECTORS All phscal uanttes n engneerng mechancs are measured usng ether scalars or vectors. Scalar. A scalar s an postve or negatve phscal uantt that can be completel specfed b ts magntude.
More informationElectrical double layer: revisit based on boundary conditions
Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer
More informationSIMPLE LINEAR REGRESSION
Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two
More informationTHE EFFECT OF BEAM TO COLUMN CONNECTION IN ARC PORTAL FRAME
THE EFFECT OF BEAM TO COLUMN CONNECTON N ARC PORTAL FRAME Asko Keronen Rakenteden Mekankka, Vol. 26 No 2 1993, ss. 35-5 SUMMARY A full scale rc (renforced concrete) portal frame has been bult n order to
More informationχ 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
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationσ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review
Chapter Four States of Stress Part Three When makng your choce n lfe, do not neglect to lve. Samuel Johnson Revew When we use matrx notaton to show the stresses on an element The rows represent the axs
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationCHAPTER 9 CONCLUSIONS
78 CHAPTER 9 CONCLUSIONS uctlty and structural ntegrty are essentally requred for structures subjected to suddenly appled dynamc loads such as shock loads. Renforced Concrete (RC), the most wdely used
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationEffect of loading frequency on the settlement of granular layer
Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationTurbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH
Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant
More informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationSecond Order Analysis
Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to
More informationStudy on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component
Study on Non-Lnear Dynamc Characterstc of Vehcle Suspenson Rubber Component Zhan Wenzhang Ln Y Sh GuobaoJln Unversty of TechnologyChangchun, Chna Wang Lgong (MDI, Chna [Abstract] The dynamc characterstc
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationSimulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM
Smulaton of 2D Elastc Bodes wth Randomly Dstrbuted Crcular Inclusons Usng the BEM Zhenhan Yao, Fanzhong Kong 2, Xaopng Zheng Department of Engneerng Mechancs 2 State Key Lab of Automotve Safety and Energy
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationPES 1120 Spring 2014, Spendier Lecture 6/Page 1
PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura
More informationVisco-Rubber Elastic Model for Pressure Sensitive Adhesive
Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA
More informationLifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model
ECCMR, Prague, Czech Republc; September 3 th, 2015 Lfetme predcton of EP and NBR rubber seal by thermos-vscoelastc model Kotaro KOBAYASHI, Takahro ISOZAKI, Akhro MATSUDA Unversty of Tsukuba, Japan Yoshnobu
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationSo far: simple (planar) geometries
Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationAPPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam
APPENDIX F A DISPACEMENT-BASED BEAM EEMENT WITH SHEAR DEFORMATIONS Never use a Cubc Functon Approxmaton for a Non-Prsmatc Beam F. INTRODUCTION { XE "Shearng Deformatons" }In ths appendx a unque development
More informationEVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES
EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge
More informationˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)
7.. = = 3 = 4 = 5. The electrc feld s constant everywhere between the plates. Ths s ndcated by the electrc feld vectors, whch are all the same length and n the same drecton. 7.5. Model: The dstances to
More information16 Circular Footing on a Semi-infinite Elastic Medium
Crcular Footng on a Sem-nfnte Elastc Medum 16-1 16 Crcular Footng on a Sem-nfnte Elastc Medum 16.1 Problem Statement Ths verfcaton problem nvolves a rgd crcular footng restng on an elastc half-space. In
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationStatistics II Final Exam 26/6/18
Statstcs II Fnal Exam 26/6/18 Academc Year 2017/18 Solutons Exam duraton: 2 h 30 mn 1. (3 ponts) A town hall s conductng a study to determne the amount of leftover food produced by the restaurants n the
More informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationPlease review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.
ME 270 Fall 2013 Fnal Exam NME (Last, Frst): Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: INSTRUCTIONS
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationANOVA. The Observations y ij
ANOVA Stands for ANalyss Of VArance But t s a test of dfferences n means The dea: The Observatons y j Treatment group = 1 = 2 = k y 11 y 21 y k,1 y 12 y 22 y k,2 y 1, n1 y 2, n2 y k, nk means: m 1 m 2
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationMECHANICS OF MATERIALS
Fourth Edton CHTER MECHNICS OF MTERIS Ferdnand. Beer E. Russell Johnston, Jr. John T. DeWolf ecture Notes: J. Walt Oler Texas Tech Unversty Stress and Stran xal oadng Contents Stress & Stran: xal oadng
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationACTM State Calculus Competition Saturday April 30, 2011
ACTM State Calculus Competton Saturday Aprl 30, 2011 ACTM State Calculus Competton Sprng 2011 Page 1 Instructons: For questons 1 through 25, mark the best answer choce on the answer sheet provde Afterward
More informationUniformity of Deformation in Element Testing
Woousng Km, Marc Loen, ruce hadbourn and Joseph Labu Unformt of Deformaton n Element Testng Abstract Unform deformaton s a basc assumpton n element testng, where aal stran tpcall s determned from dsplacement
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
More informationStatistics Chapter 4
Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment
More informationI have not received unauthorized aid in the completion of this exam.
ME 270 Sprng 2013 Fnal Examnaton Please read and respond to the followng statement, I have not receved unauthorzed ad n the completon of ths exam. Agree Dsagree Sgnature INSTRUCTIONS Begn each problem
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationPrinciples of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer
Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationStatistical Evaluation of WATFLOOD
tatstcal Evaluaton of WATFLD By: Angela MacLean, Dept. of Cvl & Envronmental Engneerng, Unversty of Waterloo, n. ctober, 005 The statstcs program assocated wth WATFLD uses spl.csv fle that s produced wth
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationOn Pressure Distributions of Drum Brakes
Yuan Mao Huang Professor e-mal: ymhuang@ccmsntuedutw J S Shyr Research Assstant Department of Mechancal Engneerng, Natonal Tawan Unversty, Tape, Tawan, Republc of hna On Pressure Dstrbutons of Drum Brakes
More informationANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION
VOL. 6, NO., NOVEMBER ISSN 89-668 6- Asan Research Publshng Network (ARPN). All rghts reserved. ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION Adel A.
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationEffect of anisotropy on laminated composite plates containing circular holes
Indan Journal of ngneerng & Materals Scences Vol. 1, June 005, pp. 07-13 ffect of ansotropy on lamnated composte plates contanng crcular holes H Murat Arslan Cukurova Unversty, Cvl ngneerng Department,
More informationAssignment 5. Simulation for Logistics. Monti, N.E. Yunita, T.
Assgnment 5 Smulaton for Logstcs Mont, N.E. Yunta, T. November 26, 2007 1. Smulaton Desgn The frst objectve of ths assgnment s to derve a 90% two-sded Confdence Interval (CI) for the average watng tme
More informationTHE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES
THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES W. C. Lao Department of Cvl Engneerng, Feng Cha Unverst 00 Wen Hwa Rd, Tachung, Tawan SUMMARY: The ndentaton etween clndrcal ndentor
More informationAnswers Problem Set 2 Chem 314A Williamsen Spring 2000
Answers Problem Set Chem 314A Wllamsen Sprng 000 1) Gve me the followng crtcal values from the statstcal tables. a) z-statstc,-sded test, 99.7% confdence lmt ±3 b) t-statstc (Case I), 1-sded test, 95%
More informationGEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam
GEO-SLOPE Internatonal Ltd, Calgary, Alberta, Canada www.geo-slope.com Introducton Vbratng Beam Ths example looks at the dynamc response of a cantlever beam n response to a cyclc force at the free end.
More informationSpring 2002 Lecture #13
44-50 Sprng 00 ecture # Dr. Jaehoon Yu. Rotatonal Energy. Computaton of oments of nerta. Parallel-as Theorem 4. Torque & Angular Acceleraton 5. Work, Power, & Energy of Rotatonal otons Remember the md-term
More informationPhysics 111: Mechanics Lecture 11
Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton
More informationLIMIT ANALYSIS IN LARGE DISPLACEMENTS OF MASONRY ARCHES SUBJECTED TO VERTICAL AND HORIZONTAL LOADS
SAC2014 9 th Internatonal Conference on Structural Analyss of storcal Constructons F. Peña & M. Chávez (eds.) Mexco Cty, Mexco, 14 17 October 2014 LIMIT ANALYSIS IN LARGE DISPLACEMENTS OF MASONRY ARCES
More informationPlease review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.
ME 270 Sprng 2014 Fnal Exam NME (Last, Frst): Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: INSTRUCTIONS
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationDESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS
Munch, Germany, 26-30 th June 2016 1 DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS Q.T. Guo 1*, Z.Y. L 1, T. Ohor 1 and J. Takahash 1 1 Department of Systems Innovaton, School
More informationStatistics for Economics & Business
Statstcs for Economcs & Busness Smple Lnear Regresson Learnng Objectves In ths chapter, you learn: How to use regresson analyss to predct the value of a dependent varable based on an ndependent varable
More informationMEASUREMENT OF MOMENT OF INERTIA
1. measurement MESUREMENT OF MOMENT OF INERTI The am of ths measurement s to determne the moment of nerta of the rotor of an electrc motor. 1. General relatons Rotatng moton and moment of nerta Let us
More information