Identification of Axial Forces on Statically Indeterminate Pin-Jointed Trusses by a Nondestructive Mechanical Test

Size: px
Start display at page:

Download "Identification of Axial Forces on Statically Indeterminate Pin-Jointed Trusses by a Nondestructive Mechanical Test"

Transcription

1 Send Orders of Reprints t reprints@benthmscience.net 50 The Open Civi Engineering Journ, 0, 7, Open Access Identifiction of Axi Forces on Stticy Indeterminte Pin-Jointed Trusses by Nondestructive Mechnic Test Emiio Turco* DADU, Università degi Studi di Sssri, 0704 Aghero, Sssri, Iti Abstrct: This work concerns the stress-stte identifiction, produced by ods or disoctions, of stticy indeterminte structures. In prticur, this pper des with the cse of pin-jointed trusses, since they re combintion of simpicity of the structur mode with very interesting technic ppictions. We present n pproch tht foows the min guideines of the fexibiity method to rete n ddition od nd the consequent nod dispcements. When sufficient number of nod dispcements re mesured, they produce system of equtions tht when soved furnishes compete reconstruction of the stress-stte to identify. In order to highight the potentiities nd the imits of the proposed pproch nd so to deinete its min chrcteristics, some simpe tests re discussed nd nyzed. Keywords: Inverse probems, stress-stte identifiction, nondestructive tests.. INTRODUCTION Knowedge of the stress-stte of mechnic prt, such s n engine, whee or ger, is significnt for the mechnic industry since the excess of stresses coud give rise to premture nd sudden brek. In order to design possibe reiforcements nd risk ssessments, n in-depth knowedge of the ctu stress-stte on structure is necessry in the fied of civi engineering. Another reserch fied where the identifiction of the ctu stress hs some reevnt ppictions is biomechnics: n enightening exmpe is reted to the determintion of rtery hypertension, s reported by []. These brief discussions contribute in highighting the importnce of these issues in the scientific community. It is known tht compete estimtion of the stress-stte is required for residu stresses, i.e. stresses corresponding to nu ods []. In the technic iterture, vrious pproches hve been reported in retion to the probem of identifying the ctu stte of stress in structure. These methods cn be subdivided into destructive nd nondestructive tests. The first css of methods incudes the tests requiring n hoe or cut on the seected prt. This cuses rextion of the stress fied nd consequent dispcements; the mesurement of such dispcements gives the dt required to reconstruct the pre-existing stress fied. In fct, it is enough to impose the opposite of the mesured dispcement fied nd derive the consequent stress fied, for exmpe the ppers of [] nd [4]. Destructive tests obviousy hve the disdvntge to prtiy or toty compromise the eement to test. On the other hnd, nondestructive tests to some extent deduce the stress fied by mesuring the effects induced by prescribed sign. Among others, the method bsed on the *Address correspondence to this uthor t the DADU, Università degi Studi di Sssri, 0704 Aghero, Sssri, Iti; Te: ; Fx: ; E-mi: emiio.turco@uniss.it X-ry diffrction is quite widespred but imited by the poor cpcity of penetrtion. Another we-estbished pproch is the utrsonic one, bsed upon the dependence between the wve veocity nd stresses. For both these methods, the book of [5] contins brief introduction nd sever references. Aso the rtice of [6] introduces the use of mgnetophotoesticity in this sme frmework. Both in the sttic nd dynmic fieds signs nd mesurements of mechnic type, forces nd dispcements, hve been considered, both in sttic nd dynmic fied. In the technic iterture [] is n ttempt to frme the probem furnishing some usefu nd firy gener guideines [7]. shows tht residu stresses re defined on the boundry when the Dirichet to Neumnn mp is known [8]. gives some usefu comments on the modeing of the residu stress suggesting the use of mode bsed on the Hrtig s w. These suggestions re so prtiy described in [9], which improves the mode redy presented in [7] showing the iniectivity of the mp. Another interesting work contined in [0] concerns the identifiction of pre-stressed cbe structures by dynmic tests. This work hs initiy portryed nondestructive stressstte identifiction strtegy ppicbe to re structur probems. The work provide detis bout simpicity nd technic ppiction of the pin-jointed trusses, which ow to expore the essenti spects of the probem without invoving hevy computtions. However, the formution deveoped for pin-jointed trusses cn be esiy extended to ny probem discretized by finite eements by using the mtrix nguge. The pper shows the significnce of the retion between n ddition sttic od nd the consequent dispcements on the exmined structure. Moving into the frmework of the force method, this retion is written using stndrd esticity hypotheses such s Hooke's w nd sm dispcements / 0 Benthm Open

