Study of key algorithms in topology optimization
|
|
- Priscilla Quinn
- 5 years ago
- Views:
Transcription
1 Int J Adv Manuf Technol (2007) 32: DOI /s ORIGINAL ARTICLE Kong-Tan Zuo L-Png Chen Yun-Qng Zhang Jngzhou Yang Study of key algorthms n topology optmzaton Receved: 30 Aprl 2004 / Accepted: 22 Aprl 2005 / Publshed onlne: 21 March 2006 # Sprnger-Verlag London Lmted 2006 Abstract The theory of topology optmzaton based on the sold sotropc materal wth penalzaton model (SIMP) method s thoroughly analyzed n ths paper In order to solve complcated topology optmzaton problems, a hybrd soluton algorthm based on the method of movng asymptotes (MMA) approach and the globally convergent verson of the method of movng asymptotes (GCMMA) approach s proposed The numercal nstablty, whch always leads to a non-manufacturng result n topology optmzaton, s analyzed, along wth current methods to control t To elmnate the numercal nstablty of topology results, a convoluton ntegral factor method s ntroduced Meanwhle, an teraton procedure based on the hybrd soluton algorthm and a method to elmnate numercal nstablty are developed The proposed algorthms are verfed wth llustratve examples The effect and functon of the hybrd soluton algorthm and the convoluton radus n optmzaton are also dscussed Keywords Topology optmzaton MMA seres algorthms Hybrd soluton algorthm Numercal nstablty Convoluton ntegral factor method 1 Introducton Topology optmzaton s consdered to be one of the most challengng felds n structural optmzaton Topology optmzaton has been an nterestng research topc snce Bendsoe and Kkuch [1] ntroduced the homogenzaton K-T Zuo L-P Chen Y-Q Zhang Center for Computer-Aded Desgn, Huazhong Unversty of Scence & Technology, Wuhan, Hube, , People s Republc of Chna e-mal: zuokt@hustcadcom J Yang (*) Center for Computer-Aded Desgn, The Unversty of Iowa, Iowa Cty, IA 52246, USA e-mal: yang@engneernguowaedu Tel: Fax: method nto the feld After more than 20 years of development, researchers have made a great deal of progress n theores as well as n real engneerng applcatons Recent research n topology optmzaton has focused on such topcs as: materal nterpolaton methods, optmal algorthms sutable for topology optmzaton, methods to elmnate numercal nstabltes, and the applcaton of topology optmzaton n engneerng Materal nterpolaton used n topology optmzaton manly nclude the followng two types of methods: the homogenzaton method [1, 2], prmarly used n theory dervatons, such as the exstence of solutons, nvestgatng numercal nstabltes, etc; and the materal densty method, sometmes named the sold sotropc materal wth penalzaton model (SIMP) method, proposed by Mlenek and Schrrmacher [3], Sgmund [4], and Bendsoe and Sgmund [5] The SIMP method, whch s used n ths paper, has been accepted by most researchers and engneers, and has been successfully used n many engneerng felds At present, two types of optmzaton algorthms are used n topology optmzaton: the optmalty crtera (OC) approach and the method of movng asymptotes (MMA) seres approach The former [6, 2] s usually deduced by a Lagrange functon composed of obectve and constrant functons accordng to the Kuhn-Tucker condton Ths approach has good convergence because t s based on heurstc formulaton and does not need the dervatves of functons However, ts dsadvantage s that t has a narrow restrcton n the sngle-constrant condton and s dffcult to use n problems wth multple constrants The MMA approach [7 9, 12] can ft sngle and multple constrants condtons, so t has wder engneerng applcatons, such as the desgn of a complant mechansm, mcroelectromechancal systems (MEMS), and multdscplne optmzaton The dsadvantage of the MMA approach exsts n the bad convergence of calculaton The globally convergent verson of the method of movng asymptotes (GCMMA) has a better convergence as compared to the MMA approach However, the calculaton speed of GCMMA s not perfect A new hybrd algorthm based on MMA and GCMMA s proposed Ths hybrd algorthm combnes the
2 788 advantages of both MMA and GCMMA together and can be wdely used n complcated engneerng problems One mportant focus of topology optmzaton s developng a method to elmnate numercal nstablty, snce t may result n a non-manufacturable structure n engneerng Many methods, such as the permeter control method [16], densty slope control [17], and meshndependent flter [18], have been proposed to deal wth ths problem But most of them have some shortcomngs Ths paper frst revews the SIMP method n topology optmzaton Then, the MMA approach and GCMMA approach based on SIMP s appled to topology optmzaton, and a hybrd soluton algorthm based on both MMA and GCMMA s proposed To elmnate numercal nstablty n the topology optmzaton, a convoluton ntegral factor method s developed, and an teraton procedure ncludng the hybrd soluton algorthm and the convoluton ntegral factor method s presented Fnally, an llustratve example s used to dscuss the effect and functon of the hybrd algorthm and the convoluton ntegral factor method 2 SIMP method n topology optmzaton The general topology optmzaton problem of mnmal complance can be defned as: mnmze : C ¼ F T U subect to : V ¼ f V 0 ¼ XN x e v e V (1a) (1b) and E are the elastc moduluses of an element before and after optmzaton, respectvely, then E ¼ ðx e Þ p E 0 : If k 0 and k e are the element s ntal stffness and the after-optmzaton stffness matrces, respectvely, the followng relatonshp exsts: k e ¼ ðx e Þ p k 0 : The parameter p s a penalty factor, and t s mportant to penalze the mddle densty n order to decrease the number of mddle-densty elements and ensure that most element denstes are close to zero or one Gven the above precondtons, every element has only one desgn varable Compared to the homogenzaton method, the SIMP approach makes excellent progress on decreasng the number of desgn varables Another advantage of ths SIMP method s that the materal characterstc after the change can be wrtten as the exponent functon of the ntal element densty and the ntal materal characterstc, so that ths approach greatly smplfes the soluton of topology optmzaton Snce a relatonshp exsts n the dscrete element: V ¼ f V 0 ¼ PN x e v e ; where v e s the element volume after optmzaton, the topology optmzaton formulaton of the mnmal complance problem based on the SIMP approach wll be: mnmze : C ¼ U T KU ¼ XN subect to : V ¼ f V 0 ¼ XN F ¼ KU u e k e u e ¼ XN x e v e V ðx e Þ p u e k 0 u e F ¼ KU (1c) 0 < x mn x e x max (2) 0 x e 1 (1d) where C s the complance of the structure, F s the force vector, U s the dsplacement vector, K s the stffness matrx of the structure, V 0 s the ntal volume of the structure, V s the structure s volume after optmzaton, and f s the rato of the volume after optmzaton wth the ntal volume The volume constrant and equlbrum equaton of the structure are ncluded In the SIMP approach, the precondtons nclude: 1 The materal characterstcs, such as the elastc modulus, n a dscrete element are constant The desgn varable s the densty of the element, represented wth x e Ifρ 0 s the ntal element densty and ρ s the element densty after optmzaton, then ρ ¼ x e ρ 0 exsts 2 The materal characterstcs n an element should be changed wth the exponent of element densty If E 0 where xmn s the lower bound of the densty, whch s ntroduced to prevent sngularty of the equlbrum problem, and x max s the upper bound of the densty, u e s the element dsplacement, and N s the total number of dscrete elements 3 A new hybrd soluton algorthm for topology optmzaton 31 MMA and GCMMA approach based on the SIMP method The MMA approach, whch was ntally proposed by Svanberg [7, 8], s based on the frst-order Taylor seres expanson of the obectve and constrant functons Wth ths method, an explct convex subproblem s generated to approxmate the mplct nonlnear problem Because the subproblem s separable and convex, a dual approach or a prmal-dual nteror-pont method can be used to solve t
3 Construct the MMA/GCMMA form of ntal problem Choose a ntal calculaton pont Calculaton functon value and ther derves of obectve and constrants Construct MMA/GCMMA subproblem Solve MMA/GCMMA subproblem and fnd the approxmate soluton of ntal problem where: eðkþ f ðxþ ¼f x þ ðkþ Xn Xn ¼1 ð u ðkþ ¼1 p ðkþ ¼ 1; ; m x ðkþ p ðkþ þ q! ðkþ u ðkþ x x l ðkþ þ x ðkþ qðkþ l ðkþ Þ 789 (4b) Update the bound of desgn varables and non-monotonous parameters No [10, 13, 14] The solutons of a sequence of subproblems can converge towards the orgnal problem A general formulaton of nonlnear optmzaton can be wrtten as: mnmze : f 0 ðxþþa 0 z þ XM ¼1 subect to : f ðxþ a z y 0 x mn x x max y 0; z 0 ðc y þ 1 2 d y 2 Þ ¼ 1; ; m ¼ 1; ; n ¼ 1; ; m where the desgn varable x ¼ ðx 1 ; ; x n Þ T 2< n ; y ¼ ðy 1 ; ; y m Þ T 2< m ; z 2<; functons f 0, f 1,, f m are contnuous and dfferentable, x mn and x max are real numbers and satsfy x mn x max, and the coeffcents a 0, a, c 0, and d are postve real numbers The obectve and constrant functons are lnearly expanded at pont 1 u x or 1 x l to construct a subproblem The subproblem s then solved to approxmate the ntal problem The subproblem of Eq 3 can be expressed as: mnmze : e f ðkþ 0 ðxþþa 0 z þ XM ¼1 ðc y þ 1 2 d y 2 Þ subect to : e f ðkþ ðxþ a z y 0 ¼ 1; ; m α ðkþ x β ðkþ y 0; z 0 Convergence? Yes Output result The end Fg 1 Flow chart of the MMA/GCMMA approach ¼ 1; ; n ¼ 1; ; m (3) (4a) 2 x þ x ðkþ þk ðkþ ¼ x ðþ k l ðþ k (4c) x ðþ k þ k ðþ p ðkþ ¼ u ðkþ q ðþ k n o ¼ max x mn ; 0:9l ðkþ þ 0:1x ðkþ ; n o (4d) ¼ mn x max ; 0:9u ðþ k þ 0:1x ðþ k α ðkþ β ðþ þ x ¼ max x x ¼ max (4e) x (k) (k) The upper asymptotes u and lower asymptotes l should be updated wth teraton as follows: l ðkþ ¼ x ðkþ s ðkþ x ðk 1Þ l ðk 1Þ ; (4f) u ðkþ ¼ x ðkþ þ s ðkþ u ðk 1Þ x ðk 1Þ Accordng to the ðx ðkþ Þ; only one of p (k) and q (k) s not zero at the same tme Therefore, the approxmaton of the subproblem to the orgnal problem n MMA s monotonous The GCMMA approach, an extenson of the MMA approach, was proposed by Svanberg [9, 11] In GCMMA, p (k) and q (k) are not zero smultaneously, and a nonmonotonous parameter ρ (k) s ntroduced, so that the followng relatons exst: p ðkþ ¼ u ðkþ x ¼ x ðþ k l ðþ q ðþ k! þ! x ðkþ þ ρðkþ u ðkþ l ðkþ 2! þ ρ ðþ k! x ðkþ 2 u ðþ k ð Þ l k (5a)
4 790 The non-monotonous parameter ρ (k) can be updated as follows: ρ ðkþ1þ ¼ 2ρ ðkþ f g ðkþ x ðkþ1þ < f x ðkþ1þ (5b) ρ ðkþ1þ ¼ ρ ðkþ f g ðkþ x ðkþ1þ f x ðkþ1þ The startng pont of ρ (0) can be selected as: ρ ð0þ ¼ 0 f ρ ð0þ 0 ¼ 1 X x ð0þ =@x 5n x max x mn ρ ð0þ ¼1 f ρ ð0þ > 0 (5c) Therefore, the GCMMA s a non-monotonous approxmaton of the ntal problem The soluton procedure of the MMA/GCMMA approach s shown n Fg 1 For mnmal the complance problem of topology optmzaton, we can construct the general formula n MMA/GCMMA as follows: mnmze : subect to : X N X N ðx e Þ p u T e k 0u e þ z þ 1000y 1 x e V 0 f y x e 1 e ¼ 1; ; N y 1 0; z 0 It s obvous that the topology optmzaton problem s a specal case of the nonlnear programmng problem n Eq 3 Appylng the followng relaton to Eq 3: f 0 ¼ XN (6) ðx e Þ p u T e k 0u e ; f 1 ¼ XN x e V 0 f (7) and a 0 =1, c 1 =1000, d 1 =0, and a 1 =1, yelds the problem n Eq 6 The obectve and constrant functon values are then calculated, as well as ther dervatves In the SIMP approach, the dervatves can be calculated as 0 ¼ pðx e Þ p 1 u T k 0u e ; 2 f 2 ¼ pðp 1Þðx e Þ p 2 u T e k 0u e 1 ¼ f 1 ¼ 2 e (8a) (8b) We can then construct the MMA/GCMMA subproblem based on Eq 6 and solve the subproblem to obtan the densty dstrbuton of the structure topology A sequence of solutons of the MMA/GCMMA subproblems gradually converges toward the global optmal of the densty dstrbuton 32 A new hybrd MMA GCMMA algorthm Some numercal tests [9] show that the GCMMA method s globally convergent, but ths has been proved only n sze optmzaton Applcatons and verfcatons for convergence and calculaton effcency of the GCMMA approach for large-scale topology optmzaton problems have not been reported It can be concluded mathematcally [15] that the optmzaton process can converge very quckly when a monotonous method, such as a one-order dervatve approxmaton, s used to approxmate a complcated problem at the desgn pont far away from the optmum, whereas the optmzaton process converges slowly f monotonous approxmaton s used at the desgn pont near the optmum On the other hand, when a nonmonotonous method, such as a second-order dervatve approxmaton, s used to approxmate a complcated problem, a stable and convergent soluton can be obtaned f the ntal desgn pont s properly selected, whereas calculaton may be slow f a bad ntal desgn pont far away from the optmum s used The topology optmzaton problem often exhbts hgh non-lnearty and a non-monotonous structural behavor The characterstc of the obectve and constrant functons s very complcated, and usually conssts of a non-convex and mplct form of desgn varables The number of desgn varables s very large MMA nvolves monotonous approxmaton to solve the problem, whereas GCMMA nvolves non-monotonous approxmaton Some research [12] has shown that t s mportant to select an approprate approxmaton scheme for calculatons The convergence process can be accelerated f a monotonous approxmaton s used for a monotonous structural response functon or f a non-monotonous approxmaton s used for a non-monotonous functon Oppostely, the convergence process may be slow f a monotonous approxmaton s used for a nonmonotonous structural response functon or a non-monotonous approxmaton s used for a monotonous functon Therefore, a non-monotonous GCMMA approach should be used to approxmate non-monotonous topology optmzaton At the same tme, the fast convergence characterstc of the MMA approach at a desgn pont far away from the optmum ought to be fully utlzed Some calculatons show that, when a smplex MMA approach s used n topology optmzaton, the calculaton can converge quckly and the obectve functons decrease rapdly at the frst few teratons However, the convergence and calculaton speed slow down, and, usually, some numercal oscllaton s exhbted when the desgn pont s near the optmum; fnally, the problem does not converge A smplex GCMMA approach used n topology optmzaton
5 791 can lead to a stable soluton, but the whole convergence process s very slow To deal wth large-scale and complex topology optmzaton, we propose a hybrd algorthm based on MMA and GCMMA, and provde a strategy for swtchng from the monotonous MMA approach to the non-monotonous GCMMA In our hybrd algorthm, the MMA approach s used when the desgn pont s far away from the optmum, where the quck convergence characterstc of the monotonous MMA approach can be fully utlzed After several teratons, the GCMMA approach s used when the desgn pont s near the optmum, where the non-monotonous characterstc of GCMMA can be fully utlzed In short, the hybrd algorthm has the combned advantages of MMA and GCMMA, and ensures calculaton effcency and convergence The flow chart of ths hybrd algorthm s llustrated n Fg 2 An mportant step n the hybrd algorthm s decdng when to swtch the optmzer from MMA to GCMMA As dscussed above, some numercal oscllaton s usually exhbted n topology optmzaton wth a smplex MMA approach when the desgn pont s near the optmum Therefore, the oscllaton value of the obectve functon can be used as a swtch condton The optmzer s swtched from MMA to GCMMA f the oscllaton s large enough The followng condton s checked as a threshold for swtchng from MMA to GCMMA The obectve nformaton of the current step and prevous two steps s used: f 0 x ðk 2Þ f0 x ðk 1Þ ðf 0 ðx ðk 2Þ Þþf 0 ðx ðk 1Þ ÞÞ=2 f 0 x ðk 1Þ f0 x ðkþ ðf 0 ðx ðk 1Þ Þþf 0 ðx ðkþ ÞÞ=2 < δ where δ s a small postve real number Convergence? Yes Output result The end No Defne the ntal desgn varables Fnte element analyss Update desgn varables wth MMA algorthm Numercal oscllaton? Yes Fnte element analyss Yes No Update desgn varables wth GCMMA algorthm Convergence? Output result The end Fg 2 Flow chart of hybrd MMA GCMMA algorthm No (9) 4 A method to elmnate numercal nstablty n topology optmzaton 41 Numercal nstablty n topology optmzaton Some numercal nstablty, such as checkerboard, mesh dependency, and porous materal, are often ncluded n topology optmzaton results, not only for the SIMP method but also for the homogenzaton method These phenomena often lead to a meanngless results because structures wth such flaws are dffcult to manufacture Porous materal means that many elements wth mddle densty are ncluded n the results, as shown n Fg 3a Results wth porous elements cannot be manufactured and are meanngless for engneerng purposes Checkerboard pattern, shown n Fg 3b, ndcates that the elements wth densty of ether 0 or 1 perodcally exst n the results Optmzaton results wth checkerboard pattern are also dffcult to manufacture and are meanngless from an engneerng standpont Mesh dependency means that the optmzaton results are correlated to the mesh densty A dfferent mesh densty may lead to a dfferent topology dstrbuton Normally, the fner mesh densty, then the greater the number of small structures ncluded n the results, as shown n Fg 3c e Results wth mesh dependency are also meanngless for engneerng applcatons 42 Current methods to elmnate numercal nstablty and ther shortcomngs Some of the reasons for numercal nstablty n topology optmzaton have been analyzed [20, 21] A few approaches [16, 17, 19] have been proposed to remove such nstablty Porous densty can be prohbted by ntroducng a penalty factor p nto the SIMP method Checkerboard pattern has some potental relaton to mesh dependency; a method desgned to avod checkerboard can often avod mesh dependency The research of Daz and Sgmund [20] demonstrates that checkerboard pattern has some connecton wth the element type A hgh-order-type element or non-conformng element can decrease or elmnate checkerboard pattern However, as the downsde, the calculaton cost ncreases wth a hgh-order element Some other methods have been used to elmnate checkerboard pattern and mesh dependency These technques nclude the permeter control method, densty slope control, and mesh-ndependent flter, etc These methods can be generally classfed nto two types The frst one s the global control approach, whch ncludes the permeter control method and the global gradent constrant method Ths approach can control and adust the global functon values and ther dervatves, but has lttle local effect Due to ths property, there are some locally superfluous structures n some engneerng problem results, such as very thn bars Furthermore, t s dffcult to determne the upper bound n permeter control and densty
6 792 a Porous materal b Checkerboard pattern c ( elements) d ( elements) e ( elements) Fg 3 a Porous materal b Checkerboard pattern Mesh dependency: c elements; d elements; e elements slope control methods If the permeter bound s too tght, there may be no soluton of the optmzaton problem, whle the method shows no strength f the constrant s too loose In addton, some extra constrant should be added to the optmzaton problem f the global approach s used, whch may lead to some convergence and numercal problems, although some advanced large-scale mathematcal programmng algorthms can solve problems wth extra constrants The second type s the local control approach, whch ncludes the densty slope control method and mesh-ndependent flter, n whch a local threshold s defned to flter the local varable The local approach may lack some global nformaton that needs to be fltered The mesh-ndependent flter s easly appled and has a good convergence wthout extra constrants added nto the optmzaton problem However, ths method s based on heurstc crtera The densty slope control method ntroduces a large number of extra constrants nto optmzaton, and, thus, t s computatonally prohbtve for practcal applcatons 43 A new method to guarantee numercal stablty As mentoned above, some shortcomngs exst n current methods to erase numercal nstablty However, It s very mportant to choose an approprate method to ensure that the results are manufacturable for engneerng purposes In theory, checkerboard pattern and porous materal can be seem as local sngulartes of numercal calculaton To elmnate these sngular numercal problems, a specal convoluton ntegral factor H (2r 1) (2r 1), n whch r, named the convoluton radus, s ntroduced to convolute the desgn varable and smooth or even the local sngularty n the desgn doman Ths process can elmnate checkerboard pattern and mesh dependency effectvely wth an approprate convoluton radus n matrx H The convoluton ntegral factor can be defned as follows: H ð2r 1Þð2r 1Þ ¼ 3 1 r 1 r 1 r r 1 r r 1 1 4r 3 r 1 r 1 r 1 r (10a) 2 3 r ¼ 2 yelds H 33 ¼ (10b) r ¼ 3 yelds H 55 ¼ (10c)
7 793 The desgn varable x k e (m, n) can be treated wth the convoluton ntegral factor as: Defne desgn doman, non-desgn doman, load and constrant Mesh the structure bx e k ð; Þ ¼ X þr 1 þr 1 X m¼ rþ1 n¼ rþ1 x e k ðm; nþ Hðm þ r; n þ rþ (11) Intalze element desgn varables Calculate stff matrx and dsplacement Calculate obectve and senstvty The dervatve ð; Þ can also be treated e wth the convoluton ntegral factor as: e ð; Þ ¼ X þr 1 X m¼ rþ1 e ðm; nþ Hðm þ r; n þ rþ (12) No extra constrants are ntroduced nto the optmzaton wth convoluton ntegral factor method Therefore, ths convoluton ntegral factor method s convenent for calculatng and programmng Ths method can be used for a varety of problems 5 Soluton procedure of topology optmzaton The soluton procedure of the topology optmzaton based on the SIMP method, hybrd soluton algorthm, and the convoluton ntegral factor method can be wrtten as follows: 1 Defne desgn doman, non-desgn doman, constrants, loads, etc Elements relatve denstes n the desgn doman should vary wth teraton procedure, whle those n the non-desgn doman are of fxed value, whch are set as a fxed value of ether 0 or 1 2 Mesh the structure, calculate the elements stffness matrces when the relatve densty of all elements s set as 1 3 Intalze desgn varables, set every element n the desgn doman to a predetermned relatve densty value 4 Calculate materal characterstcs of every element, such as Young s modulus, calculate the elements stffness matrces, assemble the stffness matrces of the entre structure, and calculate the dsplacement of nodes 5 Calculate the obectve and constrant functon value, and calculate ther dervatves Treat the desgn varables wth the convoluton ntegral factor method 6 Update the desgn varables usng the hybrd soluton algorthm 7 Check the convergence of the result Return to step 4 f the result does not converge, whle contnue to step 8 and fnsh the teraton f the result converges The relatve dfference of the desgn varable n Eq 13 or that of the obectve n Eq 14 can be used as the udgement condton of a convergent result: max ðx kþ1 Þ max ðx k Þ max ðx k Þ <" (13) C kþ1 C k <" (14) C k 8 Output the desgn varables, obectve value and constrant value, and fnsh the calculaton A flow chart of the above calculaton s shown n Fg 4 6 Numercal example Update desgn varables wth hybrd soluton algorthm Treat desgn varables wth convoluton ntegral factor method Convergence? Yes Output result The end Fg 4 Flow chart of topology optmzaton The followng example llustrates the effect and functon of the hybrd soluton algorthm and the convoluton ntegral factor method proposed n ths paper A two-dmensonal smply supported beam wth unt thckness s clamped at the rght end and supported at left end subect to a downward unt force P at the rght center, and can be approxmated as a plane-stress problem, as shown n Fg 5a The formula n Eq 14 s used to udge the convergence The beam s meshed nto quadrangle elements wth four nodes The obectve s to mnmze the complance subect to a volume fracton constrant of 05 wth an equlbrum equaton constrant 61 Effect and functon of the hybrd algorthm Results wth the MMA, GCMMA, and hybrd algorthms are compared and lsted n Fg 5b d,e g, and h, No
8 794 respectvely The convoluton ntegral factor method s used to elmnate numercal nstablty and the convoluton radus s set as 2 Fgure 5b d are the results wth the smplex MMA approach on the 100th, 400th, and 800th teratons Some checkerboard and porous elements are exhbted n the results, and the calculaton cannot converge by the end of the 800th teraton Ths means that the smplex MMA approach s not helpful to obtan a manufacturable result Fgure 5e g show the results wth the smplex GCMMA approach on the 100th, 400th, and 800th teratons The calculaton s stable and converges, but very slowly, so that a perfect and clear topology structure cannot be obtaned by the end of 800th teraton a P=1 b c d e f g h Fg 5 a Smply supported beam under sngle load b d Results of MMA on 100th, 400th, and 800th teratons, respectvely e g Results of GCMMA on 100th, 400th, and 800th teratons, respectvely h Results of hybrd algorthm on 100th, 200th, and 400th teratons, respectvely k Curves for complance versus teraton numbers
9 795 Fg 5 (contnued) Fgure 5h are the results wth the hybrd algorthm on the 100th, 200th, and 400th teraton The calculaton s stable and can converge very quckly A perfect and manufacturable topology result can be obtaned by as early as the 400th teraton Fgure 5k shows the convergence curves of obectve functon about the three dfferent algorthms It shows that the obectve functon wth MMA decreases very quckly at the frst several teratons, whle some numercal oscllaton s exhbted on the subsequent teratons, and the soluton s dffcult to converge The calculaton wth GCMMA s stable, but the obectve functon decreases very slowly The optmzaton wth the hybrd algorthm converges very quckly and exhbts no numercal oscllaton 62 Effect and functon of the convoluton radus The convoluton radus s ndependent of mesh sze An approprate convoluton radus s helpful to elmnate checkerboard pattern and porous materal n the result Fgure 6a d shows the effect of dfferent convoluton raduses wth the hybrd soluton algorthm used n calculaton Fgure 6a shows the topology result wth convoluton radus r=1, and many checkerboard patterns can be seen Fgure 6b,c are the topology results wth r=2 and 3, respectvely As can be seen, the checkerboard patterns have dsappeared and a perfect topology s found Fgure 6d shows the topology result wth r=4, and some porous materal can be seen Fgure 6 allows us to conclude that an approprate convoluton radus can beneft n the producton of a manufacturable result, so, a proper range of convoluton radus 2 r 3 s suggested 7 Concluson Topology optmzaton wth the sold sotropc materal wth penalzaton (SIMP) method has strong applcatons n engneerng A hybrd soluton algorthm based on the method of movng asymptotes (MMA) and the globally convergent verson of the method of movng asymptotes (GCMMA) can be successfully used n topology optmzaton, and calculaton effcency and convergence can be effectvely mproved The convoluton ntegral factor method s an effectve way to treat desgn varables An approprate convoluton radus can lead to a manufacturable result The work deeply dscussng the calculaton effcency and convergence of optmalty crtera (OC), MMA, GCMMA, and hybrd algorthms n topology optmzaton s currently beng undertaken a r=1 b r=2 c r=3 d r=4 Fg 6 a r=1 b r=2 c r=3 d r=4
10 796 Acknowledgments The authors are thankful for the fnancal support from the Natonal 973 Proect and the Natonal 863 Proect The assstance on the MMA/GCMMA algorthms from Professor K Svanberg s gratefully acknowledged References 1 Bendsoe MP, Kkuch N (1988) Generatng optmal topologes n structural desgn usng a homogenzaton method Comput Methods Appl Mech Eng 71(2): Hassan B, Hnton E (1999) Homogenzaton and structural topology optmzaton: theory, practce and software Sprnger, London 3 Mlenek HP, Schrrmacher R (1993) An engneerng approach to optmal materal dstrbuton and shape fndng Comput Methods Appl Mech Eng 106: Sgmund O (1994) Desgn of materal structures usng topology optmzaton PhD thess, Department of Sold Mechancs, Techncal Unversty of Denmark, Denmark 5 Bendsoe MP, Sgmund O (1999) Materal nterpolatons n topology optmzaton Arch Appl Mech 69(9 10): Zhou M, Rozvany GIN (1991) The COC algorthm, part II: topologcal, geometry and generalzed shape optmzaton Comput Methods Appl Mech Eng 89: Svanberg K (1987) The method of movng asymptotes a new method for structural optmzaton Int J Numer Methods Eng 24: Svanberg K (1999) The MMA for modelng and solvng optmzaton problems In: Proceedngs of the Thrd World Congress of Structural and Multdscplnary Optmzaton, Buffalo, New York, May Svanberg K (1999) A new globally convergent verson of the method of movng asymptotes Techncal report, Department of Mathematcs, Royal Insttute of Technology, Stockholm, Sweden 10 Svanberg K (1998) Two prmal-dual nteror-pont methods for the MMA subproblems Techncal report TRITA-MAT S12, Department of Mathematcs, Royal Insttute of Technology, Stockholm, Sweden 11 Svanberg K (1995) A globally convergent verson of MMA wthout lnesearch In: Rozvany GIN, Olhoff N (eds) Proceedngs of the Frst World Congress of Structural and Multdscplnary Optmzaton, Goslar, Germany, May/June 1995, pp Bruyneel M, Duysnx P, Fleury C (2002) A famly of MMA approxmatons for structural optmzaton Struct Multdscpl Optm 24(4): Fleury C (1979) Structural weght optmzaton by dual methods of convex programmng Int J Numer Methods Eng 14: Fleury C, Brabant V (1986) Structural optmzaton: a new dual methods usng mxed varables Int J Numer Methods Eng 23(3): Lu WX (1994) Mechancal optmal desgn Tsnghua Unversty Press, Beng, Chna 16 Haber RB, Jog CS, Bendsøe MP (1996) A new approach to varable-topology shape desgn usng a constrant on permeter Struct Optm 11: Petersson J, Sgmund O (1998) Slope constraned topology optmzaton Int J Numer Method Eng 41(8): Sgmund O, Petersson J (1998) Numercal nstabltes n topology optmzaton: a survey on procedures dealng wth checkerboards, mesh dependences and local mnma Struct Optm 16(1): Zhou M, Shyy YK, Thomas HL (2001) Checkerboard and mnmum member sze control n topology optmzaton Struct Multdscpl Optm 21(2): Daz AR, Sgmund O (1995) Checkerboard patterns n layout optmzaton Struct Optm 10(1): Jog CS, Haber RB (1996) Stablty of fnte element models for dstrbuted parameter optmzaton and topology desgn Comput Methods Appl Mech Eng 130(3 4):
Solving two-person zero-sum game by Matlab
Appled Mechancs and Materals Onlne: 2011-02-02 ISSN: 1662-7482, Vols. 50-51, pp 262-265 do:10.4028/www.scentfc.net/amm.50-51.262 2011 Trans Tech Publcatons, Swtzerland Solvng two-person zero-sum game by
More informationPlate/shell topological optimization subjected to linear buckling constraints by adopting composite exponential filtering function
Acta Mech. Sn. 2016) 324):649 658 DOI 10.1007/s10409-015-0531-5 RESEARCH PAPER Plate/shell topologcal optmzaton subjected to lnear bucklng constrants by adoptng composte exponental flterng functon Hong-Lng
More informationNUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS
ARPN Journal of Engneerng and Appled Scences 006-017 Asan Research Publshng Network (ARPN). All rghts reserved. NUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS Igor Grgoryev, Svetlana
More informationDesign of Structure Optimization with APDL
Desgn of Structure Optmzaton wth APDL Yanyun School of Cvl Engneerng and Archtecture, East Chna Jaotong Unversty Nanchang 330013 Chna Abstract In ths paper, the desgn process of structure optmzaton wth
More informationAn Application of the Dulmage-Mendelsohn Decomposition to Sparse Null Space Bases of Full Row Rank Matrices
Internatonal Mathematcal Forum, Vol 7, 2012, no 52, 2549-2554 An Applcaton of the Dulmage-Mendelsohn Decomposton to Sparse Null Space Bases of Full Row Rank Matrces Mostafa Khorramzadeh Department of Mathematcal
More informationCluster Analysis of Electrical Behavior
Journal of Computer and Communcatons, 205, 3, 88-93 Publshed Onlne May 205 n ScRes. http://www.scrp.org/ournal/cc http://dx.do.org/0.4236/cc.205.350 Cluster Analyss of Electrcal Behavor Ln Lu Ln Lu, School
More informationFinite Element Analysis of Rubber Sealing Ring Resilience Behavior Qu Jia 1,a, Chen Geng 1,b and Yang Yuwei 2,c
Advanced Materals Research Onlne: 03-06-3 ISSN: 66-8985, Vol. 705, pp 40-44 do:0.408/www.scentfc.net/amr.705.40 03 Trans Tech Publcatons, Swtzerland Fnte Element Analyss of Rubber Sealng Rng Reslence Behavor
More informationComplexity Control in the Topology Optimization of Continuum Structures
E. L. Cardoso Grupo de Mecânca Aplcada Departamento de Engenhara Mecânca UFRGS Rua Sarmento Lete, 425 95-7 Porto Alegre, RS. Brazl eduardo_lenz@yahoo.com.br J. S. O. Fonseca Grupo de Mecânca Aplcada Departamento
More informationTopology optimization considering the requirements of deep-drawn sheet metals
th World Congress on Structural and Multdscplnary Optmsaton 7 th - th, June 5, Sydney Australa Topology optmzaton consderng the requrements of deep-drawn sheet metals Robert Denemann, Axel Schumacher,
More informationLS-TaSC Version 2.1. Willem Roux Livermore Software Technology Corporation, Livermore, CA, USA. Abstract
12 th Internatonal LS-DYNA Users Conference Optmzaton(1) LS-TaSC Verson 2.1 Wllem Roux Lvermore Software Technology Corporaton, Lvermore, CA, USA Abstract Ths paper gves an overvew of LS-TaSC verson 2.1,
More informationKent State University CS 4/ Design and Analysis of Algorithms. Dept. of Math & Computer Science LECT-16. Dynamic Programming
CS 4/560 Desgn and Analyss of Algorthms Kent State Unversty Dept. of Math & Computer Scence LECT-6 Dynamc Programmng 2 Dynamc Programmng Dynamc Programmng, lke the dvde-and-conquer method, solves problems
More informationIntra-Parametric Analysis of a Fuzzy MOLP
Intra-Parametrc Analyss of a Fuzzy MOLP a MIAO-LING WANG a Department of Industral Engneerng and Management a Mnghsn Insttute of Technology and Hsnchu Tawan, ROC b HSIAO-FAN WANG b Insttute of Industral
More informationTopology Design using LS-TaSC Version 2 and LS-DYNA
Topology Desgn usng LS-TaSC Verson 2 and LS-DYNA Wllem Roux Lvermore Software Technology Corporaton, Lvermore, CA, USA Abstract Ths paper gves an overvew of LS-TaSC verson 2, a topology optmzaton tool
More informationQuality Improvement Algorithm for Tetrahedral Mesh Based on Optimal Delaunay Triangulation
Intellgent Informaton Management, 013, 5, 191-195 Publshed Onlne November 013 (http://www.scrp.org/journal/m) http://dx.do.org/10.36/m.013.5601 Qualty Improvement Algorthm for Tetrahedral Mesh Based on
More informationR s s f. m y s. SPH3UW Unit 7.3 Spherical Concave Mirrors Page 1 of 12. Notes
SPH3UW Unt 7.3 Sphercal Concave Mrrors Page 1 of 1 Notes Physcs Tool box Concave Mrror If the reflectng surface takes place on the nner surface of the sphercal shape so that the centre of the mrror bulges
More informationSolitary and Traveling Wave Solutions to a Model. of Long Range Diffusion Involving Flux with. Stability Analysis
Internatonal Mathematcal Forum, Vol. 6,, no. 7, 8 Soltary and Travelng Wave Solutons to a Model of Long Range ffuson Involvng Flux wth Stablty Analyss Manar A. Al-Qudah Math epartment, Rabgh Faculty of
More informationModule 6: FEM for Plates and Shells Lecture 6: Finite Element Analysis of Shell
Module 6: FEM for Plates and Shells Lecture 6: Fnte Element Analyss of Shell 3 6.6. Introducton A shell s a curved surface, whch by vrtue of ther shape can wthstand both membrane and bendng forces. A shell
More informationAn Optimal Algorithm for Prufer Codes *
J. Software Engneerng & Applcatons, 2009, 2: 111-115 do:10.4236/jsea.2009.22016 Publshed Onlne July 2009 (www.scrp.org/journal/jsea) An Optmal Algorthm for Prufer Codes * Xaodong Wang 1, 2, Le Wang 3,
More informationInvestigations of Topology and Shape of Multi-material Optimum Design of Structures
Advanced Scence and Tecnology Letters Vol.141 (GST 2016), pp.241-245 ttp://dx.do.org/10.14257/astl.2016.141.52 Investgatons of Topology and Sape of Mult-materal Optmum Desgn of Structures Quoc Hoan Doan
More informationAnalysis of Continuous Beams in General
Analyss of Contnuous Beams n General Contnuous beams consdered here are prsmatc, rgdly connected to each beam segment and supported at varous ponts along the beam. onts are selected at ponts of support,
More informationMathematics 256 a course in differential equations for engineering students
Mathematcs 56 a course n dfferental equatons for engneerng students Chapter 5. More effcent methods of numercal soluton Euler s method s qute neffcent. Because the error s essentally proportonal to the
More informationS1 Note. Basis functions.
S1 Note. Bass functons. Contents Types of bass functons...1 The Fourer bass...2 B-splne bass...3 Power and type I error rates wth dfferent numbers of bass functons...4 Table S1. Smulaton results of type
More informationGSLM Operations Research II Fall 13/14
GSLM 58 Operatons Research II Fall /4 6. Separable Programmng Consder a general NLP mn f(x) s.t. g j (x) b j j =. m. Defnton 6.. The NLP s a separable program f ts objectve functon and all constrants are
More informationParallelism for Nested Loops with Non-uniform and Flow Dependences
Parallelsm for Nested Loops wth Non-unform and Flow Dependences Sam-Jn Jeong Dept. of Informaton & Communcaton Engneerng, Cheonan Unversty, 5, Anseo-dong, Cheonan, Chungnam, 330-80, Korea. seong@cheonan.ac.kr
More informationSimulation Based Analysis of FAST TCP using OMNET++
Smulaton Based Analyss of FAST TCP usng OMNET++ Umar ul Hassan 04030038@lums.edu.pk Md Term Report CS678 Topcs n Internet Research Sprng, 2006 Introducton Internet traffc s doublng roughly every 3 months
More informationReview of approximation techniques
CHAPTER 2 Revew of appromaton technques 2. Introducton Optmzaton problems n engneerng desgn are characterzed by the followng assocated features: the objectve functon and constrants are mplct functons evaluated
More informationAn Iterative Solution Approach to Process Plant Layout using Mixed Integer Optimisation
17 th European Symposum on Computer Aded Process Engneerng ESCAPE17 V. Plesu and P.S. Agach (Edtors) 2007 Elsever B.V. All rghts reserved. 1 An Iteratve Soluton Approach to Process Plant Layout usng Mxed
More informationSupport Vector Machines
/9/207 MIST.6060 Busness Intellgence and Data Mnng What are Support Vector Machnes? Support Vector Machnes Support Vector Machnes (SVMs) are supervsed learnng technques that analyze data and recognze patterns.
More informationA mathematical programming approach to the analysis, design and scheduling of offshore oilfields
17 th European Symposum on Computer Aded Process Engneerng ESCAPE17 V. Plesu and P.S. Agach (Edtors) 2007 Elsever B.V. All rghts reserved. 1 A mathematcal programmng approach to the analyss, desgn and
More informationVery simple computational domains can be discretized using boundary-fitted structured meshes (also called grids)
Structured meshes Very smple computatonal domans can be dscretzed usng boundary-ftted structured meshes (also called grds) The grd lnes of a Cartesan mesh are parallel to one another Structured meshes
More informationA Saturation Binary Neural Network for Crossbar Switching Problem
A Saturaton Bnary Neural Network for Crossbar Swtchng Problem Cu Zhang 1, L-Qng Zhao 2, and Rong-Long Wang 2 1 Department of Autocontrol, Laonng Insttute of Scence and Technology, Benx, Chna bxlkyzhangcu@163.com
More informationOptimization Methods: Integer Programming Integer Linear Programming 1. Module 7 Lecture Notes 1. Integer Linear Programming
Optzaton Methods: Integer Prograng Integer Lnear Prograng Module Lecture Notes Integer Lnear Prograng Introducton In all the prevous lectures n lnear prograng dscussed so far, the desgn varables consdered
More informationImprovement of Spatial Resolution Using BlockMatching Based Motion Estimation and Frame. Integration
Improvement of Spatal Resoluton Usng BlockMatchng Based Moton Estmaton and Frame Integraton Danya Suga and Takayuk Hamamoto Graduate School of Engneerng, Tokyo Unversty of Scence, 6-3-1, Nuku, Katsuska-ku,
More informationRange images. Range image registration. Examples of sampling patterns. Range images and range surfaces
Range mages For many structured lght scanners, the range data forms a hghly regular pattern known as a range mage. he samplng pattern s determned by the specfc scanner. Range mage regstraton 1 Examples
More information5 The Primal-Dual Method
5 The Prmal-Dual Method Orgnally desgned as a method for solvng lnear programs, where t reduces weghted optmzaton problems to smpler combnatoral ones, the prmal-dual method (PDM) has receved much attenton
More informationClassification / Regression Support Vector Machines
Classfcaton / Regresson Support Vector Machnes Jeff Howbert Introducton to Machne Learnng Wnter 04 Topcs SVM classfers for lnearly separable classes SVM classfers for non-lnearly separable classes SVM
More informationAPPLICATION OF MULTIVARIATE LOSS FUNCTION FOR ASSESSMENT OF THE QUALITY OF TECHNOLOGICAL PROCESS MANAGEMENT
3. - 5. 5., Brno, Czech Republc, EU APPLICATION OF MULTIVARIATE LOSS FUNCTION FOR ASSESSMENT OF THE QUALITY OF TECHNOLOGICAL PROCESS MANAGEMENT Abstract Josef TOŠENOVSKÝ ) Lenka MONSPORTOVÁ ) Flp TOŠENOVSKÝ
More informationStructural Optimization Using OPTIMIZER Program
SprngerLnk - Book Chapter http://www.sprngerlnk.com/content/m28478j4372qh274/?prnt=true ق.ظ 1 of 2 2009/03/12 11:30 Book Chapter large verson Structural Optmzaton Usng OPTIMIZER Program Book III European
More informationLECTURE NOTES Duality Theory, Sensitivity Analysis, and Parametric Programming
CEE 60 Davd Rosenberg p. LECTURE NOTES Dualty Theory, Senstvty Analyss, and Parametrc Programmng Learnng Objectves. Revew the prmal LP model formulaton 2. Formulate the Dual Problem of an LP problem (TUES)
More informationSmoothing Spline ANOVA for variable screening
Smoothng Splne ANOVA for varable screenng a useful tool for metamodels tranng and mult-objectve optmzaton L. Rcco, E. Rgon, A. Turco Outlne RSM Introducton Possble couplng Test case MOO MOO wth Game Theory
More informationType-2 Fuzzy Non-uniform Rational B-spline Model with Type-2 Fuzzy Data
Malaysan Journal of Mathematcal Scences 11(S) Aprl : 35 46 (2017) Specal Issue: The 2nd Internatonal Conference and Workshop on Mathematcal Analyss (ICWOMA 2016) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES
More informationTopology Optimization - Practical Aspects for Industrial Applications
9 th World Congress on Structural and Multdscplnary Optmzaton June 13-17, 2011, Shzuoka, Japan Topology Optmzaton - Practcal Aspects for Industral Applcatons Mng zhou 1, Raphael Fleury 1, Sam Patten 2,
More informationSHAPE OPTIMIZATION OF STRUCTURES BY MODIFIED HARMONY SEARCH
INTERNATIONAL JOURNAL OF OPTIMIZATION IN CIVIL ENGINEERING Int. J. Optm. Cvl Eng., 2011; 3:485-494 SHAPE OPTIMIZATION OF STRUCTURES BY MODIFIED HARMONY SEARCH S. Gholzadeh *,, A. Barzegar and Ch. Gheyratmand
More informationA New Approach For the Ranking of Fuzzy Sets With Different Heights
New pproach For the ankng of Fuzzy Sets Wth Dfferent Heghts Pushpnder Sngh School of Mathematcs Computer pplcatons Thapar Unversty, Patala-7 00 Inda pushpndersnl@gmalcom STCT ankng of fuzzy sets plays
More informationThe Greedy Method. Outline and Reading. Change Money Problem. Greedy Algorithms. Applications of the Greedy Strategy. The Greedy Method Technique
//00 :0 AM Outlne and Readng The Greedy Method The Greedy Method Technque (secton.) Fractonal Knapsack Problem (secton..) Task Schedulng (secton..) Mnmum Spannng Trees (secton.) Change Money Problem Greedy
More informationClassifier Selection Based on Data Complexity Measures *
Classfer Selecton Based on Data Complexty Measures * Edth Hernández-Reyes, J.A. Carrasco-Ochoa, and J.Fco. Martínez-Trndad Natonal Insttute for Astrophyscs, Optcs and Electroncs, Lus Enrque Erro No.1 Sta.
More informationThe Codesign Challenge
ECE 4530 Codesgn Challenge Fall 2007 Hardware/Software Codesgn The Codesgn Challenge Objectves In the codesgn challenge, your task s to accelerate a gven software reference mplementaton as fast as possble.
More informationHigh-Boost Mesh Filtering for 3-D Shape Enhancement
Hgh-Boost Mesh Flterng for 3-D Shape Enhancement Hrokazu Yagou Λ Alexander Belyaev y Damng We z Λ y z ; ; Shape Modelng Laboratory, Unversty of Azu, Azu-Wakamatsu 965-8580 Japan y Computer Graphcs Group,
More informationProper Choice of Data Used for the Estimation of Datum Transformation Parameters
Proper Choce of Data Used for the Estmaton of Datum Transformaton Parameters Hakan S. KUTOGLU, Turkey Key words: Coordnate systems; transformaton; estmaton, relablty. SUMMARY Advances n technologes and
More informationAn Accurate Evaluation of Integrals in Convex and Non convex Polygonal Domain by Twelve Node Quadrilateral Finite Element Method
Internatonal Journal of Computatonal and Appled Mathematcs. ISSN 89-4966 Volume, Number (07), pp. 33-4 Research Inda Publcatons http://www.rpublcaton.com An Accurate Evaluaton of Integrals n Convex and
More informationLearning the Kernel Parameters in Kernel Minimum Distance Classifier
Learnng the Kernel Parameters n Kernel Mnmum Dstance Classfer Daoqang Zhang 1,, Songcan Chen and Zh-Hua Zhou 1* 1 Natonal Laboratory for Novel Software Technology Nanjng Unversty, Nanjng 193, Chna Department
More informationVirtual Machine Migration based on Trust Measurement of Computer Node
Appled Mechancs and Materals Onlne: 2014-04-04 ISSN: 1662-7482, Vols. 536-537, pp 678-682 do:10.4028/www.scentfc.net/amm.536-537.678 2014 Trans Tech Publcatons, Swtzerland Vrtual Machne Mgraton based on
More informationProblem Definitions and Evaluation Criteria for Computational Expensive Optimization
Problem efntons and Evaluaton Crtera for Computatonal Expensve Optmzaton B. Lu 1, Q. Chen and Q. Zhang 3, J. J. Lang 4, P. N. Suganthan, B. Y. Qu 6 1 epartment of Computng, Glyndwr Unversty, UK Faclty
More informationA Five-Point Subdivision Scheme with Two Parameters and a Four-Point Shape-Preserving Scheme
Mathematcal and Computatonal Applcatons Artcle A Fve-Pont Subdvson Scheme wth Two Parameters and a Four-Pont Shape-Preservng Scheme Jeqng Tan,2, Bo Wang, * and Jun Sh School of Mathematcs, Hefe Unversty
More informationComplex Numbers. Now we also saw that if a and b were both positive then ab = a b. For a second let s forget that restriction and do the following.
Complex Numbers The last topc n ths secton s not really related to most of what we ve done n ths chapter, although t s somewhat related to the radcals secton as we wll see. We also won t need the materal
More informationBoundary Condition Simulation for Structural Local Refined Modeling Using Genetic Algorithm
2016 Internatonal Conference on Artfcal Intellgence: Technques and Applcatons (AITA 2016) ISBN: 978-1-60595-389-2 Boundary Condton Smulaton for Structural Local Refned Modelng Usng Genetc Algorthm Zhong
More informationNAG Fortran Library Chapter Introduction. G10 Smoothing in Statistics
Introducton G10 NAG Fortran Lbrary Chapter Introducton G10 Smoothng n Statstcs Contents 1 Scope of the Chapter... 2 2 Background to the Problems... 2 2.1 Smoothng Methods... 2 2.2 Smoothng Splnes and Regresson
More informationA Fast Visual Tracking Algorithm Based on Circle Pixels Matching
A Fast Vsual Trackng Algorthm Based on Crcle Pxels Matchng Zhqang Hou hou_zhq@sohu.com Chongzhao Han czhan@mal.xjtu.edu.cn Ln Zheng Abstract: A fast vsual trackng algorthm based on crcle pxels matchng
More informationActive Contours/Snakes
Actve Contours/Snakes Erkut Erdem Acknowledgement: The sldes are adapted from the sldes prepared by K. Grauman of Unversty of Texas at Austn Fttng: Edges vs. boundares Edges useful sgnal to ndcate occludng
More informationMultiple optimum values
1.204 Lecture 22 Unconstraned nonlnear optmzaton: Amoeba BFGS Lnear programmng: Glpk Multple optmum values A B C G E Z X F Y D X 1 X 2 Fgure by MIT OpenCourseWare. Heurstcs to deal wth multple optma: Start
More informationSENSITIVITY ANALYSIS WITH UNSTRUCTURED FREE MESH GENERATORS IN 2-D AND 3-D SHAPE OPTIMIZATION.
SENSITIVITY ANALYSIS WITH UNSTRUCTURED FREE MESH GENERATORS IN 2-D AND 3-D SHAPE OPTIMIZATION. P. Duysnx, W.H. Zhang, C. Fleury. Aerospace Laboratory, LTAS, Unversty of Lège B-4000 LIEGE, BELGIUM. ABSTRACT.
More informationSupport Vector Machines
Support Vector Machnes Decson surface s a hyperplane (lne n 2D) n feature space (smlar to the Perceptron) Arguably, the most mportant recent dscovery n machne learnng In a nutshell: map the data to a predetermned
More informationSum of Linear and Fractional Multiobjective Programming Problem under Fuzzy Rules Constraints
Australan Journal of Basc and Appled Scences, 2(4): 1204-1208, 2008 ISSN 1991-8178 Sum of Lnear and Fractonal Multobjectve Programmng Problem under Fuzzy Rules Constrants 1 2 Sanjay Jan and Kalash Lachhwan
More informationLecture 4: Principal components
/3/6 Lecture 4: Prncpal components 3..6 Multvarate lnear regresson MLR s optmal for the estmaton data...but poor for handlng collnear data Covarance matrx s not nvertble (large condton number) Robustness
More informationS.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION?
S.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION? Célne GALLET ENSICA 1 place Emle Bloun 31056 TOULOUSE CEDEX e-mal :cgallet@ensca.fr Jean Luc LACOME DYNALIS Immeuble AEROPOLE - Bat 1 5, Avenue Albert
More informationLECTURE : MANIFOLD LEARNING
LECTURE : MANIFOLD LEARNING Rta Osadchy Some sldes are due to L.Saul, V. C. Raykar, N. Verma Topcs PCA MDS IsoMap LLE EgenMaps Done! Dmensonalty Reducton Data representaton Inputs are real-valued vectors
More informationSupport Vector Machines. CS534 - Machine Learning
Support Vector Machnes CS534 - Machne Learnng Perceptron Revsted: Lnear Separators Bnar classfcaton can be veed as the task of separatng classes n feature space: b > 0 b 0 b < 0 f() sgn( b) Lnear Separators
More informationA Binarization Algorithm specialized on Document Images and Photos
A Bnarzaton Algorthm specalzed on Document mages and Photos Ergna Kavalleratou Dept. of nformaton and Communcaton Systems Engneerng Unversty of the Aegean kavalleratou@aegean.gr Abstract n ths paper, a
More informationWishing you all a Total Quality New Year!
Total Qualty Management and Sx Sgma Post Graduate Program 214-15 Sesson 4 Vnay Kumar Kalakband Assstant Professor Operatons & Systems Area 1 Wshng you all a Total Qualty New Year! Hope you acheve Sx sgma
More informationc 2014 Xiaojia(Shelly) Zhang
c 2014 Xaoja(Shelly) Zhang MACRO-ELEMENT APPROACH FOR TOPOLOGY OPTIMIZATION OF TRUSSES USING A GROUND STRUCTURE METHOD BY XIAOJIA(SHELLY) ZHANG THESIS Submtted n partal fulfllment of the requrements for
More information陳申岳 S-Y. Chen, 2007, Gradient-Based Structural and CFD Global Shape Optimization with SmartDO and the Response Smoothing Technology, Proceeding of
陳申岳 S-Y. Chen, 2007, Gradent-Based Structural and CFD Global Shape Optmzaton wth SmartDO and the Response Smoothng Technology, Proceedng of the 7 th World Congress of Structural and Multdscplnary Optmzaton
More informationX- Chart Using ANOM Approach
ISSN 1684-8403 Journal of Statstcs Volume 17, 010, pp. 3-3 Abstract X- Chart Usng ANOM Approach Gullapall Chakravarth 1 and Chaluvad Venkateswara Rao Control lmts for ndvdual measurements (X) chart are
More informationSome Advanced SPC Tools 1. Cumulative Sum Control (Cusum) Chart For the data shown in Table 9-1, the x chart can be generated.
Some Advanced SP Tools 1. umulatve Sum ontrol (usum) hart For the data shown n Table 9-1, the x chart can be generated. However, the shft taken place at sample #21 s not apparent. 92 For ths set samples,
More informationA Weight Balanced Multi-Objective Topology Optimization for Automotive Development
A Weght Balanced Mult-Objectve Topology Optmzaton for Automotve Development Nkola Aulg 1, Emly Nutwell 2, Stefan Menzel 1, Duane Detwler 3 1 Honda Research Insttute Europe GmbH, Offenbach/Man, Germany
More information3D vector computer graphics
3D vector computer graphcs Paolo Varagnolo: freelance engneer Padova Aprl 2016 Prvate Practce ----------------------------------- 1. Introducton Vector 3D model representaton n computer graphcs requres
More informationPositive Semi-definite Programming Localization in Wireless Sensor Networks
Postve Sem-defnte Programmng Localzaton n Wreless Sensor etworks Shengdong Xe 1,, Jn Wang, Aqun Hu 1, Yunl Gu, Jang Xu, 1 School of Informaton Scence and Engneerng, Southeast Unversty, 10096, anjng Computer
More informationDetermining the Optimal Bandwidth Based on Multi-criterion Fusion
Proceedngs of 01 4th Internatonal Conference on Machne Learnng and Computng IPCSIT vol. 5 (01) (01) IACSIT Press, Sngapore Determnng the Optmal Bandwdth Based on Mult-crteron Fuson Ha-L Lang 1+, Xan-Mn
More informationEngineering Structures
Engneerng Structures 37 (2012) 106 124 Contents lsts avalable at ScVerse ScenceDrect Engneerng Structures journal homepage: www.elsever.com/locate/engstruct Topology optmzaton for braced frames: Combnng
More informationREFRACTIVE INDEX SELECTION FOR POWDER MIXTURES
REFRACTIVE INDEX SELECTION FOR POWDER MIXTURES Laser dffracton s one of the most wdely used methods for partcle sze analyss of mcron and submcron sze powders and dspersons. It s quck and easy and provdes
More informationOptimum Synthesis of Mechanisms For Path Generation Using a New Curvature Based Deflection Based Objective Function
Proceedngs of the 6th WSEAS Internatonal Conference on Smulaton, Modellng and Optmzaton, Lsbon, Portugal, September -4, 6 67 Optmum Synthess of Mechansms For Path Generaton Usng a ew Curvature Based Deflecton
More informationCourse Introduction. Algorithm 8/31/2017. COSC 320 Advanced Data Structures and Algorithms. COSC 320 Advanced Data Structures and Algorithms
Course Introducton Course Topcs Exams, abs, Proects A quc loo at a few algorthms 1 Advanced Data Structures and Algorthms Descrpton: We are gong to dscuss algorthm complexty analyss, algorthm desgn technques
More informationComputational layout design optimization of frame structures
Interfaces: archtecture.engneerng.scence 5-8th September, 017, Hamburg, Germany Annette Bögle, Manfred Grohmann (eds.) Computatonal layout desgn optmzaton of frame structures Jun Ye*, Paul Shepherd a,
More informationBiostatistics 615/815
The E-M Algorthm Bostatstcs 615/815 Lecture 17 Last Lecture: The Smplex Method General method for optmzaton Makes few assumptons about functon Crawls towards mnmum Some recommendatons Multple startng ponts
More informationBioTechnology. An Indian Journal FULL PAPER. Trade Science Inc.
[Type text] [Type text] [Type text] ISSN : 0974-74 Volume 0 Issue BoTechnology 04 An Indan Journal FULL PAPER BTAIJ 0() 04 [684-689] Revew on Chna s sports ndustry fnancng market based on market -orented
More informationA MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS
Proceedngs of the Wnter Smulaton Conference M E Kuhl, N M Steger, F B Armstrong, and J A Jones, eds A MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS Mark W Brantley Chun-Hung
More informationOptimal Design of Trusses With Geometric Imperfections: Accounting for Global Instability
Cleveland State Unversty EngagedScholarshp@CSU Cvl and Envronmental Engneerng Faculty Publcatons Cvl and Envronmental Engneerng -5- Optmal Desgn of Trusses Wth Geometrc Imperfectons: Accountng for Global
More informationMeta-heuristics for Multidimensional Knapsack Problems
2012 4th Internatonal Conference on Computer Research and Development IPCSIT vol.39 (2012) (2012) IACSIT Press, Sngapore Meta-heurstcs for Multdmensonal Knapsack Problems Zhbao Man + Computer Scence Department,
More informationChapter 6 Programmng the fnte element method Inow turn to the man subject of ths book: The mplementaton of the fnte element algorthm n computer programs. In order to make my dscusson as straghtforward
More informationBOTTOM STRUCTURE FOR DutchEVO CAR: FORMULATION OF THE PROBLEM AND THE ADJUSTMENT OF THE OPTIMIZATION SYSTEM
INTERNATIONAL DESIGN CONFERENCE - DESIGN 2002 Dubrovn, May 14-17, 2002. BOTTOM STRUCTURE FOR DutchEVO CAR: FORMULATION OF THE PROBLEM AND THE ADJUSTMENT OF THE OPTIMIZATION SYSTEM Natala S. Ermolaeva,
More informationModeling, Manipulating, and Visualizing Continuous Volumetric Data: A Novel Spline-based Approach
Modelng, Manpulatng, and Vsualzng Contnuous Volumetrc Data: A Novel Splne-based Approach Jng Hua Center for Vsual Computng, Department of Computer Scence SUNY at Stony Brook Talk Outlne Introducton and
More informationRouting on Switch Matrix Multi-FPGA Systems
Routng on Swtch Matrx Mult-FPGA Systems Abdel Enou and N. Ranganathan Center for Mcroelectroncs Research Department of Computer Scence and Engneerng Unversty of South Florda Tampa, FL 33620 Abstract In
More informationThe Shortest Path of Touring Lines given in the Plane
Send Orders for Reprnts to reprnts@benthamscence.ae 262 The Open Cybernetcs & Systemcs Journal, 2015, 9, 262-267 The Shortest Path of Tourng Lnes gven n the Plane Open Access Ljuan Wang 1,2, Dandan He
More informationCable optimization of a long span cable stayed bridge in La Coruña (Spain)
Computer Aded Optmum Desgn n Engneerng XI 107 Cable optmzaton of a long span cable stayed brdge n La Coruña (Span) A. Baldomr & S. Hernández School of Cvl Engneerng, Unversty of Coruña, La Coruña, Span
More informationContent Based Image Retrieval Using 2-D Discrete Wavelet with Texture Feature with Different Classifiers
IOSR Journal of Electroncs and Communcaton Engneerng (IOSR-JECE) e-issn: 78-834,p- ISSN: 78-8735.Volume 9, Issue, Ver. IV (Mar - Apr. 04), PP 0-07 Content Based Image Retreval Usng -D Dscrete Wavelet wth
More informationPrivate Information Retrieval (PIR)
2 Levente Buttyán Problem formulaton Alce wants to obtan nformaton from a database, but she does not want the database to learn whch nformaton she wanted e.g., Alce s an nvestor queryng a stock-market
More informationAn Application of Network Simplex Method for Minimum Cost Flow Problems
BALKANJM 0 (0) -0 Contents lsts avalable at BALKANJM BALKAN JOURNAL OF MATHEMATICS journal homepage: www.balkanjm.com An Applcaton of Network Smplex Method for Mnmum Cost Flow Problems Ergun EROGLU *a
More informationTsinghua University at TAC 2009: Summarizing Multi-documents by Information Distance
Tsnghua Unversty at TAC 2009: Summarzng Mult-documents by Informaton Dstance Chong Long, Mnle Huang, Xaoyan Zhu State Key Laboratory of Intellgent Technology and Systems, Tsnghua Natonal Laboratory for
More informationAVO Modeling of Monochromatic Spherical Waves: Comparison to Band-Limited Waves
AVO Modelng of Monochromatc Sphercal Waves: Comparson to Band-Lmted Waves Charles Ursenbach* Unversty of Calgary, Calgary, AB, Canada ursenbach@crewes.org and Arnm Haase Unversty of Calgary, Calgary, AB,
More informationFitting: Deformable contours April 26 th, 2018
4/6/08 Fttng: Deformable contours Aprl 6 th, 08 Yong Jae Lee UC Davs Recap so far: Groupng and Fttng Goal: move from array of pxel values (or flter outputs) to a collecton of regons, objects, and shapes.
More informationA COMPARISON OF TWO METHODS FOR FITTING HIGH DIMENSIONAL RESPONSE SURFACES
Mam, Florda, U.S.A., Aprl 6-8, 7 A COMPARISON OF TWO METHODS FOR FITTING HIGH DIMENSIONAL RESPONSE SURFACES Marcelo J. Colaço Department of Mechancal and Materals Eng. Mltary Insttute of Engneerng Ro de
More information