NII-Electronic Library Service
Cloth
Shape
Design
b(ERO
7'ifl
y
Umberto
CUGINI
PolitecnicodiMilano
9)iA<)L
ts
ti
iy
S2"-=RiYZNX\
Thispaperpresentsaphysics-basedmodel forvirtual cloth
simulation speGifically targetedforCAD applications inthe
clothing sector,A di$crete
particlebased
model inNewtonianformulation
is
considered, inwhich particlegridsdefinethefabric
structure,Gridparticles
interactwith each other throughforces
describing
Iocaltractionand compression, bendingandshear
properties
dependingon thefabricmaterial, lnteractionswith thesurrounding environment are describedinterms of external forces,imposed kinematicgeometricalconstraints, and collisionswith obstacles. Collisiondetectionis
based
on bounding box hierarchiesand region
decomposition.
The
model isimplementedintheframeworkof a 3D
graphical
en-vironment forvirtual garment
design,
andincludes
operatorsformonitoring thewhole 3D garmentshape
definition
process,
e.g. panelsewing, buttonldartinsertion,mu[ti-layered fabric composition, and other garment
finishings.
Asapplications, male and femalecloth shapes simulated on virtualmanne-quinsareshown.
Results
have
been
assessed intheframe-work of Europeanand [talian
projects
forvirtualclothdesign.1.Introduction
Enhancements inresearch
issues
on shape modellinghave
ledtothedevelopmentof methods thatallow an
increasingly
versatile representation of
3D
complex scenes, fromgenera-tionlmodificationof 2D-3D static template geometries up
todynamic simulatjon of complex-shaped objects moving
inspace, A particularly$trong
demand
is
addressed tonon-rigid body modelling. Shapes stop beingtime-invariant and
become propertieschanging withtime,whose staticand
dy-namic behaviourcan berepresented, analysed and predicted,
Handlingdeformableor moving shapes becomes a strategic
capability torvisualization and simu]ation purposes both in
computer graphicsand computer aided designforindustrial
purposes.
Modelling
deformableshapes isrequired inmanyapplica-tions.Industrialdesignand manufacturing tasksrequire,
for
instance,tomodel flexiblepartsthathaveto bemoulded byNC tools,or movedltransferred byautomated arms of
indus-trialrobots. Examplesare packscontaining fluidsand granular materials, sheet metals, wires and cables, and so on, An in-creasing demand of new soft modelling toolscomes alsofrom
thebio-medicalcontext fororganic tissuemanipulation, with virtualsimulation forinspection,surgicaltraining,
laparoscopic
surgery, etc, Besides,the fieldof computer graphicsgreatly
demands
shape animation toolsfor
movies, cartoons, and vid-eogames.
Not
surprisingly,cloth motion has recentlybecome atopicof
large
investigation.
Virtual
garments
orvirtual clothdrapingeffects are,infact,recurrent visualization elements neces-sary
for
computergraphics'
animated scenes, Beyond this,a strong impulsecomes fromindustrialsectors, e.g,for
ap-parel
and home fabricfurnitureproduction:
inthisarea, infact, modelling toolsare widely demanded toassist theprocessofclothdesignand make itfastenUntilnow, indeed,the skilled-labourdependentnature of apparel and upholstery designdid
nothave encouraged toalargeextent of automation and use
ofcomputer programs.AIthougha
quite
well-establishedtech-nology already exists inCADICAM modules for2D
pattern
editing,nesting and NC cuttinglsewing, textilecompanies still
comp[ain about thelackof more complete clothing-oriented systems, withmodelling toolsspecializedinclothshape draft-ing/designand, particular]y,virtual cloth analysis. Infact,if
such an analysis, bothat afunctionaland aesthetic level,
becomes
po$$ibleat an early stage of apparel design,thisshould aliow toreduce thenumber of intermediate
garment
prototypesand,definitively,
save timeand cost,Thiscommon interesttowardsclothmodelling on behalfof
both
computer graphicsresearchers and clothmanufacturesis
dictated
by
different
needs. Theformer
need torepresent clothbehaviour
for
visualizationand animation purposes,inwhichresults should bevisually acceptable, generallyperformedwith
lowcomputational cost.Atthislevel,geometricalinformation
is
sufficient,as noimportance
is
giventothephysicalreliabil-ityofresults.
On
theother hand,clothmanufacturers demandTTifo\mxwle
spedalissueotjapanesesocietyforthescienceofdeslgn
vol.15-4 no.60 2aOB
NII-Electronic Mbra
11
CAD toolsprovidinga more accurate and realisticc[othshape
predictionforsimulation and qualjtyevaluation tasks.
This
deeperlevelof shape representation requires aphysics-basedmodelling approach
in
which a physicalframework
supports the geometricalsubstrate, and constitutive materjal propertiesare considered tocharacterize theshape
behaviour.
Section
2 presentsa survey on some ofsuch
geometry-
orphysics-basedmethods used forclothanimation and simulation.
Regardlessofthecontinuous or discretemathematical
rep-resentation used forcloth,what emerges from
geometry-
orphysics-based
modelsproposed
bythescientific community, isthattheydescribecloth eitheras a geometricalentity or a mechanical system, butvery rarely theyhandlec[oth as a "constructionprocess".
Thisisakeyissue,inordertomake acloth mode] usable forCAD applications intheclothing indus-try.Althoughwhen themodelling core isinsertedwithin graph-icaluser interfacesor advanced virtual reality environments, most of such systems donot havesufficiently complete CAD functionalities,capable todrivethegarmentdesignand manu-facturingprocessas itisconceived
by
theapparel industryforproductiontasks.
With the intentof gettingcloser tocloth manufacturers' needs, we here presenta physics-basedmodel forvirtual
clothshape simulation,specifically oriented totext"eCAD
ap-plieat[ons.
Such a model hasbeen implementedwithinasys-tem,named SoftWorld2,O,developedbytheKAEMaRT group of UniversitbdiParma and PolitecnjcodiMilano,ltaly.
2.
CIoth
Modelling
:State-of-the-art
Therehavebeen several attempts tomodel clothand
simu-late
its
staticand dynamic response while interactingwith an externa] environment. The most naturalapproach ishandlingfabrics
by
thesamegeometricaYmechanical
standards as other highlyflexiblemateria[ sheets, such as rubber films,pa-perlayers,etc.Accordingtohighermathematical approaches,
2D models forthin
flexible
objects can be regarded inturnasdegenerategeometriesdirectlyderivedfrom
3D
volume mod-elling.Ibgether
withthesegeneralapproaches,further
modelshave
been
investigated,specificallyoriented tomodel thedis-crete
fiber
structureoffabrics.
Theoreticalstudies about cloth behaviourstarted about sixtyyearsago,
from
E:
Peirce's
first
work publishedin
1
937
[1]
tothemost promisingmethod of modelling cloth behaviour introducedinthe`70s
byS.
DeJong
and R.Postle[2].
127ifl)#{IIKe
specialissueefjapanesesocletyferthescleneeefdesignvol.15-4 no.60 200B
Fig.1
.
Weil'sclothmodel bycatenary curvesThe researchers and work developed inthoseyearshad been mainly fundedbythetextileindustry;interestjngly,
how-ever, studies on cloth modelling have been abandoned by thetextileindustryforseveral years.Aroundthe halfE80s,
as
suffjcientlymature computer resources begintobeavailable forcomputer graphicsand CAD design,the prob[emof non-rigid bodymodelling becomes aresearch subject ofincreasing
interest.Computer scientists facetheunanswered questions
proposingcomputable models and
focusing
ondiscretization
methods aspects fornumerical simulation.
Among geometry-basedrecent approaches, we
here
men-tionthe wel]-known model from
J,
Weilin1986
[3],
based
on thedescriptionofdrapes
through suspended catenary curves(Figure
1).Anothercontribution comes fromT,Agui:swork in '90[4],
modelling theeffect of abended elbow. ]n'90-'92 B.K. Hindsand J.McCartney[5]
realize one of thefirstCAD tools specialized forclothdesign,basedon surface modelling.Physics-basedapproaches conceive objects intheir
two-foldgeometricaland physicalnature. Physicsenters into
geometry through mathematical formulationsderivedfrom
dynamic laws,structuralmechanics, elasticity theories,flu[d dynamics, or other physicalcontexts. Objectsare dealtas dynamicsystems characterized bymass properties,with
spe-cificstructuralcharacteristicsdependingon material, They are
subjected tointernalinteractions,and interactwithan external environment throughtheactionofexternalforceslstressesand
theresponse tocollisionswith obstacles. Although
computa-tionallymore hesigent,thesemethods are themost suitable
toperformreliablecloth simulations
for
functional
evaluationtasks,withmore accurate
quantitative
results.The underlying mathematical models associa ±ed tophysics-based
formula-tionsaregenerally
systems of ordinary andpartial
differen-tialequations(ODE/PDE,
Iinear
or non-linear),which Ganbe
solved byrobust and well-known numerical techniques. These
NII-Electronic Library Service
timediscretizat[onof PDE models fordynamic problems,as well as
finite
element(FEM),
finite
difference
(FDM)
orspec-tralmethods
(SM)
forspace discretizationinPDE models that arise from bothstatic and dynamicproblems.
Basically,
wedistinguishtwomain different
physicsbased
approaches:con-tinuousand discrete.
(Continuous
models) These interpretthecon$idered objectas a continuous medium subjected tothelawsof
elastic-ity.They giveraise toa PDE model generallysolved by FEM, which are accurate methods forsimulations occurring instatic conditions, butare computationally more expensive. The most classical of thecontinuous mode[s ispropo$ed bythe
group
of.
D.
lerzopou[os,
J.
Platt,A.H.Barrand K.Flieischerinthe years'86-'88[6],
as an extension of thestructuralmechani-cal approach inwhich the surface deformationisdescribed
by
means ofdisplacement
functions
with respect toequilib-rium
positions.
Another
continuous model isproposed byJ.Eischen
in
'96[7],
based
on a non-linear shelltheory.
Figure2
shows Eischen:s shellconfigura±iondefiningpositions¢ offabric
shellpointsatadistance
E
along adirection
t,intermsof
fabric
midsurface:s referencepoints
th
.
M,
Aono'smethod in'90[8]
simulates insteadtheeffectof wrinkle anddrape
propagation,based on theelasticitytheory and D'Alembert principle,lnthiscase, theresultingPDE modelis
a modifiedversion of thewave equation
propagating
in
a continuous elasticmedium. L Li,M.Damoran and R.K.Gay
ini93-:96
[9]
describe
themodel ofa clo±h
immersed ina quasi-stationary viscous fluidby
combining Navier-Stokesequations andler-zopoulos' continuous model.
(Discrete
models) These represent objects as discretesystems composed of a finitenumber of simple mechanical elements interactingwith each other and with theexterior, subjected to Newton:s or Lagrange'slaws.Inthe$o-called
particle-basedmodel, inparticulactheconsMutive mechanical
elements are directlyparticles,havingmass, subjected to ln-ternaland external forces
{e.g.
mass-spring models).fypicaliy,
two main approaches are used among thediscretemode[s:
forcebasedor energy-based. Inthefirstcase, thecontribution of internaland external forcesisestimated foreach
particle,
with successive numerical solution of thecorresponding
ODE
system, derivedfrom NewtontsIaw.Inthesecond case,
suit-able energy functionsare introduced,and theoverallsystem
dynamics isderivedfrom energy conservation lawsor
La-grangianformulations.
,L.II
r- t.v
Fig.2,Eischen'sshellcontiguration
nLaul
vveft
-Fig.3,Particle-basedmodel otaweven fabric
C.R.Feynman in'86
[1
O]isthefirsttouse a particlebasedmodel torepresent thediscretemicroscopic structureof
fabric
threadswhose yarniswoven according totheorthogonal
di-rections of weft and warp
(Figure
3).
D,E.Breen
andD.
House
since early ]90s
[]
l]have been investigatingindetailwhat they call the"textilemechanism'', byuslng energy-based
par-ticlemode[s inwhich theinternalinteractionsare describedin
terms of stretchinglrepelling, bendingand trellisingenergies.X, Provotin]95-'97
[1
2]
uses a Newtonian mass-spring model inwhich theelasticforcesactnot onlyalong the weft and warp directions,butalsoalong thetwodiagonal
directions
ofeach cellinthetextile
grid.
Sincethefirst'90s untilnowadays, thecontribution of
Mi-raLab [aboratoryin
Geneve,
i.e.N.Magnenat-Thalmann, D.Thalmann, R
Volino
et al,[13,
141
hasbeen
fundamental,lnitiallyinfluenced
by
thecontinuous lerzopoulos]approach,successively
{and
stillnow) they havebeenworking mainly on discreteparticle-basedmodels, more efficient tostudynon-lineardeformationsin
dynamic
cases, using non-$tructured triangularmeshes and several techniquesof collision man-agement. Besides,importantcontributions come from M.RGascuel-Cani,
C.
Puech,
and M.Desbrunat IMAG inGrenoble[1
5],
working on structured thindeformablebodiessuch asfabrics,
and more generalmodelling aspects re[ated to2D/3D
deformation,
Recently,
aiming at improvingthecomputational efficiencyT-tf{)-'fimvewfie
specialissueo"apaneseseeletyforthescieneeofdesign vel.15-4 no.60 20e8
NII-Electronic Mbra
13
of cloth modelling methods, new techniques
for
thetimedis-cretizationof ODE systems forparticlebased
grids
havebeen particularlyinvestigated.After
a largeuse ofexplicitlow-ordermethods
(e.g,
explicit EuleOinthepreviousdecade,theinter-est hasbeen recently addressed towardsimplicitschemes, as
proposed
by
theprecursorwork ofA,Witkinand D,Baraffin' 98[1
6].Successfulcontributions intheareaof implicit meth-ods forODE
systems comefrom
YM. Kang at al.in'OO[1
7], and B,Ebherardt,O.
Etzmuss,M.Hauthand J,Gross[1
8]in' O1.
Thelatter
have
recentlyshown how discretemodels canbe derivedfromcontinuous models when applying space
dis-cretizationby FEM orFDM totheoriginal PDE problems
[19].
Asemi-implicit BDF techniquehasbeen proposedthisyearby
K.J.
Choi
and H.S.Ko[20],
handlingpost-bucklinginstabi[ityand showing significantlyrealistic results of cloth draping.
The
problem
of collisionmanagement has beenconsjderedin
some oftheabove-mentioned models, inordertotakeintoaccount interactionswith external objects and self-collisions.
two
aspects are involved:collision detectionand collision re-sponse. The formecinparticular,
playsa significantpartin
the totalcomputational timeof simulations, asproxjmity
queries
and check$ forpossiblecollisions are requiredbetween
pairs ofelementary object entities(points,
edges, triangles).fech-niques
for
collisiondetectien
optimization havebeen studied, such as voxel subdMsion, octree subdivision, bounding box hierarchy,proximitytracking,and curvature-based methods. Detailsabout collisiondetectiontechniquescan befoundin[13,12,14].
foconclude, hybridmodels
propose
a compromise betweenthecomputational speed ofpurely
geometric
representations with the higheraccuracy of physics-basedmethods applied incertain object sub-regions. Among them,we mention I.J.Rudomin's work
[21]
usinggeometrical
approximations formodelling cloth draping,and T.L Kuni'sspring-based model
[22]
introducingmetric and curvature energies, both from'90. We
also mention E
laillefeurts
study about horizontalandvertical drapes
[23],
N.fsopelas'
garment models builtfrom thinpipessubjected togravity[24],
both
fromt91,
and S.G. Dhande'stechniquesin'93
[25]
for
modeiling clothassweep-ingsurfaces subjected toelastic forces.Theaspirationtoreal timecloth animation hasrecently renewed theinteresttowards hybridapproaches. In'O1, M.
Oshita
etal.[26]
modelfabrics
throughvery coarse particiemeshes subjected todynamic
laws,
derivingglobalcloth shape motionsby
smoothing and14T-if{y\ffsckfig
specfalissueetjapanesesocietyferthesclefieeotdesignve].15-4 ne.6D 2008
(a)
(b)
Fig,4,Threadpatterns:(a)knitted;(b)woven
interpolation
techniques.A similar ideaisproposed in'Ol
byYM.
Kang etal.[27]
tofastenaparticle-basedmodel viacu-bjcsplinejnterpolationwith wrinkle generation,and in'02 by Rudomin
[28]
where mesh particlesmove accordingly tosets ofellipsoidsdefiningmannequins. Foradetailed
overview onphyslcs-based
modelling, see[29,
1
4,
30].
3.
The
Particle-basedCloth
Model
3.1.
FabricProperties
Fabrics
come ina variety of types,depending on thetype of constitutivefibers,theway fibersare combjned, spun andtwistedtogetherto
generate
threads,thetype of thread pat-tern,and thedegreeof loosenesslwidthof threadplots.Also
possiblechemical or mechanical treatments
(such
as ironing,starching, reactiontoenvironment conditions, and so on} con-tributetodjfferentbehaviouroftextiles.
It
is
clear thata propercloth model should take into ac-count allthese factors.Fabricisnota continuous medium, as threadsareinterlaced
witheach other at afinitedistance,itis highlyflexible,and presentsanisotropic,non-lineanyhysteretic and timedependent
characteristics. Differentlyfrom knittedfabrics,which presentmore complicated threacipatterns,
woven textilesare characterized bya regulac structured plot composed of twosetsofyarns,thewarp and theweft, as they are interlacedbya weaving loom
(Figure
4}.The most typical woven structure ischaracterizedby
two orthogonal weft andwarp directions.As most garments are composed of woven
textUes,thiswiH bethestructure mainly considered inour
NII-Electronic Library Service
3.2.
Particleand lnternalForceCharacterizationlncontlnuous formulations fabricsare handled as
'
formablethinshells or p[atebeams, according tostructural mechanics approaches
(Section
2).Though supportedby
nu-merically robust techniques such as finiteelement methods, nevertheiess, continuous models are more appropriate for
ma-terials$ubjected to
lower
deformability
and smallerdisplace-ments, The
discrete
nature of woven textilesas interlacedsets ofthreadsand the
highly
flexible
behaviour
seem, on thecontrary,tobebetter
de$cribed
interms of discrete represen-tationssuch as particle-basedmodels. FollowingtheideaofBreen,clothisinterpretedas a
"mechanism"
of discreteparts{point
masses, or "particles'') interactingwith each other and with an external environment. DifferentlyfromBreen'senergy-basedapproach, more suitableforstaticconditions,we here use a
force-based
representation,in
order toinclude
dynamic
problems,
e.g,tovalidate the cloth model when tested on moving virtual mannequins.Fabricsareassumed tohavea negligiblethickness, i.e.a2D topologyinwhich flatshapes arecut
from
open connected andbounded figuresF [ R2with
piecewise
regularboundary
OF(e.g.
a closed Ioopof linearorcurved edges). SimilarlytoBreen'sand Provot'smodels, thechosen
particle
meshas-sociated tofabricpanefisisas±ructured
2D
gridwhosecoor-dinatelinesare definedfromwarp and weft directions.
Interior
particlescorrespond togridnodes, locatedat warplweftthreadintersections,while
boundary
particlesaredefined
fromintersectionof
grid
lineswiththefabric
borden
Ttiangularele-ments are
derived
from
originalrectangu[ar cells byadding diagonals,according toFigure5.Panel"sgridtopology characterizes the internaldiscrete forcedistribution,
defined
bylinearor torsionalspringscom-putedfromjnteFparticleconnections. lakingintoaccount the
woven structure ofthreads,
forces
have
been classified into threemain ca±egories(Figure
5):
-
stretch[nglrepellingforces,
actingto
keep particlesat rest
distance
(modelled
as Kelvinvisco-elastic springs directedalong weftand warp) ;
-
bending
forces,acting out-of-plane tokeep objects flat
(derived
by
torsionalmoments normal totheirsupporttace)l
-
trelfising(or
sheari forces,acting tocontrast any possibledeformationof therectangular ce[ls
(modelled
again throughtorsionalmoments normal tothecells).
Fig.5.Particlegridassociated totabricpanelswith lnternaltorce characterization
The values ofIocalspring, bendingand trellisingconstant$ associated toany
group
oftwo,three,orfourparticleshavebeen derivedfromglobalmechanical dataestimated bythe
Kawabata EvaluationSystem
(KES).
Thisisa well-known sys-tem measuring structuralpropedies
offabrics
depending
on the type[31].
Itowes itsformulationtoS.
Kawabata
in1975,and isbased on experimental stress and
deformation
testsperformedon rectanguiar fabricsamples with fixed
dimen-sions. Theseare
(gf'=
gram-force):1
)
geometryand mass tests,tomeasure thethickness[mm]
of the sample at pressureO.5gf7cm!,and itsmass density
[nag/cmL'];
2) tensHetests,tomeasure thelinearity
[-],
expressed as average percentageofelongationin
thesample, itsmation energy
[
gf・cm/cm2]
and tensileresilience[%];
3) compression test$,measuring thelinearity
[-],
or averagepercentage
of shortening inthe sample, itsdeformation energy[
gf-cm/cm2]
and compression resHience[%l
;4) bendingtests,which measure thebendingrigidity
[gye・cm2/
cm
]
and itshysteresis[
gt'・cmi/cm];
5
)
shear tests,which measure the shear stiffness[
gflttm greel]andtheshearhysteresis
at O.5eand 5"respectively
[gLlcm];
6) superficial
property
tests,to measure theaverage frictioncoefficlent
[-],
the
corresponding standard deviation[-]
and the superficial rugosity
[um].
Other
systems exist, measuring similar structuralfabric
propertieswith differentunit measures(e.g,
EASZ or UNI-ISOstandard$}, Tnyifiy#ffxksee specialissueofjapanesesocietyforthescieneeofdesign vol,15-4 no,60 200S NII-Electronic Mbrar
15
3.3.
From
FabricstoGarment ModelsGarment
shape isthecomplex result ofassembled compo-nents, Basically,itlstheresult ofone- or multHayered fabricassembly, withfurtherfunctionaland aestheticfinishings.
Re-sortingtoa continuous geometricalrepresentation, a full
gar-ment shape can
be
formally
approximated as anon-manifoldentity
G=,.L.
!..,F"U,.Y.,,"i'
(])
where Fl・are
fabric
layers
(whose
supports are 2D manifolds inR3)representing thedeformedspatialconfiguratjon oforigi-nal fabricpanelsF,g RU,fori=1,2,..., nF, and A, are
domains
inR3 forrjgid orsoftaccessorjes, fori=1,2,,.., ,z,.When
exist-ing,the
latter
aregenerallyfewsmall entities, such asbuttons,
hooks,zips,paddings,etc. FabricIayersare sewed witheach
other along
portions
of theirboundary,or canbe
partiallyortotally
attachedin
interiorsub-regions with o±her
layers(as
itoccurs inmultj-layeredfabricfor
jacket
creation,by
usinglin-ings,
strengthening adhesives, etc).
Under
thesepremjses,
afterparticlegriddiscretization
of singlefabriclayers,furthersteps have tobe
considered, togenerate
a complete particle-basedmodelfor
afull
garment
shape:a)
sewing of
fabric
panels;b)dartinsertion; c) fabric
layer
overlapping ;d)buttoninsertion;
e} placementofthe
full
garmentmodel ina 3D space ration{e.g.
on a mannequin).Algorithm
(a)
dealswith2D
panelsthathavetobe sewed togetheralong assigned borders,Tb
performaseam between existing grids,particlescoming fromthetwopanels,
locatedalong each sewing edge, need tobe properlymerged along theborders.Suppose apanelA hastobesewed with a panel B,Foreach $ewing edge,
let
nbe
thenumber of original par-ticlesof theedge inA and m thenumber of originaiparticles ofthecorresponding edge inB.If
n<m, m-n paniclesare added topanel
A bylocatingthem inintermediatepositions
(analo-gous[y,
intheother case). When n=m, each i-thparticleofA ismerged withthecorresponding i-thparticleofB,fori=1,2,...,n,
by
movingparticle
pairstoan intermedjateposition,nparticles belongingtoone otthetwoborder
edges are thenremoved,Correspondingly,mesh topology around each sewed edge is
updated. Internaledges are displaced,and new
internal
edges,16ftfly#MzakKe
speciaHssueofjapamesesocietyferthes'elenceofdesignvoL15-4 ne.60 200B
edge eiA edge eie
Fig.6.Sewingprocessbetween twofabricpanels /
(1)
$eam definitionas correspondance 1to1between vertices ofpanelborders;
(2)
detaiiofthesewing processover asingle edge pajnDashedlines$how theoriginalparticlegridbeforebeingmodified bythe
sewlng process.
forceljnksand trianglesaround added
particles
areadded(case
nt m). Internaledges,force
links
and triangiesaround removedparticles
are removed. Figure6
shows a sewing process be-tween two particie-basedmodels offabricpanels.
Case
(b}
issimilar to(a}
becau$e dartscan beregarded asspecialseams along a fixedsequence ofedges, inwhich
pan-eiB iscoincident with panelrf,
ln
other words, one(or
more)vertice(s)is
(are)
sewed with another(otheO
vertice(s)belong-ing
toboundatyloopsof thesame panel.Case
(c)
models thepossib[epresenceof several fabriclayers
thatmay have tobeplacedontop
ofeach other,This
isuseful incase of
possible
paddings,linings,etc(e.g.,
tostrengthen
jackets).
The ideaisto change physics-based properties(masses,
springs, bendingand trellisingforces}in-sidea
fabric
region,substitutjng them with propervalues that includealleffectsofsingleadded layersthroughtheeffectofa unique equivalentmaterial.Ca$e
{d)
includestheeffect of buttonswhen two panelsto beconnected are sufficientlyclose toeach other. Each virtual buttonis
simulatedby
connecting thetwo gridpointsclosesttotherequired
button
positjon,
imposingsignificantly hard spring forceeffectsand massjncrement.
Final[y,step e)
defines
the3D configuration of the as-sembled cloth]s partjclegridafterdefinitionof mannequins orother rigid objects inthe scene and specificationof mapping laws2D
-,3D.
NII-Electronic Library Service
4.
DynamicSimulation
4.1.
Newtonian
Constrained
Dynamics
Thecomplete
particle-based
model ofagarment
isdefined
byasystem S=
{P,
:i= 1,,,, IV}ofIVparticles,
eachcharacter-izedbymass m,,
position
r,, velocityy,=r,
and accelerationa---r,, fori=1 ,2,..,IVLParticles'configurationatan initialtime4,is definedbyknown
positions
{r.}
and velocities{v,,},
for
i=1,.,, N. Tbpertorm
a dynamic analysis ofgarment
movements inspace, we consider Newton's lawsforparticlesystems,
The
analysis iscarried out intheframeworkofconstrained dynam-ics,totakeintoaccount thepossibleexistence of imposedconstraints,i.e.specific
geometric
and klnematicconditionsthatobjects must respect.
The dynamic behaviourof thesystem isthusobtained as solution ofthetollowingNewton'sproblem
Vi
--
1,2,..,N,・
ri<to)=rio
(2)
Vi(to)= Vio miii = Fi
(ri
,..rrv
,"i ,..iN,t)
whereE.
= FiCi"`}+F}(ext)(3}
is
thetotalforce
actingon thei-thparticle,sum of aH internaland externalcontributions.Eq.
{2)-{3}
isaCauchy problemforasystem of3N ordinary
differential
equations ofsecond order,equivalent toa system of 6Nequations of firstorder, Solutions are trajectories
{
r,(t)},
as well as velocities{
v,(t)},
forany par-ticlei= 1,2,,.,ACfortE[4, T].The motion of thesystem isuniquely determinedunder
hy-potheses
of sufficientlyregulardataF,as functionsofparticle
positions,velocities and time,allinfactsatisfied inour linear forceestimation method
[32].
The firstaddend inEq.(3),
in fact,iscomputed fromFiCust).FiCsiretcit)+FJrep}+FiCbe,td)+Fi{shettr),
(4)
where
E("'e`'h)
andECmp)
aretheresultantstretching andrepel-lingforcesacting on
particle
i,F,U'""diistheresultantbendingforce
on i,and Fi('Se"Dtheshearingforce
on i.The
approxima-tionmethod forinternalforcecomputation depends,evidently, on thechosen particlediscretizationcriterion
($ection
3.2).
Thesecond term is
given
byE{ext}=E(a)
÷IilC'},(5)
where F,C")isthe sum of active external forceson i,due to
gravity,wind, viscosity, orimposed constant stresses, and
F,(")thesum of allreactive forcescompensating
possible
con-straints imposed on i,Among theseveral methods for con-straintmanagement
(e,g.,
penaitymethod, Lagrangemultipli-ers, rate-controlled constraints, and dynamic constraints
[33,
34,
35]),
the dynamicconstraint method hasbeenhereimple-mented, since
it
permitstoapply multiple constraints tothesame
particle
and ensures the respect of alltheconstraints ateach time$tep ofthesimulation,
We
hererefer toequality(or
bilateral)constraints,i.e.havingthe generaiformC(rl,.,rA7,fr,,,.tN,t)=
±
O,{6)
beingC a known sufficientlyregularscalaror vector function. Differently
from
themostgeneral
rheonomous constraints,scleronomous constraints do not
depend
explicitelyon time.Inour model we have imposed,
for
instance,
constraints suchas fixed
positions,
constant distancesbetweenpoints,
area orvolume conservation, as well as
fixed
trajectoriesfor
pointsorparticles.
lermsF,`",tobeincludedintheRHS ofEq,(2)
for each particleiinvolvedinsome constraint,arederived
from
soiutionofacertainlinearalgebraicsystem[32].
The solution of Eq.
(2)
can benumerically estimated byus-ingproperexplicitiimp[icit timediscretizationschemes on a
discretesequence to,ti,.., t.oftimesteps.1nour case, a mod-ifiedtwo-stepsEulerscheme hasbeen utilized
[I2],
in
which particles'positionsat k-thtimesteparederivedfromvelocitiesalready updated atk-thtimestep and corrected with
post-collisionresponse. The modified Euler'smethod turns out to
bestabler thanexplicit
Euler's
and remains computationally cheap, Tbincreasethespeed ofconvergence and manage stiffequations, however,higher-orderimplicitor semiimplicit time
discretizationschemes are currently under implementationfor thefutureversion ofour
ODE
solver,according to ideassug-gested
by
[1
6,
l7,
18,
20].
4.2.
Collision
Management
fo
complete theanalysis of interactlonwith the surround-ingenvironment also theeffect of possiblecolfiisions withob-stacles
has
tobe
taken intoaccount. Co]iisionsoccur, forin-stance, when partsofflexibleobjects
(e,g.,
garments)hitsomerigidobjects
{e.g,,
mannequins) or penetratetowards each other{selfcollisions).
InFigure7,forinstance,collision effects were dominantinsimulating cloth positioningon a tableor asphere, orclothfallon thefloor.
The two aspects of coliision detectionand response are
.-ift)\ffsckese
specialissueofjapamesesocletyforthescieneeofdesign
vol,15-4 mo.60 2eOSNII-Electronic
Mbra
17
Fig.7.Simulationoffabricson rigjdobjects :
(1)
tableclothon atable;
(2)
carpet on asphere ;(3)tableclothfallingon the tloor;(4)
ribbontallingonthefloor,
analysed jnseparate algorithms. Ateach timestep, possible
collisionsagainst rigidobjects
{e.g.
garmentson mannequin$,fabricson rigidframes)or other softparts
(self-collisions)
aredetected,
by
doing
a locaisearch on elementary grid enti-ties(points
orparticles,
edges, triangles).Then,velocities are modified onlyforcoHiding particles.From a computational pointof view, collisjon detection constitutes thebottleneckof simulations, as checks haveto
be
done
in
aIarge
numbe4 between pairsof elementaryenti-ties.Among optimization techniques
proposed
inSection2, aligned axjsbounding
box
(AABB)
hierarchieshavebeen here adopted, as a good compromise between simplicityandef-ficiency.
Accordingtothistechnique, objects are
grouped
hierarchi-caily ciependingon proximityrules,and thedetec±ioniscarriedout byexploring bounding
box
intersect[onsinthehierarchy, Furthermore,since garmentfabricscould self-collide,objectshave been split intosub-regions, each
in
turnassociated to an AABB, Otheroptimization techniquesare currentlyunderimplementationforfutureversions of our model.
When
twosub-regions are foundtobe intersecting,i.e.candidate for
col-lision,checks are finally
detailed
atthe
Ievel
ofgridelementary entities.Collisionsare detectedwhen particlesorpoints
are near asurface or edges are clese tootheredges(i.e.
their dis-tanceisIessthana thresholdvalue},lnthecollision response phase,the velocities ofthe
par-ticlesbelongingtocolliding entities are modified toavoid compenetration. The new velocity vectors are computed
de-18THvrf
v#divekscespecialissueofiapanesesecietyforthescienceofdesign
vol.15-4 no.60 2eOS
pending
on thecoefficjent of restitution, byimposingthe con-servationoflinearmomentum. AIIcollisions together giveraisetoa ljnearsystem
{usually
verysparse) whose $olutiongives
thevelocity variations to beapplied tothe model.
More precisely,lett=
(
r,, r, ,.., h,)be
thesystemconfigura-tionat a certain time step, and sbethenumber ofcollisions
detected,havingtheformof unsatisfiedscalarunilateral con-straints
d,(r,t)<O, i=I,2,".s.
(7)
WedefineD(r,t)=(4{gt},.,,4(gt})and
an(r,t)
J:=
o,'{8)
the sx31V
-sized
Jacobian
matrjxof D with respect tor.First,we impose theconservation oflinearmomentum. Defining
Av asthe
3N
-sized
vectorofvelocity variations, M as the31Vx 3Ndiagonal
matrix whose elements are particlemasses(3
vaiues foreach partjcle),andJ' as thetrasposeof
J,
we canthus
impose
M・Av+J{A=O,
(9)
where A isas
-sized
auxiliaryvector definedasA=-F,
{10)
beingI"the
(unknown)
vector ofsLagrangemultipljers as$oci-ated tothe constraints(7),
The
first
termin(9}
isthe momen-tum of theparticles.Inthesecond term,matrixJ'
describes how thevelocity variations mustbe
distributedamong the threevelocity components of each collidingparticle.Inotherwords,
J'
describes
both
thedirection
ofvelocityvariation vec-torsand the proportionalityamong thejrmoduli. Ifaparticle
is not involvedinany collisionthereis
no velocityvariationandwe can imposedirectly
M・Ay=O.
(ll)
Each collision invoivinga particley'adds a contribution
IA,,
which impliesavariation ofthelinearmomentum.Furthe4we havetoimposetheamount ofvelocityvariation,
dependingon theelasticityof thecol[idingobjects.
A
constraint violationmeans thatD{g
4
changes froma posjtiveor zero va[ue toanegative value, thereforeD(;
4
must benegative.fo
pushthesystem toan acceptable state we must vary D(g
4
insuch away thatD
tr,
4
results positiveor zero. Thus,ifEi$theNII-Electronic Library Service
AD(r,t)=-(1+s)D(r,t).
(1
2)SubstitutingAD
(g
4
=J・Av, we obtain
J・Ay=-(1+a)D(r,t).
(13)
Eqs.
(9)-(1
1}
and(1
3)
constitutea non-singularlinear
system with3N+ sequations and 3AJ+sunknowns, i.e.[Y
el
l"AY]=:[-a+e9・D(r,t>]'
(1
4) whose solution{e,g,
bya conjugate gradientmethod) returnsthevelocity variations thatavoid compenetration. Pre-collision velocity y
{
4}
foreach particleJ'involvedincollisions during thetimestep 4can befinallycorrected with thepost-collision velocity y・(
tn)
+ Ay・A particularsituation fora collis]onsolver arises when there isacontinued contact
(e.g.,
a cloth resting on a rigid object) :our sy$tem has
been
tested
alsoforthesecase$ with a good response bothinprecisionand stability.
Generally,
thesimulationstepis
so smallthatitisnot neces-sary tosearch fornew col]isions at each time step. DuringtheODE
step,the veiocity ofthe fa$testparticleiscomputed. Inthe worst case, thisparticleisclose tocollision with atriangle, the
distance
between
thetwobeing
onlya littlelargerthanthethresholdvalue, say
6.
By
computing thetime4neededforthisparticletocover halfofO,we can
be
surethatno un-detectedcollisions willleadtoa penetrationinthenext 4 sec-onds. Itispossibletomake a more conservative computation of t.,basedon a quarteror lessof dor tobemore aggressive using biggerfractions.Accordingtothisideaithasbeen pos-sible toreduce thecomputational timeforcollision detection.5.
Implementation lssuesThe
particle-based
model describedinSections3 and 4 hasbeenimplementedina system named Softworld2.0,run-ning on Windows, Unix/Linux,SGI-IRIXplatforms.Thesystem iscomposed oftwo fundamentalmodules :
・
2D13D
Modeller,creating theparticle-based
modelassoci-ated to3D
garment
models ;-
3D
Simulator,
genera
±ing
simulationframes
atdiscrete
timesteps4,ti,..,t.,
The
systemis
providedwith aGraphical
User
lnterface
tocreate a unique environment
from
whichboth
the Modellerand the
Simulator
are executed. Figure8
displaystheoveraliFig.8.lable,1Ba$icSoftworld:sarchitecture
,
Datastructure ofthe2D13DModellerlable.2.Mainalgorithm ofthe3D Simulator
Softworld
architecture,lable
1displays
inputand output dataofthe
2D/3D
Modellec
whilelable
2describesthemain algo-rithm inthe3DSimulator.
Paretleiversion.
Some
partsof thesimuiation algorithm al-[owdatato besplitintoindependentsets, suitable forparallel processing.Beowulfclusters, networks ot workstations usingrifify#ffxtsseg
speclalissueefjapanesesocietytorthescienceotdesign
val.15-4 ne.60 2e08
NII-ElectronicMbra
19
MPI
libraries,
have
been used with thatintent.Ibdecreasethecomputational time,in
particulany
thecollisiondetection
algo-rithmhas been
parallelized.
The approach basedonbounding
box hierarchywith sub-region subdivision,
in
fact,
allows man-agement ofeach region as an independentdata$et.When two regions are checked to
find
co[lisionstheirdata
are sent toa separate processorthat
performs
thecheck,Each processorcan work
in
parallelusing only itslocaldata and itsresults(a
listofdata
representing detectedcollisions) can beeasily merged withtheones fromotherprocessors.6.
Applications
Several
testshave
been
done
toeva[uate Softworld2.0,Fjgure7 already showed
frames
of simulated fabricsfallingoverrigidobjects such as atable,a sphere, or thefloor. Actu-ally,forour
design
intents,theseexamples arenot sosignffi-cant because theyhand[esimp[e-shaped one-layered fabric
wjthout any ofthosedesigncomplexities
(seams,
complexborders,
djfferentmaterials,etc,)encountered while designingand manufacturing real cloth. Forour purposes,testson
com-plete
garment
shapes havebeen more interesting,
The
nextSections
show how thesystem hasbeenapplied inthe design processofmen's and women]s garments.Results have
been
obtainedin
thecontext of Europeanand ltalian re-search projects,incooperation with industrialpartners
fromthe clothing industry.Althoughbasedon thesame modelling and simulation algorithms, the two ways of proceedingwhile definingma[e and
female
clothshavebeendifferent,6.1.
WomenCloth
Design
Theactivity on women garmentdesignhasbeencarriedout
byKAEMaRT group
in
theframework
ofMASCOT
Brjte-EuramProject.Here,we havestarted from
3D
dressmodels, directly chosen froma libraryof templatecloth shapes(skirt,
blouse,
etc.), thatare modified by3D surface model[ing toolsbased
on NURBS representation, acting on characteristic lines.
Seams, pences,and other textilecharacterizatjons aredefined on the3D model, fromwhich 2D panelsare extracted with allnecessary intormation
(e,g.,
sewing lines)necessary tothesimulator fortheircorrect locationon virtual mannequins,
The modelling algorithm generatesthepadicle-basedgrid
of each 2D
panel
composing the cloth, associating proper KES datacorresponding tothechosen fabricmateriai, Byus-inginformation
provided
bythemodeller, 2D panelsare semtedand located
properly
on themannequin forthe finalsimulation207ifl)\Mvek#g
speclalissueofjapenesesocietyforthesclemceotdesignvol.15-4 no.60 200S
Fig,9,Particle-basedmodels otwomenis garments placedon
quinsbeforesimulation
Fig,1O,Simulationof women's garments:initialconfiguration and simulated resultofaskirtand tunic
of itsbehaviourinstatic conditions
[36].
Here,
themodelling and sjmu[ation algorithms developedbyKAEMaRT grouphave
been accompanied by a CAD intenfacefor3D
geometricalmodeiling,
provided
bytheUniversityofValenciennes.In
Fig-ure 9,we displayparticle-basedmodeis of a skirt, atop anda woman's dressbeforesimulation. Figure1O presentsthe simulated resultsofaskirtand a linentunic.Detailsabout the modelling and simulation phases,results on testmodels, and
a
performance
analysisarereported in[37].
6.2.
Men
Cloth
Designlaskstormen's
garment
designhave been developed inthe
framework
oftheNationaiproject
IA2000[38],
Differentlyfromfemale
fashion,
theproduction
of men:sgarments
NII-Electronic Library Service
thatonly slightlydifferfromone season and theother,Here, therefore,we havechosen towork bymodification of existing appare[articles.
The
3D
modeller has been implemented byu$ing MayaDeformers(htrp:/C'tvww.afiaswavof-ont.comfezarebroductimaya/1'ndex. shtm4. An interface
{modification
module) iscreated inMayatsowner languageallowing interactivemodification of sizes in
jacket
parts(sleeves,
shoulders, breast,waist, etc.}, with free form deformationdriven
by
characteri$tic lines.An exportmodule automatically generatesafifieinCEX
{Cloth
EXchange} format,containing informationabout themodified 3D mode[, tobeused as inputfor2D pane]extraction. Generationof 2D panelsisdone bymeans of an owner textileCAD
software used at Cornelianiclothing industry.A program,called cexlblVcreates the IVformats
for
the3D
clothmodel and its2D
pan-els.
Softworld's
modellergenerates
theparticie-basedgrids
ofthepanels,thenassembled togetheron thevirtualmannequin thanks ±othe
2D13D
mapping rulesstoredfrom
the3D
mod-eller.Once
theparticle-basedmodelhas
been
createdin
3D,
simulationcan
be
performed.
Figure
1
1
displays
themain modelling $±eps. Figure12 shows severalpoints
ofview ofthe simulated effect of thejacket
leaned
on thebust
ofa mannequinin
staticconditions,The man
jacket
example was computationally expensive, due toa fineparticle
discretization
necessaryfor
guaranteeinggood
accuracy foracomplex model composed ofatleast12panels.Detailsabout the
jacket
model and simulationresu[ts,as well as performance
(for
bothserial and parallelSoftworld:s versions), can be read again in[37],
Parallelsimulations have been executed on thecluster Feronia,40 APIbi-processor nodes, locatedatEnea-Casaccia,Rome.7.
Conclusions
and FutureWorkPhysics-basedmodels show theirpotentialinrealistically
simulating cloth shapes, as the geometricalsubstrate is en-riched with informationderivedfromstructural fabric proper-tiesand
dynamic
laws.
In
thiswork,the
physics-basedcloth shape simulation isderived
from a discreteparticle-basedapproach, with
dynamic
analysisexpressedby
Newton]s
law,
Fig.11.Virtuaidesignofmale garments.Main modelling steps :
(1)
Smoothed geometricalmodel fromdigitalizationofareal
et,
{2)
ModMcation module.(3)
Extractionof2D patterns,(4)
Particlebasedmodel ofeach panel,
(5)
Panelassembly withfinalparticle-basedmodel readytor 3D simulation.
Fig.12.Some views otasimulated jacket
7on)7Mverekg specialissueo"apanesesocietyforthescienceofdesigm vol.15-4 no.60 200S NII-Electronic Mbra
21
yinvolving
internal
and external forcecomputation, constraint management and collisiondetectionand response. Thesys-tem includesspecificalgorithms fortextileoperations such as
sewing, buttonand dartinsertion,etc. and has
been
included
inCAD environments forcloth design,Applicationsand
ex-amples on simulated testcases were shown formen:s and
women's garmentdesign.
Modellingand simu[ation module$ of
Softworld
have beenassessed byend-users participatingto European projects,
in-volving bothacademic and jndus±rialpartners,The modellingmodule has
been
evaluated positivelyby
end-users,because
itdoes not require specific
knowledge
ongeometric
model-lingor professionalskillson apparel
design,
On
theotherhand,simulation re$ults arereliable.A significantaccuracy is reached forsimple- or complexshaped
fabric
models, when designedthrough simple seams, one-layered fabric,e.g.fortunics and nonaccessoried garments composed of Iight fab-ric. Models such as male
jackets
providegood simulationresults as well,although theseturnout tobe computation-al[ymore
hesigent,
due thepresence
of multi-layered fabrics,abundance ofseams between constitutive
panels,
and use ofheaviermaterials.
New featureswillbe added, to expand thesystem
capabili-tiesforvirtual cloth
design
and generateshapes thatareclos-er totheones conceived by stylists.These are,
for
instance, parametricsewing algorithms givingspecial aesthetic effects(collars,
puffedsleves, closefitting sewings,jacket
rollers, etc,),multi-layered fabricmodelling
(such
as liningsand paddings), ormodelljng effectsdue
to
manufacturing proces$es,e.g.,producingpermanent deformations
(as
jtoccursby
ironing,starching, etc.), Alltheseissuescorrespond togeometrical tricksapplied on associated particlemeshes.
The most critical issueisthecalcu]ation time,particulariy duetothecollisiondetectionphase,as possiblecollisions at
each timestep are searched among pairsof elementary
enti-ties
{particies
orpoints,
edges and triangles).New optimization techniquescan beintroducedtoreduce the number ofneces-sarychecks. Higher-orderimplicitand semi-implicitODE
solv-ersfortheassocia ±e Newtonian problem and adaptive time
discretization
criteriaare also currentiy under development toreduce thecomputational time.Besides,furthereffortsareunder way
for
parallelimplementation
ofthecollision detection phase,and othersimulationmodules.Many
questionsstillremain open.Cloth
modelling issome-22T-if{>\ffftksug
speclaiissueetjapamesesecietyforthescienceotdesignvei.15-4 no.60 200e
how a challenge forcomputer graphic$community. At
pres-ent,our group
is
working on thesemodelling issues:
-
geometrical
modelling aspects(e.g.,
grid
generation
forcomplex-shaped objects,
grid
assembly, 2D-3D
ping)i
-
materialmechanjcs issues(e.g,,
fabricmaterial properties,extension of KES data,otherfabricmeasure systems);
-
specifictextileoperators(e,g.,
sewing, multi-layeredrics,small accessory insertion,manufacturing processes),
Concerningthe3D Simulator,arguments come from:
-
dynamic analysis problems(such
as internalforcetationand collision detectionand response);
-
numerical aspeots{e.g.,
stabi[ity and convergence ofODE
solvers,ill-conditioningof algebraic solvers}.
Despitetheseopen issues,simulation results obtained on complete garmentshapes are realistic and encouraging.
They
pointout potentialitiesof physics-basedapproaches in
model-linga complex matter such as cloth,
for
CAD
applicationsin
clothing, home fabricfurnitureand upholstety
industry,
Acknowledgements
Thiswork has
been
carried out intheframework oftheBrite-EuramProjectMASCOZ fundedbythe European
com-mission, and theItalianProjectIA2000, fundedbyMIUR.The
authors are
gratefui
toDrs.A.GalimbertiandG.
Frugoli,
for
contributions with implementationcode, and students from theDjpartimentodilngegnerialndustriale,Universitadi
Parma,developingtheGUIinterface.
References
1ET.
Peirce,"Thegeometry
ofclo±hstructure",Journaloffextile
lnstitute,
Vol.28,]937,pp.45-97.
2
S.
DeJong
and R.Postle,"An energy analysis of wovenfabric
mechanics by means of optimal control theory.Parts1and ll",Journalof lextilelnsMute,Vol.68,1977,
pp.350-369.
3
J,
Weil,
"The synthesis ofCloth
Objects",
Computer
Graphics
(Proc.
SIGGRAPH
'86),
Vdl.20(4},1986, 54.4T Agui,Y Nagao and M.Nakajma,
"An
Expressionod of Cy[indrical
Cloth
Objects-An
Expression
ofFolds
ofa Sleeveusing Computer Graphics",frans,
Soc.
ofNII-Electronic Library Service
1990,
pp,1
095-1097,5 B.K.Hinds,J.McCartney and
G.
Woods,
"Patternvelopment for3D surfaces",
Computer
Aided
Design,Vol,23(8},1991,
pp.583-592.
6 D.ferzopoulos,J.Platt,A.Barrand K.Fleischer,
cally Deformable Models",Computer
Graphics
(Proc.
SIGGRAPH ・87}, Vol,21
{4),
1987,
pp,205-21
4.
7 J.W.Eischen,S,Deng,and 11G.Clapp,"Finite ElementModelingand Controlof FIexibleFabricParts",Computer
Graphicsinlextilesand Apparel,IEEEComputer
ics,Vol.l6(5),l996, pp,71
-80.
8 M.Aono,'`A WrinklePropagationModel forCloth",Proc.
Computer
Graphicslnternational,1990, pp.95-105.9 L.Li,M. Damoran and R,K.Gay, [`Aerodynamic force
models foranimating cloth motion inairflow'',The Visual
Computer,
Vol.12,
1
996,
pp.84-104.1
O
C,R.
Feynman,
"Modelling
theAppearanceof Cloth",MScthesis,
Dept,
ofElectrical
Engineeringand Computerence, Massachussetts]n$tituteoffechnology,US,1986,
11
D.H.
House
andD,E.
Breen,
"Representationof Woven
Fabrics",
SIGGRAPH
"98
Course
Notes,N.31,l998. I2 ×,
Provot,"CoHision and self-collisionhandlinginclothmodel dedicatedto
design
garments",Proc.Graphicsterface'97,1997,pp.177-189,
13
R Volinoand N,Magnenat
Thalmann,
"CoHislon andCo[lisionDetectlon:Robustand
Efficient
fechniques
for
Highly
Deformable
Surfacesi',
EurographicsWorkshop onAnimationand Simulation,1995.
14 R Volinoand N,Magnenat-Thalmann, "Virtual Clothing,
Theoryand Practise",Springer-Verlag,Heidelberg,
many, 2000.
15 M.Desbrun,R Schroederand A.H.Barr,`tlnteractive
mation of structured deformableobjects", Graphicsjntee
face,1999,pp.1-8.
16 D.Baraffand A.Witkin,"Large Steps inCloth tion",Computer Graphics
{Proc.
SIGGRAPH '98), 1998,pp,43-54.
I7 M Kang,J.Choiand C.RH.Cho,"Fast and stable tionof cloth with an approximated implicitmethod:',
Proc.
Computer Graphicslnternational
(CGI
'OO),2000.
18 M.Hauthand O.EtzmuB,"A HighPerformanceSolver
fortheAnimationof Deformable
Objects
usingAdvanced
merical Methods",Proc,Eurographics
'Ol
, 200I.
19 O.EtzmuB,J.Grossand W.
StraBer,
`'Deriving
a ParticleSystem
from
Continuum
Mechanics
fortheAnimationofDeformableObjects",IEEE
fransactions
on Visualizationand Computer Graphics,
2002.
20 K-J.Choiand H-S. Ko,"Stable
but
ResponsiveC)oth",Computer Graphics
(Proc,
SIGGRAPH:02),
2002,
21 l.J.Rudomin, "Simulating clothusing a mixed
physicalmethod", Ph.DThesis,Dept.of
Computer
ence, Universityof Pennsylvania,US,
1990,
22 T.L Kuniiand H,Godota, "Modelling and animation of garment wrinkle formationprocesses",Proc.
Computer
Animation,
Springer-Verlag,
Berlin,1990,pp,1
31
-1
47,23
F.Taiilefer,"Mixed
modeling", Proc.Compugraphics,
1991,pp.467-478.
24
N.isopelas,"Animating
thecrump[ing behaviourofments'',
Proc,
2nd
Eurographics
Workshop on Animationand
Simulation,
Blackwekk
UK, 1991 , pp.11-24.
25 S.G,Dhande, RNLM, Rao,
S.
lavakkoiiand C.L Moore.,"Geometric modeling ofdrapedfabricsurface", Proc.IFIP
lnt.ConferenceOn Computer
Graphics,
North-Ho]land, Amsterdam, 1993, pp.349-356.26 M.Oshitaand A.Makinouchi,
"Real-Time
Cloth
tionmtithSparseParticles",SIGGRAPH
:OlSketches
AndApplications,2001.
27 Y,M. Kang,J.H.
Choi,
H.
G,
Cho
andD.H.
Lee,"Anficientanimation ofwrinkled clothwith approximate
plicitintegration",The VisualComputer,Vol.17(3),2001, pp,147-157,
28 1.J,Rudomin and J.
Castilio,
"Real-Time
Clothing: etry and Physics",WSCG
Posters
'02, 2002.29 Eds.IEEE
Computer
Graphicsand Applications,specialissueon
Computer
graphicsinfextilesand Apparel,IEEE ComputerGraphics,
Vbl.16(5),1996,30
Eds.
D,H,
House
and D.E.Breen,"Clothmodelling and
Animationi',
A KPetersLtd.,Natick,MA(US),
2000.31 S.Kawabata,
``The
Standardizationand Analysisof HandEva[uation"',ThefextileMachinerySocietyofJapan,
ka,
1980,32
G.
Frugoli,
A.Galimberti,
C,Rizziand M.Bordegoni,crete element approach fornon-rigid materlal modelling",
Robot
Manipulationof DeformableObjects,Eds.D,
rich &H,Worn,Advanced ManufacturingSeries,Springer,
2000, pp.29-42.
33 K.Fleisher,A.Witkinand A.Barr,"Energy constraints on parametrizedmodeis",
Computer
Graphics,
Vo[.21
,7if{y#anvekEt
specialissueetjapanesesocietyferthescienceatdesign vel.15-4 no.60 2008
NII-Electronic Mbrar
l987,
pp.225-232.
34
J.Plattand A.H. Barr,"ConstraintMethods
forible
Models",
ComputerGraphjcs,
Vol,22(4),1988,
pp.279-288.
35
A.Witkin,"lnteractiveDynamic$",
Computer Graphics,Voi.24,Proc.
Symposium
on3D
lnteractiveGraphics,1990,pp.l1-20
36 r Bonte,
A.
Galimberti
andC.
Rizzi,[`A 3Dgraphic ronment to design
garments",
From GeometricModeljng toShape
Modeling,
Eds.
U.Cugini& M.Wozny, Kluwer Academic Publishers,2000,
pp,137-1
50{also
Proc,7thIFIP
Workshop
on GeometricModeling :Fundamentalsand Appljcatjons
-
GE07,
Parma,2-4 ottobre 2000}, 37 A.Galimberti,
"Sviluppo djun sistema disimulazione per
]a
progettazjonedicapi diabbigliamento", Ph.D.Thesis,
Dipartlmentodilngegnerialndustriale,Universitadi
ma,
ltaly,
2001,
38
U.Cugini
and O.Rizzi,"3D Designand Simulationofmen
garments",