2 Identifiction of Axi Forces on Stticy Indeterminte Pin-Jointed Trusses The Open Civi Engineering Journ, 0, Voume 7 5 Some brief considertion retive to the use of the Hrtig's w re so presented. Here, the choices in the design of the proposed gorithm were simpe in order to be esiy impemented on n ctu probem. This, for exmpe, ws the motivtion to use the ddition sttic ods with respect to those dynmic. After this brief introduction, Section contins brief discussion on the sttement of the underying identifiction probem nd discussion on the min toos used to buid the identifiction strtegy for stticy indeterminte pin-jointed trusses. Aso some remrks on the possibe use of Hrtig's w re briefy sketched. Section comprises of some resuts on simpe pin-jointed truss structures usefu to iustrte the min chrcteristics of the proposed pproch. Concuding remrks re eborted in Section 4.. IDENTIFICATION PROBLEM: STATEMENT AND TOOLS.. Probem Sketch It is considered n ssigned configurtion of pin-jointed truss, denoted s current, the structure is so considered in equiibrium under the extern ods. Axi forces on the truss re coected in the vector t nd cn be produced by ods nd disoctions such s therm vrition or ck of fit between the brs. The identifiction of vector t shoud be done by simpe procedure possiby dptbe to in situ tests. In cse the system is stticy determined the probem cn be soved through the soution of equiibrium equtions, obviousy. Conversey, if the truss is stticy undetermined it is necessry n d hoc identifiction procedure. Beginning from the current configurtion, n ddition od ws ppied nd the truss moves to new configurtion denoted s test configurtion. The increment of xi forces nd of nod dispcements were denoted by nd respectivey. The question is wheter it is possibe to identify the xi force vector by ssuming the dequte number of the dispcement increments of the vector s known? It hs been ssumed tht identifiction probem cn be treted s iner nd estic by considering the position of ech node in the current configurtion, nd so the ength of ech br, the mteri constnt, such s the Young's moduus, nd, finy, the cross-section re of ech br s known. There re, bsicy, two terntive wys to sove the identifiction probem s mentioned erier. The first foows the guideines outined by the stiffness method tht substntiy ssumes the nod dispcements s primry unknows. On the other hnd, the second method is the fexibiity method which ssumes the xi forces on the truss s primry unknowns. Here, the unknowns fctors re xi forces nd so the preferences wi be given to the second method. The potentiities of the fexibiity method ws expored to buid n gorithm be to sove our probem. In the next section, the identifiction probem wi be defined in more form wy by buiding too which wi generte the system of equtions nd soving those equtions the xi force vector wi be obtined... Toos Foowing the wy trced in [] nd more recenty in [], n ssemby of b brs nd j joints on which c kinemtic constrints prevent some dispcements hs been considered. Extern forces cting on uncostrined joints re coected in the -dimension vector where d is equ to or for two- nd three-dimension probem respectivey. The equiibrium mtrix inks the xi forces on the brs coected in the b-dimension vector with the od vector nd, consequenty, hs dimensions. Kinemtic comptibiity is ensured by the mtrix, see [], which inks the -dimension vector of nod dispcements to the b-dimension vector binding br eongtions () The retion between br eongtions e nd xi forces t cn be written by mens of the xi fexibiity mtrix F nd the vector of the imposed br eongtion coects br eongtions deriving for exmpe from therm vrition or ck of fit between brs. We cn write this constitutive retion s () Fexibiity mtrix is squre, with b-dimension, nd digon with the i-th term, in the cse we use Hooke's w, which is equ to nd depends on the ength, the cross-section re nd the Young's moduus of the i-th br. It hs been emphsized tht the bove retion is deduced from the St. Vennt's probem nd therefore is referred to the initi nd stress-free configurtion. Here the cse of stticy indeterminte nd kinemticcy determined pin-jointed trusses s defined in [] re considered. A structure is supposed to beong to this css if the equiibrium system of equtions hs infinite soutions for ny right hnd side for the equiibrium system (), whie the comptibiity system of equtions () hs unique soution for some prticur right hnd side, but otherwise hs no soution. This choice derives from the technic importnce of this css of structures. Systems of equtions (), () nd () cn be rrnged in vrious wys to find the unknown vector, nd of the so-ced direct probem, giving the weknown soution gorithm, such s the force method nd the dispcement method. For the stress stte identifiction probem defined bove the force method ppers to be the most promising. The rod trced by the force method works by ssuming the br xi forces s primriy unknowns. This method cn be synthesized strting from the representtion of br xi forces which stisfy, priori, the equiibrium equtions: (5) where soves the system of equtions nd the k-th coumn of, tht is, soves the homogeneous system of equtions. Using the constitutive equtions (), the br eongtion vector cn be expressed s (6) () (4)

3 5 The Open Civi Engineering Journ, 0, Voume 7 Emiio Turco nd, finy, the imposition of the kinemtic comptibiity in the form proposed in [] (7) produces the fin system of equtions (8) The sme resut cn be obtined by imposing the kinemtic comptibiity through the use of the virtu work theorem, precisey the virtu force theorem, see for exmpe []. The theorem of virtu forces cn so be used to evute nod dispcements by using the retion (9) where nd re, respectivey, the xi force vector nd the nod dispcement vector deriving from the soution of (8), whie is soution of the probem. It is enough to choose s s vector with the ony nonzero component equ to unity in correspondence with the dispcement to ccute nd obtin the formu (0) where corresponds through the equiibrium mtrix, to the prticur choice of. Writing (0) for two distinct soutions of (8) beed with nd we obtin () where to sim the nottion we hve used the positions nd. In the i-th digon term of the mtrix, see (4), there is the initi br ength tht is unknown in our probem whie the current br ength is known. For generic br nd two distinct configurtion nd, the eongtion increment is worth whie the br ength referred to stte is Combining the st two retions we obtin () () (4) nd so the corresponding digon term of the fexibiity mtrix referred to the current configurtion is (5) Finy, nod dispcement increments coected in the vector cn be written s (6) where is formed simpy by binding the corresponding to the considered nod dispcements. With these simpe toos cn be identified the xi forces in n ssigned configurtion, for exmpe in the current, which cn be reched under the ction of ods nd imposed eongtions. By imposing n ddition od on the structure produces both the increment of xi forces nd of nod dispcements. In gener, ony prt of this dispcement increments is mesurbe. We denote with the subset of mesured. Axi forces to identify nd their increments produced by the ddition od, cn be represented in such wy tht they stisfy the equiibrium equtions: (7) where the unknown quntities re coected in nd. Retion (6) cn be seen s system of equtions contining the unknowns nd whie is the known term... Remrks on the constitutive w We hve discussed bove simpe strtegy for the identifiction of the xi forces in pin-jointed trusses bsed on stndrd hypotheses of iner esticity. Here, s suggested in [8], some possibe modifictions to the constitutive w wi be considered briefy in order to modify (5) in ccordnce with the Hrtig's w insted of the cssic Hooke's w. Hrtig's w, see [4], cn simpy be written s (8) where the Young's moduus, defined s increment rtio between the stress nd the engineering strin, depends inery on through two constitutive coefficients:, id est the Young's moduus for, nd the dimension coefficient. The stress cn expicity be written s (9) from which the engineering strin cn esiy be deduced (0) Considering the bove two different sttes, denoted with nd, for generic br of the truss, simpe gebric mniputions permit to write the eongtion increment s () nd by using (0) () where the xi force nd its increment on the considered br re used. This retion inks the eongtion increments nd xi force increments for the considered br nd consequenty furnish the definition of the fexibiity coefficient for the br. The fexibiity coefficient for generic br cn so be derived by inerizing Hrtig's w s reported in [8]. Strting from the ineriztion () the combintion of () nd (0) gives (4) tht is the ink between nd which gin defines gin the fexibiity coefficient of the considered br. Hrtig's w or its inerized version produces ony different evution of the fexibiity coefficients nd its use in stress-stte identifiction gorithm shoud be evuted by compring costs nd benefits.. SOME RESULTS Here some resuts referred to simpe pin-jointed trusses re being discussed. The im is to expore the min fetures of the proposed identifiction strtegy nd to render the princip drwbcks evidenty. In the foowing, trusses formed by brs with the sme Young's moduus nd cross-section re hve been considered. Where necessry for numeric ccutions,

