THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES
|
|
- Raymond Byrd
- 5 years ago
- Views:
Transcription
1 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 and angle-pl compostes s nvestgated. An exact Green s functon for the surface dsplacement of an angle-pl eam s derved usng Jones s dsplacement formulaton. B matchng the dsplacements of the top surface of the angle-pl lamnate and those of the ndentor wthn the contact regon, we can fnd the contact stress dstruton, the magntude of contact force and the ndentaton. In order to speed up the formulaton of fundamental deflecton mode shape of Green s functon, method of ntal functon (MIF) s adopted to fnd the dsplacements n the top and ottom surface of an angle-pl lamnate. In the MIF, the order of sstem equaton s onl 3 3, whch wll reduce the computaton tme. Numercal results for angle-pl lamnates sujected to clndrcal ndentaton are evaluated to verf the applcalt of the proposed method. KEYWORDS: angle-pl lamnates, ndentaton, pont matchng method, method of ntal functons INTRODUCTION The contact prolems of composte lamnate are mportant n the analss of foregn oject mpact and flexural measurement of composte materals. Several authors have studed the contact ehavors etween a clndrcal ndentor and orthotropc eams[,,3]. Sanar solved the
2 contact prolem of orthotropc eam the superposton of half space soluton and eam soluton[]. B applng the exact dsplacement soluton for clndrcal endng of orthotropc eams, the exact Green s functon s formulated n solvng the contact prolems[3]. However the ndentaton etween clndrcal ndentor and angle-pl compostes s seldom nvestgated. In ths stud, an exact Green s functon for the surface dsplacement of an angle-pl eam s derved usng Jones s[4] dsplacement formulaton. Ths Green s functon satsfes the dsplacement and stress contnut condton at the nterface, and the prescred stress condton at the top and ottom surfaces. B matchng the dsplacements of the top surface of the angle-pl lamnate and those of the ndentor wthn the contact regon, we can fnd the contact stress dstruton, the magntude of contact force and the ndentaton. In constructng the exact Green s functon, we have to solve a matrx of order 6N 6N, wth M tmes, where N s the pl numer and M s the numer of asc flexural modes. If the pl numer ncreases, the computaton tme wll ncrease drastcall. In order to speed up the calculaton, method of ntal functon (MIF)[5] s adopted to fnd the dsplacements n the top and ottom surface of an angle-pl lamnate. In the MIF, the order of sstem equaton s onl 3 3, whch wll reduce the computaton tme. Numercal examples are llustrated to demonstrate the contact ehavors of angle-pl lamnates. CYLINDRICAL BENDING OF ANGLE-PLY LAMINATES In order to stud the contact ehavor of angle-pl lamnates, the clndrcal endng prolems of an ansotropc lamnate s nvestgated frst whch wll provde the Green s functon for contact prolem. Assumng plane stan condton n the x-drecton, an ansotropc eam has wdth and thcness h n the -drecton s shown n Fgure. The stran feld s defned as u v w ε x 0, ε, ε x v w u u γ +, γ x, γ x () The consttutve relaton for a generall ansotropc materal s xx c c c c c c c c ε c c c c ε () c c γ x c c 0 γ x x c c c 0 0 c γ x
3 In whch c j s the elastc constant for each laer. The eam s smpl supported along 0 and, wth the followng oundar condtons w(0,)w(,)0 (3a) ( 0, ) (, ) 0 (3) x ( 0, ) (, ) 0 (3c) x (, 0) (, 0) (, 0) 0 (3d) x ( h, ) q ( ), ( h, ) ( h, ) 0 (3e) x If the nerta force s neglected, the equlrum equaton can e wrtten as + 0 (4a) x, x, + 0 (4),, + 0 (4c),, The Dsplacement felds for ths prolem can e expressed as nπ u (, ) Un ( ) cos nπ v(, ) Vn ( ) cos nπ w (, ) Wn ( ) sn The aove canddates satsf Eqns. (3a), (3) and (3c). When susttute Eqn. (5) nto Eqn.(4) and Eqn.(), the followng equatons can e derved. '' '' pc66 U( ) + c55 U ( ) pc6 V( ) + c45v ( ) + pc ( 36 + c45) W' 0 (6) '' '' pc6 U( ) + c45 U ( ) pc V( ) + c44 V ( ) + pc ( 3 + c44 ) W' 0 (7) '' pc ( 45 + c36) U'( ) pc ( 44 + c3) V'( ) + c33 W( ) pc44 W ( ) 0 (8) In whch pnπ /, the suscrpt n of U, V and W are omtted for smplct. In order to solve the U (), V () and W (), we let [4] U ( ) U V ( ) V 0 0 W ( ) W e 0 e λ λ e λ a c55, a c44, a3 c33 p c66, p c, 3 p c44 (0) c c45, d pc6, e pc ( 36 + c45), e pc ( 3 + c44 ) B susttutng Eqns. (9) and (0) nto Eqns. (6)-(8), the nontrval solutons exst f the followng relaton holds λ a + λ c d λe det λ c d λ a + λe 0 () λe λe λ a (5) (9)
4 After solvng the 6th order characterstc equaton of λ, the fnal form of U n (), V n () and W n () can e elded as 6 * λn n U ne 6 * λn n ( ) V ne () 6 * λn n ( ) W ne U ( ) V W * * * In whch U, V and W are the correspondng egenvectors when solvng Eqn.(). The onl unnowns are n. For a fxed p, there are sx n s n each laer. These constants can e solved from the nterface contnut and the tracton condtons at the top and ottom surfaces. If the lamnate s made of N laers, the contnut of dsplacements and stresses at the nterface requres that + u (, h ) u (, 0) + v ( h, ) v (, 0) + w (, h ) w (, 0) + ( h, ) (, 0 ) (3) + ( h, ) (, 0) + x ( h, ) x (, 0),,,..., N In whch h s the thcness of the -th laer. The tracton condtons on the top and ottom surface mpose that N N (, 0) (, 0) x (, 0) ( h, N) x( h, N ) 0 N ( h, N ) q( ) (4) If the normal stress on the top surface s q( ) o sn p (5) then the coeffcents n can e solved through a set of matrx of order 6N 6N. The vertcal dsplacements w(,h) and w(,0) are of great mportance n the contact prolem. CYLINDRICAL INDENTATION OF SIMPLY SUPPORTED ANGLE-PLY LAMINATES Consder a smpl supported angle-pl lamnate sujected to the ndentaton of a rgd clnder as shown n Fg.. The radus of the ndentor s R. The frcton etween the ndentor and lamnate s neglected. The ndentaton s assumed to e smmetrc aout the mdspan of the eam. The contact length s c. In the contact prolems the dstruton of contact pressure has to e determned. The pont matchng method s adopted to fnd the contact stress dstrutons.
5 From prevous secton, the dsplacement at the top surface of the eam due to a normal stress o sn( nπ/ ) appled at h s shown to e nπ wh (, ) Wn ( h) sn (6) For an artrar normal stress q() at the top surface, the dsplacement of the eam wll e nπ wh (, ) qw n n( h) sn (7) In whch nπ qn q d ( )sn 0 (8) + c nπ q d ( )sn c From Eqns. (7) and (8), we have nπξ nπ w(, h) Wn ( h) q( ξ )sn dξsn Ω Ω G( ξ) q( ξ ) dξ Where Ω s contact regon, and the Green s functon s defned as[3] nπ nπξ G( ξ ) Wn ( h) sn sn (0) The relatve vertcal dsplacements of an pont w(,h) nsde the contact regon and the center pont w(0.5,h) can e expressed as[] 0.5 ( 0.5) w(, h) w(0.5, h) R R R () [ G( ξ) G(0.5 ξ )] q( ξ ) dξ Ω The contact regon s dscreted nto several segments, wthn each segment the contact stress s assumed to e constant. B the least square method and pont matchng method[,3], we can fnd the contact stress dstrutons. (9) METHOD OF INITIAL FUNCTIONS In constructng the Green s functon for the contact prolems, the vertcal deflecton at top surface W n (h) due to a snusodal stress s needed. If the laer numer ncreases, t wll e tme consumng to fnd the asc deflecton mode. The method of ntal functon[5] s appled to smplf the calculaton of (h). From the stran defntons, the stress-stran relaton and W n equlrum equatons, the followng equaton can e derved[5]:
6 U V Z 0 A X B 0 Y W In whch c44 / d c45 / d 0 A c / / 45 d c55 d 0 0 U V Z X Y W () (3a) * c66 * B c6 c36 c33 c c * 6 * c c 3 33 c36 c 33 c3 c 33 / c33 (3) d c44c55 c45,,, c c * 3 j3 cj cj,, j,, 3, 6 (4) c33 X, Y, Z x The materal constants for the -th laer should e susttuted n Equatons () and (3). The soluton of Eqn. () s shown to e R( ) e Where 0 B A 0 R(0) R( ) T { U( ), V( ), Z( ), X( ), Y( ), W( )} (6) The state vector R() n the -th laer at artrar heght can e expressed n terms of the state vector at the ntal reference heght of that laer. B a successve multplcaton of the state vector at each laer, W n (h) can e found solvng a set of 6 6 smultaneous equaton. (5)
7 NUMERICAL EXAMPLES A [± 45] S composte lamnate made of Hexcel graphte/epox s analed as an llustratve example. The eam has a wdth of 0 cm and thcness of cm. Each lamna assumes the followng propertes, E 05.7GPa, E E3 8.57GPa, υ υ3 0.37, υ G G 4.39GPa, G 3.05GPa 3 3 The radus R of the rgd ndentor s 7.6cm. The through-the-thcness stress dstruton s calculated for contact length cmm. Fgure 3 shows the thcness-wse dstruton of the normal stress. The results at three dfferent postons (the contact center, the contact edge and end of eam) are compared. It seems that the normal stress wll have maxmum value at the top surface of the contact center. The nterlamnar stress dstrutons of and x are demonstrated n Fgs. 4 and 5. The solutons satsf the tracton free condtons at oth top and ottom surfaces. Orthotropc solutons are also compared, n whch the effectve elastc constants of ths angle-pl lamnate s calculated usng a 3-D effectve modul theor[6]. The orthotropc results are comparale to the laer-wse solutons. Fgure 6 shows the n-plane stress dstruton of. The relaton of ndentaton and contact force s evaluated n Fg. 7. The ndentaton s xx defned as the dfference of vertcal dsplacement of w(0.5, h) and w(0.5, 0). Fgure 8 s the contact pressure dstruton for dfferent contact length. In order to verf the applcalt of method of ntal functon, results of W n (h) for dfferent n are revealed n Fg. 9. It seems that oth methods agree well for ths case. In some other case, f n s large, the successve multplcaton of state vectors wll cause some numercal nstalt. If ths nstalt can e overcome, the method of ntal functon wll e an effcent tool for contact analss. CONCLUSIONS A method for the analss of ndentaton of angle-pl lamnates sujected to rgd clndrcal ndentor s proposed. Numercal results show that ths method satsfes the gven oundar condtons and tracton condtons. The contact pressure dstruton, magntude of contact force and ndentaton can e evaluated. Ths analtcal results wll compare to the expermental results n the near future.
8 REFERENCES. Sanar, B. V., Smooth Indentaton of Orthotropc Beams, Compostes Scence and Technolog, 34, (989) Sun, C. T. and Sanar, B. V., Smooth Indentaton of an ntall stressed orthotropc eam, Int. J. Solds and Structures,, (985), Wu, E., Chao, J. C. and Yen, C. S., Smooth Contact of Orthotropc Lamnates Rgd Clnders, AIAA Journal, 3, 0, (993), Jones, A. T., Exact Natural Frequences and Modal Functons for a Thc Off-Axs Lamna, J. of Composte Materals, 5, Oct. (97), Sundara Raja Iengar, K. T. and Panda, S. K., Analss of Orthotropc Rectangular Thc Plates, Fre Scence and Technolog, 8 (983), Sun, C. T. and Lee, S., Three-dmensonal Effectve Elastc Constants for Thc Lamnates, J. of Composte Materals, 988,,
Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
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 informationApplication to Plane (rigid) frame structure
Advanced Computatonal echancs 18 Chapter 4 Applcaton to Plane rgd frame structure 1. Dscusson on degrees of freedom In case of truss structures, t was enough that the element force equaton provdes onl
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 informationCOMPOSITE 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 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 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 informationUNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS
OCD0 UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 07/08 FINITE ELEMENT AND DIFFERENCE SOLUTIONS MODULE NO. AME6006 Date: Wednesda 0 Ma 08 Tme: 0:00
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 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 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 informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
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 informationCHAPTER 4. Vector Spaces
man 2007/2/16 page 234 CHAPTER 4 Vector Spaces To crtcze mathematcs for ts abstracton s to mss the pont entrel. Abstracton s what makes mathematcs work. Ian Stewart The man am of ths tet s to stud lnear
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 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 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 informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s a geometrcal concept arsng from parallel forces. Tus, onl parallel forces possess a centrod. Centrod s tougt of as te pont were te wole wegt of a pscal od or sstem of partcles
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 informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
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 informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s a geometrcal concept arsng from parallel forces. Tus, onl parallel forces possess a centrod. Centrod s tougt of as te pont were te wole wegt of a pscal od or sstem of partcles
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 informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
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 informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
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 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 informationNumerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method
Appled Mathematcs, 6, 7, 5-4 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton
More informationINDETERMINATE STRUCTURES METHOD OF CONSISTENT DEFORMATIONS (FORCE METHOD)
INETNTE STUTUES ETHO OF ONSISTENT EFOTIONS (FOE ETHO) If all the support reactons and nternal forces (, Q, and N) can not be determned by usng equlbrum equatons only, the structure wll be referred as STTIY
More informationThe Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites
7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT
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 informationCHAPTER 4 MAX-MIN AVERAGE COMPOSITION METHOD FOR DECISION MAKING USING INTUITIONISTIC FUZZY SETS
56 CHAPER 4 MAX-MIN AVERAGE COMPOSIION MEHOD FOR DECISION MAKING USING INUIIONISIC FUZZY SES 4.1 INRODUCION Intutonstc fuzz max-mn average composton method s proposed to construct the decson makng for
More informationCHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS COURSE CODES: FFR 35, FIM 72 GU, PhD Tme: Place: Teachers: Allowed materal: Not allowed: January 2, 28, at 8 3 2 3 SB
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
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 informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
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.
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: Instructor s Name and Secton: (Crcle Your Secton) Sectons:
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationEffects of Boundary Conditions on Cross-Ply Laminated Composite Beams
Internatonal Journal of Engneerng Research And Advanced Technology (IJERAT) DOI: http://dx.do.org/0.734/ijerat.344 E-ISSN : 454-635 Vol.3 (0) Oct -07 Effects of Boundary Condtons on Cross-Ply Lamnated
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationThe Study of Teaching-learning-based Optimization Algorithm
Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute
More informationPlate Theories for Classical and Laminated plates Weak Formulation and Element Calculations
Plate heores for Classcal and Lamnated plates Weak Formulaton and Element Calculatons PM Mohte Department of Aerospace Engneerng Indan Insttute of echnolog Kanpur EQIP School on Computatonal Methods n
More informationNumerical analysis of buried pipes subjected to lateral soil movements
Numercal analss of bured ppes subjected to lateral sol movements P. Casamchele, M. Mauger & E. Motta Department of Cvl and Envronmental Engneerng, Catana Unverst, Ital Abstract he response of bured ppes
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationSMOOTHED PARTICLE HYDRODYNAMICS METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM
The 5th Internatonal Symposum on Computatonal Scences (ISCS) 28 May 2012, Yogyakarta, Indonesa SMOOTHED PARTICLE HYDRODYNAMICS METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM Dede Tarwd 1,2 1 Graduate School
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
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 informationNovember 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project
SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationBuckling analysis of single-layered FG nanoplates on elastic substrate with uneven porosities and various boundary conditions
IOSR Journal of Mechancal and Cvl Engneerng (IOSR-JMCE) e-issn: 78-1684,p-ISSN: 30-334X, Volume 15, Issue 5 Ver. IV (Sep. - Oct. 018), PP 41-46 www.osrjournals.org Bucklng analyss of sngle-layered FG nanoplates
More informationImportant Instructions to the Examiners:
Summer 0 Examnaton Subject & Code: asc Maths (70) Model Answer Page No: / Important Instructons to the Examners: ) The Answers should be examned by key words and not as word-to-word as gven n the model
More informationA new method of boundary treatment in Heat conduction problems with finite element method
Avalale onlne at www.scencedrect.com Proceda Engneerng 6 (11) 675 68 Frst Internatonal Smposum on Mne Safet Scence and Engneerng A new method of oundar treatment n Heat conducton prolems wth fnte element
More informationPreliminary Design of Moment-Resisting Frames
Prelmnary Desgn of Moment-Resstng Frames Preprnt Aamer Haque Abstract A smple method s developed for prelmnary desgn of moment-resstng frames. Preprnt submtted to Elsever August 27, 2017 1. Introducton
More informationBallot Paths Avoiding Depth Zero Patterns
Ballot Paths Avodng Depth Zero Patterns Henrch Nederhausen and Shaun Sullvan Florda Atlantc Unversty, Boca Raton, Florda nederha@fauedu, ssull21@fauedu 1 Introducton In a paper by Sapounaks, Tasoulas,
More informationNote: Please use the actual date you accessed this material in your citation.
MIT OpenCourseWare http://ocw.mt.edu 6.13/ESD.13J Electromagnetcs and Applcatons, Fall 5 Please use the followng ctaton format: Markus Zahn, Erch Ippen, and Davd Staeln, 6.13/ESD.13J Electromagnetcs and
More information8.592J: Solutions for Assignment 7 Spring 2005
8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationA boundary element method with analytical integration for deformation of inhomogeneous elastic materials
Journal of Physcs: Conference Seres PAPER OPEN ACCESS A boundary element method wth analytcal ntegraton for deformaton of nhomogeneous elastc materals To cte ths artcle: Moh. Ivan Azs et al 2018 J. Phys.:
More informationResource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis
Resource Allocaton and Decson Analss (ECON 800) Sprng 04 Foundatons of Regresson Analss Readng: Regresson Analss (ECON 800 Coursepak, Page 3) Defntons and Concepts: Regresson Analss statstcal technques
More informationAPPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS
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
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 informationCHAPTER-5 INFORMATION MEASURE OF FUZZY MATRIX AND FUZZY BINARY RELATION
CAPTER- INFORMATION MEASURE OF FUZZY MATRI AN FUZZY BINARY RELATION Introducton The basc concept of the fuzz matr theor s ver smple and can be appled to socal and natural stuatons A branch of fuzz matr
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 informationMODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava
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
More informationORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp.
Clck to Vew Mathcad Document 2011 Knovel Corp. Buldng Structural Desgn. homas P. Magner, P.E. 2011 Parametrc echnology Corp. Chapter 3: Renforced Concrete Slabs and Beams 3.2 Renforced Concrete Beams -
More informationNumerical Simulation of Wave Propagation Using the Shallow Water Equations
umercal Smulaton of Wave Propagaton Usng the Shallow Water Equatons Junbo Par Harve udd College 6th Aprl 007 Abstract The shallow water equatons SWE were used to model water wave propagaton n one dmenson
More informationPHYS 1443 Section 003 Lecture #17
PHYS 144 Secton 00 ecture #17 Wednesda, Oct. 9, 00 1. Rollng oton of a Rgd od. Torque. oment of Inerta 4. Rotatonal Knetc Energ 5. Torque and Vector Products Remember the nd term eam (ch 6 11), onda, Nov.!
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationx yi In chapter 14, we want to perform inference (i.e. calculate confidence intervals and perform tests of significance) in this setting.
The Practce of Statstcs, nd ed. Chapter 14 Inference for Regresson Introducton In chapter 3 we used a least-squares regresson lne (LSRL) to represent a lnear relatonshp etween two quanttatve explanator
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 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 informationLecture 8 Modal Analysis
Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal
More informationComputers and Mathematics with Applications
Computers and Mathematcs wth Applcatons 56 (2008 3204 3220 Contents lsts avalable at ScenceDrect Computers and Mathematcs wth Applcatons journal homepage: www.elsever.com/locate/camwa An nnovatve egenvalue
More informationDynamic Programming. Preview. Dynamic Programming. Dynamic Programming. Dynamic Programming (Example: Fibonacci Sequence)
/24/27 Prevew Fbonacc Sequence Longest Common Subsequence Dynamc programmng s a method for solvng complex problems by breakng them down nto smpler sub-problems. It s applcable to problems exhbtng the propertes
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
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 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 informationAGC Introduction
. Introducton AGC 3 The prmary controller response to a load/generaton mbalance results n generaton adjustment so as to mantan load/generaton balance. However, due to droop, t also results n a non-zero
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 informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationSolution 1 for USTC class Physics of Quantum Information
Soluton 1 for 018 019 USTC class Physcs of Quantum Informaton Shua Zhao, Xn-Yu Xu and Ka Chen Natonal Laboratory for Physcal Scences at Mcroscale and Department of Modern Physcs, Unversty of Scence and
More informationDEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica
demo8.nb 1 DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA Obectves: - defne matrces n Mathematca - format the output of matrces - appl lnear algebra to solve a real problem - Use Mathematca to perform
More informationA comprehensive study: Boundary conditions for representative volume elements (RVE) of composites
Insttute of Structural Mechancs A comprehensve study: Boundary condtons for representatve volume elements (RVE) of compostes Srhar Kurukur A techncal report on homogenzaton technques A comprehensve study:
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 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 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 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 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 informationRichard Socher, Henning Peters Elements of Statistical Learning I E[X] = arg min. E[(X b) 2 ]
1 Prolem (10P) Show that f X s a random varale, then E[X] = arg mn E[(X ) 2 ] Thus a good predcton for X s E[X] f the squared dfference s used as the metrc. The followng rules are used n the proof: 1.
More informationΔ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions
ME67 - Handout 4 Vbratons of Contnuous Systems Axal vbratons of elastc bars The fgure shows a unform elastc bar of length and cross secton A. The bar materal propertes are ts densty ρ and elastc modulus
More informationFlow Induced Vibration
Flow Induced Vbraton Project Progress Report Date: 16 th November, 2005 Submtted by Subhrajt Bhattacharya Roll no.: 02ME101 Done under the gudance of Prof. Anrvan Dasgupta Department of Mechancal Engneerng,
More informationAssessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion
Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,
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 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 information