Meshless Elasticity Model and Contact Mechanics-based Verification Technique
|
|
- Lenard Wilson
- 6 years ago
- Views:
Transcription
1 Meshless Elastcty Model and Contact Mechancs-based Verfcaton Technque Rfat Aras 1, Yuzhong Shen 1, Mchel Audette 1, Stephane Bordas 2 1 Department of Modelng, Smulaton, and Vsualzaton Engneerng, Old Domnon Unversty, Norfolk, VA, Unted States 2 Insttute of Mechancs and Advanced Materals, Cardff Unversty, Wales, UK Abstract Mesh-based technques are well studed and establshed methods for solvng contnuum bomechancs problems. When the problem at hand nvolves extreme deformatons or artfcal dscontnutes, meshless methods provde several advantages over the mesh-based methods. Ths work dscusses the Movng Least Square approxmaton-based meshless collocaton method for smulatng deformable objects and presents a verfcaton technque that s based on the Hertzan theory of non-adhesve elastc contact. The effectveness of the Hertzan contact theory as a means for verfcaton was frst tested and proven through a well-establshed FEM code, FEBo. The meshless method was mplemented as a reusable component for the SOFA framework, an open source software lbrary for real-tme smulatons. Through expermentaton, the Hertzan theory has been tested aganst SOFA hexahedral FEM and the meshless models wthn the SOFA framework. Convergence studes and L2 error curves are provded for both models. Expermental results demonstrated the effectveness of the mplementaton of the meshless method. 1 Introducton Soft tssue models have a wde range of applcaton areas wth partcular focus on real-tme medcal smulatons. As a result, accurate nteractve modelng of soft tssue s an mportant and well-establshed research feld n contnuum bomechancs studes. A contnuum model typcally reles on an underlyng mesh structure ether n 2D or 3D dependng on the nature and the requrements of the problem. A breadth-frst classfcaton of mesh-based contnuum models ncludes mass-sprng networks [1], Fnte Element Methods [2], Fnte Volume Methods [3], and Fnte Dfference Methods [4]. Orgnally developed for molecular dynamcs problems and astrophyscal smulatons, meshless (mesh-free) methods offer an appealng alternatve to meshbased methods when the problem nvolves large deformatons and mposed ds-
2 2 contnutes such as cracks and cuts. In ths work, the meshless method s mplemented as a component under the open-source SOFA lbrary, whch focuses on real-tme nteractve smulatons wth an emphass on medcal smulatons. In conjuncton wth the nteractve smulaton requrement, we assume the lnear property of the soft tssue model wthn a specfc range of the stran-stress curve. Ths assumpton s parallel wth the small stran and lnear materal assumptons of the Hertzan contact theory. There have been numerous approaches to modelng mechancs of soft tssues. Among those s the early work of Sederberg and Parry [5] that presented a technque for deformng sold geometrc models n an ntutve free-form manner. The deformatons were based on nterpolatng trvarate Bernsten polynomals, and could be appled ether globally or locally wth volume preservaton. Free-form deformaton s an approxmate and smple method for deformng sold objects; however the lack of physcal bass s grounds for excludng t as an opton for realstc medcal smulatons. The work of Frsken-Gbson [6] modfed the tradtonal voxel based representaton of volumetrc objects and presented a lnked volume representaton that was capable of handlng nteractve object manpulatons such as carvng, cuttng, tearng, and jonng, but stll unable to produce physcally realstc results. An alternatve to volumetrc methods s to use mass-sprng models [7] and membrane based approxmatons that utlze sprng meshes [8]. A sprng mesh s composed of vertces and edges, n whch each edge s realzed as a sprng that connects vertces par-wse, and each vertex s dealzed as a pont mass. Although sprng meshes employ physcal equatons lke Hooke s law, t s dffcult to reproduce specfc elastc materal propertes even wth very careful dstrbuton of sprng stffness through the mesh and suffer from poor numercal stablty. The early work of Bro-Nelsen dscussed a fast adaptaton of fnte element modelng to satsfy speed and robustness requrements n a surgcal smulaton settng [9]. In ths work, the body part was modeled as a 3D lnear elastc sold that conssted of partcles, whch were deformed nto a new shape when forces were appled to the elastc sold. The author ncorporated a technque called condensaton n order to acheve nteractve smulaton. In the fnte element modelng context, condensaton translates nto obtanng a more compact verson of the system model by rearrangng the terms of the matrx equatons. In ths work, the author condensed the equatons by only consderng the fnte element nodes on the surface of the model. A number of recent technques have addressed the fdelty versus effcency trade-off. Another mportant work n the area s the fnte element model based on Total Lagrangan Explct Dynamcs (TLED) by Mller et al. [10]. The dfference between the TLED based fnte element model and other approaches s usng the orgnal reference confguraton of the object to calculate the stress and stran tensors durng a smulaton step. As a result of expressng computatons n the reference coordnates, the authors were able to pre-compute spatal dervatves. The pre-computaton of the spatal dervatves leads to effcency n terms of computatonal resources, whle beng capable of handlng geometrc and materal non-
3 lneartes. In ther work, the authors employed the central dfference based explct ntegraton rather than the mplct ntegraton scheme. Wth ths choce, they were able to avod solvng the set of non-lnear algebrac equatons that are requred by the mplct ntegraton at each tme step. However, the use of explct ntegraton brngs lmtaton on the tme step sze. The authors justfed ther mplementaton choce by statng that the relatvely lower stffness (Young s modulus) value of the soft tssue relaxes the tme step lmtaton consderably compared to the typcal smulatons nvolvng more stff materal lke steel or concrete. Another attempt to ncrease the computatonal effcency of the elastc model n the context of nteractve smulaton was dscussed n the work of Marchesseau et al. [11]. The authors presented a new dscretzaton method called Multplcatve Jacoban Energy Decomposton (MJED), whch allows the smulaton to assemble the stffness matrx of the system faster than the tradtonal Galerkn FEM formulaton. The authors reported up to fve tmes faster computatons for the St. Venant Krchoff materals. For valdaton purposes, the authors compared the MJED approach to the tradtonal FEM mplementaton n the SOFA framework [12], referred to as the standard FEM mplementaton. FEM technques have domnated the feld of computatonal mechancs n the past several decades. They have been wdely used for modelng physcal phenomena such as elastcty, heat transfer, and electromagnetsm and they heavly rely on the assumpton of a contnuous doman. When the problem doman no longer comples wth ths contnuum assumpton, the ratonale behnd usng an FEM based soluton dsappears. FEM s also not well suted to problems nvolvng extreme mesh dstortons that result n degenerate element shapes, movng dscontnutes that do not algn wth the element edges such as propagatng cracks, and advanced materal transformatons such as meltng of a sold or freezng. To address these ssues, sgnfcant nterest has been developed towards a dfferent class of methods for solvng PDEs, namely meshless or mesh-free methods [13, 14]. The very frst meshless method dated back to 1977 [15] and proposed a smoothed partcle hydrodynamcs (SPH) method that was used to model theoretcal astrophyscal phenomena such as galaxy formaton, star formaton, stellar collsons, and dust clouds. Its meshless Lagrangan nature allowed dverse usage areas besdes astrophyscs such as flud flow, ballstcs, volcanology, and oceanography [16]. Although the SPH method elmnates the necessty of a mesh structure and allows the soluton of unbounded problems, t also has ts lmtatons. Because of ts smplstc kernel based approxmaton scheme, t fals to reproduce even frst order polynomals, resultng n severe consstency problems [13]. To allevate ths problem, methods that utlze movng least squares (MLS) approxmatons have been developed. The frst work that used MLS approxmatons n a Galerkn method s the work of Nayroles et al. [17], whch was refned by Belytschko et al. [18] and named Element-Free Galerkn (EFG) method. Ths class of methods, dfferent from the SPH method, use shape functons n approxmatons that are essentally corrected versons of compact supported weght functons. The shape functons 3
4 4 are obtaned by frst representng the approxmaton as a product of a polynomal bass and a vector of unknown coeffcents. Then, a functonal s created by takng the weghted sum of square of the approxmaton error. By takng the dervatve of ths functonal wth respect to the unknown coeffcents and settng t to zero for mnmzng the approxmaton error, we obtan a set of equatons that are reorganzed to solve for the MLS shape functons. The order of consstency of the MLS approxmaton scheme depends on the order and completeness of the used bass functon. If the bass functon used n the approxmaton s a complete polynomal of order k, then the MLS approxmaton s sad to be k th order consstent. Ths fact makes the MLS based approxmatons more consstent than the SPH method. Another technque that has used the MLS approxmaton s the work of Mueller et al. [19] and forms the bass of the meshless method dscussed n ths paper. In ths work, the authors calculated the spatal dervatves of the deformaton gradent only at the partcle locatons smlar to the meshless collocaton methods. Ths technque s capable of smulatng a wde range of materal propertes from very stff materals to soft ones and also handles plastc deformatons as well. Horton et al. [20] proposed a new knd of meshless method named meshless total Lagrangan explct dynamcs method. In ths work, the authors extended ther prevous TLED algorthm [10] to the meshless dscretzaton methodology by precomputng the stran-dsplacement matrces. Ther method s a fully explct method, meanng not requrng mplct tme ntegraton and costly soluton of large system of equatons. Dfferent than our nodal ntegraton approach, ther proposed algorthm ntegrates over a regular background grd wth sngle ntegraton pont per cell. In ths work, we present the detals of the MLS approxmaton-based meshless collocaton method for soft-tssue deformaton. Specfc contrbutons are SOFA mplementaton of the presented meshless method, a new verfcaton technque wth the Hertzan non-adhesve contact theory, and verfcaton of the algorthm wth well-establshed FEM code and our SOFA component that mplements the presented meshless method. The rest of the paper s organzed as follows. Secton 2 descrbes the detals of the meshless elastcty algorthm, ncludng node dstrbuton, constructon of support domans, used weght functon, and calculaton of nodal smulaton values such as mass, volume, and densty. Operatons nvolved n the MLS approxmaton of the deformaton gradent are dscussed, followed by the nternal elastc force calculaton procedure. Secton 3 descrbes the mplementaton detals of the presented meshless method n SOFA and presents the verfcaton technques for the elastcty model along wth the vrtual experment setup and the results. Secton 4 concludes the paper and dscusses possble future extensons to the presented work.
5 5 2 Algorthm 2.1 Dscretzaton of the Contnuum Node dstrbuton s the frst step n the presented meshless collocaton algorthm, whch supports both regular and herarchcal dstrbuton of the nodes through the smulaton doman. In the case of a smulaton doman wth a regular geometrc shape, regular dstrbuton of the nodes s the natural choce. On the other hand, f the smulaton doman has a complex geometry, the regular dstrbuton smply becomes napplcable. In ths case, we sample the volumetrc smulaton doman bounded by the complex boundary surface by herarchcally samplng the volume. In the work of Pauly et al.[21], the authors used a balanced octree data structure to dstrbute the nodes nsde the volume. In ths work, we frst tetrahedralze the smulaton doman wth well-establshed computatonal geometry lbrares lke TetGen [22] and CGAL [23] and then use the set of vertces of the tetrahedra as the meshless node locatons. In ths way, smlar to graded fnte element meshng technques, we can have hgher node densty close to the doman boundary and fewer nodes towards the nteror of the volume where the materal s contnuous. Meshless methods represent a deformable body by a cloud of partcles wth overlappng support domans. Quanttes such as mass, volume, support sze, stran, and stress are stored and updated per partcle for the duraton of the smulaton. In ths work, the support domans of the partcles are sphercal and ther rad are computed by fndng the average dstance of the central node to ts k-nearest neghbors. For effcent neghborhood search purposes, a k-d tree data structure s used. Weght (kernel) functon n the meshless method context s an element that descrbes the way meshless nodes affect each other and how the materal values of the contnuum such as mass, volume, and densty are dstrbuted among the nodes. The neghborng partcles that fall nsde the support doman of a central partcle are weghted usng the polynomal kernel functon wth r j rj 3 rj rj, rj 1 wr ( j) 0, rj 1 j x x, where x and x j are the current locatons of the neghbor- h ng and central partcles, respectvely and h s the support radus of the central partcle. The mass and densty of a meshless node are assgned at the begnnng and kept fxed throughout the smulaton. The mass values are ntalzed wth (1)
6 6 m sr, (2) where ρ s the materal densty value, r s the average dstance of the th node to ts k-nearest neghbors, and s s a scalng factor that s chosen so that the average of the assgned denstes s close to the actual materal densty. The assgned mass value of a meshless node s spread around the node wth the kernel functon. Therefore the densty of a meshless node s calculated after the mass allocaton step by takng the weghted average of the masses of the neghborng nodes m w( r ) (3) 3 j j j In ths work, spatal ntegraton s performed through the nodal ntegraton technque. Compared to other spatal ntegraton technques that utlze a background mesh or grd wth multple ntegraton ponts per regon, nodal ntegraton s fast and effcent wth the added dsadvantage of decreased stablty. We calculate the spatal dervatves of the deformaton gradent only at the partcle locatons smlar to the meshless collocaton methods. 2.2 Movng Least Square (MLS) Approxmaton wth Taylor Seres In contnuum sold mechancs, the elastc stresses nsde a deformable body are computed based on the dsplacement vector feld, whch s typcally defned by a deformable model wth two confguratons named the reference and the current confguratons, respectvely. The coordnates of the partcles n the reference and current confguratons are represented by the materal X ( X, Y, Z) and spatal x ( xy,,z) coordnates respectvely. The dsplacement vector of a partcle s u u u u x X. therefore defned as T X Y Z The Green-St.Venant stran tensor s used to measure the stran T T u u u u, (4) where the gradent of the contnuous dsplacement vector feld u s essentally the dervatves of ( ux, uy, u z ) wth respect to ( X, Y, Z ) arranged n the Jacoban format. For an sotropc lnear-elastc materal, the stran s mapped to the stress by the C tensor that approxmates the materal propertes and s composed of two ndependent coeffcents, Young s Modulus and Posson s Rato C. (5) We need the partal dervatves of the dsplacement vector feld n order to compute the stran, stress, and the nternal elastc forces appled to the meshless partcles. Movng Least Squares (MLS) approxmaton s used to compute ths gradent ( u ).
7 7 For a central meshless partcle and ts neghbor j, the value of u at the locaton of j can be approxmated by the frst order Taylor expanson as u ux u.( X X ) X x j x j x. (6) The weghted sum of squared dfferences between the dsplacement vector and ts approxmaton obtaned from the equaton (6) gves the error measure of the MLS approxmaton 2 e ( ux u ) ( ) j x w r j j. (7) j The expanded equaton of the error measure for a partcle s therefore obtaned by ux ux ux 2 e ( ux ( )X ( ) ( ) ) ( ) j Yj Zj ux w r j j. (8) X Y Z j We want to mnmze ths error measure for some values of ( u x )/ X, ( u x )/ Y, and ( u x )/ Z, therefore we set the dervatve of the error measure e wth respect to the partal dervatves of the dsplacement vector to zero, resultng n three equatons for three unknowns T ux ( X j X)( X j X) w rj ux u ( ) j x j w rj j X X. (9) X j The coeffcent of the partal dervatve on the left-hand sde of the equaton (9) s the 3x3 matrx called the moment matrx (A). A can be nverted and premultpled wth both of the sdes of the equaton for computng the partal dervatves. 2.3 Force Calculaton The elastc body forces that are appled to the ndvdual partcles n the meshless collocaton method are calculated through the stran energy densty, whch s a functon of the partcle dsplacements. For a partcle wth volume v, stran, and stress, the stran energy densty becomes 1 U v ( ). (10) 2 The elastc force per unt volume at a meshless node s locaton s the negatve drectonal dervatve of the above stran energy densty wth respect to ths node s dsplacement. The forces appled to the partcle and ts neghbors j are then
8 8 f U v u u x f U v xj x j u u xj. (11) Eq. (11) s terated over all partcles and the total force appled to a partcle s the sum of all the forces wth the same ndex as that partcle. These force components are obtaned by usng the Green-Sant-Venant stran tensor, whch measures the lnear and shear elongaton. In case of a volume nvertng dsplacement though, ths stran becomes zero. In order to ntroduce restorng body forces n cases of volume nverson, Muller et al. [19] added another energy term to the system that penalzes devatons from a volume conservng transformaton 1 2 U kv( J 1). (12) 2 In ths energy term, J s the Jacoban of the dsplacement vector feld mappng and v k s the volume restoraton constant. Although n our experments we have not worked on cases that cause the nverson of the volume, the effect of ths term and the volume restoraton constant parameter n cases of large deformatons s yet to be examned. For each of the meshless nodes, these force components are accumulated n the force vectors and then passed to the SOFA tme ntegraton module. The modular mplementaton of the SOFA lbrary allows the use of dfferent tme ntegraton schemes wthout changng the actual mplementaton of the algorthm. In ths work, we used the 4 th order Runge-Kutta ntegraton scheme. 3 Implementaton and Verfcaton Technques In modelng and smulaton studes, verfcaton and valdaton of the model s a crucal step. For nteractve smulatons wth a focus on tranng n partcular, t s mportant to valdate the behavor of the deformable body n order to prevent false learnng outcomes. Specfcally, realstc contact representaton and force feedback are sgnfcant features for ensurng the valdty of physcally based smulators. 3.1 Meshless Method Implementaton n SOFA SOFA [12] s an open-source object-orented software lbrary that s targeted towards nteractve medcal smulatons. SOFA has a modular structure that allows users to quckly prototype smulaton scenes wth ready-to-use components. Its object-orented, modular archtecture makes t easy for developers to extend the
9 9 functonaltes of the lbrary by dervng new components from the exstng ones. The meshless elastc model descrbed n ths work has been mplemented as a force feld component, whch s easly nterchangeable wth the exstng force feld components such as hexahedral fnte elements or mass-sprng networks. The requred component nterfaces such as ntalzaton and force accumulaton were mplemented accordng to the algorthm steps descrbed n the prevous secton. Due to the nature of the presented verfcaton method, accurate contact handlng s essental to obtan precse results. Unfortunately, because of the approxmatng nature of the meshless methods, mposng Drchlet boundary condtons s a challengng task on ts own [24]. In ths work, we have followed the approaches that were orgnally adopted n the SOFA lbrary. In SOFA, constrants are flter lke components, whch cancel out the forces and dsplacements appled to ther assocated partcles. For example for fxed node boundary condtons, SOFA's fxed constrant component s used to attach a meshless partcle to a fxed pont n the current confguraton. SOFA has support for several contact handlng methods such as penalty-based and constrant-based methods. Among these, the constrant-based methods are more appealng than the former class of methods because they use Lagrange multplers to handle complex constrants and produce physcally accurate results wth the addtonal computaton cost. Lagrange multplers wth unlateral nteracton laws are used to handle complex constrants. The constrants depend on the relatve postons of the nteractng objects, whch are the meshless partcles and the sphercal rgd ndenter n our case [12]. 3.2 Contact Mechancs Theory Hertzan theory of non-adhesve elastc contact [25] defnes analytcal solutons for the nteracton of elastc half-spaces wth smple shapes n terms of appled force and object ndentaton. For example, the amount of ndentaton of an elastc half-space under a sphercal load s gven by 4 * 3 2 F E Rd, (13) 3 where F s the vertcal force appled on the sphercal load, R s the radus of the * sphercal load, d s the ndentaton amount, and E s the combned Young s E, E ) modulus of the two materals and calculated usng the Young s modul ( 1 2 and Posson s ratos ( v 1, v 2 ) of the two materals as 1 1 v 1 v. (14) * E E1 E2
10 10 The Hertzan theory assumes 1) small strans wthn the elastc materal, 2) much smaller area of contact compared to the areas of the objects n contact, and 3) contnuous and frctonless contact surfaces. The Hertzan theory s an mportant steppng stone n the feld of contact mechancs; therefore there have been numerous fnte element analyss studes about the subject that use both research and commercal fnte element code [26-29]. 3.3 Verfcaton Experments In order to verfy the usablty of the Hertzan contact theory as a means of verfcaton of soft-tssue deformaton, we frst conducted experments usng wellestablshed fnte element code. FEBo s an open-source software sute that s prmarly targeted towards bomechancs and bophyscs problems wth a specfc focus on nonlnear large deformaton problems n bosold mechancs [30]. FEBo provdes several models and optons to represent the non-adhesve Hertzan contact theory. In our experments, we selected the facet-to-facet sldng algorthm that s based on Laursen s contact formulaton [31]. In ths algorthm, the contact constrants are enforced through Lagrange multplers. The FEBo experment was setup by defnng the fxed-poston boundary condtons of the deformable block at ts bottom and sde faces, facet-to-facet sldng contact between the top face of the deformable block and the rgd sphercal ndenter, and the sphere s ndentaton amount. The smulaton was run for 10 tme steps of 0.1s each and the smulaton runtme took over 4 mnutes. The node at the mddle of the top of the deformable block was tracked for the vertcal dsplacement and the vertcal component of the contact force. We compared the obtaned load-dsplacement curve to the theoretcal soluton, whch were n very good agreement (Fgure 1), therefore verfyng the usefulness of the Hertzan contact theory as a verfcaton method. Fgure 1 Comparson of the FEBo FEM Code and the theoretcal soluton of the Hertzan non-adhesve frctonless contact theory.
11 11 The contact mechancs experment was setup as a SOFA scene. In order to assess the ground-truth performance of the contact handlng n SOFA, the rectangular deformable block was represented wth the hexahedral fnte element model n addton to the descrbed meshless method. The valdaton of the hexahedral FEM mplementaton of SOFA s studed by Marchal et al. [32]. For a gven sphere radus, smulatons were performed for varyng force values (F) appled to the sphercal load (Fgure 2). Wth the appled force, the rgd sphere comes nto contact wth the block and deforms t. The vertcal velocty of the sphere s montored and the ndentaton of the materal (d) s measured when the sphere comes to rest. Ths F-d par s compared to the theoretcal soluton. (a) Fgure 2 (a) Intal setup of the ndentaton experment for the SOFA FEM model, (b) the close-up vew of the ndented deformable materal. SOFA allows the user to track and montor smulaton values of ndexed partcles. The meshless nodes are dstrbuted unformly nsde a cubcal volume wth 2m long edge length. The ndentaton experment s repeated for several dstrbuton confguratons, whch play crtcal role especally for the MLS approxmatonbased collocaton methods. The convergence rate n the L2 (vector) error norm of the force-ndentaton pars wth respect to the theoretcal values (Fgure 3) s nvestgated. The effect of dfferent dstrbuton schemes on the accuracy, stablty, and performance of the meshless collocaton methods s yet to be examned. In our mplementaton, the number of neghborng nodes s lmted to 16 for each of the meshless nodes. (b)
12 12 Fgure 3 Error n the L2 norm wth respect to the theoretcal soluton as functons of total number of the degrees of freedom for the (a) FEM model and (b) meshless method. We compared the SOFA FEM mplementaton and the meshless collocaton method wth close accuracy (Fgure 4). For our meshless collocaton method wth nodal ntegraton, we used an explct tme ntegraton scheme wth a tme step of 0.001s wthout any stablty problem. For the SOFA FEM mplementaton, we had to use mplct ntegraton wth a tme step of 0.01s. The calculatons were performed wthn the SOFA applcaton on a sngle Intel Core 5 CPU runnng at 2.67 GHz wth 16 GB of RAM under Wndows 7 operatng system. The SOFA FEM mplementaton took 195ms of calculatons per tme step, whereas the meshless method consumed 20.11ms for calculatons per tme step. Therefore, the meshless collocaton mplementaton n SOFA (along wth other SOFA related operatons such as collson detecton) s roughly 25 tmes slower than the real-tme operaton, whch s slghtly better than the 30 tmes slower performance reported by the Meshless TLED algorthm [20]. The calculaton speed of the meshless collocaton algorthm s governed by the number of partcles and the number of neghbors assgned to each partcle. Fgure 4 Comparson of the SOFA FEM mplementaton and the meshless collocaton method wth close ndentaton accuracy and the theoretcal soluton.
13 13 We also performed a mesh convergence study for the meshless collocaton method by nvestgatng the convergence of the ndentaton amount to the theoretcal value for a fxed amount of force (Fgure 5). After around 6000 partcles, the ndentaton value converges to the theoretcal ndentaton value. Fgure 5 Convergence of the ndentaton value wth ncreasng number of meshless partcles. 4 Conclusons and Future Work The dscussed Movng Least Square approxmaton-based meshless method may be useful n contnuum problems wth extreme deformatons or movng dscontnutes such as cracks or cuts. The presented algorthm for the meshless collocaton method was mplemented as a component for the open source SOFA lbrary, whch s prmarly targeted at real-tme medcal smulaton. Contact mechancs-based verfcaton experments were conducted wth FEBo FEM code, SOFA hexahedral FEM method, and the presented meshless collocaton method. The convergence study results and L2 error norm curves are promsng for the meshless deformable model. For the meshless collocaton methods, the nfluence of the meshless node dstrbuton and the sze of the nodal support rad of the nodes are not well studed. The fndngs of our own experments also suggest further studes of these mportant aspects of the meshless collocaton methods.
14 14 References [1] L. P. Nedel and D. Thalmann, "Real tme muscle deformatons usng mass-sprng systems," n Computer Graphcs Internatonal, Proceedngs, 1998, pp [2] K.-J. Bathe, Fnte element procedures vol. 2: Prentce hall Englewood Clffs, [3] R. Eymard, et al., "Fnte volume methods," Handbook of numercal analyss, vol. 7, pp , [4] A. R. Mtchell and D. F. Grffths, "The fnte dfference method n partal dfferental equatons," Chchester, Sussex, England and New York, Wley- Interscence, p, [5] T. W. Sederberg and S. R. Parry, "Free-form deformaton of sold geometrc models," ACM Sggraph Computer Graphcs, vol. 20, pp , [6] S. F. Frsken-Gbson, "Usng lnked volumes to model object collsons, deformaton, cuttng, carvng, and jonng," Vsualzaton and Computer Graphcs, IEEE Transactons on, vol. 5, pp , [7] K. Waters and D. Terzopoulos, "A physcal model of facal tssue and muscle artculaton," n Vsualzaton n Bomedcal Computng, 1990, pp [8] A. Van Gelder, "Approxmate smulaton of elastc membranes by trangulated sprng meshes," Journal of graphcs tools, vol. 3, pp , [9] M. Bro-Nelsen, "Fnte element modelng n surgery smulaton," Proceedngs of the IEEE, vol. 86, pp , [10] K. Mller, et al., "Total Lagrangan explct dynamcs fnte element algorthm for computng soft tssue deformaton," Communcatons n numercal methods n engneerng, vol. 23, pp , [11] S. Marchesseau, et al., "Multplcatve jacoban energy decomposton method for fast porous vsco-hyperelastc soft tssue model," n Medcal Image Computng and Computer-Asssted Interventon MICCAI 2010, ed: Sprnger, 2010, pp [12] F. Faure, et al., "SOFA: a mult-model framework for nteractve physcal smulaton," n Soft Tssue Bomechancal Modelng for Computer Asssted Surgery, ed: Sprnger, 2012, pp [13] Y. Chen, et al., Meshless methods n sold mechancs: Sprnger Verlag, [14] V. P. Nguyen, et al., "Meshless methods: a revew and computer mplementaton aspects," Mathematcs and Computers n Smulaton, vol. 79, pp , [15] L. B. Lucy, "A numercal approach to the testng of the fsson hypothess," The astronomcal journal, vol. 82, pp , [16] N. Sukumar. (2002, 08/01/2012). The Natural Element Method (NEM) n Sold Mechancs. Avalable: [17] B. Nayroles, et al., "Generalzng the fnte element method: dffuse approxmaton and dffuse elements," Computatonal mechancs, vol. 10, pp , 1992.
15 [18] T. Belytschko, et al., "Element free Galerkn methods," Internatonal Journal for Numercal Methods n Engneerng, vol. 37, pp , [19] M. Müller, et al., "Pont based anmaton of elastc, plastc and meltng objects," n ACM SIGGRAPH/Eurographcs symposum on Computer anmaton, 2004, pp [20] A. Horton, et al., "A meshless total Lagrangan explct dynamcs algorthm for surgcal smulaton," Internatonal Journal for Numercal Methods n Bomedcal Engneerng, vol. 26, pp , [21] M. Pauly, et al., "Meshless anmaton of fracturng solds," ACM Trans. Graph., vol. 24, pp , [22] H. S, "TetGen: A qualty tetrahedral mesh generator and three-dmensonal delaunay trangulator," Weerstrass Insttute for Appled Analyss and Stochastc, Berln, Germany, [23] J.-D. Bossonnat, et al., "Trangulatons n CGAL," n Proceedngs of the sxteenth annual symposum on Computatonal geometry, 2000, pp [24] G.-R. Lu, Meshfree methods: movng beyond the fnte element method: CRC press, [25] H. Hertz, "{Über} de {B} erührung fester elastscher {K} örper," J. für de rene u. angew. Math., vol. 92, [26] D. Franke, et al., "The p verson of the FEM for computatonal contact mechancs," PAMM, vol. 8, pp , [27] N. Schwarzer, et al., "Comparson between analytcal and FEM calculatons for the contact problem of sphercal ndenters on layered materals," Thn sold flms, vol. 270, pp , [28] D. Franke, et al., "A comparson of the h-, p-, hp-, and rp-verson of the FEM for the soluton of the 2D Hertzan contact problem," Computatonal mechancs, vol. 45, pp , [29] F. Pennec, et al., "Verfcaton of contact modelng wth COMSOL multphyscs software," [30] S. A. Maas, et al., "FEBo: fnte elements for bomechancs," Journal of bomechancal engneerng, vol. 134, [31] T. Laursen and B. Maker, "An augmented Lagrangan quas Newton solver for constraned nonlnear fnte element applcatons," Internatonal Journal for Numercal Methods n Engneerng, vol. 38, pp , [32] M. Marchal, et al., "Towards a framework for assessng deformable models n medcal smulaton," n Bomedcal Smulaton, ed: Sprnger, 2008, pp
Modeling, 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 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 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 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 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 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 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 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 informationHermite Splines in Lie Groups as Products of Geodesics
Hermte Splnes n Le Groups as Products of Geodescs Ethan Eade Updated May 28, 2017 1 Introducton 1.1 Goal Ths document defnes a curve n the Le group G parametrzed by tme and by structural parameters n the
More informationAnalysis of 3D Cracks in an Arbitrary Geometry with Weld Residual Stress
Analyss of 3D Cracks n an Arbtrary Geometry wth Weld Resdual Stress Greg Thorwald, Ph.D. Ted L. Anderson, Ph.D. Structural Relablty Technology, Boulder, CO Abstract Materals contanng flaws lke nclusons
More informationHelsinki University Of Technology, Systems Analysis Laboratory Mat Independent research projects in applied mathematics (3 cr)
Helsnk Unversty Of Technology, Systems Analyss Laboratory Mat-2.08 Independent research projects n appled mathematcs (3 cr) "! #$&% Antt Laukkanen 506 R ajlaukka@cc.hut.f 2 Introducton...3 2 Multattrbute
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 informationSimulation of a Ship with Partially Filled Tanks Rolling in Waves by Applying Moving Particle Semi-Implicit Method
Smulaton of a Shp wth Partally Flled Tanks Rollng n Waves by Applyng Movng Partcle Sem-Implct Method Jen-Shang Kouh Department of Engneerng Scence and Ocean Engneerng, Natonal Tawan Unversty, Tape, Tawan,
More informationComplex Deformable Objects in Virtual Reality
Complex Deformable Obects n Vrtual Realty Young-Mn Kang Department of Computer Scence Pusan Natonal Unversty ymkang@pearl.cs.pusan.ac.kr Hwan-Gue Cho Department of Computer Scence Pusan Natonal Unversty
More informationADJUSTING DEFORMATION METHODS FOR VIRTUAL REALITY SYSTEMS
ADJUSTING DEFORMATION METHODS FOR VIRTUAL REALITY SYSTEMS Suene Ferrera Campos 1, Rone Marcos de Moraes 2, Llane dos Santos Machado 3 Abstract Vrtual Realty Systems for tranng have been proposed and developed
More informationSPH and ALE formulations for sloshing tank analysis
Int. Jnl. of Multphyscs Volume 9 Number 3 2015 209 SPH and ALE formulatons for sloshng tank analyss Jngxao Xu 1, Jason Wang 1 and Mhamed Soul*, 2 1 LSTC, Lvermore Software Technology Corp. Lvermore CA
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 informationWavefront Reconstructor
A Dstrbuted Smplex B-Splne Based Wavefront Reconstructor Coen de Vsser and Mchel Verhaegen 14-12-201212 2012 Delft Unversty of Technology Contents Introducton Wavefront reconstructon usng Smplex B-Splnes
More informationLecture #15 Lecture Notes
Lecture #15 Lecture Notes The ocean water column s very much a 3-D spatal entt and we need to represent that structure n an economcal way to deal wth t n calculatons. We wll dscuss one way to do so, emprcal
More informationA unified library of nonlinear solution schemes
A unfed lbrary of nonlnear soluton schemes Sofe E. Leon, Glauco H. Paulno, Anderson Perera, Ivan F. M. Menezes, Eduardo N. Lages 7/27/2011 Motvaton Nonlnear problems are prevalent n structural, flud, contnuum,
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 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 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 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 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 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 informationSimulation of Surface Mesh Deformation in Orthodontics by Mass-Spring Model
292 ECTI TRANSACTIONS ON ELECTRICAL ENG., ELECTRONICS, AND COMMUNICATIONS VOL.9, NO.2 August 2011 Smulaton of Surface Mesh Deformaton n Orthodontcs by Mass-Sprng Model Pathomphong Phannurat 1, Wcht Tharanon
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 information12/2/2009. Announcements. Parametric / Non-parametric. Case-Based Reasoning. Nearest-Neighbor on Images. Nearest-Neighbor Classification
Introducton to Artfcal Intellgence V22.0472-001 Fall 2009 Lecture 24: Nearest-Neghbors & Support Vector Machnes Rob Fergus Dept of Computer Scence, Courant Insttute, NYU Sldes from Danel Yeung, John DeNero
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 informationOverview. Basic Setup [9] Motivation and Tasks. Modularization 2008/2/20 IMPROVED COVERAGE CONTROL USING ONLY LOCAL INFORMATION
Overvew 2 IMPROVED COVERAGE CONTROL USING ONLY LOCAL INFORMATION Introducton Mult- Smulator MASIM Theoretcal Work and Smulaton Results Concluson Jay Wagenpfel, Adran Trachte Motvaton and Tasks Basc Setup
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 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 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 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 informationUser Authentication Based On Behavioral Mouse Dynamics Biometrics
User Authentcaton Based On Behavoral Mouse Dynamcs Bometrcs Chee-Hyung Yoon Danel Donghyun Km Department of Computer Scence Department of Computer Scence Stanford Unversty Stanford Unversty Stanford, CA
More informationy and the total sum of
Lnear regresson Testng for non-lnearty In analytcal chemstry, lnear regresson s commonly used n the constructon of calbraton functons requred for analytcal technques such as gas chromatography, atomc absorpton
More informationAssignment # 2. Farrukh Jabeen Algorithms 510 Assignment #2 Due Date: June 15, 2009.
Farrukh Jabeen Algorthms 51 Assgnment #2 Due Date: June 15, 29. Assgnment # 2 Chapter 3 Dscrete Fourer Transforms Implement the FFT for the DFT. Descrbed n sectons 3.1 and 3.2. Delverables: 1. Concse descrpton
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 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 informationLecture 5: Multilayer Perceptrons
Lecture 5: Multlayer Perceptrons Roger Grosse 1 Introducton So far, we ve only talked about lnear models: lnear regresson and lnear bnary classfers. We noted that there are functons that can t be represented
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 informationModeling of Fillets in Thin-walled Beams Using Shell/Plate and Beam Finite. Elements
Modelng of Fllets n Thn-walled Beams Usng Shell/Plate and Beam Fnte Elements K. He Graduate Research Assstant W.D. Zhu Professor and Correspondng Author Emal: wzhu@umbc.edu; Tel: 410-455-3394 Department
More informationIMAGE FUSION TECHNIQUES
Int. J. Chem. Sc.: 14(S3), 2016, 812-816 ISSN 0972-768X www.sadgurupublcatons.com IMAGE FUSION TECHNIQUES A Short Note P. SUBRAMANIAN *, M. SOWNDARIYA, S. SWATHI and SAINTA MONICA ECE Department, Aarupada
More informationDynamic wetting property investigation of AFM tips in micro/nanoscale
Dynamc wettng property nvestgaton of AFM tps n mcro/nanoscale The wettng propertes of AFM probe tps are of concern n AFM tp related force measurement, fabrcaton, and manpulaton technques, such as dp-pen
More informationInverse Kinematics (part 2) CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Spring 2016
Inverse Knematcs (part 2) CSE169: Computer Anmaton Instructor: Steve Rotenberg UCSD, Sprng 2016 Forward Knematcs We wll use the vector: Φ... 1 2 M to represent the array of M jont DOF values We wll also
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 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 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 informationHigh resolution 3D Tau-p transform by matching pursuit Weiping Cao* and Warren S. Ross, Shearwater GeoServices
Hgh resoluton 3D Tau-p transform by matchng pursut Wepng Cao* and Warren S. Ross, Shearwater GeoServces Summary The 3D Tau-p transform s of vtal sgnfcance for processng sesmc data acqured wth modern wde
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 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 informationReproducing Works of Calder
Reproducng Works of Calder Dongkyoo Lee*, Hee-Jung Bae*, Chang Tae Km*, Dong-Chun Lee*, Dae-Hyun Jung*, Nam-Kyung Lee*, Kyoo-Ho Lee*, Nakhoon Baek**, J. Won Lee***, Kwan Woo Ryu* and James K. Hahn*** *
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 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 informationCloth Simulation. Burak Ertekin i MSc Computer Animation and Visual Effects
Cloth Smulaton Burak Ertekn 7637445 MSc Computer Anmaton and Vsual Effects November 17, 2014 1 Abstract In ths paper, we explan how to make a Cloth Smulaton usng Mass Sprng Model. We explaned the nternal
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 informationAn inverse problem solution for post-processing of PIV data
An nverse problem soluton for post-processng of PIV data Wt Strycznewcz 1,* 1 Appled Aerodynamcs Laboratory, Insttute of Avaton, Warsaw, Poland *correspondng author: wt.strycznewcz@lot.edu.pl Abstract
More informationPhysically based morphing of point-sampled surfaces
COMPUTER ANIMATION AND VIRTUAL WORLDS Comp. Anm. Vrtual Worlds 005; 16: 509 518 Publshed onlne n Wley InterScence (www.nterscence.wley.com). DOI: 10.100/cav.100 Anmatng Geometrcal Models Physcally based
More informationGeometric Error Estimation
Geometrc Error Estmaton Houman Borouchak Project-team GAMMA 3 UTT Troyes, France Emal: houman.borouchak@utt.fr Patrck Laug Project-team GAMMA 3 INRIA Pars - Rocquencourt, France Emal: patrck.laug@nra.fr
More informationLoad-Balanced Anycast Routing
Load-Balanced Anycast Routng Chng-Yu Ln, Jung-Hua Lo, and Sy-Yen Kuo Department of Electrcal Engneerng atonal Tawan Unversty, Tape, Tawan sykuo@cc.ee.ntu.edu.tw Abstract For fault-tolerance and load-balance
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 informationHierarchical clustering for gene expression data analysis
Herarchcal clusterng for gene expresson data analyss Gorgo Valentn e-mal: valentn@ds.unm.t Clusterng of Mcroarray Data. Clusterng of gene expresson profles (rows) => dscovery of co-regulated and functonally
More informationModeling of the Absorption of the Electromagnetic Wave Energy in the Human Head Induced by Cell Phone
Journal of Appled Mathematcs and Physcs, 14,, 179-184 Publshed Onlne November 14 n ScRes. http://www.scrp.org/ournal/amp http://dx.do.org/1.436/amp.14.114 Modelng of the Absorpton of the Electromagnetc
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 informationLobachevsky State University of Nizhni Novgorod. Polyhedron. Quick Start Guide
Lobachevsky State Unversty of Nzhn Novgorod Polyhedron Quck Start Gude Nzhn Novgorod 2016 Contents Specfcaton of Polyhedron software... 3 Theoretcal background... 4 1. Interface of Polyhedron... 6 1.1.
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 informationThe Theory and Application of an Adaptive Moving Least. Squares for Non-uniform Samples
Xanpng Huang, Qng Tan, Janfe Mao, L Jang, Ronghua Lang The Theory and Applcaton of an Adaptve Movng Least Squares for Non-unform Samples Xanpng Huang, Qng Tan, Janfe Mao*, L Jang, Ronghua Lang College
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 informationAccounting for the Use of Different Length Scale Factors in x, y and z Directions
1 Accountng for the Use of Dfferent Length Scale Factors n x, y and z Drectons Taha Soch (taha.soch@kcl.ac.uk) Imagng Scences & Bomedcal Engneerng, Kng s College London, The Rayne Insttute, St Thomas Hosptal,
More informationCHAPTER 3 SEQUENTIAL MINIMAL OPTIMIZATION TRAINED SUPPORT VECTOR CLASSIFIER FOR CANCER PREDICTION
48 CHAPTER 3 SEQUENTIAL MINIMAL OPTIMIZATION TRAINED SUPPORT VECTOR CLASSIFIER FOR CANCER PREDICTION 3.1 INTRODUCTION The raw mcroarray data s bascally an mage wth dfferent colors ndcatng hybrdzaton (Xue
More informationFEATURE EXTRACTION. Dr. K.Vijayarekha. Associate Dean School of Electrical and Electronics Engineering SASTRA University, Thanjavur
FEATURE EXTRACTION Dr. K.Vjayarekha Assocate Dean School of Electrcal and Electroncs Engneerng SASTRA Unversty, Thanjavur613 41 Jont Intatve of IITs and IISc Funded by MHRD Page 1 of 8 Table of Contents
More informationAADL : about scheduling analysis
AADL : about schedulng analyss Schedulng analyss, what s t? Embedded real-tme crtcal systems have temporal constrants to meet (e.g. deadlne). Many systems are bult wth operatng systems provdng multtaskng
More informationSimplification of 3D Meshes
Smplfcaton of 3D Meshes Addy Ngan /4/00 Outlne Motvaton Taxonomy of smplfcaton methods Hoppe et al, Mesh optmzaton Hoppe, Progressve meshes Smplfcaton of 3D Meshes 1 Motvaton Hgh detaled meshes becomng
More informationSubspace clustering. Clustering. Fundamental to all clustering techniques is the choice of distance measure between data points;
Subspace clusterng Clusterng Fundamental to all clusterng technques s the choce of dstance measure between data ponts; D q ( ) ( ) 2 x x = x x, j k = 1 k jk Squared Eucldean dstance Assumpton: All features
More informationCell Count Method on a Network with SANET
CSIS Dscusson Paper No.59 Cell Count Method on a Network wth SANET Atsuyuk Okabe* and Shno Shode** Center for Spatal Informaton Scence, Unversty of Tokyo 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
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 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 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 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 informationMachine Learning 9. week
Machne Learnng 9. week Mappng Concept Radal Bass Functons (RBF) RBF Networks 1 Mappng It s probably the best scenaro for the classfcaton of two dataset s to separate them lnearly. As you see n the below
More informationD shape manipulation via topology-aware rigid grid. Introduction. By Wenwu Yang and Jieqing Feng *
COMPUTER ANIMATION AND VIRTUAL WORLDS Comp. Anm. Vrtual Worlds 2009; 20: 175 184 Publshed onlne 1 June 2009 n Wley InterScence (www.nterscence.wley.com).285 2D shape manpulaton va topology-aware rgd grd
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 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 informationSimulation: Solving Dynamic Models ABE 5646 Week 11 Chapter 2, Spring 2010
Smulaton: Solvng Dynamc Models ABE 5646 Week Chapter 2, Sprng 200 Week Descrpton Readng Materal Mar 5- Mar 9 Evaluatng [Crop] Models Comparng a model wth data - Graphcal, errors - Measures of agreement
More informationSteps for Computing the Dissimilarity, Entropy, Herfindahl-Hirschman and. Accessibility (Gravity with Competition) Indices
Steps for Computng the Dssmlarty, Entropy, Herfndahl-Hrschman and Accessblty (Gravty wth Competton) Indces I. Dssmlarty Index Measurement: The followng formula can be used to measure the evenness between
More informationReal-time Motion Capture System Using One Video Camera Based on Color and Edge Distribution
Real-tme Moton Capture System Usng One Vdeo Camera Based on Color and Edge Dstrbuton YOSHIAKI AKAZAWA, YOSHIHIRO OKADA, AND KOICHI NIIJIMA Graduate School of Informaton Scence and Electrcal Engneerng,
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 information3D Virtual Eyeglass Frames Modeling from Multiple Camera Image Data Based on the GFFD Deformation Method
NICOGRAPH Internatonal 2012, pp. 114-119 3D Vrtual Eyeglass Frames Modelng from Multple Camera Image Data Based on the GFFD Deformaton Method Norak Tamura, Somsangouane Sngthemphone and Katsuhro Ktama
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 informationFeature Reduction and Selection
Feature Reducton and Selecton Dr. Shuang LIANG School of Software Engneerng TongJ Unversty Fall, 2012 Today s Topcs Introducton Problems of Dmensonalty Feature Reducton Statstc methods Prncpal Components
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 informationApplication of Learning Machine Methods to 3 D Object Modeling
Applcaton of Learnng Machne Methods to 3 D Object Modelng Crstna Garca, and José Al Moreno Laboratoro de Computacón Emergente, Facultades de Cencas e Ingenería, Unversdad Central de Venezuela. Caracas,
More informationKinematics of pantograph masts
Abstract Spacecraft Mechansms Group, ISRO Satellte Centre, Arport Road, Bangalore 560 07, Emal:bpn@sac.ernet.n Flght Dynamcs Dvson, ISRO Satellte Centre, Arport Road, Bangalore 560 07 Emal:pandyan@sac.ernet.n
More informationQuerying by sketch geographical databases. Yu Han 1, a *
4th Internatonal Conference on Sensors, Measurement and Intellgent Materals (ICSMIM 2015) Queryng by sketch geographcal databases Yu Han 1, a * 1 Department of Basc Courses, Shenyang Insttute of Artllery,
More informationMultiblock method for database generation in finite element programs
Proc. of the 9th WSEAS Int. Conf. on Mathematcal Methods and Computatonal Technques n Electrcal Engneerng, Arcachon, October 13-15, 2007 53 Multblock method for database generaton n fnte element programs
More informationHarmonic Coordinates for Character Articulation PIXAR
Harmonc Coordnates for Character Artculaton PIXAR Pushkar Josh Mark Meyer Tony DeRose Bran Green Tom Sanock We have a complex source mesh nsde of a smpler cage mesh We want vertex deformatons appled to
More informationELEC 377 Operating Systems. Week 6 Class 3
ELEC 377 Operatng Systems Week 6 Class 3 Last Class Memory Management Memory Pagng Pagng Structure ELEC 377 Operatng Systems Today Pagng Szes Vrtual Memory Concept Demand Pagng ELEC 377 Operatng Systems
More informationReading. 14. Subdivision curves. Recommended:
eadng ecommended: Stollntz, Deose, and Salesn. Wavelets for Computer Graphcs: heory and Applcatons, 996, secton 6.-6., A.5. 4. Subdvson curves Note: there s an error n Stollntz, et al., secton A.5. Equaton
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 information