4 Identifiction of Axi Forces on Stticy Indeterminte Pin-Jointed Trusses The Open Civi Engineering Journ, 0, Voume 7 5 GP nd cm corresponding to n iron L- profie 40 x 5 (UNI 578-7) hve been ssumed... First exmpe: Two-Br Truss The current configurtion of the two-br pin-jointed truss sketched in Fig. () hs been considered s first exmpe. Fig. (). Two-brs truss probem: geometry of the current configurtion. If the ntur or stress-free ength of ech br is, on ech br there is n xi force. In order to identify the xi forces, strting from the current configurtion, horizont force on the centr node hs been ppied. The symmetry of the probem produces horizont dispcement which chrcterizes the test configurtion. Assuming the equiibrium of the xi forces in the current configurtion, both nd its increment, cn be represented s (5) since becuse there re no extern ods in the current configurtion. Vectors nd hve ony one component nd represent, for exmpe, the xi force nd its increment on the nd br. It is simpe to prove tht the sef-stress mtrix is equ to whie the vector is worth (6) (7) nd, finy, the mtrix ssumes the foowing form (8) Consequenty, the retion between the horizont dispcement increment of the centr node nd the stressstte, current nd its increments, is (9) nd cn be used to evute the unknown if it is observed tht for the symmetry nd tht the increment dispcement is, in our identifiction probem, mesured quntity nd therefore shoud be treted s known. It resuts (0) which, using the first of (5) to evute, produces the sme vue of the reference soution. First of the symmetry condition ws used in this cse s there ws ony one mesurbe dispcement. Secondy, n error free soution ws obtined since the position of the centr node is the sme both in the current nd ntur configurtions. In more gener exmpes beow, we hve the possibiity to evute the identifiction error when current nd ntur configurtions re different... Second Exmpe: Three-Br Truss The stress-free configurtion of the three-br truss is iustrted in Fig. (). If n imposed eongtion is ppied on the nd br, for exmpe deriving from n increment of temperture, this produces vertic dispcement of the free node nd xi forces on the brs. In order to simpify the ccutions, we choose for the imposed eongtion the vue () tht produces vertic dispcement of the free-node just equ to nd the xi force vector () The identifiction of this xi force vector foowing the proposed procedure strts from the ccution of the vector, of the sef-stress mtrix nd of the mtrix Fig. (b). It is convenient for these ccutions to consider the stticy determined schemes reported in Fig. (). More precisey, derives from the soution of the equiibrium equtions reted to the scheme reported in Fig. () nd resuts whie () derives from the soution of the equiibrium equ- e f o () (b),,,, Fig. (). Three-br truss probem: Figur geometry e. of Three-br the stress-free trussconfigurtion probem: geometry () nd current of the stress-free configurtion configu- with test od (b). f v

5 54 The Open Civi Engineering Journ, 0, Voume 7 Emiio Turco tions reted to the scheme reported in Fig. (b). (4) nd, finy, the mtrix cn be written by referring to the schemes reported in Fig. (c) nd (d) nd produces the foowing resuts (5) The dispcement increment vector retive to the ony free node cn be be ccuted s (6) being the term retive to -th br of the digon mtrix equ to (7) The use of xi force representtion (7) permits to ccute the expicit expressions of the horizont nd vertic dispcement increment of the free node. The resut is then (8) (9) nd using (7) the xi force vector to identify. It foows tht (40) where it ws noted tht for this probem it is necessry to use horizont force s test od in order to void nu vue of horizont dispcement. For the purpose of this work, reference soution for dispcement increments nd cn be ccuted by studying the probem reported in Fig. (). Using the stiffness method, the equiibrium imposed on the current configurtion, see gin Fig. (), immeditey produces the vues (4) Assuming tht both nd re mesured nd consequenty known, we cn evute the unknowns nd tht tke on the vues where the esticity coefficients of the i-th br re referred to the stress-free stte tht in this ccution is considered s known. Now, the identifiction error is defined s id ref (4) ref f o f v () (b) (c) (d) Fig. (). Three-br truss probem: uxiiry schemes for the identifiction probem. Figur e. Three-br truss probem: uxiiry schemes for the identifiction probem. Error Vrying Tbe. Three-br Test: Identifiction Rtio

6 Identifiction of Axi Forces on Stticy Indeterminte Pin-Jointed Trusses The Open Civi Engineering Journ, 0, Voume 7 55 where s is the ssumed vue of (), simpe ccutions show tht this error is ony dependent on the rtio. Tbe shows the identifiction error s the rtio increses in the cse of m, N nd N. It ws noted tht, in the rnge considered for the rtio, the identifiction error depends prcticy on iner wy from in such wy tht mkes the proposed identifiction strtegy ccurte... Third Exmpe: Six-Br Truss In the third test we consider the stress-free stte of the six-br pin-joined truss ws considered s reported in Fig. (4). An equiibrted sef-stress fied cn be chieved by incresing temperture on the -rd br. This produces dispcement of the 4-th node, depicted in Fig. (4b), denoted by the prmeter. It is very simpe to ccute the dispcement prmeter by using the fexibiity or the stiffness method. It ssumes the vue (4) where is the thermic dittion coefficient chrcteristic of the considered mteri. The corresponding vector of xi forces which needs to be identified is equ to (44) Our identifiction process uses the forces nd on the 4-th node s mechnic sign, see Fig. (4c), nd s mesured dispcement increments tht do not ssume the zero vue (48) Mtrix cn be ccuted by mking use of the uxiiry schemes reported in Fig. (6). It tkes the vue (49) In order to ccute (6) it is necessry to dd the expression of the digon eement of the fexibiity mtrix corresponding to the i-th br (50) Assuming the known the dispcements coected in the soution of (6) cn be chieved in the foowing two steps (5) This is necessry since the probem defined by (6) is noniner in the unknowns coected in nd so convenient for the prticur structure of mtrix. For the first system of (5), simpe ccution provides (5) (45) where nd re the horizont nd vertic dispcement of the k-th node, respectivey. As in the previous exmpes, both the vector of xi forces nd its increment cn be represented s (46) whie there re no conditions on nd The prticur structure of permits the foowing rerrngement In this cse, nd coect the xi forces nd their increments produced by the ddition forces, nd, on rd nd 4th br. Vector nd mtrix cn esiy be ccuted by referring to the schemes reported in Fig. (5) nd Fig. (5,b-c) respectivey. Using the prmeters nd, they ssume the vues which, when soved, gives the unknowns (5) (54) (47) where the foowing positions re used

