Modeling of Finite Gradient Elasticity Using Optimal Transportation Method
|
|
- Claire Powers
- 5 years ago
- Views:
Transcription
1 9 th Iteratioal Cogress of Croatia Society of Mechaics September 2018 Split, Croatia Modelig of Fiite Gradiet Elasticity Usig Optimal Trasportatio Method Boris JALUŠIĆ *, Christia WEIßENFELS +, Jurica SORIĆ *, eter WRIGGERS + * Faculty of Mechaical Egieerig ad Naval Architecture, Uiversity of Zagreb, Ivaa Lučića 5, Zagreb, Croatia. s: {boris.jalusic, jurica.soric}@fsb.hr + Istitut für Kotiuumsmechaik, Leibiz Uiversität Haover, Appelstraße 11, Haover, Germay. {weissefels, wriggers}@ikm.ui-haover.de Abstract. A meshfree Optimal Trasportatio Method (OTM) is applied for modelig of gradiet hyperelasticity usig higher-order theory based o oly oe microstructural parameter. The OTM method is utilized for the purpose of structural behaviour modelig udergoig large deformatio. The idepedet displacemet variables are approximated usig iterpolatig maximum-etropy fuctios. This eables the satisfactio of the boudary coditios i a simple ad straightforward maer without the eed for the calculatio of additioal parameters. Novel stabilizatio algorithm employig oe material poit is preseted. The aalysis of umerical performace of the proposed approach is demostrated by a represetative umerical bechmark example. The obtaied results are commeted ad arise umerical problems are disscussed. 1 Itroductio At the momet, a wide variety of differet meshfree methods are beig employed for modelig of various egieerig problems. The reaso behid this applicatio are beefits ad advatages i compariso with stadard mesh-based methods, like Fiite Elemet Method (FEM) [1]. The meshless umerical approaches are able to overcome problems such as elemet distortio ad time-demadig mesh geeratio process. Nevertheless, the umerical itegratio ad calculatio of meshfree shape fuctios ad its derivatives still represet a major task due to the associated high computatioal costs [2]. I the preset cotributio the OTM formulatio [3] is adapted for modelig of deformatio of homogeeous structures employig strai gradiet elasticity theory [2]. I additio, the strai gradiet elasticity based o the theory with oly oe microstructural parameter is cosidered. The gradiet theory is utilized i order to more accurately capture the material behaviour, to remove discotiuities from strai ad stress fields ad to simulate the size effects. Size effects deped o the microstructure of the material ad ca be observed whe the size of the specime is sufficietly small, i.e. approximately the size of the microstructure. At this legth scale, specimes with similar shape but differet dimesios show differet mechaical behaviors. The solutio of fourth-order differetial equatios arisig i cosidered theories requires a high-order cotiuity of approximatio fuctios. Hece, usig the FEM for solvig this
2 type of problems is ot a wise choice sice stadard formulatios eed to possess C 1 cotiuity, which leads to complicated shape fuctios with large umber of odal degrees of freedom [4], eve if mixed elemets are utilized [5]. I compariso, the required cotiuity i the meshfree methods is obtaiable i a simple by addig more odes withi the approximatio domai [6]. The applied OTM method is based o the combiatio of the optimal trasportatio theory [7] with material-poit samplig usig the max-etropy meshless approximatio [8]. It eforces the mass trasport ad essetial boudary coditios exactly ad should be free from tesio istabilities. Furthermore, the method coserves liear ad agular mometum [3]. The ivestigated global domai is approximated with two differet sets of poits, the discretizatio odes ad material poits. At the odes the displacemets as the primary variables are computed by solvig the discretized equatio of motio. O the other had, at the material poits data associated with material is determied. This icludes desity, strai, stress, secod-order stress, etc. Furthermore, the material poits are used as itegratio poits. The algorithms are based o the explicit time itegratio. I every time step the positio of discretizatio poits ad the material poits is calculated usig updated Lagragia scheme. The OTM algorithm is combied with the stadard meshfree search algorithm where a arbitrary umber of odes ca be assiged to a material poit. Thus, the shape fuctios ca become highly oliear if a larger umber of odes are located withi the approximatio domai. For the itegratio of the goverig equatios oly oe material poit is used as i the origially proposed paper [3] to lower the computatioal time effort. It has bee observed that oe poit is eough whe solvig the problem usig classical cotiuum theory [9], but ot might be whe o-local (higher-order) cotiuum equatios is employed. Hece, curretly the ovel stabilizatio algorithm usig pealty method similar to oe i [9] is beig developed that could alleviate the methods stability issues. The OTM framework for the modelig of deformatio resposes of homogeeous material usig gradiet elasticity is preseted ad explaied at large i Sectio 2. Herei, the cosidered stabilizatio algorithm is also preseted. I Sectio 3. the proposed stabilizatio algorithm is tested o a problem of the truss subjected to tesile test usig uiform displacemets at both eds. I Sectio 4 cocludig remarks ad further research guidelies are give. 2 Optimal Trasportatio Method for Fiite Gradiet Elasticity 2.1 Goverig equatios for a higher-order oliear cotiuum The three dimesioal higher-order oliear cotiuum which occupies the global computatioal domai surrouded by the global outer boudary is cosidered. Accordig to [10], the goverig equatio is the weak form of the dyamic equilibrium writte as ηdiv(div ) dv ηdiv dv ηbdv ηφ d V =0, (1) V V V V 2
3 where is the desity ad φ defies the secod-order time derivative of deformatio map (positio vector φ( X, t) x i the curret cofiguratio) defied as φ 12φ φ 1 φ. (2) 2 t Furthermore, the displacemet vector is defied ad calculated as ux (, t) φ( X, t) X. I relatio (1) η deotes the admissible test fuctio, the iola stress tesor ad the double stress tesor. The tesors ad are defied as the partial derivatives of free eergy desity W W ( F) W G. The free eergy desity is depedet o the first-order deformatio gradiet F φ ad the secod-order deformatio gradiet x G φ. Thus, the secod- ad third-order stress tesors are computed as x x W, W. (3) F G The utilized free eergy desity fuctio correspods to the Neo-Hooke material used i [11] ad is take as 1 1 F F F, (4) 2 2 l : 3 2l W 2 J J W G l G G, (5) 2 where ad deote the Lame`s costats, J the determiat of F ad l the parameter associated with the size of the microstructure of the higher-order cotiuum. I order to derive the appropriate equatio to discretize the relatio (1) is firstly itegrated by parts (the first part two times), after that the divergece theorem is applied ad the decompositio of the gradiet operator is doe to properly idetify the higherorder boudary coditios [10] associated with the used cotiuum leadig to ( η)dv η t d ηt d ηbdv ηφ d V =0, (6) x x x V V V where t ad t represet the stadard tractio vector ad the higher-order tractio vector, respectively. The weak form i (6) has to fulfill the essetial boudary coditios u u, (7) u u, (8) x x 3
4 ad atural boudary coditios t =( Div ) N xn( I N N) : I : ( N N) ( x( N) ( I N N)) : I t (9) t N N t, (10) = : outer boudary. 2.2 Optimal Trasportatio Meshfree algorithm I compariso to classical meshfree methods like Elemet Free Galerki (EFG) [12] method ot oly the positio of odes but also the positios of material poits are updated durig calculatio withi the updated Lagragia framework. I this cotributio the OTM algorithm is derived usig the cetral differece time itegratio scheme (2) i the preseted weak form (6). Furthermore, the domai aroud each material poit is defied as a port (itegratio) domai usig the appropriate search algorithm. At every material poit the test fuctio η ad the deformatio map φ are approximated usig the max-etropy shape fuctios [8] ad odal values N p ηp N I x p η I, N I 1 N p φ x φ. (11) p I p I I 1 Herei N p is the umber of odes i the curret port domai for each material poit. Usig the approximatios (11) equatio (6) is trasformed to a algebraic equatio M 2 I 2 I ( t) φ F G R ( t) φ. (12) I relatio (12) the icremets used are defied as φ φ φ ad -1 φ φ φ. The stadard ad higher-order forces at the outer boudary are give as ad, equal to mp N p N ( x ) t v, (13) p1 I 1 I p p mp N p N ( x ) N t v, (14) p1 I 1 while the global mass matrix I p p M ad residual vectors R, F ad G are 4
5 mp Np Np M N ( x ) I N ( x ) m, (15) p1 I1 J1 I p J p p mp N p R N ( x ) b m, (16) p1 I 1 I p p mp N p F B ( x ) v, (17) p1 I 1 I p p mp N p G H ( x ) : v. (18) p1 I 1 I p p I relatios (13) - (18) the m p represets the mass ad v p the volume at each material poit. The matrices B I ad H I cotai the first- ad secod-order derivatives of the +1 shape fuctios of ode I. By solvig (12) the icremetal positio vector φ I is determied ad the odal positio vector is updated φ φ φ. (19) I I I Furthermore, i order to reduce the computatioal time the lumped mass cocept [13] withi the OTM framework is utilized. 2.3 Stabilizatio procedure For the purpose of preservig the same computatioal time ad icreasig the accuracy of the method the ovel stabilizatio procedure is proposed. The procedure is based o method of pealizig the egative effects resultig from uderitegratio [14]. The developed stabilizatio algorithm is based o the assumptio of oly quadratic distributio of displacemets withi the ifluece domai. This is doe by computig the ormalized error at every material poit iside the ifluece domai. The calculated error is the multiplied with the appropriate pealty parameter ad subtracted from the odal residual vector r I assembled aalogously usig the vectors R, F ad G but utilizig the material poits withi the domai of ifluece. The ormalized error criterio for every poit p withi the ifluece domai of the ode I is e Ip x x x x dx dx I p I p = d X X X I 1 p 1, (20) where the distace vectors calculated from the OTM method are defied as 5
6 1 1 : x I p p X I p p I p X I p (21) I order to eforce the error to be zero the pealty method is applied ad the stabilized odal residual vector is computed as if N I stab I I N I p Ip p1 r r x e. (22) 3 Numerical example 3.1 Truss subjected to large displacemets A three-dimesioal truss subjected to mixed boudary coditios as show i Figure 2 is cosidered. The truss i subjected to large displacemets at the frot ad back side, while the other compoets of boudary coditios are the stadard tractios ad higher-order trasctio which are take as zero. The geometry of the truss is defied by the legth L 1. For the aalysis of deformatio the bechmark material properties 6 6 of the truss are take as E 1 10, 0.3 ad The large desity is chose here i order to compute the example faster i a quasi-static case to test the stabilizatio algorithm for the higher-order cotiuum. For the purpose of comparig the solutios, firstly the deformatio of the geometrically idetical truss uder the assumpsio of the classical cotium without ay microstructural effects icluded is computed. I Figures 3 ad 5, the deformed shapes of the truss with the cotours of the displacemets i the z directio ad the Vo Mises stress are show. As see from the deformed shape i Figure 3, the truss at both eds reached the fial values of the displacemet boudary coditios. Figure1: Truss with boudary coditios Secodly, a higher-order cotiuum with a very small value of the microstructural parameter of l/ L 0.05 is applied i computatios. Here it should be metioed that by itroducig the higher-order cotiuum ad withi that the calculatio of secod- 6
7 order derivatives, o solutio ca be achieved without the applicatio of stabilizatio procedure. Hece, the uder itegratio of the weak form has a very large effect o the deformatio of the structure, umerical calculatios ca ot be performed because of very large errors arisig durig computatio process. Figure2: Truss classical cotiuum distributio of displacemet u z Figure4: Truss classical cotiuum distributio of Vo Mises stress Figure3: Truss higher-order cotiuum distributio of displacemet u z Figure5: Truss higher-order cotiuum distributio of Vo Mises stress I Figures 3 ad 5, the obtaied deformed shapes of the higher-order cotiuum with very high value of the stabilizatio parameter are portrayed. By the fial positios of the odes, it is clear that the applicatio of the preseted stabilizatio procedure causes the destructio of the displacemet boudary coditios ad deformed cotiuum forms which are ot accurate. Hece, here additioal umerical studies are eccesary i order to determie why the metioed problems persist. Furthermore, as a alterative to max-etropy approximatio fuctios differet more robust approximatio fuctios should be take i cosideratio ad the obtaied results should be compared. 4 Coclusio The meshfree OTM formulatio has bee applied for the modelig of deformatio of the homogeeous structure. The ifluece of the material microstructure is take ito accout usig oe parameter withi the free eergy desity of the utilized higher-order cotiuum. Furthermore, as a remedy to the methods stability issues coected to uderitegratio of the weak form ovel stabilizatio procedure is employed. The OTM algorithm is tested o oe bechmark truss example. Usig the preseted stabilizatio algorithm, solutios ca be obtaied, however some problems have bee idetified. The first problem is the determiatio of the value of the pealty parameter if higher-order 7
8 cotiuum is computed. It is really hard to choose the value of the pealty due to the icorrect itegratio of the secod-order derivatives of the max-etropy shape fuctios that are complex ad their values is ot easy to cotrol. The secod problem observed is that the stabilizatio procedure destroys the essetial boudary coditios. I further research the metioed issues should be addressed. The chage to the more robust approximatio fuctios should be cosidered. This could improve the smoothess of secod-order derivatives ad more accurate calculatios should be carried out. Ackowledgemet: This work has bee ported i part by Croatia Sciece Foudatio uder the project The ivestigatios are also a part of the project ported by Alexader vo Humboldt Foudatio. Refereces [1]. Wriggers, Noliear Fiite Elemet Methods, Spriger Sciece & Busiess Media2008. [2] H. Askes, E.C. Aifatis, Gradiet elasticity i statics ad dyamics: A overview of formulatios, legth scale idetificatio procedures, fiite elemet implemetatios ad ew results, Iteratioal Joural of Solids ad Structures, 48 (2011) [3] B. Li, F. Habbal, M. Ortiz, Optimal trasportatio meshfree approximatio schemes for fluid ad plastic flows, Iteratioal Joural for Numerical Methods i Egieerig, 83 (2010) [4] S.-A. apaicolopulos, A. Zervos, I. Vardoulakis, A three-dimesioal C1 fiite elemet for gradiet elasticity, Iteratioal Joural for Numerical Methods i Egieerig, 77 (2009) [5] E. Amaatidou, N. Aravas, Mixed fiite elemet formulatios of strai-gradiet elasticity problems, Computer Methods i Applied Mechaics ad Egieerig, 191 (2002) [6] Z. Tag, S. She, S.N. Atluri, Aalysis of materials with strai-gradiet effects: A Meshless Local etrov-galerki(mlg) approach, with odal displacemets oly, Cmes-Computer Modelig i Egieerig & Scieces, 4 (2003) [7] C. Villai, Topics i optimal trasportatio theory, Graduate Studies i Mathematics, America Mathematical Society, rovidece, Rhode Islad, USA, [8] M. Arroyo, M. Ortiz, Local maximum-etropy approximatio schemes: a seamless bridge betwee fiite elemets ad meshfree methods, Iteratioal Joural for Numerical Methods i Egieerig, 65 (2006) [9] C. Weißefels,. Wriggers, Stabilizatio algorithm for the optimal trasportatio meshfree approximatio scheme, Computer Methods i Applied Mechaics ad Egieerig, 329 (2018) [10] R. Suyk,. Steima, O higher gradiets i cotiuum-atomistic modellig, Iteratioal Joural of Solids ad Structures, 40 (2003) [11].E. Fischer, C1 Cotiuous Methods i Computatioal Gradiet Elasticity, Uiversitat Erlage-Nürberg, Germay, [12] T. Belytschko, Y.Y. Lu, L. Gu, Elemet-free Galerki methods, Iteratioal Joural for Numerical Methods i Egieerig, 37 (1994) [13] L. Fedeli, A. adolfi, M. Ortiz, Geometrically-exact time-itegratio mesh-free schemes for advectio-diffusio problems derived from optimal trasportatio theory ad their coectio with particle methods, [14] G.C. Gazemüller, A hourglass cotrol algorithm for Lagragia Smooth article Hydrodyamics, Computer Methods i Applied Mechaics ad Egieerig, 286 (2015)
Analysis of composites with multiple rigid-line reinforcements by the BEM
Aalysis of composites with multiple rigid-lie reiforcemets by the BEM Piotr Fedeliski* Departmet of Stregth of Materials ad Computatioal Mechaics, Silesia Uiversity of Techology ul. Koarskiego 18A, 44-100
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More informationCO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS
CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,
More informationA numerical Technique Finite Volume Method for Solving Diffusion 2D Problem
The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed
More informationMETHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS
Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationBoundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary
Boudary layer problem o coveyor belt Gabriella Bogár Uiversity of Miskolc 355 Miskolc-Egyetemváros, Hugary e-mail: matvbg@ui-miskolc.hu Abstract: A techologically importat source of the boudary layer pheomeo
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationDECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan
Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed
More informationNumerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem
Leoardo Joural of Scieces ISSN 1583-0233 Issue 9, July-December 2006 p. 25-32 Numerical Simulatio of Thermomechaical Problems i Applied Mechaics: Applicatio to Solidificatio Problem Vicet Obiajulu OGWUAGWU
More informationPAijpam.eu ON TENSOR PRODUCT DECOMPOSITION
Iteratioal Joural of Pure ad Applied Mathematics Volume 103 No 3 2015, 537-545 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://wwwijpameu doi: http://dxdoiorg/1012732/ijpamv103i314
More informationDYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS
DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS Ivaa Štimac 1, Ivica Kožar 1 M.Sc,Assistat, Ph.D. Professor 1, Faculty of Civil Egieerig, Uiverity of Rieka, Croatia INTRODUCTION The vehicle-iduced
More informationNumerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets
Australia Joural of Basic ad Applied Scieces, 5(): 98-5, ISSN 99-878 Numerical Solutio of the Two Poit Boudary Value Problems By Usig Wavelet Bases of Hermite Cubic Splie Wavelets Mehdi Yousefi, Hesam-Aldie
More informationMATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS
6 th INTERNATIONAL MULTIDISCIPLINARY CONFERENCE MATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS Istva eppler SZIE Faculty of Mechaics, H-2103 Gödöllő Páter. 1., Hugary Abstract: The mathematical
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationTHE KALMAN FILTER RAUL ROJAS
THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider
More informationNUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS
THERMAL SCIENCE, Year 07, Vol., No. 4, pp. 595-599 595 NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS by Yula WANG *, Da TIAN, ad Zhiyua LI Departmet of Mathematics,
More informationMicroscopic traffic flow modeling
Chapter 34 Microscopic traffic flow modelig 34.1 Overview Macroscopic modelig looks at traffic flow from a global perspective, whereas microscopic modelig, as the term suggests, gives attetio to the details
More informationMonte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem
Australia Joural of Basic Applied Scieces, 5(): 097-05, 0 ISSN 99-878 Mote Carlo Optimizatio to Solve a Two-Dimesioal Iverse Heat Coductio Problem M Ebrahimi Departmet of Mathematics, Karaj Brach, Islamic
More informationExact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method
Exact Solutios for a Class of Noliear Sigular Two-Poit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationDamped Vibration of a Non-prismatic Beam with a Rotational Spring
Vibratios i Physical Systems Vol.6 (0) Damped Vibratio of a No-prismatic Beam with a Rotatioal Sprig Wojciech SOCHACK stitute of Mechaics ad Fudametals of Machiery Desig Uiversity of Techology, Czestochowa,
More informationNumerical Method for Blasius Equation on an infinite Interval
Numerical Method for Blasius Equatio o a ifiite Iterval Alexader I. Zadori Omsk departmet of Sobolev Mathematics Istitute of Siberia Brach of Russia Academy of Scieces, Russia zadori@iitam.omsk.et.ru 1
More informationTHE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS
R775 Philips Res. Repts 26,414-423, 1971' THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS by H. W. HANNEMAN Abstract Usig the law of propagatio of errors, approximated
More informationA STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD
IRET: Iteratioal oural of Research i Egieerig ad Techology eissn: 39-63 pissn: 3-7308 A STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD Satish
More informationA new formulation of internal forces for non-linear hypoelastic constitutive models verifying conservation laws
A ew formulatio of iteral forces for o-liear hypoelastic costitutive models verifyig coservatio laws Ludovic NOELS 1, Lauret STAINIER 1, Jea-Philippe PONTHOT 1 ad Jérôme BONINI. 1 LTAS-Milieux Cotius et
More informationPrinciple Of Superposition
ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give
More informationResearch Article A New Second-Order Iteration Method for Solving Nonlinear Equations
Abstract ad Applied Aalysis Volume 2013, Article ID 487062, 4 pages http://dx.doi.org/10.1155/2013/487062 Research Article A New Secod-Order Iteratio Method for Solvig Noliear Equatios Shi Mi Kag, 1 Arif
More informationOrthogonal Gaussian Filters for Signal Processing
Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios
More informationA widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α
Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows
More information2C09 Design for seismic and climate changes
2C09 Desig for seismic ad climate chages Lecture 02: Dyamic respose of sigle-degree-of-freedom systems I Daiel Grecea, Politehica Uiversity of Timisoara 10/03/2014 Europea Erasmus Mudus Master Course Sustaiable
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationTrue Nature of Potential Energy of a Hydrogen Atom
True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial
More informationDiscrete Orthogonal Moment Features Using Chebyshev Polynomials
Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical
More information3 Balance equations ME338A CONTINUUM MECHANICS
ME338A CONTINUUM MECHANICS lecture otes 1 thursy, may 1, 28 Basic ideas util ow: kiematics, i.e., geeral statemets that characterize deformatio of a material body B without studyig its physical cause ow:
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationJournal of Computational Physics 149, (1999) Article ID jcph , available online at
Joural of Computatioal Physics 149, 418 422 (1999) Article ID jcph.1998.6131, available olie at http://www.idealibrary.com o NOTE Defiig Wave Amplitude i Characteristic Boudary Coditios Key Words: Euler
More informationDynamical stability and parametrical vibrations of the laminated plates with complex shape
(3) 75 88 Dyamical stability ad parametrical vibratios of the lamiated plates with comple shape Abstract The problem of oliear vibratios ad stability aalysis for the symmetric lamiated plates with comple
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationModified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations
Noliear Aalysis ad Differetial Equatios, Vol. 5, 27, o. 4, 57-7 HIKARI Ltd, www.m-hikari.com https://doi.org/.2988/ade.27.62 Modified Decompositio Method by Adomia ad Rach for Solvig Noliear Volterra Itegro-
More informationLecture 19: Convergence
Lecture 19: Covergece Asymptotic approach I statistical aalysis or iferece, a key to the success of fidig a good procedure is beig able to fid some momets ad/or distributios of various statistics. I may
More informationInformation-based Feature Selection
Iformatio-based Feature Selectio Farza Faria, Abbas Kazeroui, Afshi Babveyh Email: {faria,abbask,afshib}@staford.edu 1 Itroductio Feature selectio is a topic of great iterest i applicatios dealig with
More informationRandom Walks on Discrete and Continuous Circles. by Jeffrey S. Rosenthal School of Mathematics, University of Minnesota, Minneapolis, MN, U.S.A.
Radom Walks o Discrete ad Cotiuous Circles by Jeffrey S. Rosethal School of Mathematics, Uiversity of Miesota, Mieapolis, MN, U.S.A. 55455 (Appeared i Joural of Applied Probability 30 (1993), 780 789.)
More informationThe Analysis of the Non-linear Deflection of Non-straight Ludwick type Beams Using Lie Symmetry Groups
Proceedigs of the 3 rd Iteratioal Coferece o Cotrol, Dyamic Systems, ad Robotics CDSR 16 Ottawa, Caada May 9 10, 016 Paper No. 107 DOI: 10.11159/cdsr16.107 The Aalysis of the No-liear Deflectio of No-straight
More informationOn forward improvement iteration for stopping problems
O forward improvemet iteratio for stoppig problems Mathematical Istitute, Uiversity of Kiel, Ludewig-Mey-Str. 4, D-24098 Kiel, Germay irle@math.ui-iel.de Albrecht Irle Abstract. We cosider the optimal
More informationTime-Domain Representations of LTI Systems
2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable
More informationAER Benchmark Specification Sheet
AER Bechmark Specificatio Sheet. Test ID: AER-DYN-002 2. Short Descriptio: A three-dimesioal hexagoal dyamic bechmark describig the ejectio of a peripheral cotrol rod i a VVER-440 core with eutro kietics
More informationCS537. Numerical Analysis and Computing
CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be
More informationNumerical Study on MHD Flow And Heat Transfer With The Effect Of Microrotational Parameter In The Porous Medium
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 4 (April. 5), V PP 8-7 www.iosrje.org Numerical Study o MHD Flow Ad Heat rasfer With he Effect Of Microrotatioal Parameter
More informationUsing Spreadsheets as a Computational Tool in Teaching Mechanical. Engineering
Proceedigs of the th WSEAS Iteratioal Coferece o COMPUTERS, Vouliagmei, Athes, Greece, July 3-5, 6 (pp35-3) Usig Spreadsheets as a Computatioal Tool i Teachig Mechaical Egieerig AHMADI-BROOGHANI, ZAHRA
More informationNumerical simulation of a climate energy balance model with continents distribution.
Numerical simulatio of a climate eergy balace model with cotiets distributio. Lourdes Tello & Arturo Hidalgo Uiversidad Politécica de Madrid. Physical problem Some processes ivolved i global climate models:
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationAlgebra of Least Squares
October 19, 2018 Algebra of Least Squares Geometry of Least Squares Recall that out data is like a table [Y X] where Y collects observatios o the depedet variable Y ad X collects observatios o the k-dimesioal
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationAvoidance of numerical singularities in free vibration analysis of Euler-Bernoulli beams using Green functions
WSEAS TRASACTIOS o APPLIED ad THEORETICAL MECHAICS Goraka Štimac Ročević, Braimir Ročević, Ate Skoblar, Saji Braut Avoidace of umerical sigularities i free vibratio aalysis of Euler-Beroulli beams usig
More information8. Applications To Linear Differential Equations
8. Applicatios To Liear Differetial Equatios 8.. Itroductio 8.. Review Of Results Cocerig Liear Differetial Equatios Of First Ad Secod Orders 8.3. Eercises 8.4. Liear Differetial Equatios Of Order N 8.5.
More informationGUIDELINES ON REPRESENTATIVE SAMPLING
DRUGS WORKING GROUP VALIDATION OF THE GUIDELINES ON REPRESENTATIVE SAMPLING DOCUMENT TYPE : REF. CODE: ISSUE NO: ISSUE DATE: VALIDATION REPORT DWG-SGL-001 002 08 DECEMBER 2012 Ref code: DWG-SGL-001 Issue
More informationTHE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT
Europea Joural of Egieerig ad Techology Vol. 3 No., 5 ISSN 56-586 THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE Atif Nazir, Tahir Mahmood ad
More informationNumerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION
alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty
More informationLecture 2: Finite Difference Methods in Heat Transfer
Lecture 2: Fiite Differece Methods i Heat Trasfer V.Vuorie Aalto Uiversity School of Egieerig Heat ad Mass Trasfer Course, Autum 2016 November 1 st 2017, Otaiemi ville.vuorie@aalto.fi Overview Part 1 (
More informationPoincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains
Advaces i Pure Mathematics 23 3 72-77 http://dxdoiorg/4236/apm233a24 Published Olie Jauary 23 (http://wwwscirporg/oural/apm) Poicaré Problem for Noliear Elliptic Equatios of Secod Order i Ubouded Domais
More informationA Block Cipher Using Linear Congruences
Joural of Computer Sciece 3 (7): 556-560, 2007 ISSN 1549-3636 2007 Sciece Publicatios A Block Cipher Usig Liear Cogrueces 1 V.U.K. Sastry ad 2 V. Jaaki 1 Academic Affairs, Sreeidhi Istitute of Sciece &
More informationOn Strongly Consistent Finite Dierence Approximations
D.Michels et al. (KAUST,JINR,SSU) Strogly Cosistet Approximatios 17 April 2018 1 / 25 O Strogly Cosistet Fiite Dierece Approximatios Domiik Michels 1, Vladimir Gerdt 2, Dmitry Lyakhov 1, ad Yuri Blikov
More informationFinite Element Analysis of Rubber Bumper Used in Air-springs
Available olie at www.sciecedirect.com Procedia Egieerig 48 (0 ) 388 395 MMaMS 0 Fiite Elemet Aalysis of Rubber Bumper Used i Air-sprigs amas Makovits a *, amas Szabó b a Departmet of Mechaical Egieerig,
More informationA general finite element model for numerical simulation of structure dynamics
440 Sciece i Chia Series G: Physics, Mechaics & Astroomy 2006 Vol.49 No.4 440 450 DOI: 10.1007/s11433-006-2011-1 A geeral fiite elemet model for umerical simulatio of structure dyamics WANG Fuju 1, LI
More informationL 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!
aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (
More informationGoal. Adaptive Finite Element Methods for Non-Stationary Convection-Diffusion Problems. Outline. Differential Equation
Goal Adaptive Fiite Elemet Methods for No-Statioary Covectio-Diffusio Problems R. Verfürth Ruhr-Uiversität Bochum www.ruhr-ui-bochum.de/um1 Tübige / July 0th, 017 Preset space-time adaptive fiite elemet
More informationL = n i, i=1. dp p n 1
Exchageable sequeces ad probabilities for probabilities 1996; modified 98 5 21 to add material o mutual iformatio; modified 98 7 21 to add Heath-Sudderth proof of de Fietti represetatio; modified 99 11
More informationInelastic spherical caps with defects
Ielastic spherical caps with defects JAAN LLLP RNST TUNGL Istitute of athematics Uiversity of Tartu Liivi str Tartu 549 STONIA aalellep@utee Tartu College Talli Uiversity of Techology Puiestee 78 58 Tartu
More informationANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION
Molecular ad Quatum Acoustics vol. 7, (6) 79 ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION Jerzy FILIPIAK 1, Lech SOLARZ, Korad ZUBKO 1 Istitute of Electroic ad Cotrol Systems, Techical Uiversity of Czestochowa,
More informationLinear Elliptic PDE s Elliptic partial differential equations frequently arise out of conservation statements of the form
Liear Elliptic PDE s Elliptic partial differetial equatios frequetly arise out of coservatio statemets of the form B F d B Sdx B cotaied i bouded ope set U R. Here F, S deote respectively, the flux desity
More informationChapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics:
Chapter 6 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals (which is what most studets
More informationA New Solution Method for the Finite-Horizon Discrete-Time EOQ Problem
This is the Pre-Published Versio. A New Solutio Method for the Fiite-Horizo Discrete-Time EOQ Problem Chug-Lu Li Departmet of Logistics The Hog Kog Polytechic Uiversity Hug Hom, Kowloo, Hog Kog Phoe: +852-2766-7410
More informationChapter 2: Numerical Methods
Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,
More information) +m 0 c2 β K Ψ k (4)
i ħ Ψ t = c ħ i α ( The Nature of the Dirac Equatio by evi Gibso May 18, 211 Itroductio The Dirac Equatio 1 is a staple of relativistic quatum theory ad is widely applied to objects such as electros ad
More informationFinite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw
Fiite Differece Derivatios for Spreadsheet Modelig Joh C. Walto Modified: November 15, 2007 jcw Figure 1. Suset with 11 swas o Little Platte Lake, Michiga. Page 1 Modificatio Date: November 15, 2007 Review
More informationApply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.
Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α
More informationNew Version of the Rayleigh Schrödinger Perturbation Theory: Examples
New Versio of the Rayleigh Schrödiger Perturbatio Theory: Examples MILOŠ KALHOUS, 1 L. SKÁLA, 1 J. ZAMASTIL, 1 J. ČÍŽEK 2 1 Charles Uiversity, Faculty of Mathematics Physics, Ke Karlovu 3, 12116 Prague
More informationPOWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS
Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity
More informationSequences of Definite Integrals, Factorials and Double Factorials
47 6 Joural of Iteger Sequeces, Vol. 8 (5), Article 5.4.6 Sequeces of Defiite Itegrals, Factorials ad Double Factorials Thierry Daa-Picard Departmet of Applied Mathematics Jerusalem College of Techology
More informationSenvion SE Franz-Lenz-Str. 1, Osnabrück, Germany
Iteratioal Wid Egieerig Coferece IWEC 014 WAVE INDUCED FATIGUE LOADS ON MONOPILES - NEW APPROACHES FOR LUMPING OF SCATTER TABLES AND SITE SPECIFIC INTERPOLATION OF FATIGUE LOADS M. SEIDEL Sevio SE Fraz-Lez-Str.
More informationDISTRIBUTION LAW Okunev I.V.
1 DISTRIBUTION LAW Okuev I.V. Distributio law belogs to a umber of the most complicated theoretical laws of mathematics. But it is also a very importat practical law. Nothig ca help uderstad complicated
More informationChapter 7: Numerical Series
Chapter 7: Numerical Series Chapter 7 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals
More informationA statistical method to determine sample size to estimate characteristic value of soil parameters
A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig
More informationKolmogorov-Smirnov type Tests for Local Gaussianity in High-Frequency Data
Proceedigs 59th ISI World Statistics Cogress, 5-30 August 013, Hog Kog (Sessio STS046) p.09 Kolmogorov-Smirov type Tests for Local Gaussiaity i High-Frequecy Data George Tauche, Duke Uiversity Viktor Todorov,
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationTR Use of the Implicit HHT-I3 and the Explicit ADAMS Methods with the Absolute Nodal Coordinate Formulation
R-007-05 Use of the Implicit HH-I3 ad the Explicit ADAMS Methods with the Absolute Nodal Coordiate Formulatio Bassam Hussei, Da Negrut, Ahmed A. Shabaa August 007 Abstract his ivestigatio is cocered with
More informationSECTION 2 Electrostatics
SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial
More informationUniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations
Global Joural of Sciece Frotier Research Mathematics ad Decisio Scieces Volume 3 Issue Versio 0 Year 03 Type : Double Blid Peer Reviewed Iteratioal Research Joural Publisher: Global Jourals Ic (USA Olie
More informationAnalysis and Numerical Approximation of a Singular Boundary Value Problem for the One-dimensional p-laplacian
Proceedigs of the 4th WSEAS Iteratioal Coferece o Fluid Mechaics Gold Coast Queeslad Australia Jauary 7-9 007 4 Aalysis ad Numerical Approximatio of a Sigular Boudary Value Problem for the Oe-dimesioal
More informationGoodness-Of-Fit For The Generalized Exponential Distribution. Abstract
Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationPRELIM PROBLEM SOLUTIONS
PRELIM PROBLEM SOLUTIONS THE GRAD STUDENTS + KEN Cotets. Complex Aalysis Practice Problems 2. 2. Real Aalysis Practice Problems 2. 4 3. Algebra Practice Problems 2. 8. Complex Aalysis Practice Problems
More informationStudy on Coal Consumption Curve Fitting of the Thermal Power Based on Genetic Algorithm
Joural of ad Eergy Egieerig, 05, 3, 43-437 Published Olie April 05 i SciRes. http://www.scirp.org/joural/jpee http://dx.doi.org/0.436/jpee.05.34058 Study o Coal Cosumptio Curve Fittig of the Thermal Based
More informationFinally, we show how to determine the moments of an impulse response based on the example of the dispersion model.
5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).
More informationNumerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials
IOSR Joural of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 11, Issue 6 Ver. IV (Nov. - Dec. 15), PP 1-11 www.iosrjourals.org Numerical Solutios of Secod Order Boudary Value Problems
More informationThe Numerical Solution of Singular Fredholm Integral Equations of the Second Kind
WDS' Proceedigs of Cotributed Papers, Part I, 57 64, 2. ISBN 978-8-7378-39-2 MATFYZPRESS The Numerical Solutio of Sigular Fredholm Itegral Equatios of the Secod Kid J. Rak Charles Uiversity, Faculty of
More informationSequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018
CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical
More informationOn an Application of Bayesian Estimation
O a Applicatio of ayesia Estimatio KIYOHARU TANAKA School of Sciece ad Egieerig, Kiki Uiversity, Kowakae, Higashi-Osaka, JAPAN Email: ktaaka@ifokidaiacjp EVGENIY GRECHNIKOV Departmet of Mathematics, auma
More informationTaylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH
Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )
More information