Journal of Mechanical Science and Technology 26 (10) (2012) 3231~ DOI /s
|
|
- Penelope Allison
- 5 years ago
- Views:
Transcription
1 Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~34 DOI /s A DXDR large deflectio aalysis of uiformly loaded square, circular ad elliptical orotropic plates usig o-uiform rectagular fiite-differeces M. Kadkhodaya 1,*, A. Erfai Moghadam 1, G.J. Turvey ad J. Alamatia 3 1 Departmet of Mechaical Egieerig, The Uiversity of Ferdowsi, , Mashhad, Ira Egieerig Departmet, Lacaster Uiversity, Bailrigg, Lacaster LA1 4YR, U.K. 3 Departmet of Civil Egieerig, Azad Uiversity, Mashhad, Ira (Mauscript Received October 17, 011; Revised Jauary 7, 01; Accepted April 30, 01) Abstract A fiite-differece aalysis of e large deflectio respose of uiformly loaded square, circular ad elliptical clamped ad simplysupported orotropic plates is preseted. Several types of o-uiform (graded) mesh are ivestigated ad a mesh suited to e curved boudary of e orotropic circular ad elliptical plate is idetified. The DXDR meod a variat of e DR (dyamic relaatio) meod is used to solve e fiite-differece forms of e goverig orotropic plate equatios. The DXDR meod ad irregular rectiliear mesh are combied alog wi e Cartesia coordiates to treat all types of boudaries ad to aalyze e large deformatio of o-isotropic circular/elliptical plates. The results obtaied from plate aalyses demostrate e potetial of e o-uiform meshes employed ad it is show at ey are i good agreemet wi oer results for square, circular ad elliptical isotropic ad orotropic clamped ad simply-supported plates i bo fied ad movable cases subjected to trasverse pressure loadig. Keywords: Orotropic circular ad elliptical plates; No-uiform rectagular mesh; Large deflectio; Dyamic relaatio Itroductio The icreasig use of composite materials i i-walled structures, which may be regarded as assemblies of plate ad shell elemets, provides bo e catalyst ad justificatio for e cotiual developmet ad refiemet of umerical procedures for e solutio of orotropic ad oer plate fleure problems. The solutios of particular iterest here are for e large deflectio respose of clamped ad simply-supported circular ad elliptical plates wi rectagular Cartesia orotropy. Such plates commoly arise i practice as ispectio covers ad oer types of ed closures. Bo e small ad large deflectio resposes of orotropic plates have bee ivestigated by a umber of researchers durig e course of e past half-cetury. For eample, i e 1970 s Chia [1, ] studied e oliear fleural behavior of rectagular orotropic plates usig a approimate aalytical solutio. Dalaei ad Kerr [3] ivestigated e small deflectio respose of clamped rectagular orotropic plates subjected to uiform trasverse pressure usig e eteded Katorovich meod. I later studies, Mbakogu ad Pavlovica [4] ad Tabalov et al. [5] also ivestigated e small deflectio * Correspodig auor. Tel.: , Fa.: address: kadkhoda@um.ac.ir Recommeded by Associate Editor Heug Soo Kim KSME & Spriger 01 respose of clamped rectagular orotropic plates. Mistou et al. [6] ivestigated e behavior of clamped lamiated rectagular orotropic plates subjected to uiform trasverse pressure bo eperimetally ad umerically usig a Ritz approimate aalysis. Bhaskar ad Kaushik [7] used a double Fourier series approach to aalyze usymmetric cross-ply lamiated clamped rectagular plates. Recetly, Salehi ad Sobhai [8] used Midli plate eory i cojuctio wi e DR meod to ivestigate bo e small ad large deflectio respose symmetrically lamiated fiber-reiforced sector plates. I e majority of e aforemetioed papers e plate was eier of rectagular or circular geometry. Therefore, for umerical mesh-based solutios, a uiform square or rectagular mesh would be adequate for rectagular plates wi rectagular Cartesia orotropy. Similarly, for circular plates wi polar orotropy, a orogoal polar mesh would suffice but would ot readily accommodate rectagular Cartesia orotropy. Here, however, e focus of iterest is o e use of a sigle o-uiform rectagular mesh at may be used to aalyze square, circular ad elliptical plates wi rectagular Cartesia orotropy. It is e mai objective of e curret study to show if e DXDR meod ad irregular rectiliear mesh are combied alog wi e Cartesia coordiates oe ca treat all types of
2 33 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~34 Fig. 1. Square plate wi its pricipal aes of orotropy parallel to e ad y aes: pla view; cross-sectio. boudaries which is also etedable to e large deformatio of o-isotropic plates. The mai applicatio of e curret meod ca be e aalysis of orotropic materials which are made iitially i rectagular dies ad e fibers are iserted i rectiliear dictios ad e are cut ito circular ad elliptical forms. The large deflectio respose of clamped ad simplysupported square, circular ad elliptical orotropic plates subjected to uiform trasverse pressure is ivestigated usig several types of o-uiform rectagular fiite-differece mesh. Cosequetly, e same goverig equatios may be used for square, circular ad elliptical plates. Moreover, e mesh could be used to aalyze orotropic plates wi more geeral boudary shapes. A variat of e DR meod, amely e DXDR meod [9] is employed to solve e fiite-differece forms of e large deflectio orotropic plate equatios. The results are preseted for bo deflectios ad stress couples i order to demostrate e accuracy of e preset o-uiform mesh large deflectio solutios by compariso wi couterpart solutios obtaied wi uiform meshes ad by oer approimate aalysis techiques.. Mesh geeratio Cosider a square orotropic plate wi sides of leg a ad ickess h, as show i Fig. 1. For coveiece, e origi of e co-ordiate system, y, z is located at e ceter of e mid-plae of e plate. I order to carry out a fiitedifferece aalysis of a square or rectagular plate it is geerally sufficiet to use a uiform square or rectagular mesh etedig over e plate domai. However, for plates wi curved boudaries it is ot usually possible to use a uiform square or regular rectagular mesh because e odes of e mesh do ot coicide wi curved boudary of e plate. Differet types of o-uiform rectagular meshes may be geerated i eier or bo of e co-ordiate directios ad y. The mesh may be uiform i oe directio ad ouiform i e oer directio. I is case e mesh is deoted as a Type 1 mesh. Fig. shows a eample of a Type 1 mesh, i which e mesh is uiform i e y directio ad o-uiform i e directio. Here, e mesh size i e directio reduces as e odal distace from e ceter of e plate icreases (it is assumed at symmetry is eploited ad oly e positive quadrat of e plate is modeled). Clearly, e mesh may be arraged to be o-uiform wi respect to bo co-ordiate directios. Oe eample of ese o-uiform meshes, desigated Type, is illustrated i Fig.. Such o-uiform bi-directioal meshes may be geerated i a variety of ways. Here, simple power laws are employed. The co-ordiates of e odes of e mesh i e ad y directios are defied as follows: i = a k Fig.. Eamples of o-uiform meshes: uiform i e y directio y = 1 ad o-uiform i e directio 1 [Type 1]; o- = y = 0.5 [Type ]. uiform i bo ad y directios ( ) (1a)
3 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~ j y= a k y. (1b) I Eq. (1) a is e half side leg of e square plate or e radius of e circular plate, i ad j are e ode umbers i e ad y directios, respectively, where 0 i, j k, k+ 1 is e total umber of odes i e ad y directios, respectively, ad ad y are bo positive costats. If oe of e costats or y is equal to oe ad e oer is ot equal to oe, e a Type 1 mesh is defied, i.e. = 1 ad y 1 or 1 ad y = 1. O e oer had, e mesh will be o-uiform i bo ad y directios whe = 1 ad y 1. Moreover, whe = y= 1, a uiform mesh wi respect to bo e ad y aes is defied. For plates wi a curved boudary, ad y caot be selected idepedetly. I oer words, ese quatities must be determied so at e mesh odes coicide wi curved boudary of e plate. For eample, e boudary odes of a circular plate must satisfy followig equatio: + y = a () where a is e radius of plate. Substitutig Eqs. (1a) ad (1b) ito Eq. () leads to e followig result: y i j + = 1. k k For a Type 1 mesh (uiform i oe directio ad ouiform i e oer directio), e ukow power or y i e o-uiform directio may be determied by solvig Eq. (3). I oer words, (or y ) is determied so at e mesh lies itersect e circle boudary. It is clear at (or y ) will have a specific value for each mesh lie. However, a Type mesh is defied by = y= 0.5, where e mesh has e same degree of o-uiformity i bo e ad y directios. Oer values of ad y, betwee 0 ad 1 or greater a 1, are ulikely to be of iterest because ey oly icrease e rate of directioal o-uiformity of e mesh. I such cases, ad y would have differet values for each lie of e mesh for e oer values, i.e. whe y 0.5. Nevereless, for e sake of completeess, ey are ivestigated to a limited etet ad discussed later i e paper. Similar remarks apply to oer plates wi curviliear boudaries, e.g. elliptical plates. I order to ivestigate e accuracy of ese meshes, e circular ad elliptical forms of Eq. (3) are used first to geerate meshes for square ad rectagular plates, respectively. This mesh geeratio meod may be eplaied by a simple eample. Cosider a square plate of side leg a. To defie a Type 1 mesh over a square plate, a circle of radius a wi its ceter at e ceter of e square plate is defied, as show i Fig.. The, a Type 1 mesh is geerated i oe directio ( or y ) over e (3) circle by, for eample, settig = 1 ad determiig y from Eq. (3). Here, e total umber of mesh odes i e ad y directios is assumed to be 13 (k = 1). The mesh lies over circle are e eteded to e sides of e square to produce a Type 1 mesh over e square plate. Oer mesh types (for eample, Type mesh which is show i Fig. ), may be geerated i a similar maer. For rectagular plates, a iscribed ellipse should be used. To demostrate e mesh distributio i e curret study, e used mesh for e circular ad elliptical plates are show i Figs. 3 ad 3, respectively. 3. Goverig orotropic plate equatios The Cartesia forms of e equilibrium equatios for a orotropic plate may be writte i e followig form [10]: N N y + = 0 y N y Ny + = 0 y M M y M y w N y y Fig. 3. The geeral form of mesh distributio of: circular plate; elliptical plate. w w Ny + N y + q= 0. y y (4) (5) (6)
4 334 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~34 It should be appreciated at Eqs. (4) ad (5) are e iplae equilibrium equatios ad Eq. (6) is e out-of-plae equilibrium equatio. I ese equatios, N, N y ad N y are e stress resultats, M, M y ad M y are e stress couples, w is e deflectio ad q is e trasverse pressure. The stress resultats ad stress couples may be epressed as: (, y, y ) h( σ, σ y, τ y ) h N N N dz = (7a) ( y y ) h ( σ σ y τ y ) h M, M, M z,, dz. = (7b) I Eq. (7) σ, σ y ad τ y are e direct ad shear stress compoets. The stresses uder e itegrals o e right-had sides of Eq. (7) are give i Ref. [11] as: σ Q11 Q 1 σ y = Q1 Q τ 0 0 y 0 ε 0 ε y. Q66 γ y I Eq. (8), ε, ε y ad γ y are e direct ad shear strais ad e terms Qij( i, j = 1,, 6) are e reduced stiffess which may be epressed as: Q11= E11 ( 1 υ1υ1) υ υ 1 E 1 E Q 11 1= = Q= Q ( 1 υ1υ1) ( 1 υ1υ1) E ( 1 υ1υ1) (8) (9a) (9b) (9c) = G (9d) I Eq. (9), E11 ad E are e elastic moduli wi respect to e major ad mior aes of orotropy, respectively, ad G 1 is e shear modulus wi respect to e same aes. For orotropic plates, which are cosidered here, e pricipal aes of orotropy coicide wi e ad y aes. The terms υ 1 ad υ 1 deote e major ad mior Poisso s ratios. The strais of poits at a distace z above e mid-plae of e plate are give as: 0 0 k, k y ad 0 k y are e curva- of e plate mid-plae ad tures ad twist of e plate mid-plae. The relatioships betwee e mid-plae strais ad displacemets are give as: 0 u 1 w ε = + 0 v 1 w ε = + y y 0 u v w w γ y = + +. y y (11a) (11b) (11c) I Eq. (11) u ad v are e i-plae displacemets i e ad y directios, respectively. Moreover, it should also be appreciated at geometric oliearity, due to stretchig of e plate mid-plae, is reflected by e presece of e quadratic terms i Eq. (11). Similar relatioships eist betwee e mid-plae curvatures ad deflectio of e plate as follows: 0 w k = (1a) 0 w k y y = (1b) 0 w k y. y = (1c) The boudary coditios i e plate aalyses for e square plate are: Fied clamped edge Alog =± a : w u= 0, v= 0, w= 0, = 0. Alog y =± a : w u= 0, v= 0, w= 0, = 0. y Movable clamped edge Alog =± a : (13a) (13b) 0 0 ε ε k 0 0 εy = ε y + z k y. γ 0 0 y γ y k y (10) w v= 0, w= 0, N = 0, = 0. Alog y =± a : (13c) I Eq. (10) 0 0 y ε, ε ad 0 γ y are e direct ad shear strai w u= 0, w= 0, N y = 0, = 0. y (13d)
5 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~ (c) Fied simply-supported edge Alog =± a : u= 0, v= 0, w= 0, M = 0. (13e) ( σ ) = 0 & ε = 0 membrae u u u = 0 cosθ + siθ = 0. y t (13k) Alog y =± a : u= 0, v= 0, w= 0, M y = 0. (13f) (d) Movable simply-supported edge Alog =± a : v= 0, w= 0, M = 0, N = 0. (13g) Alog y =± a : u= 0, w= 0, M = 0, N = 0. (13h) y y Because bo e plate geometry ad e loadig are symmetric about e ad y aes oly oe quarter of e plate has to be aalyzed. It is clear at w, N, N y, M ad M y are symmetric about e ad y aes, M y is atisymmetric about e ad y aes, u is symmetric about e ais ad atisymmetric about e y ais ad υ is atisymmetric about e ais ad symmetric about e y ais. Hece, e foregoig symmetry/atisymmetry coditios have to be eforced alog e ad y aes. I e case of e circular ad elliptical plates, e same symmetry/atisymmetry coditios have to be eforced alog e ad y aes. However, alog e circumferetial boudary e coditios are somewhat more complicated. For eample i e fied clamped boudary coditio, at all odal poits aroud e circumferece of e plate u, v ad w are set to zero. I additio, e slope ormal to e circumferece, i.e. i e directio (see Fig. ), has also to be set to zero. This boudary coditio is defied as: w w w y w = + = cosθ+ w si θ. y y (13i) I Eq. (13i) e agle θ is defied as show i Fig.. It is easy to show at e coditio for simply-supported cases (fied ad movable) would be: w w M = 0 = 0 cos θ + w w siθ cosθ + si θ = 0. y y (13j) Moreover, e coditio for movable cases (simplysupported ad clamped) is: 4. DXDR solutio of e goverig large deflectio orotropic plate equatios The DXDR meod is a modified versio of e origial DR meod. I is modified meod e startig vector ad dampig factor are selected more deliberately compared to e origial DR [1]. It is a iterative time-steppig solutio techique for modelig e respose of damped dyamic systems. However, it has maily bee used to obtai e static respose by artificially dampig out e oscillatios of equivalet quasidyamic systems. The startig poit of ay DXDR aalysis to determie e static structural respose is e formulatio of e quasi-dyamic system of goverig equatios. Thus, e secod-order quasi-dyamic plate equatios describig e large deflectio respose of e square, circular ad elliptical orotropic plates, i.e. Eqs. (4)-(6), may be epressed as follows for e time icremet [1]: m D&& c D& f p (14) ii i + ii i + i = i i i where D & ad 3 1 vectors of velocity ad acceleratio, fi ad p i are vectors of e iteral forces (correspodig pricipally to e LHS s of Eqs. (4)-(6)) ad e applied pressure q at e i ode ad m ii ad c ii are ( 3 3) diagoal mass ad dampig matrices for e ree displacemet compoets u, v ad w at ode i durig e iteratio, respectively. Now, because oly e static large deflectio respose of e orotropic plate is of iterest, Eq. (14) may be trasformed ito time-steppig iitial value format for e plate velocities [1]. The equatio e becomes: D && are ( ) ( τ ςi ) τ D ( + τ ςi ) ( + τ ςi ) r D i i = i +. mii & & (15) I Eq. (15), τ is e fictitious time step, ς are e odal dampig factors ad ri = ( pi fi ) is e force vector at ode i. 1 The odal displacemets D i + are determied from e velocities calculated i Eq. (15) by meas of e followig simple itegratio rule: Di Di τ Di. = + & (16) Furermore, it is coveiet to relate e dampig ad mass matrices as follows [1]: i
6 336 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~34 ii i ii c = ς m. (17) I e DXDR meod, e mass matri, odal dampig factors ad e time icremet should be defied i such a way as to guaratee e stability ad covergece of e iterative procedure. The most commo meod of achievig is objective is to determie m ii by meas of Gerschgöri eorem. Accordig to is eorem, e followig iequality must be satisfied i order to guaratee e stability of e iteratios [13]: Table 1. Elastic costats for e ree types of orotropic materials used i e large deflectio plate aalyses. Elastic costats Glass-epoy Graphite-epoy Boro-epoy E11/E G1/E υ ii ( τ ) k S, m ij T. 4 (18) j= 1 I Eq. (18) k is e umber of degrees of freedom of e structure ad S ij, T is e ij elemet of e taget stiffess matri durig e iteratio ad is give by: f S i ij, T. D j = (19) Alteratively, by applyig Rayleigh's Priciple [14] at each ode, e istataeous critical dampig factor for ode i durig e iteratio [1] may be epressed as: ς 1 T ( Di ) fi Di i =. T ( Di ) mii D i (0) Oer approaches have also bee suggested [13]. The fictitious time icremet of e iteratio is usually assumed to be costat ad equal to uity. However, some researchers have proposed at is parameter should also be based o e Rayleigh quotiet [15]. Techiques have also bee proposed to determie e iitial displacemet vector i order to reduce e computatioal time [1, 16, 17]. 5. Large deflectio results for square, circular ad elliptical plates usig o-uiform meshes Differet results are preseted i is sectio. The total umber of odes i e ad y directios, is equal to 13 (k = 1). The first set of results has bee computed for a square clamped isotropic plate subjected to uiform trasverse pressure loadig. The results are iteded to demostrate at coverged ad accurate deflectios ca be obtaied wi ouiform fiite-differece meshes. I Fig. 4 e plate ceter deflectios wi Type 1 ad Type o-uiform meshes are compared wi deflectios obtaied from Chia s [18] approimate aalytical solutio. It is clear at e ceter deflectios obtaied for bo meshes are i reasoably good agreemet wi e deflectios predicted by Chia s approimate aalytical solutio. Fig. 4. Compariso of load-ceter deflectio resposes of clamped square isotropic plates. A secod set of compariso aalyses has bee udertake for clamped square orotropic plates subjected to uiform trasverse pressure i order to demostrate at e DXDR aalysis produces coverged ad accurate results for orotropic materials. For ese comparisos it has bee ecessary to defie e elastic moduli of e orotropic plate materials to be used i e aalyses. Three materials were selected wi properties represetative of uidirectioal glass, carbo ad boro fiber reiforced polymers. The elastic costats, o-dimesioalized wi respect to e mior elastic modulus, are preseted i Table 1. The pressure-ceter deflectio resposes of glass-epoy, carbo-epoy ad boro-epoy plates are preseted i Figs. 5, 5 ad 5(c), respectively. Agai, Chia s [18] approimate aalytical solutio is used as e bechmark solutio for e compariso. It is evidet at e ceter deflectios predicted wi e Type 1 mesh (o-uiform i e directio) are i closer agreemet wi Chia s [18] values a e ceter deflectios predicted wi e Type mesh (ouiform i bo directios). It is presumed at e additioal o-uiformity of e mesh is resposible for e uderestimatio of e ceter deflectio respose by e Type mesh. It has to be poited out here at i orotropic plates e depedecy of results to e type of mesh is more severe a at for e isotropic plates. However, usig e Type 1 mesh is always more reliable i all cases. I order to eamie furer e differeces betwee e results obtaied wi bo uiform ad o-uiform meshes, e stress couple distributio alog e ais is computed for a square glass-epoy clamped plate whe e major orotropic ais is ormal to e ais. It is clear at ere is very good agreemet betwee e stress couples for bo uiform ad Type 1 o-uiform meshes. The Type o-uiform mesh appears to predict sigificatly smaller (i e egative sese) values for
7 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~ Fig. 6. Compariso of stress couple profiles alog e ais of a square glass-epoy plate subjected to uiform trasverse pressure (computed usig uiform ad Type 1 ad meshes) pricipal ais of orotropy parallel to e y ais. (c) Fig. 5. Compariso of load-ceter deflectio resposes of square orotropic clamped plates subjected to uiform trasverse pressure (computed usig Type 1 ad meshes): glass-epoy plate; carbo-epoy plate; (c) boro-epoy plate. Fig. 7. Deflectio profiles alog e ad y aes of square orotropic clamped plates subjected to uiform trasverse pressure computed usig uiform ad Type 1 meshes: major ais of orotropy i e directio; major ais of orotropy i e y directio. e stress couple whe / a 0.7, as show i Fig. 6. It is perhaps wor emphasizig at i e Type 1 mesh aalysis, show i Fig. 6, e major orotropic ais was parallel to y ais (e uiform mesh directio). This choice of directio for e major orotropic ais may be justified by comparig e results obtaied for e deflectio profiles alog e ad y aes. The clamped square orotropic plate was re-aalyzed twice wi bo e uiform ad ouiform Type 1 meshes, i.e. oce wi e major ais of orotropy parallel to e y ais ad oce wi major ais parallel to e ais. It is clear from e deflectio profiles alog e ad y aes, show i Fig. 7, at e results obtaied wi e major ais of orotropy parallel to e y ais are more accurate. It is well kow at for orotropic plates subjected to uiform trasverse pressure over e cetral regio of e plate e deflectio profile is flatter trasverse to e pricipal orotropic ais a it is parallel to is ais [18]. Moreover, for e Type 1 mesh, wi e pricipal ais of orotropy i e y directio, e mesh i e directio, i.e. trasverse to e major orotropic ais, is coarser earer to e ceter ad fier earer e edge of e plate. Hece, because e mesh is coarser i e flat regio of e plate, e error is smaller a if e coarse mesh is used i e directio of e major orotropic ais (e y directio). Oer o-uiform meshes may be geerated by usig differet values of e epoet y i Eq. (1b). Square glassepoy clamped plates subjected to uiform trasverse pressure
8 338 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~34 Table. The values correspodig to each mesh lie for y = 1., 1.3 ad 1.4 (k = 1). i j y = 1. y = 1.3 y = Fig. 9. A compariso betwee e covergece rate of two differet meshes: Type 1 mesh; mesh wi y = 1.4. y Fig. 8. Deflectio profiles alog e y ais of a square orotropic clamped plate subjected to uiform trasverse pressure computed usig uiform ad a rage of o-uiform meshes: pricipal ais of orotropy i e directio; pricipal ais of orotropy i e y directio. have bee aalyzed usig uiform ad o-uiform meshes defied by y = 1., 1.3 ad 1.4. Here, 1 ad it is calculated from Eq. (3) for each mesh lie (k = 1). Table shows e values for y = 1., 1.3 ad 1.4. The deflectio profiles alog e y ais obtaied wi e uiform ad o-uiform meshes are show i Fig. 8. It is clear at as e epoet icreases e o-uiformity of e mesh icreases ad e calculated deflectios are over-estimated compared to e values obtaied from e uiform mesh aalyses. For eample, whe y = 1.4 e mid-spa deflectio is about twice of at predicted wi e uiform mesh. Geerally e mesh used to aalyze e orotropic curved boudaries should have adequate accuracy as well as satisfactory covergece rate. As aforemetioed, e Type 1 mesh could provide sufficiet accuracy i almost all cases. Furermore, is type of mesh could also give higher covergece rate. The Figs. 9 ad 9 show e variatio of residual eergy wi iteratio umber for Type 1 mesh ad for e mesh wi y = 1.4, respectively. The o-uiform mesh may ot just be used to aalyze square isotropic ad orotropic plates. It may also be used to aalyze isotropic ad orotropic circular ad elliptical plates. It is assumed at e aes of orotropy are parallel to e ad y aes. I order to demostrate e fleibility ad accuracy of e Type 1 mesh e large deflectio respose of a clamped isotropic circular plate subjected to uiform trasverse pressure is aalyzed. The ceter deflectio-pressure respose is show i Fig. 10 ad compared wi e results give i Ref. [18] ad also wi e results obtaied by usig a DXDR polar co-ordiate fiite-differece aalysis of e problem. I order to demostrate e ability of e Type 1 mesh for bo e small ad large deflectio resposes of clamped orotropic circular plates, two furer results comparisos are preseted for a glass-epoy plate. A small deflectio aalysis of e load-deflectio respose for e case of uiform tras-
9 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~ Fig. 10. Load-ceter-deflectio respose of a clamped isotropic circular plate subjected to uiform trasverse pressure, compariso of Type 1 mesh ( y = 1, 1 ) ad oer approimate resposes. Fig. 13. Large deflectio of e load-ceter deflectio respose of a movable clamped isotropic circular plate subjected to uiform trasverse pressure, compariso of Type 1 mesh ( y = 1, 1 ) ad approimate resposes. Fig. 11. Small deflectio aalysis of e load-ceter deflectio respose of a clamped orotropic circular plate subjected to uiform trasverse pressure, compariso of Type 1 mesh ( y = 1, 1 ) ad eact resposes. Fig. 1. Deflectio profiles alog e ad y aes of a clamped orotropic circular plate subjected to uiform trasverse pressure, profiles computed wi a Type 1 mesh ( y = 1, 1 ) wi e major ais of orotropy parallel to e y ais. verse pressure loadig is show i Fig. 11. As epected, e load-deflectio respose is liear ad e Type 1 mesh results are i good agreemet wi e eact results give i Ref. [19]. Large deflectio aalysis results for e Type 1 mesh are show i Fig. 1. The deflectio profiles i e ad y directios are preseted for e pricipal aes of orotropy parallel to e y directio wi e o-uiformity of e mesh i e directio, Type 1 mesh. The Type oe mesh is also capable to ivestigate e circular ad elliptical plates wi movable boudary coditio i bo fied ad simply-supported cases. Fig. 13 shows e large deflectio aalysis of a movable clamped isotropic circular plate subjected to uiform trasverse pressure. As it is see, ere is a good agreemet betwee e results obtaied usig e Type 1 mesh ( y = 1, 1 ) ad e Baerjee results [0]. Moreover, e employed mesh ca also be used to aalyze e isotropic elliptical plates. The results obtaied by e aid of curret study ad ose of Baerjee [1] are compared to each oer i Figs for movable clamped ad simplysupported edges, respectively. Figs ad 14(c)- 15(c) show e profiles alog e ad y aes ad e stress couple profiles alog e ais for e elliptical plate for e metioed boudary coditios. Furermore, e meod ca be eteded to aalyze e large orotropic circular ad elliptical plates wi clamped ad simply-supported movable edges. Figs. 16 ad 17 show e results obtaied from curret study for glass-epoy orotropic plates. Geerally, by e aid of proposed mesh ad meod eplaied above it is ot oly possible to aalyze plates wi differet boudaries but e adequate covergece ad accuracy are also achievable. For istace, e Table 3 shows e accuracy of e results which is attaiable i is meod compared to oer published data. A quick glace at Figs displays e capability of e proposed meod ad at it could properly treat e circular ad elliptical boudaries while Cartesia aes are used. The employed proposed irregular mesh may be applied to ay oer curved boudaries wi more complicated shapes which are especially useful to aalyze e orotropic plates whe e reiforced fibers are placed i rectiliear directios. I fact, it quite happes at e reiforced sheets are made iitially i large scale ad e e plates are cut out wi a arbitrary shape ad boudary form.
10 340 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~34 (c) Fig. 14. Large deflectio aalysis of a movable clamped isotropic elliptical plate subjected to uiform trasverse pressure (a/b = ): loaddeflectio curve ad compariso of Type 1 mesh ( y= 1, 1 ) ad approimate resposes; profiles alog e ad y aes; (c) stress couple profiles alog e ais. (c) Fig. 15. Large deflectio aalysis of a movable simply-supported isotropic elliptical plate subjected to uiform trasverse pressure (a/b = ): load-deflectio curve ad compariso of Type 1 mesh ( y = 1, 1) ad approimate resposes; profiles alog e ad y aes, (c) stress couple profiles alog e ais.
11 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~ Table 3. A geeral compariso betwee e curret results wi ose of oers published data. Large deflectio aalysis Curret result Referece result Average differece (%) Square orotropic clamped plates (Glass-epoy) Fig. (5) Chia [18].74 Clamped isotropic circular plate Fig. (10) Chia [18].6 Movable clamped isotropic circular plate Fig. (13) Baerjee [0].5 Movable clamped isotropic elliptical plate Fig. (14) Baerjee [1] 1.96 Movable simply-supported isotropic elliptical plate Fig. (15) Baerjee [1] 0.68 best results whe e o-uiformity is trasverse to e major ais of orotropy. Fig. 16. Large deflectio aalysis of e load-ceter deflectio respose of a movable simply-supported ad clamped glass-epoy orotropic circular plate subjected to uiform trasverse pressure usig Type 1 mesh ( y = 1, 1 ). Fig. 17. Large deflectio aalysis of e load-ceter deflectio respose of a movable simply-supported ad clamped glass-epoy orotropic elliptical plate subjected to uiform trasverse pressure usig compariso of Type 1 mesh ( y = 1, 1 ). 6. Cocludig remarks Clamped ad simply-supported square, circular ad elliptical plates wi Cartesia orotropy subjected to uiform trasverse pressure have bee aalyzed usig e DXDR techique i cojuctio wi bo uiform ad o-uiform Cartesia fiite-differece meshes. Bo fied ad movable boudary coditios have bee ivestigated. It has bee show at e use of o-uiform fiite-differece meshes permits square, circular ad elliptical plate geometries to be aalyzed usig e same Cartesia mesh. Moreover, bo e small ad large deflectio resposes of e plates may be predicted wi good accuracy provided e o-uiformity of e mesh is ot ecessive. It appears at, for orotropic plates, e Type 1 mesh (o-uiform i oe directio oly) gives Ackowledgmet The research reported i is paper was carried out durig part of e first auor s sabbatical leave spet i e Egieerig Departmet at Lacaster Uiversity. The first auor wishes to record his appreciatio to e Iraia Miistry of Higher Educatio ad Ferdowsi Uiversity for providig fiacial support durig his sabbatical leave. Auors wish to record eir appreciatio to e Egieerig Departmet for supportig is research collaboratio. Refereces [1] C. Y. Chia, Large deflectio of rectagular orotropic plates, ASCE Proceedigs, Joural of e Egieerig Mechaics Divisio, 98 (5) (197) [] C. Y. Chia ad M. K. Prabhakara, Noliear aalysis of orotropic plates, proceedigs of e istitutio of mechaical egieers, Joural of Mechaical Egieerig Sciece, 17 (3) (1975) [3] M. Dalaei ad A. D. Kerr, Aalysis of clamped rectagular orotropic plates subjected to a uiform lateral load, Iteratioal Joural of Mechaical Scieces, 37 (5) (1995) [4] F. C. Mbakogu ad M. N. Pavlovic, Bedig of clamped orotropic rectagular plates: a variatioal symbolic solutio, Computers ad Structures, 77 () (000) [5] P. Tabakov, V. Verijeko ad B. Verijeko, Refied eory for e aalysis of lamiated orotropic structures, Composite Structures, 6 (3-4) (003) [6] S. Mistou, M. Karama ad M. Sabarots, Eperimetal ad umerical simulatios of e static behavior of orotropic plates, Measuremet, 30 (3) (001) [7] K. Bhaskar ad B. Kaushik, Aalysis of clamped usymmetric cross-ply rectagular plates by superpositio of simple eact double fourier series solutios, Composite Structures, 68 (3) (005) [8] M. Salehi ad A. R. Sobhai, Elastic liear ad o-liear aalysis of fibers-reiforced symmetrically lamiated sector midli plate, Composite Structures, 65 (1) (004) [9] M. Kadkhodaya, L. C. Zhag ad R. Sowerby, Aalysis of wriklig ad bucklig of elastic plates by DXDR meod,
12 34 M. Kadkhodaya et al. / Joural of Mechaical Sciece ad Techology 6 (10) (01) 331~34 Computers ad Structures, 65 (4) (1997) [10] M. Kadkhodaya, Usig iterlacig ad irregular meshes i plate aalysis, Iraia Joural of Sciece ad Techology, 6 (B3) (00) [11] J. M. Whitey, Structural aalysis of lamiated aisotropic plates, Techomic Publishig Compay, New York (1987). [1] L. C. Zhag, M. Kadkhodaya ad Y. W. Mai, Developmet of e madr meod, Computers ad Structures, 5 (1) (1994) 1-8. [13] P. Uderwood, Dyamic Relaatio, i T. Belytshko ad T.J.R. Hughes, eds., Computatioal Meod for Trasiet Aalysis, Chapter 5, Elsevier, Amsterdam (1983) [14] W. T. Thomso, Vibratio eory ad applicatios, Pretice-Hall, 5 Editio, New York (1997). [15] S. Qiag, A adaptive dyamic relaatio meod for oliear problems, Computers ad Structures, 30 (4) (1988) [16] L. C. Zhag ad T. X. Yu, Modified adaptive dyamic relaatio meod ad its applicatio to elastic-plastic bedig ad wriklig of circular plates, Computers ad Structures, 33 () (1989) [17] S. E. Ha ad K. S. Lee, A study of e stabilizig process of ustable structures by dyamic relaatio meod, Com- puters ad Structures, 81 (17) (003) [18] C. Y. Chia, Noliear Aalysis of Plates, McGraw-Hill, New York (1980). [19] S. G. Lekhitskii, Aisotropic Plates, Gordo ad Breach Sciece, New York (1968). [0] B. Baerjee, Large deflectios of circular plates of variable ickess, Iteratioal Joural of Solids ad Structures, 19 () (1983) [1] G. C. Siharay ad B. Baerjee, A modified approach to large deflectio aalysis of i elastic plates, Iteratioal Joural of Mechaical Scieces, 8 (3) (1986) Mehra Kadkhodaya received his B.Sc. ad M.Sc. degrees from Tehra Uiversity, Ira i 1987 ad Ph.D degree from e Uiversity of Sydey, Australia i sheet metal formig area i Curretly as a full professor he is workig i Ferdowsi Uiversity of Mashhad, Ira. He has published about 36 joural papers ad more a 50 coferece papers.
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 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 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 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 informationPILOT STUDY ON THE HORIZONTAL SHEAR BEHAVIOUR OF FRP RUBBER ISOLATORS
Asia-Pacific Coferece o FRP i Structures (APFIS 2007) S.T. Smith (ed) 2007 Iteratioal Istitute for FRP i Costructio PILOT STUDY ON THE HORIZONTAL SHEAR BEHAVIOUR OF FRP RUBBER ISOLATORS T.B. Peg *, J.Z.
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 informationFAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES
LECTURE Third Editio FAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES A. J. Clark School of Egieerig Departmet of Civil ad Evirometal Egieerig Chapter 7.4 b Dr. Ibrahim A. Assakkaf SPRING 3 ENES
More informationANALYSIS OF THERMOELASTIC STRESSES IN LAYERED PLATES
6t Iteratioal DAAAM Baltic Coferece "INDUSTRIAL ENGINEERING" 4-6 April 8, Talli, Estoia ANALYSIS OF THERMOELASTIC STRESSES IN LAYERED PLATES Kõo, J. & Valgur, J. Abstract: Tis paper deals wit termoelastic
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 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 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 informationFREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES
Abdullah : FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS USING FOURIER SERIES SAWA MUBARAK ABDUAH Assistat ecturer Uiversity of Mosul Abstract Fourier series will be
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 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 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 informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
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 informationApproximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square
WSEAS TRANSACTONS o BUSNESS ad ECONOMCS S. Khotama, S. Boothiem, W. Klogdee Approimatig the rui probability of fiite-time surplus process with Adaptive Movig Total Epoetial Least Square S. KHOTAMA, S.
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 informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationU8L1: Sec Equations of Lines in R 2
MCVU Thursda Ma, Review of Equatios of a Straight Lie (-D) U8L Sec. 8.9. Equatios of Lies i R Cosider the lie passig through A (-,) with slope, as show i the diagram below. I poit slope form, the equatio
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 informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
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 informationU8L1: Sec Equations of Lines in R 2
MCVU U8L: Sec. 8.9. Equatios of Lies i R Review of Equatios of a Straight Lie (-D) Cosider the lie passig through A (-,) with slope, as show i the diagram below. I poit slope form, the equatio of the lie
More informationCALCULUS BASIC SUMMER REVIEW
CALCULUS BASIC SUMMER REVIEW NAME rise y y y Slope of a o vertical lie: m ru Poit Slope Equatio: y y m( ) The slope is m ad a poit o your lie is, ). ( y Slope-Itercept Equatio: y m b slope= m y-itercept=
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 informationEVALUATION OF GLASS FIBER/EPOXY INTERFACIAL STRENGTH BY THE CRUCIFORM SPECIMEN METHOD
EVALUATION OF GLASS FIBER/EPOX INTERFACIAL STRENGTH B THE CRUCIFORM SPECIMEN METHOD Ju KOANAGI, Hajime KATO, Akihiro KASHIMA, uichi IGARASHI, Keichi WATANABE 3, Ichiro UENO 4 ad Shiji OGIHARA 4 Istitute
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 information1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y
Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every o-zero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:
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 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 information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationHonors Calculus Homework 13 Solutions, due 12/8/5
Hoors Calculus Homework Solutios, due /8/5 Questio Let a regio R i the plae be bouded by the curves y = 5 ad = 5y y. Sketch the regio R. The two curves meet where both equatios hold at oce, so where: y
More informationMTH Assignment 1 : Real Numbers, Sequences
MTH -26 Assigmet : Real Numbers, Sequeces. Fid the supremum of the set { m m+ : N, m Z}. 2. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a
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 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 informationwavelet collocation method for solving integro-differential equation.
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
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 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 7: Polar representation of complex numbers
Lecture 7: Polar represetatio of comple umbers See FLAP Module M3.1 Sectio.7 ad M3. Sectios 1 ad. 7.1 The Argad diagram I two dimesioal Cartesia coordiates (,), we are used to plottig the fuctio ( ) with
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More informationMULTIPLE TIME SCALES SOLUTION OF AN EQUATION WITH QUADRATIC AND CUBIC NONLINEARITIES HAVING FRAC- TIONAL-ORDER DERIVATIVE
Mathematical ad Computatioal Applicatios, Vol. 6, No., pp. 3-38,. Associatio for Scietific Research MULIPLE IME SCALES SOLUION OF AN EQUAION WIH QUADRAIC AND CUBIC NONLINEARIIES HAVING FRAC- IONAL-ORDER
More informationANALYSIS OF EXPERIMENTAL ERRORS
ANALYSIS OF EXPERIMENTAL ERRORS All physical measuremets ecoutered i the verificatio of physics theories ad cocepts are subject to ucertaities that deped o the measurig istrumets used ad the coditios uder
More informationNumerical Solution of the First-Order Hyperbolic Partial Differential Equation with Point-Wise Advance
Iteratioal oural of Sciece ad Research (ISR) ISSN (Olie): 39-74 Ide Copericus Value (3): 4 Impact Factor (3): 4438 Numerical Solutio of the First-Order Hyperbolic Partial Differetial Equatio with Poit-Wise
More informationTHE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
More informationNBHM QUESTION 2007 Section 1 : Algebra Q1. Let G be a group of order n. Which of the following conditions imply that G is abelian?
NBHM QUESTION 7 NBHM QUESTION 7 NBHM QUESTION 7 Sectio : Algebra Q Let G be a group of order Which of the followig coditios imply that G is abelia? 5 36 Q Which of the followig subgroups are ecesarily
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationAN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS
http://www.paper.edu.c Iteratioal Joural of Bifurcatio ad Chaos, Vol. 1, No. 5 () 119 15 c World Scietific Publishig Compay AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC
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 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 informationμ are complex parameters. Other
A New Numerical Itegrator for the Solutio of Iitial Value Problems i Ordiary Differetial Equatios. J. Suday * ad M.R. Odekule Departmet of Mathematical Scieces, Adamawa State Uiversity, Mubi, Nigeria.
More informationTHEORETICAL RESEARCH REGARDING ANY STABILITY THEOREMS WITH APPLICATIONS. Marcel Migdalovici 1 and Daniela Baran 2
ICSV4 Cairs Australia 9- July, 007 THEORETICAL RESEARCH REGARDING ANY STABILITY THEOREMS WITH APPLICATIONS Marcel Migdalovici ad Daiela Bara Istitute of Solid Mechaics, INCAS Elie Carafoli, 5 C-ti Mille
More informationSection 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations
Differece Equatios to Differetial Equatios Sectio. Calculus: Areas Ad Tagets The study of calculus begis with questios about chage. What happes to the velocity of a swigig pedulum as its positio chages?
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 informationFinite Difference Approximation for Transport Equation with Shifts Arising in Neuronal Variability
Iteratioal Joural of Sciece ad Research (IJSR) ISSN (Olie): 39-764 Ide Copericus Value (3): 64 Impact Factor (3): 4438 Fiite Differece Approimatio for Trasport Equatio with Shifts Arisig i Neuroal Variability
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 informationResearch Article Two Expanding Integrable Models of the Geng-Cao Hierarchy
Abstract ad Applied Aalysis Volume 214, Article ID 86935, 7 pages http://d.doi.org/1.1155/214/86935 Research Article Two Epadig Itegrable Models of the Geg-Cao Hierarchy Xiurog Guo, 1 Yufeg Zhag, 2 ad
More information(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is
Calculus BC Fial Review Name: Revised 7 EXAM Date: Tuesday, May 9 Remiders:. Put ew batteries i your calculator. Make sure your calculator is i RADIAN mode.. Get a good ight s sleep. Eat breakfast. Brig:
More informationNewton Homotopy Solution for Nonlinear Equations Using Maple14. Abstract
Joural of Sciece ad Techology ISSN 9-860 Vol. No. December 0 Newto Homotopy Solutio for Noliear Equatios Usig Maple Nor Haim Abd. Rahma, Arsmah Ibrahim, Mohd Idris Jayes Faculty of Computer ad Mathematical
More informationFundamental Concepts: Surfaces and Curves
UNDAMENTAL CONCEPTS: SURACES AND CURVES CHAPTER udametal Cocepts: Surfaces ad Curves. INTRODUCTION This chapter describes two geometrical objects, vi., surfaces ad curves because the pla a ver importat
More information( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.
Joh Riley 8 Jue 03 ANSWERS TO EXERCISES IN APPENDIX B Sectio B VECTORS AND SETS Exercise B-: Covex sets (a) Let 0 x, x X, X, hece 0 x, x X ad 0 x, x X Sice X ad X are covex, x X ad x X The x X X, which
More informationExact analysis of resonance frequency and mode shapes of isotropic and laminated composite cylindrical shells; Part I: analytical studies
Joural of Mechaical Sciece ad Techology ( (0 ~ www.sprigerlik.com/cotet/- DOI 0.00/s0-0-00- Eact aalysis of resoace frequecy ad mode shapes of isotropic ad lamiated composite cylidrical shells; Part I:
More informationModified Logistic Maps for Cryptographic Application
Applied Mathematics, 25, 6, 773-782 Published Olie May 25 i SciRes. http://www.scirp.org/joural/am http://dx.doi.org/.4236/am.25.6573 Modified Logistic Maps for Cryptographic Applicatio Shahram Etemadi
More informationStability Analysis of the Euler Discretization for SIR Epidemic Model
Stability Aalysis of the Euler Discretizatio for SIR Epidemic Model Agus Suryato Departmet of Mathematics, Faculty of Scieces, Brawijaya Uiversity, Jl Vetera Malag 6545 Idoesia Abstract I this paper we
More informationFinite Difference Approximation for First- Order Hyperbolic Partial Differential Equation Arising in Neuronal Variability with Shifts
Iteratioal Joural of Scietific Egieerig ad Research (IJSER) wwwiseri ISSN (Olie): 347-3878, Impact Factor (4): 35 Fiite Differece Approimatio for First- Order Hyperbolic Partial Differetial Equatio Arisig
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 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 information577. Estimation of surface roughness using high frequency vibrations
577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia
More informationApplication of Dynamic Relaxation in Thermo-Elastic Structural Analysis of Highway Pavement Structures
9 th Iteratioal LS-DYNA Users Coferece Simulatio Techology () Applicatio of Dyamic Relaxatio i Thermo-Elastic Structural Aalysis of Highway Pavemet Structures Samir N. Shoukry, Gergis W. William, Mourad
More informationTeaching Mathematics Concepts via Computer Algebra Systems
Iteratioal Joural of Mathematics ad Statistics Ivetio (IJMSI) E-ISSN: 4767 P-ISSN: - 4759 Volume 4 Issue 7 September. 6 PP-- Teachig Mathematics Cocepts via Computer Algebra Systems Osama Ajami Rashaw,
More informationSome New Iterative Methods for Solving Nonlinear Equations
World Applied Scieces Joural 0 (6): 870-874, 01 ISSN 1818-495 IDOSI Publicatios, 01 DOI: 10.589/idosi.wasj.01.0.06.830 Some New Iterative Methods for Solvig Noliear Equatios Muhammad Aslam Noor, Khalida
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 informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationStopping oscillations of a simple harmonic oscillator using an impulse force
It. J. Adv. Appl. Math. ad Mech. 5() (207) 6 (ISSN: 2347-2529) IJAAMM Joural homepage: www.ijaamm.com Iteratioal Joural of Advaces i Applied Mathematics ad Mechaics Stoppig oscillatios of a simple harmoic
More informationCS321. Numerical Analysis and Computing
CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca
More informationA COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS
A COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS Y. S. PARK a adw.t.ang b Computatioal ad Mathematical Scieces Group Faculty of Iformatio
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationComparison of Minimum Initial Capital with Investment and Non-investment Discrete Time Surplus Processes
The 22 d Aual Meetig i Mathematics (AMM 207) Departmet of Mathematics, Faculty of Sciece Chiag Mai Uiversity, Chiag Mai, Thailad Compariso of Miimum Iitial Capital with Ivestmet ad -ivestmet Discrete Time
More informationAn Analysis of a Certain Linear First Order. Partial Differential Equation + f ( x, by Means of Topology
Iteratioal Mathematical Forum 2 2007 o. 66 3241-3267 A Aalysis of a Certai Liear First Order Partial Differetial Equatio + f ( x y) = 0 z x by Meas of Topology z y T. Oepomo Sciece Egieerig ad Mathematics
More information(a) (b) All real numbers. (c) All real numbers. (d) None. to show the. (a) 3. (b) [ 7, 1) (c) ( 7, 1) (d) At x = 7. (a) (b)
Chapter 0 Review 597. E; a ( + )( + ) + + S S + S + + + + + + S lim + l. D; a diverges by the Itegral l k Test sice d lim [(l ) ], so k l ( ) does ot coverge absolutely. But it coverges by the Alteratig
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 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 informationSubject: Differential Equations & Mathematical Modeling -III. Lesson: Power series solutions of Differential Equations. about ordinary points
Power series solutio of Differetial equatios about ordiary poits Subject: Differetial Equatios & Mathematical Modelig -III Lesso: Power series solutios of Differetial Equatios about ordiary poits Lesso
More informationCHAPTER 8 SYSTEMS OF PARTICLES
CHAPTER 8 SYSTES OF PARTICLES CHAPTER 8 COLLISIONS 45 8. CENTER OF ASS The ceter of mass of a system of particles or a rigid body is the poit at which all of the mass are cosidered to be cocetrated there
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 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 informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
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 informationSample Size Estimation in the Proportional Hazards Model for K-sample or Regression Settings Scott S. Emerson, M.D., Ph.D.
ample ie Estimatio i the Proportioal Haards Model for K-sample or Regressio ettigs cott. Emerso, M.D., Ph.D. ample ie Formula for a Normally Distributed tatistic uppose a statistic is kow to be ormally
More informationMathematics Extension 2
009 HIGHER SCHOOL CERTIFICATE EXAMINATION Mathematics Etesio Geeral Istructios Readig time 5 miutes Workig time hours Write usig black or blue pe Board-approved calculators may be used A table of stadard
More informationCommutativity in Permutation Groups
Commutativity i Permutatio Groups Richard Wito, PhD Abstract I the group Sym(S) of permutatios o a oempty set S, fixed poits ad trasiet poits are defied Prelimiary results o fixed ad trasiet poits are
More informationSolving a Nonlinear Equation Using a New Two-Step Derivative Free Iterative Methods
Applied ad Computatioal Mathematics 07; 6(6): 38-4 http://www.sciecepublishiggroup.com/j/acm doi: 0.648/j.acm.070606. ISSN: 38-5605 (Prit); ISSN: 38-563 (Olie) Solvig a Noliear Equatio Usig a New Two-Step
More informationSize, shape and temperature effect on nanomaterials
Idia Joural of Pure & Applied Physics Vol. 53, November 2015, pp. 768-775 Size, shape ad temperature effect o aomaterials G Sharma, S Bhatt, R Kumar & M Kumar* Departmet of Physics, G.B. Pat Uiversity
More informationIN many scientific and engineering applications, one often
INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS, VOL 3, NO, FEBRUARY 07 5 Secod Degree Refiemet Jacobi Iteratio Method for Solvig System of Liear Equatio Tesfaye Kebede Abstract Several
More informationChapter 1. Complex Numbers. Dr. Pulak Sahoo
Chapter 1 Complex Numbers BY Dr. Pulak Sahoo Assistat Professor Departmet of Mathematics Uiversity Of Kalyai West Begal, Idia E-mail : sahoopulak1@gmail.com 1 Module-2: Stereographic Projectio 1 Euler
More informationFailure analysis of single-bolted joint for lightweight composite laminates and metal plate
IOP Coferece Series: Materials Sciece ad Egieerig PAPER OPEN ACCESS Failure aalysis of sigle-bolted joit for lightweight composite lamiates ad metal plate To cite this article: Lijie Li et al 018 IOP Cof.
More informationA beam nite element based on layerwise trigonometric shear deformation theory
A beam ite elemet based o layerwise trigoometric shear deformatio theory R.P. Shimpi *, A.V. Aiapure Departmet of Aerospace Egieerig, Idia Istitute of Techology, Bombay, Powai, Mumbai 400 076, Idia Abstract
More informationCUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE
CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE Ryutaro SEGAWA 1, Shizuo YAMAMOTO, Akira SONE 3 Ad Arata MASUDA 4 SUMMARY Durig a strog earthquake, the respose of a structure
More information