7 56 The Open Civi Engineering Journ, 0, Voume 7 Emiio Turco (55) As in the exmpes iustrted bove, in order to compre the obtined resuts it is necessry to estimte the nod dispcements coected in. This cn be done by evuting the dispcements of the scheme reported in Fig. (4c) cused by the forces nd. By referring to the current configurtion nd in the frmework of the fexibiity method, xi force increments on the rd nd 4th br, see (8), cn be computed from the foowing system tht when it is soved it produces (56) (57) from which, using (5), cn be ccuted. Successivey, the requested dispcements cn be ccuted by mens of virtu work (9). Resuts deriving from the proposed identifiction pproch depends gin ony on the rtio. In the rnge it vries in prcticy iner fshion from. It is so interesting to study how the proposed identifiction pproch is ffected by the errors tht certiny exist in the considered identifiction probem. For exmpe, nod dispcements re obtined from mesurements nd thus re, certiny, ffected by errors. In order to simpy evute the effects of mesurement errors, we indicte with nd the re nd the mesured dispcements respectivey. In this cse the dt error prmeter cn impicity be defined s (58) Tbe reports for the identifiction error incresing dt error prmeter. In ddition to dt reported in the beginning of this section we hve ssumed m, N nd N. These resuts show, in n unmbiguous wy, the necessity of speci techniques to fiter the errors tht re wys present in the dt. A fitering strtegy such s tht proposed by Tikhonov [5] coud be usefu to sove ctu probem. 4. CONCLUDING REMARKS This pper hs sketched n identifiction procedure be to reconstruct ctu xi forces on pin-jointed truss. The proposed procedure is bsed on the resuts gthered through simpe mechnic test tht uses prescribed od nd the mesured dispcement induced. In the frmework of the fexibiity method, the ink between the xi forces to be identified nd the structur response to the test od is described in the context of stndrd iner esticity. In Section simpe numeric tests showed the min chrcteristics of the proposed procedure nd so highighted its performnces. When mesurement errors re tken into ccount, quick view of the identifiction of xi forces suggests tht n d hoc procedure is strongy necessry to fiter such errors s it is usu in the identifiction procedures. Tbe. Six-Br Truss Probem: Identifiction Error Vrying Dt Error Prmeter for. IDENTIFICATION OF AXIAL FORCES ON STATICALLY INDETERMINATE PIN-JOINTED TRUSSES f v f o 4 5 () (b) (c) Fig. (4). Six-br truss probem: stress-free configurtion (), current configurtion (b) nd test od (c). Figur e 4. Six-br truss probem: stress-free configurtion (), current configurtion f v (b) nd test od (c). f o Fig. (5). Six-br truss probem: uxiiry schemes for the ccutions of (t) nd (s) () (b) (c)

8 Identifiction of Axi Forces on Stticy Indeterminte Pin-Jointed Trusses The Open Civi Engineering Journ, 0, Voume 7 57 Fig. (6). Six-br truss probem: uxiiry schemes for the ccution of mtrix. Figur e 6. Six-br truss probem: uxiiry schemes for the ccution ACKNOWLEDGEMENTS Some future deveopments re syntheticy isted beow:. In this work we hve considered ony stticy indeterminte nd kinemticy determined pin-jointed trusses, it ppers interesting to extend this pproch to prestressed mechnisms.. Ony simpe probems were considered here. Probems of medium nd rge dimensions, speciy from the technic point of view, require prticur cre nd specific strtegies to be computtiony efficient.. Athough the fexibiity method ppers to be more promising s frmework to identify the stress-stte, the cpbiities of the stiffness method shoud be dequtey expored. 4. Impementtions of more sophisticted constitutive ws s briefy described in Section, hve to be expored to understnd computtion costs nd dvntges. 5. Pin-jointed trusses were considered importnt due to their technic ppictions nd so for the simpicity of the mthemtic mode which used to study them. Viewing the procedure presented here in finite eement frmework, sm effort woud be dequte for the ntur extension of the resuts presented here. 6. As it is we known, the soution of n identifiction probem necessriy requires specific strtegies cpbe of fitering errors which cnnot be eiminted, e.g. those on the dt. Such errors occur from vrious cuses for e.g. inccurcy of the dt, pproximtion of the mthemtic mode nd sometimes due to the numeric mode nd rough ccurcy of the soution gorithm. Athough there re numerous gener purpose strtegies tht cn be used to reduce the forementioned errors but specificy designed procedure hs its own benefits s it cn de with specific probem. 7. Finy, it woud be interesting to expore the resuts obtined from geometricy noniner mode. CONFLICT OF INTEREST The uthors confirm tht this rtice content hs no conficts of interest. () (b) (c) (d) None decred. REFERENCES [] Y. C. Fung, A First Course in Continuum Mechnics, Prentice- H: USA, 994. [] A. Hoger, On the determintion of residu stress in n estic body, Journ of Esticity, vo. 6, pp. 0-4, 986. [] B. C. P. Burke, S. O. Kim nd K. S. Kim, Prti por decomposition inverse method ppied to determintion of intern stresses in n estic compex structures, Interntion Journ of Soids nd Structures, vo. 44, pp.00-00, 007. [4] F. Zhng, A. J. Kssb nd D. W. Nichoson, A boundry eement soution of n inverse esticity probem nd ppiction to determining residu stress nd contct stress, Interntion Journ of Soids nd Structures, vo. 4, no. 6, pp , 997. [5] H. D. Bui, Inverse Probems in the Mechnics of Mteris: An Introduction, CRC Press: USA, 994. [6] A. Puro, Mgnetophotoesticity s prmetric tensor fied tomogrphy, Inverse Probems, vo. 4, pp. 5-0, 998. [7] R. L. Robertson, Boundry identifibiity of residu stress vi the Dirichet to Neumnn mp, Inverse Probems, vo., pp. 07-9, 997. [8] C. S. Mn, Hrtig's w nd iner esticity with initi stress, Inverse Probems, vo. 4, pp. -9, 998. [9] R. L. Robertson, Determining residu stress from boundry mesurements, Journ of Esticity, vo. 5, pp. 6-7, 998. [0] G. E. B. Tn nd S. Peegrino, Non-iner dynmic identifiction: An ppiction to prestressed cbe structures, Journ of Sound nd Vibrtion, vo. 08, no., pp. -45, 997. [] S.Peegrino, Anysis of prestressed mechnism, Interntion Journ of Soids nd Structures, vo.6. no., 9-50, 990. [] S.Guest, The stiffness of prestressed frmeworks: A unifying pproch, Interntion Journ of Soids nd Structures, vo. 4, pp , 006. [] R. K. Livesey, Mtrix method of structur nysis, nd ed. Pergmon Press: Oxford, 975. [4] J. F. Be, The experiment foundtions of soid mechnichs, vo. I, Springer-Verg: Berin, 97. [5] A. N. Tikhonov nd V. Y. Arsenin, Soution of i-posed probems, John Wiey & Sons: New York, 977. Received: December 0, 0 Revised: December 8, 0 Accepted: December 8, 0 Emiio Turco; Licensee Benthm Open. This is n open ccess rtice icensed under the terms of the Cretive Commons Attribution Non-Commerci License ( by-nc/.0/) which permits unrestricted, non-commerci use, distribution nd reproduction in ny medium, provided the work is propery cited.

1. The vibrating string problem revisited.

1. The vibrating string problem revisited. Weeks 7 8: S eprtion of Vribes In the pst few weeks we hve expored the possibiity of soving first nd second order PDEs by trnsforming them into simper forms ( method of chrcteristics. Unfortuntey, this

More information

Introduction to statically indeterminate structures

Introduction to statically indeterminate structures Sttics of Buiding Structures I., EASUS Introduction to stticy indeterminte structures Deprtment of Structur echnics Fcuty of Civi Engineering, VŠB-Technic University of Ostrv Outine of Lecture Stticy indeterminte

More information

Math 124B January 24, 2012

Math 124B January 24, 2012 Mth 24B Jnury 24, 22 Viktor Grigoryn 5 Convergence of Fourier series Strting from the method of seprtion of vribes for the homogeneous Dirichet nd Neumnn boundry vue probems, we studied the eigenvue probem

More information

Fluid Flow through a Tube

Fluid Flow through a Tube . Theory through Tube In this experiment we wi determine how we physic retionship (so ced w ), nmey Poiseue s eqution, ppies. In the suppementry reding mteri this eqution ws derived s p Q 8 where Q is

More information

Statistical Physics. Solutions Sheet 5.

Statistical Physics. Solutions Sheet 5. Sttistic Physics. Soutions Sheet 5. Exercise. HS 04 Prof. Mnfred Sigrist Ide fermionic quntum gs in hrmonic trp In this exercise we study the fermionic spiness ide gs confined in three-dimension hrmonic

More information

MAGIC058 & MATH64062: Partial Differential Equations 1

MAGIC058 & MATH64062: Partial Differential Equations 1 MAGIC58 & MATH646: Prti Differenti Equtions 1 Section 4 Fourier series 4.1 Preiminry definitions Definition: Periodic function A function f( is sid to be periodic, with period p if, for, f( + p = f( where

More information

Development of the Sinc Method for Nonlinear Integro-Differential Eequations

Development of the Sinc Method for Nonlinear Integro-Differential Eequations Austrin Journ of Bsic nd Appied Sciences, 4(): 558-555, ISS 99-878 Deveopment of the Sinc Method for oniner Integro-Differenti Eequtions K. Jei, M. Zrebni, 3 M. Mirzee Chi,3 Ismic Azd University Brnch

More information

P R A S I C ' 02. Simpozionul na.ional cu participare interna.ional/ PRoiectarea ASIstat/ de Calculator

P R A S I C ' 02. Simpozionul na.ional cu participare interna.ional/ PRoiectarea ASIstat/ de Calculator UNIVERSITTE TRNSILVNI DIN R0OV tedr Design de Produs i Robotic Simpozionu n.ion cu prticipre intern.ion/ PRoiectre SIstt/ de cutor P R S I ' 0 Vo. I Mecnisme i Triboogie 7-8 Noiembrie rov Români ISN 97-65-064-9

More information

Complete Description of the Thelen2003Muscle Model

Complete Description of the Thelen2003Muscle Model Compete Description o the he23usce ode Chnd John One o the stndrd musce modes used in OpenSim is the he23usce ctutor Unortuntey, to my knowedge, no other pper or document, incuding the he, 23 pper describing

More information

arxiv: v1 [math.co] 5 Jun 2015

arxiv: v1 [math.co] 5 Jun 2015 First non-trivi upper bound on the circur chromtic number of the pne. Konstnty Junosz-Szniwski, Fcuty of Mthemtics nd Informtion Science, Wrsw University of Technoogy, Pond Abstrct rxiv:1506.01886v1 [mth.co]

More information

In this appendix, we evaluate the derivative of Eq. 9 in the main text, i.e., we need to calculate

In this appendix, we evaluate the derivative of Eq. 9 in the main text, i.e., we need to calculate Supporting Tet Evoution of the Averge Synptic Updte Rue In this ppendi e evute the derivtive of Eq. 9 in the min tet i.e. e need to ccute Py ( ) Py ( Y ) og γ og. [] P( y Y ) P% ( y Y ) Before e strt et

More information

GEOMETRIC PROBABILITY MODELS TO ANALYZE STRATEGIES FOR FINDING A BURIED CABLE

GEOMETRIC PROBABILITY MODELS TO ANALYZE STRATEGIES FOR FINDING A BURIED CABLE Journ of the Opertions Reserch Society of Jpn Vo. 6, No. 3, Juy 17, pp. 4 417 c The Opertions Reserch Society of Jpn GEOMETRIC PROBABILITY MODELS TO ANALYZE STRATEGIES FOR FINDING A BURIED CABLE Ken-ichi

More information

Stability condition of finite difference solution for viscoelastic wave equations

Stability condition of finite difference solution for viscoelastic wave equations Erthq Sci (009): 79 85 79 Doi: 0007/s589-009-079- Stbiity condition of finite difference soution for viscoestic wve equtions Chengyu Sun, Yunfei Xio Xingyo Yin nd Hongcho Peng Chin University of Petroeum,

More information

Combined RKDG and LDG Method for the Simulation of the Bipolar Charge Transport in Solid Dielectrics

Combined RKDG and LDG Method for the Simulation of the Bipolar Charge Transport in Solid Dielectrics PIERS ONLINE, VOL. 5, NO. 1, 009 51 Combined RKDG nd LDG Method for the Simution of the Bipor Chrge Trnsport in Soid Dieectrics Jihun Tin, Jun Zou, nd Jinsheng Yun Deprtment of Eectric Engineering, Tsinghu

More information

A Slipping and Buried Strike-Slip Fault in a Multi-Layered Elastic Model

A Slipping and Buried Strike-Slip Fault in a Multi-Layered Elastic Model Geosciences 7, 7(): 68-76 DOI:.59/j.geo.77. A Sipping nd Buried Strike-Sip Fut in Muti-Lyered Estic Mode Asish Krmkr,*, Snjy Sen Udirmpur Pisree Sikshytn (H.S.), Udirmpur, P.O. Knyngr, Pin, Indi Deprtment

More information

Generalized Uncertainty Principle, Extra-dimensions and Holography

Generalized Uncertainty Principle, Extra-dimensions and Holography 470 Brziin Journ of Physics, vo. 35, no. 2B, June, 2005 Generized Uncertinty Principe, Extr-dimensions nd Hoogrphy Fbio Scrdigi CENTRA - Deprtmento de Fisic, Instituto Superior Tecnico, Av. Rovisco Pis

More information

Game Characterizations of Process Equivalences

Game Characterizations of Process Equivalences Gme Chrcteriztions of Process Equivences Xin Chen nd Yuxin Deng Deprtment of Computer Science nd Engineering, Shnghi Jio Tong University, Chin Abstrct In this pper we propose hierrchy of gmes tht ows us

More information

Linear Systems COS 323

Linear Systems COS 323 Liner Systems COS 33 Soving Liner Systems of Equtions Define iner system Singurities in iner systems Gussin Eimintion: A gener purpose method Nïve Guss Guss with pivoting Asymptotic nysis Tringur systems

More information

Bending Analysis of Castellated Beams. By Sahar Elaiwi Boksun Kim Long-Yuan Li

Bending Analysis of Castellated Beams. By Sahar Elaiwi Boksun Kim Long-Yuan Li Athens Journ of Technoogy nd Engineering X Y Bending Anysis of Csteted Bems By Shr Eiwi Boksun Kim Long-Yun Li Existing studies hve shown tht the od-crrying cpcity of csteted bems cn be infuenced by the

More information

Energy Balance of Solar Collector

Energy Balance of Solar Collector Renewbe Energy Grou Gret Ides Grow Better Beow Zero! Wecome! Energy Bnce of Sor Coector Mohmd Khrseh E-mi:m.Khrseh@gmi.com Renewbe Energy Grou Gret Ides Grow Better Beow Zero! Liuid Ft Pte Coectors. Het

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com M Dynmics - Dmped nd forced hrmonic motion. A P α B A ight estic spring hs ntur ength nd moduus of esticity mg. One end of the spring is ttched to point A on pne tht is incined to the horizont t n nge

More information

ANALYZING EFFECT OF TEMPERATURE VARIATION ON SKIN SUB-LAYERS USING FINITE ELEMENT APPROACH

ANALYZING EFFECT OF TEMPERATURE VARIATION ON SKIN SUB-LAYERS USING FINITE ELEMENT APPROACH Origin Reserch Artice Aied sciences Interntion Journ of Phrm nd Bio Sciences ISSN 0975-699 ANALYZING EFFECT OF TEMPERATURE VARIATION ON SKIN SUB-LAYERS USING FINITE ELEMENT APPROACH MONA DWIVEDI 1 * AND

More information

Homework Assignment #5 Solutions

Homework Assignment #5 Solutions Physics 506 Winter 008 Textbook probems: Ch. 9: 9., 9., 9.4 Ch. 10: 10.1 Homework Assignment #5 Soutions 9. A spheric hoe of rdius in conducting medium cn serve s n eectromgnetic resonnt cvity. Assuming

More information

Mutipy by r sin RT P to get sin R r r R + T sin (sin T )+ P P = (7) ffi So we hve P P ffi = m (8) choose m re so tht P is sinusoi. If we put this in b

Mutipy by r sin RT P to get sin R r r R + T sin (sin T )+ P P = (7) ffi So we hve P P ffi = m (8) choose m re so tht P is sinusoi. If we put this in b Topic 4: Lpce Eqution in Spheric Co-orintes n Mutipoe Expnsion Reing Assignment: Jckson Chpter 3.-3.5. Lpce Eqution in Spheric Coorintes Review of spheric por coorintes: x = r sin cos ffi y = r sin sin

More information

Dynamic Analysis of the Turnout Diverging Track for HSR with Variable Curvature Sections

Dynamic Analysis of the Turnout Diverging Track for HSR with Variable Curvature Sections Word Journ of Engineering nd Technoogy, 07, 5, 4-57 http://www.scirp.org/journ/wjet ISSN Onine: -449 ISSN Print: -4 Dynmic Anysis of the Turnout Diverging Trck for HSR with Vribe Curvture Sections Wdysw

More information

FRACTURE OF PIEZOELECTRIC MATERIALS

FRACTURE OF PIEZOELECTRIC MATERIALS CF000OR FRACUR OF PZOCRC MARAS ong-yi Zhng* nd Mingho Zho eprtment of Mechnic ngineering Hong ong University of Science nd echnoogy Cer Wter By owoon Hong ong Chin * -mi: mezhngt@ust.hk ABSRAC he present

More information

Where does Oxygen Extinction Occur in a Soil Profile?

Where does Oxygen Extinction Occur in a Soil Profile? 2st Interntion Congress on Modeing nd Simution, God Cost, Austri, 29 Nov to 4 ec 25 www.mssnz.org.u/modsim25 Where does Oxygen Extinction Occur in Soi Profie? Freemn J Cook nd John H Knightb Freemn Cook

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

L v. removal. elastic. the body is. Hooke s force. [ M L 1 T 2 ] (1) (2) (3) Normal Stress Tensile Stress. stress. parallel. Shearing.

L v. removal. elastic. the body is. Hooke s force. [ M L 1 T 2 ] (1) (2) (3) Normal Stress Tensile Stress. stress. parallel. Shearing. Esticity Definition Esticity: A wire is cmped t one end nd oded t its free end. It is found tht the ength of the wire chnges. The force is known s deforming force, the chnges known s deformtion. If fter

More information

SIMULATION OF WOOD DRYING STRESSES USING CVFEM

SIMULATION OF WOOD DRYING STRESSES USING CVFEM Ltin Americn Appied Reserch 4:2- (2) SIMULATION OF WOOD DRYING STRESSES USING CVFEM C. H. SALINAS, C. A. CHAVEZ, Y. GATICA nd R.A. ANANIAS Deprtmento de Ingenierí Mecánic, Universidd de Bío-Bío, Av. Coo

More information

Coordinate Geometry. Coordinate Geometry. Curriculum Ready ACMNA: 178, 214, 294.

Coordinate Geometry. Coordinate Geometry. Curriculum Ready ACMNA: 178, 214, 294. Coordinte Geometr Coordinte Geometr Curricuum Red ACMNA: 78, 4, 94 www.mthetics.com Coordinte COORDINATE Geometr GEOMETRY Shpes ou ve seen in geometr re put onto es nd nsed using gebr. Epect bit of both

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

SUMMER KNOWHOW STUDY AND LEARNING CENTRE SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18

More information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

arxiv: v1 [math.na] 28 Dec 2017

arxiv: v1 [math.na] 28 Dec 2017 An dptive hp-refinement strtegy with computbe gurnteed bound on the error reduction fctor Ptrik Dnie Aexndre Ern Iin Smers Mrtin Vohrík December 9, 07 rxiv:7.098v [mth.na] 8 Dec 07 Abstrct We propose new

More information

Lecture 6 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell

Lecture 6 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell Lecture 6 Notes, Eectrognetic Theory I Dr. Christopher S. Bird University of Msschusetts Lowe. Associted Legendre Poynois - We now return to soving the Lpce eqution in spheric coordintes when there is

More information

Spanning tree congestion of some product graphs

Spanning tree congestion of some product graphs Spnning tree congestion of some product grphs Hiu-Fi Lw Mthemticl Institute Oxford University 4-9 St Giles Oxford, OX1 3LB, United Kingdom e-mil: lwh@mths.ox.c.uk nd Mikhil I. Ostrovskii Deprtment of Mthemtics

More information

Figure XX.1.1 Plane truss structure

Figure XX.1.1 Plane truss structure Truss Eements Formution. TRUSS ELEMENT.1 INTRODUTION ne truss struture is ste struture on the sis of tringe, s shown in Fig..1.1. The end of memer is pin juntion whih does not trnsmit moment. As for the

More information

Suggested Solution to Assignment 5

Suggested Solution to Assignment 5 MATH 4 (5-6) prti diferenti equtions Suggested Soution to Assignment 5 Exercise 5.. () (b) A m = A m = = ( )m+ mπ x sin mπx dx = x mπ cos mπx + + 4( )m 4 m π. 4x cos mπx dx mπ x cos mπxdx = x mπ sin mπx

More information

Dynamic analysis of an underwater suspension system in ocean currents with method of averaging

Dynamic analysis of an underwater suspension system in ocean currents with method of averaging Indin Journ of Geo Mrine Sciences Vo. (),November 7, pp. -7 Dynmic nysis of n underwter suspension system in ocen currents with method of verging Xing Li,, Jifeng ui,, Min Zho,,* & Tong Ge, Stte Key Lbortory

More information

A Fully-coupled Thermo-Hydro-Mechanical Model for the Description of the Behavior of Swelling Porous Media

A Fully-coupled Thermo-Hydro-Mechanical Model for the Description of the Behavior of Swelling Porous Media Act Poytechnic Hungric Vo. 8, No. 4, 211 A Fuy-couped Thermo-Hydro-Mechnic Mode for the Description of the Behvior of Sweing Porous Medi Hnifi Missoum, Ndi Lredj, rim Bendni, Mustph Miki Construction,

More information

CBE 291b - Computation And Optimization For Engineers

CBE 291b - Computation And Optimization For Engineers The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn

More information

SPECTRAL APPROXIMATION OF TIME-HARMONIC MAXWELL EQUATIONS IN THREE-DIMENSIONAL EXTERIOR DOMAINS. 1. Introduction

SPECTRAL APPROXIMATION OF TIME-HARMONIC MAXWELL EQUATIONS IN THREE-DIMENSIONAL EXTERIOR DOMAINS. 1. Introduction PECTRAL APPROXIMATION OF TIME-HARMONIC MAXWELL EQUATION IN THREE-DIMENIONAL EXTERIOR DOMAIN LINA MA 1, JIE HEN 1 AND LI-LIAN WANG 2 Abstrct. We deveop in this pper n efficient nd robust spectr-gerkin method

More information

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE Determintion RevAdvMterSci of mechnicl 0(009) -7 properties of nnostructures with complex crystl lttice using DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING

More information

Chapter 5 Bending Moments and Shear Force Diagrams for Beams

Chapter 5 Bending Moments and Shear Force Diagrams for Beams Chpter 5 ending Moments nd Sher Force Digrms for ems n ddition to illy loded brs/rods (e.g. truss) nd torsionl shfts, the structurl members my eperience some lods perpendiculr to the is of the bem nd will

More information

Fourier Series in Complex notation. cos(x) = eix + e ix 2. A n cos + B n sin l. i 2 B n. e inx=l + A n + ib n 2 8. ( 0 m 6= n. C n = 1 2l.

Fourier Series in Complex notation. cos(x) = eix + e ix 2. A n cos + B n sin l. i 2 B n. e inx=l + A n + ib n 2 8. ( 0 m 6= n. C n = 1 2l. Fourier Series in Compex nottion sin(x) = eix e ix i = i eix e ix cos(x) = eix + e ix So So '(x) = A 1 0 + nx nx A n cos + B n sin = A 1 0 + e inx= + e inx= A n = A 0 + 1 = C n = C n 1 n= 1 A n ib n C

More information

ELASTICITY AND PLASTICITY I

ELASTICITY AND PLASTICITY I ELASTICITY AD PLASTICITY I Litertur in Czech Ing. Lenk Lusová LPH 407/1 te. 59 73 136 enk.usov@vs.cz http://fst10.vs.cz/usov Literture: Higgeer Sttics nd mechnics of mteri, USA 199 Beer, Johnston, DeWof,

More information

Predict Global Earth Temperature using Linier Regression

Predict Global Earth Temperature using Linier Regression Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

Student Activity 3: Single Factor ANOVA

Student Activity 3: Single Factor ANOVA MATH 40 Student Activity 3: Single Fctor ANOVA Some Bsic Concepts In designed experiment, two or more tretments, or combintions of tretments, is pplied to experimentl units The number of tretments, whether

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

Scientific notation is a way of expressing really big numbers or really small numbers.

Scientific notation is a way of expressing really big numbers or really small numbers. Scientific Nottion (Stndrd form) Scientific nottion is wy of expressing relly big numbers or relly smll numbers. It is most often used in scientific clcultions where the nlysis must be very precise. Scientific

More information

Improper Integrals, and Differential Equations

Improper Integrals, and Differential Equations Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted

More information

Chapters 4 & 5 Integrals & Applications

Chapters 4 & 5 Integrals & Applications Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions

More information

In-Class Problems 2 and 3: Projectile Motion Solutions. In-Class Problem 2: Throwing a Stone Down a Hill

In-Class Problems 2 and 3: Projectile Motion Solutions. In-Class Problem 2: Throwing a Stone Down a Hill MASSACHUSETTS INSTITUTE OF TECHNOLOGY Deprtment of Physics Physics 8T Fll Term 4 In-Clss Problems nd 3: Projectile Motion Solutions We would like ech group to pply the problem solving strtegy with the

More information

KINEMATICS OF RIGID BODIES

KINEMATICS OF RIGID BODIES KINEMTICS OF RIGID ODIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. Description

More information

Acceptance Sampling by Attributes

Acceptance Sampling by Attributes Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

Numerical Integration

Numerical Integration Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

On the Generalized Weighted Quasi-Arithmetic Integral Mean 1

On the Generalized Weighted Quasi-Arithmetic Integral Mean 1 Int. Journl of Mth. Anlysis, Vol. 7, 2013, no. 41, 2039-2048 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/10.12988/ijm.2013.3499 On the Generlized Weighted Qusi-Arithmetic Integrl Men 1 Hui Sun School

More information

221B Lecture Notes WKB Method

221B Lecture Notes WKB Method Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using

More information

Finite Element Determination of Critical Zones in Composite Structures

Finite Element Determination of Critical Zones in Composite Structures Finite Element Determintion of Criticl Zones in Composite Structures Alexey I. Borovkov Dmitriy V. Klimshin Denis V. Shevchenko Computtionl Mechnics Lb., St. Petersburg Stte Polytechnicl University, Russi

More information

An approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95

An approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95 An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further

More information

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b

More information

Math& 152 Section Integration by Parts

Math& 152 Section Integration by Parts Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible

More information

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), ) Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s

More information

V. More General Eigenfunction Series. A. The Sturm-Liouville Family of Problems

V. More General Eigenfunction Series. A. The Sturm-Liouville Family of Problems V. More Gener Eigenfunction Series c 214, Phiip D. Loewen A. The Sturm-Liouvie Fmiy of Probems Given n re interv, b], functions p(x), q(x), r(x), nd constnts c, c 1, d, d 1, consider (ODE) (p(x)y (x))

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

1 1D heat and wave equations on a finite interval

1 1D heat and wave equations on a finite interval 1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion

More information

VI. The Heat Equation. A. The Diffusion Equation Conservation in Action

VI. The Heat Equation. A. The Diffusion Equation Conservation in Action VI. The Het Eqution c 214, Phiip D. Loewen A. The Diffusion Eqution Conservtion in Action The simpest PDE describing diffusion (of het, perfume, poutnt) ong the x-xis is u t = α 2 u xx. ( ) Theme: Concvity

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

Deteriorating Inventory Model for Waiting. Time Partial Backlogging

Deteriorating Inventory Model for Waiting. Time Partial Backlogging Applied Mthemticl Sciences, Vol. 3, 2009, no. 9, 42-428 Deteriorting Inventory Model for Witing Time Prtil Bcklogging Nit H. Shh nd 2 Kunl T. Shukl Deprtment of Mthemtics, Gujrt university, Ahmedbd. 2

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

Arithmetic & Algebra. NCTM National Conference, 2017

Arithmetic & Algebra. NCTM National Conference, 2017 NCTM Ntionl Conference, 2017 Arithmetic & Algebr Hether Dlls, UCLA Mthemtics & The Curtis Center Roger Howe, Yle Mthemtics & Texs A & M School of Eduction Relted Common Core Stndrds First instnce of vrible

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

Riemann is the Mann! (But Lebesgue may besgue to differ.)

Riemann is the Mann! (But Lebesgue may besgue to differ.) Riemnn is the Mnn! (But Lebesgue my besgue to differ.) Leo Livshits My 2, 2008 1 For finite intervls in R We hve seen in clss tht every continuous function f : [, b] R hs the property tht for every ɛ >

More information

LECTURE 14. Dr. Teresa D. Golden University of North Texas Department of Chemistry

LECTURE 14. Dr. Teresa D. Golden University of North Texas Department of Chemistry LECTURE 14 Dr. Teres D. Golden University of North Texs Deprtment of Chemistry Quntittive Methods A. Quntittive Phse Anlysis Qulittive D phses by comprison with stndrd ptterns. Estimte of proportions of

More information

Matrices, Moments and Quadrature, cont d

Matrices, Moments and Quadrature, cont d Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl

More information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

Continuous Random Variables

Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht

More information

TRANSFUSION OF TIME-DOMAIN AND FREQUENCY-DOMAIN EDDY

TRANSFUSION OF TIME-DOMAIN AND FREQUENCY-DOMAIN EDDY 4th Midde Est ND Conference nd Ehibition, Kingdom of Bhrin, Dec 007 For ppers of this pubiction cick: www.ndt.net/serch/docs.php3?minsource57 RANSFUSON OF ME-DOMAN AND FREQUENCY-DOMAN EDDY CURREN SGNAS

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

7.2 The Definite Integral

7.2 The Definite Integral 7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where

More information

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below. Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

Heat flux and total heat

Heat flux and total heat Het flux nd totl het John McCun Mrch 14, 2017 1 Introduction Yesterdy (if I remember correctly) Ms. Prsd sked me question bout the condition of insulted boundry for the 1D het eqution, nd (bsed on glnce

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

Lecture Notes PH 411/511 ECE 598 A. La Rosa Portland State University INTRODUCTION TO QUANTUM MECHANICS

Lecture Notes PH 411/511 ECE 598 A. La Rosa Portland State University INTRODUCTION TO QUANTUM MECHANICS Lecture Notes PH 4/5 ECE 598. L Ros Portlnd Stte University INTRODUCTION TO QUNTUM MECHNICS Underlying subject of the PROJECT ssignment: QUNTUM ENTNGLEMENT Fundmentls: EPR s view on the completeness of

More information

arxiv:math/ v2 [math.ho] 16 Dec 2003

arxiv:math/ v2 [math.ho] 16 Dec 2003 rxiv:mth/0312293v2 [mth.ho] 16 Dec 2003 Clssicl Lebesgue Integrtion Theorems for the Riemnn Integrl Josh Isrlowitz 244 Ridge Rd. Rutherford, NJ 07070 jbi2@njit.edu Februry 1, 2008 Abstrct In this pper,

More information

Module 6: LINEAR TRANSFORMATIONS

Module 6: LINEAR TRANSFORMATIONS Module 6: LINEAR TRANSFORMATIONS. Trnsformtions nd mtrices Trnsformtions re generliztions of functions. A vector x in some set S n is mpped into m nother vector y T( x). A trnsformtion is liner if, for

More information

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses

More information

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,

More information

For convenience, we rewrite m2 s m2 = m m m ; where m is repeted m times. Since xyz = m m m nd jxyj»m, we hve tht the string y is substring of the fir

For convenience, we rewrite m2 s m2 = m m m ; where m is repeted m times. Since xyz = m m m nd jxyj»m, we hve tht the string y is substring of the fir CSCI 2400 Models of Computtion, Section 3 Solutions to Homework 4 Problem 1. ll the solutions below refer to the Pumping Lemm of Theorem 4.8, pge 119. () L = f n b l k : k n + lg Let's ssume for contrdiction

More information

HPLP310 Biblio_35 Fissures radial intern in a thick cylinder under pressure and thermal loading

HPLP310 Biblio_35 Fissures radial intern in a thick cylinder under pressure and thermal loading defult Titre : HPLP310 - Biblio_35 Fissure rdile interne dns u[...] Dte : 0/09/011 Pge : 1/15 Responsble : TRAN Vn Xun Clé : V7.0.310 Révision : HPLP310 Biblio_35 Fissures rdil intern in thick cylinder

More information

A Study on Prediction of Cavitation for Centrifugal Pump Myung Jin Kim, Hyun Bae Jin, and Wui Jun Chung

A Study on Prediction of Cavitation for Centrifugal Pump Myung Jin Kim, Hyun Bae Jin, and Wui Jun Chung Word Acdemy of Science, Engineering nd Technoogy Interntion Journ of Mechnic nd Mechtronics Engineering Vo:6, No:12, 2012 A Study on Prediction of Cvittion for Centrifug Pump Myung Jin Kim, Hyun Be Jin,

More information

arxiv: v1 [math.ra] 1 Nov 2014

arxiv: v1 [math.ra] 1 Nov 2014 CLASSIFICATION OF COMPLEX CYCLIC LEIBNIZ ALGEBRAS DANIEL SCOFIELD AND S MCKAY SULLIVAN rxiv:14110170v1 [mthra] 1 Nov 2014 Abstrct Since Leibniz lgebrs were introduced by Lody s generliztion of Lie lgebrs,

More information