10.1.1.61.7862

8/11/2019 10.1.1.61.7862

http://slidepdf.com/reader/full/1011617862 1/10

I t e r a t i v e M e s h G e n e r a t i o n

f o r F E - C o m p u t a t i o n s o n F r e e F o r m S u r f a c e s

L e i f K o b b e l t

T o r s t e n H e s s e

y

H a r t m u t P r a u t z s c h

z

K a r l S c h w e i z e r h o f

y

M a y 3 0 , 1 9 9 6

A b s t r a c t

W e p r e s e n t a n i n t e r p o l a t o r y s u b d i v i s i o n s c h e m e

t o g e n e r a t e a d a p t i e l y r e n e d q u a d r i l a t -

e r a l m e s h e s w h i c h a p p r o x i m a t e a s m o o t h s u r f a c e

o f a r b i t r a r y t o p o l o g y . T h e d e s c r i b e d m e t h o d s i g -

n i c a n t l y d i e r s f r o m c l a s s i c a l m e s h g e n e r a t i o n

t e c h n i q u e s b a s e d o n s p l i n e s u r f a c e s o r i m p l i c i t

r e p r e s e n t a t i o n s s i n c e n o e x p l i c i t d e s c r i p t i o n o f

t h e l i m i t s u r f a c e i s u s e d . I n s t e a d , s i m p l e a n e

c o m b i n a t i o n s a r e a p p l i e d t o c o m p u t e n e w v e r -

t i c e s i f a f a c e o f t h e n e t i s s p l i t . T h e s e r u l e s

a r e d e s i g n e d t o g u a r a n t e e a s y m p t o t i c s m o o t h -

n e s s , i . e . , t h e s e q u e n c e o f r e n e d n e t s c o n v e r g e s

t o a s m o o t h l i m i t s u r f a c e . S u b d i v i s i o n t e c h n i q u e s

a r e u s e f u l m a i n l y i n a p p l i c a t i o n s w h e r e a g i v e n

q u a d r i l a t e r a l n e t i s a c o a r s e a p p r o x i m a t i o n o f a

s u r f a c e a n d p o i n t s o n a r e n e d g r i d h a v e t o b e

e s t i m a t e d . T o e v a l u a t e o u r a p p r o a c h , w e s h o w

e x a m p l e s f o r F E - c o m p u t a t i o n s o n s u r f a c e s g e n -

e r a t e d b y t h i s a l g o r i t h m .

I n t r o d u c t i o n

T o p e r f o r m n u m e r i c a l c o m p u t a t i o n s o n s u r f a c e s ,

o n e u s u a l l y s t a r t s b y d i s c r e t i z i n g t h e p r o b l e m ,

i . e . , a s m o o t h m a n i f o l d i s a p p r o x i m a t e d b y a t r i -

a n g u l a r o r q u a d r i l a t e r a l n e t . T h e f a c e s o f s u c h

n e t s c a n b e l l e d b y l i n e a r o r b i l i n e a r e l e m e n t s

r e s p e c t i v e l y . T h e s i z e o f t h e f a c e s c o n t r o l s t h e

C o m p u t e r S c i e n c e s D e p a r t m e n t , U n i v e r s i t y o f W i s -

c o n s i n | M a d i s o n , 1 2 1 0 W e s t D a y t o n S t r e e t , M a d i s o n ,

W I 5 3 7 0 6 - 1 6 8 5 , U S A , k o b b e l t @ c s . w i s c . e d u

y

U n i v e r s i t a t K a r l s r u h e , I n s t i t u t f u r M e c h a n i k , 7 6 1 2 8

K a r l s r u h e , G e r m a n y

z

U n i v e r s i t a t K a r l s r u h e , I n s t i t u t f u r B e t r i e b s - u n d D i -

a l o g s y s t e m e , 7 6 1 2 8 K a r l s r u h e , G e r m a n y

q u a l i t y o f t h e a p p r o x i m a t i o n a n d t h u s t h e r e -

l i a b i l i t y o f t h e c o m p u t a t i o n . T h e t r a d e - o b e -

t w e e n c o m p u t a t i o n a l c o m p l e x i t y a n d e r r o r t o l e r -

a n c e m a k e s i t n e c e s s a r y t o a d a p t t h e g r i d s i z e o f

t h e d i s c r e t i z a t i o n d u r i n g t h e c o m p u t a t i o n . U s u -

a l l y o n e s t a r t s w i t h a r a t h e r c o a r s e a p p r o x i m a -

t i o n a n d t h e n i t e r a t i v e l y r e n e s t h e n e t u n t i l t h e

r e q u i r e d p r e c i s i o n i s o b t a i n e d .

I n m a n y a p p l i c a t i o n s , t h e u n d e r l y i n g s u r f a c e i t -

s e l f i s g i v e n b y a s e t o f d a t a p o i n t s , e . g . , s c a n n e d

f r o m a r e a l o b j e c t . I n a p r e - p r o c e s s i n g s t e p , t h e

t o p o l o g y h a s t o b e d e t e r m i n e d , i . e . , t h e p o i n t s

h a v e t o b e c o n n e c t e d t o b u i l d a n i n i t i a l m e s h

( r e v e r s e e n g i n e e r i n g ) . T h e n o n e c a n a p p l y a p -

p r o x i m a t i o n o r i n t e r p o l a t i o n t e c h n i q u e s t o c o n -

s t r u c t a s m o o t h s u r f a c e t t i n g t h e d a t a . H o w -

e v e r , t o p e r f o r m n i t e e l e m e n t c o m p u t a t i o n s ,

t h i s s u r f a c e a g a i n h a s t o b e s a m p l e d o n s o m e

( n e r ) g r i d t h u s r e t u r n i n g t o a d i s c r e t e r e p r e -

s e n t a t i o n ( c f . F i g . 1 ) . T h e s t a n d a r d m e t h o d s d e -

s c r i b e s u c h c o n s t r u c t i o n s u s u a l l y b y m e a n s o f

s p l i n e s u r f a c e s H L 9 3 ] o r i m p l i c i t r e p r e s e n t a t i o n s

B C X 9 4 a , B C X 9 4 b ] .

I f w e a r e i n t e r e s t e d o n l y i n t h e m e s h e s t h e m -

s e l v e s , i t m a k e s s e n s e t o a v o i d t h e e x p l i c i t c o n -

s t r u c t i o n o f a n i n t e r p o l a t i n g s u r f a c e a n d t o d e -

r i v e t h e r e n e d n e t d i r e c t l y f r o m t h e g i v e n o n e .

T h i s t e c h n i q u e i s c a l l e d i n t e r p o l a t o r y s u b d i v i -

s i o n o r r e n e m e n t . A r e n e m e n t s c h e m e i s c o m -

p l e t e l y d e n e d b y a s e t o f r u l e s h o w t o c o m -

p u t e n e w v e r t i c e s . T h e s e r u l e s a r e u s u a l l y s i m p l e

a n e c o m b i n a t i o n s o f p o i n t s f r o m t h e u n r e n e d

n e t . I f t h e c o e c i e n t s o f t h e s e c o m b i n a t i o n s s a t -

i s f y c e r t a i n c o n d i t i o n s R e i 9 5 , P r a 9 5 ] t h e n t h e

i t e r a t i v e r e n e m e n t g e n e r a t e s a s e q u e n c e o f n e r

a n d n e r n e t s w h i c h c o n v e r g e t o a t a n g e n t - p l a n e

c o n t i n u o u s s u r f a c e .

T h e j u s t i c a t i o n f o r t h i s a p p r o a c h i s t h a t o n e

1

8/11/2019 10.1.1.61.7862


Interpolation

Subdivision

Discretization

F i g u r e 1 : G e n e r a t i o n o f i n t e r p o l a t i n g s u r f a c e s

a v o i d s t h e m a t h e m a t i c a l l y i n v o l v e d c o n s t r u c t i o n

o f s m o o t h p a r a m e t r i c s u r f a c e s b u t u s e s t h e i n i -

t i a l m e s h o b t a i n e d d u r i n g t h e r e v e r s e e n g i n e e r -

i n g p r o c e s s d i r e c t l y o r a f t e r s o m e r e g u l a r i z a t i o n .

I f a r e n e m e n t o f t h e m e s h i s n e e d e d t o a c h i e v e

a p r e s c r i b e d t o l e r a n c e , t h e n e t i s r e n e d b y s i m -

p l e x e d r u l e s t o c o m p u t e t h e l o c a t i o n s o f t h e

n e w p o i n t s .

I n t h e l i t e r a t u r e t h e t e r m s u b d i v i s i o n u s u a l l y

r e f e r s t o a m o r e g e n e r a l c l a s s o f a l g o r i t h m s

w h i c h g e n e r a t e s e q u e n c e s o f p o l y h e d r a l n e t s

D S 7 8 , C C 7 8 , L o o 8 7 , L o o 9 4 , K o b 9 6 ] . M o s t o f

t h e s e s c h e m e s a r e d e r i v e d f r o m k n o t - i n s e r t i o n a l -

g o r i t h m s o f p i e c e w i s e p o l y n o m i a l t e n s o r - p r o d u c t

s p l i n e s . I n t h i s p a p e r w e r e s t r i c t o u r s e l v e s t o t h e

s p e c i a l c a s e o f i n t e r p o l a t o r y s u b d i v i s i o n w h e r e

t h e n e t r e n e m e n t i s a c h i e v e d b y i n s e r t i n g n e w

v e r t i c e s a n d k e e p i n g t h e e x i s t i n g v e r t i c e s u n -

c h a n g e d .

A f a i r l y w e l l - k n o w n a l g o r i t h m o f t h i s s u b c l a s s

i s t h e b u t t e r y - a l g o r i t h m o f D G L 9 0 , Z S S 9 6 ] .

H o w e v e r , t h i s a l g o r i t h m o p e r a t e s o n t r i a n -

g u l a r n e t s o n l y . F o r o u r a p p l i c a t i o n o f F E -

c o m p u t a t i o n s o n s u r f a c e s , q u a d r i l a t e r a l ( b i l i n -

e a r ) e l e m e n t s t u r n o u t t o b e m o r e a p p r o p r i a t e

s i n c e , e s p e c i a l l y f o r v a r y i n g g r a d i e n t s , t h e g e o m -

e t r y i s a p p r o x i m a t e d b e t t e r w i t h t h e s a m e n u m -

b e r o f d e g r e e s o f f r e e d o m a n d l o c k i n g p r o b l e m s

a r e o v e r c o m e B D 8 6 , S F 8 9 a , S F 8 9 b ] .

I n t h e f o l l o w i n g , w e r s t d e s c r i b e a n i n t e r p o l a -

t o r y s u b d i v i s i o n s c h e m e f o r c l o s e d q u a d r i l a t e r a l

n e t s w h i c h , i n t h e l i m i t , g e n e r a t e s s m o o t h ( C

1

)

s u r f a c e s o f a r b i t r a r y t o p o l o g i c a l t y p e . A s l i g h t

m o d i c a t i o n o f t h e s c h e m e a l l o w s u s t o r e n e

o p e n n e t s i n w h i c h c a s e s m o o t h b o u n d a r y c u r v e s

i n t e r p o l a t e t h e g i v e n b o u n d a r y p o l y g o n .

S i n c e t h e c o m p u t a t i o n o f n e w v e r t i c e s o n l y r e -

q u i r e s l o c a l i n f o r m a t i o n f r o m t h e p r e c e e d i n g r e -

n e m e n t l e v e l , l o c a l r e n e m e n t s t r a t e g i e s c a n

b e a p p l i e d . H o w e v e r , t o k e e p t h e m e s h C

0

-

c o n s i s t e n t t h e l o c a l l y r e n e d n e t h a s t o r e m a i n

b a l a n c e d , i . e . , t h e r e n e m e n t l e v e l o f t w o f a c e s

h a v i n g a t l e a s t o n e v e r t e x i n c o m m o n , m u s t

n o t d i e r b y m o r e t h a n 1 . W e r e f e r t o w e l l -

k n o w n t e c h n i q u e s t o h a n d l e b a l a n c e d q u a d r i l a t -

e r a l m e s h e s P S B 9 0 ] .

W e d e r i v e t h e a p p r o x i m a t i o n o r d e r , i . e . , t h e f a c -

t o r b y w h i c h t h e a p p r o x i m a t i o n t o a r e a l s u r -

f a c e S i s i m p r o v e d i f w e s t a r t t h e i t e r a t i v e r e n e -

m e n t o n a n e r i n i t i a l m e s h ( w i t h p o i n t s l y i n g

e x a c t l y o n S ) . F u r t h e r , w e c o m p a r e t h e r e s u l t s

o f F E - c o m p u t a t i o n s o n m e s h e s g e n e r a t e d b y o u r

r e n e m e n t s c h e m e t o r e s u l t s o b t a i n e d b y o t h e r

m e s h r e n e m e n t t e c h n i q u e s . T h e e x p e r i m e n t a l

r e s u l t s d e m o n s t r a t e t h e p r a c t i c a l u s e f u l n e s s o f

t h e s c h e m e .

S u b d i v i s i o n R u l e s

T h e f o l l o w i n g d e s c r i p t i o n i s i n t e n d e d t o b e a s

c o n c i s e a s p o s s i b l e . A m o r e d e t a i l e d p r e s e n t a t i o n

c a n b e f o u n d i n K o b 9 5 , K o b 9 6 ] .

U n i f o r m r e n e m e n t o f a q u a d r i l a t e r a l m e s h

m e a n s t o c o m p u t e o n e n e w v e r t e x f o r e a c h e d g e

a n d o n e f o r e a c h f a c e . W e c a l l t h e s e n e w v e r -

t i c e s e d g e - v e r t i c e s ( E ) a n d f a c e - v e r t i c e s ( F ) r e -

s p e c t i v e l y ( c f . F i g . 2 ) . F o r b o t h t y p e s o f n e w

v e r t i c e s , w e d e n e a s e p a r a t e r e n e m e n t r u l e t o

d e t e r m i n e t h e i r p o s i t i o n .

W e l a b e l t h e v e r t i c e s a s s h o w n i n F i g . 3 . L e t e b e

t h e v e r t e x a s s o c i a t e d w i t h t h e e d g e p

0

q

0

t h e n

e : =

4 + w

8

( p

0

+ q

0

) ;

w

8

( p

+ q

) ( 1 )

w h e r e p

, q

a n d w a r e d e n e d a s f o l l o w s : T h e

v a l u e w i s a t e n s i o n p a r a m e t e r i n u e n c i n g h o w

c l o s e l y t h e r e s u l t i n g s u r f a c e f o l l o w s t h e i n i t i a l

n e t . F o r w = 0 , e h a p p e n s t o b e t h e m i d p o i n t

o f t h e e d g e p

0

q

0

. I n o r d e r t o g u a r a n t e e t h e

s m o o t h n e s s o f t h e l i m i t i n g s u r f a c e t h e v a l u e w

s h o u l d b e i n t h e i n t e r v a l ( 0

p

5 ; 1 ] . B e s t r e s u l t s

a r e o b t a i n e d b y c h o o s i n g w c l o s e t o t h e s t a n d a r d

2

8/11/2019 10.1.1.61.7862


E

E

E

E F

F i g u r e 2 : S p l i t t i n g a q u a d r i l a t e r a l f a c e b y E -

v e r t i c e s a n d F - v e r t i c e s .

v a l u e w = 1 D G L 8 7 ] . T h e p o i n t q

i s g i v e n b y

q

: =

4

n

n ; 1

X

i = 0

p

i

; ( p

n ; 1

+ p

0

+ p

1

)

+

w

4 + w

( k

n ; 2

+ k

n ; 1

+ k

0

+ k

1

)

;

4 w

( 4 + w ) n

n ; 1

X

i = 0

k

i

!

( 2 )

a n d t h e p o i n t p

i s d e n e d a n a l o g o u s l y . H e r e , n

d e n o t e s t h e v a l e n c e o f t h e v e r t e x q

0

. I n t h e r e g u -

l a r c a s e , n = 4 , t h e f o r m u l a r e d u c e s t o q

= p

2

.

S i n c e t h e n e w v e r t i c e s w h i c h a r e i n s e r t e d d u r i n g

t h e r e n e m e n t p r o c e s s a l w a y s h a v e v a l e n c e 4 , t h e

g e n e r a l f o r m u l a ( 2 ) i s o n l y n e e d e d t o s u b d i v i d e

t h e e d g e s d i r e c t l y e m a n a t i n g f r o m a s i n g u l a r v e r -

t e x ( v a l e n c e 6= 4 ) . I n t h e r e g u l a r r e g i o n s w h i c h

c o v e r m o s t o f t h e n e t , n o f u r t h e r c o m p u t a t i o n i s

n e c e s s a r y ( c f . F i g . 4 ) .

A f a c e - p o i n t f i s c o m p u t e d b y a p p l y i n g ( 1 ) t o

f o u r s u c c e e d i n g e d g e - p o i n t s :

f : =

4 + w

8

( b + c ) ;

w

8

( a + d ) : ( 3 )

T h e r e a r e t w o p o s s i b i l i t i e s f o r c h o o s i n g t h i s s e t

o f e d g e - p o i n t s a , b , c , a n d d . H o w e v e r , r u l e ( 2 )

i s e s p e c i a l l y d e s i g n e d t o m a k e b o t h a l t e r n a t i v e s

e q u i v a l e n t . H e n c e , i t d o e s n o t m a t t e r i n w h i c h

d i r e c t i o n r u l e ( 3 ) i s a p p l i e d s i n c e b o t h l e a d t o

t h e s a m e r e s u l t .

I f a r e n e m e n t s c h e m e i s t o b e u s e d i n p r a c -

t i c a l m o d e l i n g o r r e c o n s t r u c t i o n a p p l i c a t i o n s , i t

m u s t p r o v i d e f e a t u r e s t h a t a l l o w t h e d e n i t i o n o f

c r e a s e s a n d c u s p s H o p 9 4 ] . T h e s e r e q u i r e m e n t s

p

2

n−2

0k 1

1k

2

k

n−2k

n−1n−1k

q

a

bcd f

e

p

p

p

0p 0

F i g u r e 3 : L a b e l i n g t h e v e r t i c e s i n t h e v i c i n i t y o f

t h e e d g e p

0

q

0

.

c a n b e s a t i s e d i f t h e s c h e m e i n c l u d e s s p e c i a l

r u l e s f o r t h e r e n e m e n t o f o p e n n e t s w h i c h y i e l d

w e l l - b e h a v e d b o u n d a r y c u r v e s t h a t i n t e r p o l a t e

t h e b o u n d a r y p o l y g o n s o f t h e g i v e n n e t . H a v -

i n g s u c h a s c h e m e , c r e a s e s c a n b e m o d e l e d b y

j o i n i n g t w o s e p a r a t e s u b d i v i s i o n s u r f a c e s a l o n g a

c o m m o n b o u n d a r y c u r v e a n d c u s p s r e s u l t f r o m

a t o p o l o g i c a l h o l e i n t h e i n i t i a l n e t w h i c h g e o -

m e t r i c a l l y s h r i n k s t o a s i n g l e p o i n t ( c f . F i g . 5 ) .

I n o r d e r t o r e n e o p e n n e t s , w e h a v e t o d e -

n e a d d i t i o n a l r u l e s f o r r e n i n g e d g e s o n o r

n e a r t h e b o u n d a r y . C

0

- j o i n s o f s u b d i v i s i o n s u r -

f a c e s w h o s e i n i t i a l n e t s h a v e a c o m m o n b o u n d -

a r y p o l y g o n , r e q u i r e t h a t t h e r u l e s f o r t h e s u b -

d i v i s i o n o f b o u n d a r y e d g e s d o e s n o t d e p e n d o n

i n n e r v e r t i c e s . H e n c e , w e c o m p u t e n e w e d g e v e r -

t i c e s o n t h e b o u n d a r y b y a p p l y i n g r u l e ( 3 ) w h e r e

i n t h i s c a s e a , b , c a n d d r e p r e s e n t f o u r n e i g h -

b o r i n g b o u n d a r y v e r t i c e s .

I n c a s e s w h e r e t h r e e o r m o r e s u b d i v i s i o n s u r -

f a c e s m e e t a t a c o m m o n p o i n t c , i t i s n e c e s s a r y

t o a l l o w p i e c e w i s e s m o o t h b o u n d a r y c u r v e s . L e t

a , b , c , a n d d b e f o u r n e i g h b o r i n g b o u n d a r y

v e r t i c e s o f o n e o f t h e n e t s . S i n c e c b e l o n g s t o

m o r e t h a n t h r e e a d j a c e n t s u r f a c e s , i t h a s t o b e

a b r e a k p o i n t i n t h e p i e c e w i s e s m o o t h b o u n d a r y

c u r v e s . W e d e n e t h e c o r r e s p o n d i n g r e n e m e n t

r u l e b y r e p l a c i n g d b y d

: = 2 c ; b i n ( 3 ) .

W e l l - d e n e d r u l e s f o r t h e s u b d i v i s i o n o f i n n e r

e d g e s e m a n a t i n g f r o m a b o u n d a r y v e r t e x c a n b e

o b t a i n e d b y l i n e a r l y e x t r a p o l a t i n g t h e b o u n d a r y

f a c e s . T h i s g i v e s a n a d d i t i o n a l l a y e r o f f a c e s a t

t h e b o u n d a r i e s o f t h e n e t , m a k i n g f o r m e r b o u n d -

3

8/11/2019 10.1.1.61.7862


F i g u r e 5 : M o d e l i n g s h a r p f e a t u r e s ( p i e c e w i s e s m o o t h b o u n d a r y , c r e a s e , c u s p )

F i g u r e 4 : S i n g u l a r i t i e s i n t h e r e n e d n e t .

a r y v e r t i c e s t o i n n e r v e r t i c e s s u c h t h a t t h e a b o v e

r u l e s f o r i n n e r e d g e s c a n b e a p p l i e d . M o r e p r e -

c i s e l y , i f q

1

: : : q

n

a r e t h e i n n e r v e r t i c e s c o n -

n e c t e d t o t h e b o u n d a r y v e r t e x p t h e n w e c a n

e x t r a p o l a t e b y p

: = 2 p ;

1

n

P

i

q

i

.

I n t h e c a s e o f a b o u n d a r y f a c e

f p q r s

g w h e r e

t h e b o u n d a r y v e r t e x p i s n o t c o n n e c t e d t o a n y

i n n e r v e r t e x w e h a v e t o d i s t i n g u i s h t w o c o n s t e l -

l a t i o n s . I f p i s a b r e a k p o i n t i n t h e p i e c e w i s e

s m o o t h b o u n d a r y c u r v e , w e a d d t h e e x t r a p o -

l a t e d f a c e f p q

r

s

g w i t h q

: = 2 p ; q ,

s

: = 2 p ; s a n d r

= q

+ s

; p . I f p i s

n o b r e a k p o i n t , w e a d d t h e v e r t e x p

: = 2 p ; r .

C o n n e c t i n g t h e e x t r a p o l a t e d v e r t i c e s c o r r e -

s p o n d i n g t o n e i g h b o r i n g b o u n d a r y v e r t i c e s g i v e s

t h e a d d i t i o n a l l a y e r o f e x t r a p o l a t e d f a c e s . A f t e r

t h e r e n e m e n t , t h e s e f a c e s h a v e t o b e d e l e t e d .

C o n v e r g e n c e A n a l y s i s

T h e r e n e m e n t s c h e m e d e n e d i n t h e l a s t s e c -

t i o n b e l o n g s t o t h e c l a s s o f s t a t i o n a r y s u b d i -

v i s i o n s c h e m e s D y n 9 1 ] s i n c e t h e s a m e f o r m u -

l a s a r e u s e d o n e v e r y r e n e m e n t l e v e l . I n t h e

v i c i n i t y o f a n y v e r t e x , s u c h s c h e m e s c a n b e d e -

s c r i b e d i n t e r m s o f a m a t r i x w h i c h m a p s v e r t i c e s

f r o m t h e m t h l e v e l t o t h e ( m + 1 ) s t l e v e l . I n

R e i 9 5 , P r a 9 5 ] s u c i e n t c o n d i t i o n s o n t h e e i g e n

s t r u c t u r e o f t h i s m a t r i x a r e d e r i v e d w h i c h g u a r -

a n t e e t h e c o n v e r g e n c e o f t h e s e q u e n c e o f n e t s t o

a t a n g e n t - p l a n e c o n t i n u o u s l i m i t i n g s u r f a c e .

T h e s e c o n d i t i o n a r e t h a t t h e l e a d i n g e i g e n v a l -

u e s h a v e t o b e

1

= 1 a n d

2

=

3

< 1 a n d

t h e e i g e n v e c t o r s c o r r e s p o n d i n g t o

2

a n d

3

a l -

l o w t h e d e n i t i o n o f a r e g u l a r p a r a m e r i z a t i o n

w i t h r e s p e c t t o w h i c h C

1

- c o n t i n u i t y i s a c h i e v e d

R e i 9 5 ] . W e o m i t t h e d e t a i l s o f t h e s t r a i g h t f o r -

w a r d v e r i c a t i o n o f t h e s e c r i t e r i a a n d r e f e r t o

K o b 9 6 ] .

A p p r o x i m a t i o n O r d e r

T h e i n t e r p o l a t i o n s c h e m e p r e s e n t e d i n t h i s p a p e r

i s d e s i g n e d t o g e n e r a t e f r e e f o r m s u r f a c e s . S i n c e

w e d o n o t e x p l o i t a n y m e t a - k n o w l e g d e a b o u t

t h e g e o m e t r y b u t o n l y t h e i n f o r m a t i o n g i v e n b y

t h e i n i t i a l n e t , w e c a n n o t e x p e c t t o e x a c t l y r e -

p r o d u c e b a s i c s h a p e s l i k e s p h e r e s o r c y l i n d e r s .

H o w e v e r , t h i s l a c k i s t o l e r a b l e i f w e c a n g u a r -

a n t e e t h a t d e n s e r i n f o r m a t i o n a b o u t t h e c o n s i d -

e r e d o b j e c t , i . e . , a n e r i n i t i a l n e t w i t h c o r r e c t

v e r t e x p o s i t i o n s ( o n t h e s u r f a c e t o b e a p p r o x i -

m a t e d ) , s i g n i c a n t l y i m p r o v e s t h e q u a l i t y o f t h e

a p p r o x i m a t i o n .

F r o m D G L 8 7 ] w e k n o w t h a t , o n r e g u l a r t e n s o r

p r o d u c t m e s h e s , t h e a b o v e p r e s e n t e d i n t e r p o l a -

t i o n s c h e m e r e p r o d u c e s t h e p o l y n o m i a l s 1 x y

a n d x y . F u r t h e r , i f w e r e s t r i c t t h e t e n s i o n p a -

r a m e t e r t o w = 1 , t h e s c h e m e e v e n r e p r o d u c e s

4

8/11/2019 10.1.1.61.7862


a l l p o l y n o m i a l s u p t o b i - d e g r e e 3 . F o r i r r e g u -

l a r m e s h e s , w e h a v e t o a p p l y t h e r u l e ( 2 ) i n

i t s f u l l g e n e r a l i t y . S i n c e t h e s c h e m e g e n e r a t e s

C

1

- s u r f a c e s i n t h e l i m i t , w e c a n u s e t h e s p e c i a l

p a r a m e t r i z a t i o n i n d u c e d b y t h e e i g e n s t r u c t u r e

o f t h e s u b d i v i s i o n m a t r i x w i t h r e s p e c t t o w h i c h

t h e r e n e m e n t s c h e m e h a s a t l e a s t l i n e a r p r e c i -

s i o n .

K e e p i n g i n m i n d t h a t t h e p o s i t i o n o f t h e n e w

p o i n t s o n l y d e p e n d s o n a x e d n u m b e r o f n e i g h -

b o r i n g v e r t i c e s , i t i s o b v i o u s t h a t f o r s u c i e n t l y

s m o o t h s u r f a c e s , t h e a p p r o x i m a t i o n o r d e r d i -

r e c t l y f o l l o w s f o r m t h e p o l y n o m i a l p r e c i s i o n o f

t h e s c h e m e . T h i s c a n b e s e e n b y l o o k i n g a t a

l o c a l T a y l o r e x p a n s i o n o f t h e s u r f a c e . H e n c e , w e

o b t a i n a n a p p r o x i m a t i o n o r d e r o f O ( h

2

) w h i c h

i n c r e a s e s t o O ( h

4

) o v e r r e g u l a r r e g i o n s o f t h e n e t

i f t h e t e n s i o n p a r a m e t e r i s s e t t o w = 1 . H e r e ,

h d e n o t e s a n y u n i f o r m m e a s u r e o f t h e g r i d s i z e ,

e . g . , a v e r a g e o r m a x i m a l l e n g t h o f e d g e s i n t h e

i n i t i a l m e s h . I n o r d e r t o a c h i e v e t h e b e s t a p p r o x -

i m a t i o n , l e t w : = 1 f o r t h e r e s t o f t h i s s e c t i o n .

O b v i o u s l y , t h e O ( h

2

) - t e r m i n t h e v i c i n i t y o f e x -

t r a o r d i n a r y v e r t i c e s a s y m p t o t i c a l l y d o m i n a t e s

t h e e r r o r f u n c t i o n a l s u c h t h a t i n g e n e r a l , w e o n l y

g e t q u a d r a t i c c o n v e r g e n c e . H o w e v e r , i f w e c a n

e n s u r e t h a t t h e n u m b e r o f s i n g u l a r i t i e s i n t h e

i n i t i a l n e t d o e s i n c r e a s e s i g n i c a n t l y w i t h t h e

d e n s i t y o f t h e i n p u t d a t a t h e n w e a c h i e v e t h e

f u l l O ( h

4

) c o n v e r g e n c e . T h e a s s u m p t i o n t h a t t h e

n u m b e r o f s i n g u l a r v e r t i c e s i n a m e s h P

h

a p p r o x -

i m a t i n g a s u r f a c e S c a n b e b o u n d e d b y a c o n -

s t a n t i s q u i t e r e a l i s t i c s i n c e t h e a m o u n t o f s i n -

g u l a r i t i e s s o m e w h a t r e e c t s t h e t o p o l o g i c a l o r

g e o m e t r i c a l c o m p l e x i t y o f t h e s h a p e o f S w h i c h

i s i n d e p e n d e n t o f h .

S u p p o s e S : ! I R

3

i s a p a r a m e t r i z a t i o n o f

t h e s u r f a c e a n d P

h

i s a n i n t e r p o l a t i n g q u a d r i -

l a t e r a l m e s h w h o s e a v e r a g e m e s h s i z e i s m e a -

s u r e d b y t h e s t e p w i d t h h ( i n t h e d o m a i n ) .

T h e s i z e o f t h e l o c a l r e g i o n a r o u n d s i n g u l a r v e r -

t i c e s w h e r e O ( h

2

) - c o n v e r g e n c e o c c u r s , i s d e t e r -

m i n e d b y t o p o l o g i c a l d i s t a n c e o n P

h

, i . e . , b y t h e

m i n i m a l n u m b e r o f e d g e s c o n n e c t i n g s o m e v e r -

t e x t o a s i n g u l a r o n e . H e n c e t h e c o r r e s p o n d i n g

p o r t i o n o f t h e d o m a i n c o v e r s a n a r e a o f

k k = O ( h

2

) i f t h e n u m b e r o f s i n g u l a r v e r t i c e s

i s n i t e . U n d e r t h e s e c o n d i t i o n s t h e e r r o r f u n c -

t i o n a l c a n b e w r i t t e n a s

E

h

( S ) =

Z

n

O ( h

4

) +

Z

O ( h

2

)

= O ( h

4

) :

S e v e r a l a l g o r i t h m s t o r e d u c e t h e n u m b e r o f s i n -

g u l a r v e r t i c e s i n a m e s h h a v e b e e n p r o p o s e d

E c k 9 5 , S c h 9 5 ] . T h e m e s h e s o n e e n c o u n t e r s i n

p r a c t i c a l a p p l i c a t i o n s u s u a l l y h a v e a c o n s t a n t o r

v e r y s l o w l y i n c r e a s i n g n u m b e r o f s i n g u l a r i t i e s a s

t h e m e s h - s i z e h i s r e d u c e d , e . g . , 3 D - s c a n n e r s o f -

t e n m e a s u r e a n o b j e c t a l o n g p a r a l l e l p r o l e s o r

a l o n g r a d i a l m e r i d i a n s a n d t h u s g e n e r a t e t e n s o r -

p r o d u c t - t y p e d a t a .

E x a m p l e s

W e u s e t h e r e n e m e n t s c h e m e t o a p p r o x i m a t e

t w o t e s t s u r f a c e s : t h e u n i t - s p h e r e a n d t h e g r a p h

o f t h e f u n c t i o n c o s (

p

x

2

+ y

2

) o v e r t h e d o m a i n

= 0 1 ]

2

. I n t h e s p h e r e c a s e t h e n u m b e r o f s i n -

g u l a r v e r t i c e s i s c o n s t a n t , s i n c e t h e i n i t i a l n e t s

a r e g e n e r a t e d b y u s i n g n e t s a s d e p i c t e d i n F i g . 1

a n d p r o j e c t i n g t h e v e r t i c e s o n t o t h e s p h e r e . T h e

g r a p h - s u r f a c e i s a p p r o x i m a t e d b y t a k i n g e q u i d i s -

t a n t s a m p l e s i n x - a n d y - d i r e c t i o n .

h

R

j " j q

1

m a x j " j q

1

1 3 . 1 1 e - 0 4 5 . 6 4 e - 0 4

. 5 1 . 9 5 e - 0 5 1 5 . 9 3 . 5 9 e - 0 5 1 5 . 7

. 2 5 1 . 2 2 e - 0 6 1 6 . 0 2 . 2 6 e - 0 6 1 5 . 9

. 1 2 5 7 . 5 7 e - 0 8 1 6 . 1 1 . 4 1 e - 0 7 1 6 . 0

T a b l e 1 : A p p r o x i m a t i o n e r r o r s f o r t h e g r a p h o f

c o s ( r ) , r

2

= x

2

+ y

2

o v e r = 0 1 ]

2

, s t e p -

w i d t h h . W e s h o w t h e e r r o r s w i t h r e s p e c t t o t h e

L

1

- a n d t h e L

1

- n o r m s . T h e q - c o l u m n s c o n t a i n

t h e r a t i o o f s u c c e s s i v e e r r o r s i n d i c a t i n g O ( h

4

) -

c o n v e r g e n c e .

I n o r d e r t o d e m o n s t r a t e t h e d i e r e n t r a t e s o f

c o n v e r g e n c e i n t h e p a r a m e t r i c c a s e , w e c o m p u t e

t h e a p p r o x i m a t i o n e r r o r s e p a r a t e l y f o r t h e r e g u -

l a r p o r t i o n s o f t h e n e t . I n t h e f u n c t i o n a l s e t t i n g s

t h i s d i s t i n c t i o n i s n o t n e c e s s a r y s i n c e t h e n e t s

h a v e n o e x t r a o r d i n a r y v e r t i c e s .

N u m e r i c a l E x p e r i m e n t s

I n o r d e r t o e v a l u a t e t h e p r a c t i c a l u s e f u l n e s s

o f t h e p r e s e n t e d s c h e m e , w e a p p l y i t t o a s e t

o f e x a m p l e s a n d c o m p a r e o u r r e s u l t s t o t h o s e

o b t a i n e d b y s t a n d a r d t e c h n i q u e s b a s e d o n i m -

p l i c i t r e p r e s e n t a t i o n s u s e d , e . g . , i n t h e c o n t e x t

o f a d a p t i v e n i t e e l e m e n t a n a l y s i s b y S R 9 1 ,

B a u 9 5 ] . F o r b o t h t y p e s o f r e p r e s e n t a t i o n w e

5

8/11/2019 10.1.1.61.7862


h

R

r e g

j " j q

r e g

R

t o t

j " j q

t o t

. 2 5 1 . 6 4 e - 0 3 5 . 3 6 e - 0 3

. 1 2 5 1 . 2 2 e - 0 4 1 3 . 5 4 . 0 2 e - 0 4 1 3 . 3

. 0 6 2 5 8 . 0 9 e - 0 6 1 5 . 1 3 . 1 2 e - 0 5 1 2 . 9

. 0 3 1 2 5 4 . 8 8 e - 0 7 1 6 . 6 2 . 4 6 e - 0 6 1 2 . 7

T a b l e 2 : A p p r o x i m a t i o n e r r o r s f o r t h e s p h e r e .

T h e e r r o r s a r e n o r m a l i z e d t o t h e v o l u m e o f t h e

s p h e r e . T h e s t e p w i d t h h = 1 c o r r e s p o n d s t o a

c u b e a s i n i t i a l a p p r o x i m a t i o n .

c o m p a r e a p p r o x i m a t i o n s o b t a i n e d b y a d a p t i v e

u n i f o r m r e n e m e n t a n d a d a p t i v e i r r e g u l a r r e -

n e m e n t B a u 9 5 ] . B o t h s c h e m e s a r e a l s o t e s t e d

i n c o m b i n a t i o n w i t h m e s h s m o o t h i n g a s s u g -

g e s t e d b y S c h 9 3 , S c h 9 5 , R i c 9 6 ] , w h e r e t h e e l -

e m e n t s h a p e i s o p t i m i z e d a f t e r e a c h r e n e m e n t

s t e p b y u s i n g a l o c a l s m o o t h i n g t e c h n i q u e . W e

u s e t h e f o l l o w i n g a b b r e v i a t i o n s :

u n i f : : : u n i f o r m r e n e m e n t , e x a c t g e o m e t r y

a d a p : : : n o n - u n i f o r m r e n e m e n t , e x a c t g e -

o m e t r y

s m o o : : : w i t h m e s h s m o o t h i n g

o r i g : : : w i t h o u t m e s h s m o o t h i n g

i n t e : : : m e s h g e n e r a t i o n b y t h e p r o p o s e d

r e n e m e n t s c h e m e

C y l i n d r i c a l S h e l l w i t h H o l e

T h e c y l i n d r i c a l s h e l l , o r i g i n a l l y p r o p o s e d b y

S c o 6 9 ] i s r i g i d l y s u p p o r t e d a t b o t h c u r v e d

b o u n d a r i e s ( c f . F i g . 6 ) a n d l o a d e d b y a u n i f o r m l y

d i s t r i b u t e d l o a d i n v e r t i c a l d i r e c t i o n . T h e v e r t i -

c a l d i s p l a c e m e n t v a t t h e c e n t e r o f t h e f r e e e d g e s

i s t h e q u a n t i t y c o m p a r e d i n t h e a n a l y s i s .

D u e t o s y m m e t r y , o n l y a q u a r t e r o f t h e s h e l l h a s

t o b e d i s c r e t i z e d i n t h e n i t e e l e m e n t m o d e l . T h e

s t a r t i n g m e s h e s u s e d i n t h e a n a l y s i s a r e s h o w n i n

F i g . 7 . T h e y c o n s i s t o f 3 2 e l e m e n t s a n d 4 5 n o d e s .

I n o r d e r t o a c h i e v e a s u c i e n t l y g o o d a p p r o x i -

m a t i o n o f t h e g e o m e t r y , 5 n o d e s a r e c h o s e n i n

t h e l a t i t u d i n a l d i r e c t i o n . U s i n g a n i n i t i a l m e s h

h a v i n g 8 e l e m e n t s ( 1 5 n o d e s ) , i t t u r n s o u t t h a t

t h e g e o m e t r y i s n o t w e l l a p p r o x i m a t e d a n d t h e

n u m e r i c a l a n a l y s i s r e s u l t s i n a d i s p l a c e m e n t b e -

i n g 1 5 % l a r g e r t h a n t h e e x a c t v a l u e .

T h e m e s h g i v e n i n F i g . 7 b i s a s l i g h t m o d i c a -

t i o n o f t h e m e s h i n F i g . 7 a , w h i c h r e s u l t s f r o m

F i g u r e 6 : C i r c u l a r s h e l l w i t h h o l e , g e o m e t r y , m a -

t e r i a l p r o p e r t i e s , a n d l o a d i n g

m e s h s m o o t h i n g . A s t h e e l e m e n t b o u n d a r i e s a r e

n o l o n g e r a l i g n e d w i t h t h e d i r e c t i o n o f t h e m a i n

c u r v a t u r e s , t h i s l e a d s t o a s i g n i c a n t l y d i e r e n t

i n i t i a l g e o m e t r y f o r t h e i n t e r p o l a t i o n a l g o r i t h m .

T h e m o d i c a t i o n o f t h e m e s h l e a d s t o d i s p l a c e -

m e n t s d i e r i n g a b o u t 1 0 % f r o m t h e c o r r e c t v a l -

u e s .

I n F i g . 8 t h e m e s h e s f o r u n i f o r m a n d i r r e g u l a r

a d a p t i v e r e n e m e n t b a s e d o n t h e i n t e r p o l a t o r y

r e n e m e n t s c h e m e a r e g i v e n . C o n v e r g e n c e d i a -

g r a m s a r e d e p i c t e d f o r t h e d i s p l a c e m e n t v i n F i g .

9 a n d f o r t h e r e l a t i v e e r r o r i n e n e r g y i n F i g . 1 0 .

F o r c o m p a r i s o n , t h e c u r v e s f o r o t h e r r e n e m e n t

s t r a t e g i e s w h i c h a r e p e r f o r m e d o n t h e i m p l i c i t

r e p r e s e n t a t i o n o f t h e e x a c t g e o m e t r y , a r e a l s o

i n c l u d e d .

-0.36

-0.34

-0.32

-0.3

-0.28

-0.26

10 100 1000 10000 100000{Number of degrees of freedom}

inte, adapinte, unifunif, origunif, smooadap, origadap, smoo

F i g u r e 9 : C i r c u l a r s h e l l w i t h h o l e , C o n v e r g e n c e

o f m i d s i d e d i s p l a c e m e n t

D i s c u s s i o n

I f t h e s t a r t i n g m e s h i s r e a s o n a b l y n e , t h e i n -

t e r p o l a t i o n a l g o r i t h m l e a d s t o a l m o s t t h e s a m e

r e s u l t s a s t h e r e n e m e n t o n t h e i m p l i c i t l y g i v e n

6

8/11/2019 10.1.1.61.7862


F i g u r e 7 : C i r c u l a r s h e l l w i t h h o l e , a ) b a s i c , b ) s m o o t h e d s t a r t i n g m e s h

F i g u r e 8 : C i r c u l a r s h e l l w i t h h o l e , a ) u n i f o r m , b ) a d a p t i v e i r r e g u l a r n a l m e s h

3

5

10

20

30

40

50

100 1000 10000{Number of degrees of freedom}

inte, adapinte, unifunif, origunif, smooadap, origadap, smoo

F i g u r e 1 0 : C i r c u l a r s h e l l w i t h h o l e , C o n v e r g e n c e

o f e n e r g y e r r o r

g e o m e t r y f o r u n i f o r m r e n e m e n t . F o r a d a p t i v e

i r r e g u l a r r e n e m e n t t h e r e s u l t s c a n b e w o r s e i f

t h e s t a r t i n g m e s h i s t o o c o a r s e b e c a u s e i n t h i s

c a s e t h e g e o m e t r y i t s e l f i s n o t a p p r o x i m a t e d c o r -

r e c t l y . A c c o r d i n g t o t h e o b s e r v a t i o n s i n t h e l a s t

s e c t i o n , t h e a p p r o x i m a t i o n o r d e r o f t h e r e n e -

m e n t s c h e m e g u a r a n t e e s t h a t a h i g h e r d e n s i t y

o f t h e i n i t i a l m e s h s t r o n g l y i m p r o v e s t h e q u a l i t y

o f t h e a p p r o x i m a t i o n o f t h e g e o m e t r y .

A c o m b i n a t i o n o f t h e r e n e m e n t s c h e m e w i t h

s m o o t h i n g t e c h n i q u e s S c h 9 3 ] i s n o t p o s s i b l e , a s

t h e i n t e r p o l a t i o n n o d e s a r e x e d i n t h e c u r r e n t

f o r m o f t h e i n t e r p o l a t i o n a l g o r i t h m . H o w e v e r ,

t h e c o n v e r g e n c e o f t h e d i s p l a c e m e n t s i s v e r y

g o o d | o n l y t h e a l g o r i t h m s w h i c h u s e s m o o t h -

i n g t e c h n i q u e s l e a d t o b e t t e r r e s u l t s . I t t u r n s o u t

t h a t t h e a l g o r i t h m i s w e l l s u i t e d f o r p r a c t i c a l e n -

g i n e e r i n g a p p l i c a t i o n s .

C o m p o s i t e C y l i n d e r - P l a t e s t r u c t u r e

T h e s e c o n d e x a m p l e i s t h e i n t e r s e c t i o n o f a c y l i n -

d e r a n d a p l a t e b y a n a n g l e o f 4 5 d e g r e e s . T h e

w i d t h o f t h e q u a d r a t i c p l a t e i s w = 1 5 a n d t h e

c y l i n d e r o f r a d i u s r = 6 h a s a n a v e r a g e l e n g t h o f

l = 1 0 . T h e t h i c k n e s s o f t h e c y l i n d e r a n d o f t h e

p l a t e a r e b o t h t = 1 . T h e m o d u l u s o f E l a s t i c i t y

i s E = 1 0

6

a n d t h e P o i s s o n r a t i o i s s e t t o z e r o .

T h e p l a t e i s s i m p l y s u p p o r t e d . A u n i f o r m l y d i s -

t r i b u t e d l o a d i s a p p l i e d a t t h e f r e e e d g e o f t h e

c y l i n d e r a n d i s d i r e c t e d a l o n g t h e a x i s o f t h e

c y l i n d e r .

A s w e e x p l a i n e d e a r l i e r , t h e c r e a s e a l o n g t h e i n -

t e r s e c t i o n c u r v e c a n b e m o d e l e d b y j o i n i n g t w o

s e p a r a t e p a r t s , t h e c y l i n d e r a n d t h e t h e p l a t e ,

h a v i n g a c o m m o n b o u n d a r y p o l y g o n . H e n c e ,

c o m m o n n o d e s h a v e t o b e s t o r e d t w i c e s i n c e o t h -

e r w i s e t h e i n t e r p o l a t o r y r e n e m e n t w o u l d g e n -

e r a t e a s m o o t h j o i n .

7

8/11/2019 10.1.1.61.7862


F i g u r e 1 1 : C o m p o s i t e s t r u c t u r e ( c o m m o n

b o u n d a r y i n d i c a t e d ) , s t a r t i n g m e s h

T h e i n i t i a l a n d t h e n a l m e s h e s o f t h e a d a p t i v e

a n a l y s i s a r e s h o w n i n F i g . 1 1 a n d 1 2 . F o r a n e x -

a c t i n t e r p o l a t i o n o f t h e p l a t e i t i s s u c i e n t t o

d e n e 4 c o r n e r v e r t i c e s ( l i n e a r p r e c i s i o n ) . T h e

c y l i n d e r i s d e n e d b y t h e b o u n d a r y p o l y g o n s a t

t h e f r e e e d g e a n d a t t h e i n t e r s e c t i o n . U s i n g 8

n o d e s a l o n g t h e s e b o u n d a r i e s g u a r a n t e e s a s u f -

c i e n t l y a c c u r a t e r e p r o d u c t i o n .

D i s c u s s i o n

T h e c o n v e r g e n c e d i a g r a m f o r t h e e n e r g y e r r o r

i n F i g . 1 3 s h o w s a s l i g h t l y r e d u c e d c o n v e r g e n c e

f o r t h e i n t e r p o l a t i o n a l g o r i t h m i n c o m b i n a t i o n

w i t h u n i f o r m r e n e m e n t c o m p a r e d t o t h e r e -

s u l t s o b t a i n e d b y u s i n g t h e o t h e r a l g o r i t h m s .

W i t h a d a p t i v e i r r e g u l a r r e n e m e n t t h e r e s u l t s

o b t a i n e d w i t h o u r s c h e m e a r e a l m o s t i d e n t i c a l

t o a l l o t h e r r e s u l t s .

C o n c l u s i o n

T h e p r e s e n t e d s c h e m e t o g e n e r a t e a d a p t i v e l y

r e n e d q u a d r i l a t e r a l m e s h e s a l l o w s t o a p p r o x i -

m a t e a r b i t r a r y s u r f a c e s w i t h h o l e s a n d s h a r p l o -

c a l f e a t u r e s . T h e q u a l i t y o f t h e a p p r o x i m a t i o n i s

s u c i e n t l y w e l l , e v e n i f t h e s c h e m e i s a p p l i e d t o

m o d e r a t e l y d e t a i l e d i n i t i a l m e s h e s . D u e t o t h e

s y m m e t r y o f t h e r e n e m e n t r u l e s f o r t h e e d g e s

a r o u n d a v e r t e x , t h e l o c a l d i s t o r t i o n i s r e d u c e d ,

i . e . , i f n e d g e s m e e t a t a c o m m o n v e r t e x , t h e

a n g l e s b e t w e e n t h e m t e n d t o

2

n

. T h i s r e g u l a r -

i z a t i o n e e c t c a u s e s t h e s h a p e s o f a d j a c e n t f a c e s

t o b e c o m e a l i k e .

T h e s l i g h t l y l e s s a c c u r a t e r e s u l t s o b t a i n e d b y u s -

i n g t h e i n t e r p o l a t o r y r e n e m e n t s c h e m e i n s t e a d

o f a n e x a c t r e p r e s e n t a t i o n o f t h e g i v e n o b j e c t s

( c f . l a s t S e c t i o n ) s t e m s f r o m t h e f a c t t h a t t h e

s u b d i v i s i o n s c h e m e c o m p u t e s n e w v e r t i c e s w i t h -

o u t a n y m e t a - k n o w l e d g e a b o u t t h e s h a p e ( i n

t h i s c a s e t h a t t h e s u r f a c e s a r e c y l i n d r i c a l ) . H o w -

e v e r , t h e s c h e m e i s m a i n l y d e s i g n e d f o r c o m -

p u t a t i o n s o n f r e e f o r m s u r f a c e s o f w h i c h o n l y

d i s c r e t e p o i n t s a r e g i v e n a n d a n e x a c t r e p r e s e n -

t a t i o n i s n o t k n o w n . T h e r e f o r e a l s o o t h e r a p -

p r o a c h e s b a s e d o n s p l i n e s w o u l d h a v e t o g u e s s

t h e a c t u a l s h a p e .

T h e p o w e r o f t h e s c h e m e i s d e m o n s t r a t e d b e s t

b y a p p l y i n g i t t o a m o r e s o p h i s t i c a t e d m o d e l .

T h e s u r f a c e i n F i g . 1 4 i s g e n e r a t e d b y u s i n g t h e

c o a r s e a p p r o x i m a t i o n a s i n p u t d a t a a n d r e n i n g

t w i c e . N o a d d i t i o n a l i n f o r m a t i o n i s n e c e s s a r y .

C o n s t r u c t i n g a p i e c e w i s e s m o o t h s p l i n e s u r f a c e

i n t e r p o l a t i n g t h i s i r r e g u l a r n e t w o u l d b e r a t h e r

c o m p l i c a t e d .

T h e m a j o r a d v a n t a g e o f t h i s s c h e m e c o m p a r e d

t o o t h e r a p p r o a c h e s t h e r e f o r e i s i t s s t r i k i n g s i m -

p l i c i t y . O n l y l o c a l a n e c o m b i n a t i o n s o f p o i n t s

h a v e t o b e c o m p u t e d i n o r d e r t o r e n e a g i v e n

n e t . A l m o s t n o s p e c i a l c a s e s c o n c e r n i n g t h e l o -

c a l t o p o l o g y o f t h e m e s h h a v e t o b e c o n s i d e r e d .

T h e e x a m p l e s h a v e b e e n g e n e r a t e d b y o u r i m -

p l e m e n t a t i o n w h i c h p r o v i d e s a l l t h e p r e s e n t e d

f e a t u r e s i n c l u d i n g a d a p t i v e i r r e g u l a r r e n e m e n t

w i t h C

0

- c o n s i s t e n c y p r e s e r v a t i o n . ( p l e a s e c o n -

t a c t : k o b b e l t @ c s . w i s c . e d u )

R e f e r e n c e s

B C X 9 4 a ] C . B a j a j / J . C h e n / G . X u , M o d -

e l i n g w i t h C u b i c A - P a t c h e s , A C M

T r a n s . G r a p h . 1 9 9 4

B C X 9 4 b ] C . B a j a j / J . C h e n / G . X u , S m o o t h

L o w D e g r e e A p p r o x i m a t i o n s o f P o l y -

h e d r a , p r e p r i n t 1 9 9 4

B D 8 6 ] K . B a t h e / E . D v o r k i n , A f o r m u l a -

t i o n f o G e n e r a l S h e l l E l e m e n t s - T h e

U s e o f M i x e d I n t e r p o l a t i o n o f T e n s o -

r i a l C o m p o n e n t s , I n t . J . N u m . M e t h .

E n g . 2 2 ( 1 9 8 6 )

B a u 9 5 ] M . B a u m a n n / K . S c h w e i z e r h o f ,

A d a p t i v e m e s h g e n e r a t i o n o n a r b i -

t r a r i l y c u r v e d s h e l l s t r u c t u r e s , i n

M . P a p a d r a k a k i s a n d B . T o p p i n g ,

e d s , P r o c e e d i n g s o f 2 n d I n t . C o n f .

o n C o m p u t . S t r u c t u r e s T e c h n . i n

8

8/11/2019 10.1.1.61.7862


F i g u r e 1 2 : C o m p o s i t e s t r u c t u r e , n a l m e s h e s a ) u n i f o r m b ) i r r e g u l a r r e n e m e n t

10

100

1000 10000{Number of degrees of freedom}

inte, adapinte, unifadap, smooadap, origunif, origunif, smoo

F i g u r e 1 3 : C o m p o s i t e s t r u c t u r e , c o n v e r g e n c e d i a g r a m f o r e n e r g y e r r o r

" ' A d v a n c e s i n F i n i t e E l e m e n t T e c h -

n i q u e s " ' , C i v i l - C o m p P r e s s E d i n -

b u r g h , 1 9 9 5 .

C C 7 8 ] E . C a t m u l l / J . C l a r k , R e c u r s i v e l y

G e n e r a t e d B - S p l i n e S u r f a c e s o n A r -

b i t r a r y T o p o l o g i c a l M e s h e s , C A D 1 0

( 1 9 7 8 ) , p p . 3 5 0 { 3 5 5

D G L 8 7 ] N . D y n / J . G r e g o r y / D . L e v i n ,

A 4 - P o i n t I n t e r p o l a t o r y S u b d i v i s i o n

S c h e m e f o r C u r v e D e s i g n , C A G D

4 ( 1 9 8 7 ) , p p . 2 5 7 { 2 6 8

D G L 9 0 ] N . D y n / J . G r e g o r y / D . L e v i n ,

A B u t t e r y S u b d i v i s i o n S c h e m e

f o r S u r f a c e I n t e r p o l a t i o n w i t h T e n -

s i o n C o n t r o l l , A C M T r a n s . G r a p h . 9

( 1 9 9 0 ) , p p . 1 6 0 { 1 6 9

D S 7 8 ] D . D o o / M . S a b i n , B e h a v i o u r o f R e -

c u r s i v e D i v i s i o n S u r f a c e s N e a r E x -

t r a o r d i n a r y P o i n t s , C A D 1 0 ( 1 9 7 8 ) ,

p p . 3 5 6 { 3 6 0

D y n 9 1 ] N . D y n , S u b d i v i s i o n S c h e m e s i n

C o m p u t e r A i d e d G e o m e t r i c D e s i g n ,

A d v a n c e s i n N u m e r i c a l A n a l y s i s I I ,

W a v e l e t s , S u b d i v i s i o n s a n d R a d i a l

F u n c t i o n s , W . A . L i g h t e d . , O x f o r d

U n i v e r s i t y P r e s s , 1 9 9 1 , p p . 3 6 { 1 0 4 .

E c k 9 5 ] E c k / D e R o s e / D u c h a m p / H o p p e

/ L o u n s b e r y / S t u e t z l e , M u l t i r e s o l u -

t i o n A n a l y s i s o f A r b i t r a r y M e s h e s ,

R e p o r t , U n i v e r s i t y o f W a s h i n g t o n ,

U W - C S E - 9 5 - 0 1 - 0 2

H o p 9 4 ] H . H o p p e , S u r f a c e R e c o n s t r u c t i o n

f r o m u n o r g a n i z e d p o i n t s , T h e s i s ,

U n i v e r s i t y o f W a s h i n g t o n , 1 9 9 4

H L 9 3 ] J . H o s c h e k / D . L a s s e r , F u n d a m e n -

t a l s o f c o m p u t e r - a i d e d g e o m e t r i c d e -

s i g n , A . K . P e t e r s , 1 9 9 3

K o b 9 5 ] L . K o b b e l t , I t e r a t i v e E r z e u g u n g

g l a t t e r I n t e r p o l a n t e n , T h e s i s , U n i -

v e r s i t a t K a r l s r u h e , 1 9 9 5

9

8/11/2019 10.1.1.61.7862


F i g u r e 1 4 : G e n e r a t i n g c o m p l e x s h a p e s b y i n t e r p o l a t o r y s u b d i v i s i o n

K o b 9 6 ] L . K o b b e l t , I n t e r p o l a t o r y S u b d i v i -

s i o n o n O p e n Q u a d r i l a t e r a l N e t s

w i t h A r b i t r a r y T o p o l o g y , t o a p p e a r

i n C o m p u t e r G r a p h i c s F o r u m

L o o 8 7 ] C . L o o p , S m o o t h S u b d i v i s i o n f o r

S u r f a c e s B a s e d o n T r i a n g l e s , M a s -

t e r T h e s i s , U n i v e r s i t y o f U t a h , 1 9 8 7

L o o 9 4 ] C . L o o p , A G

1

t r i a n g u l a r s p l i n e

s u r f a c e o f a r b i t r a r y t o p o l o g i c a l t y p e ,

C A G D 1 1 ( 1 9 9 4 ) , p p . 3 0 3 { 3 3 0

P r a 9 5 ] H . P r a u t z s c h , C

k

- a n a l y s i s o f s u b -

d i v i s i o n s u r f a c e s a t e x t r a o r d i n a r y

p o i n t s , s u b m i t t e d , 1 9 9 5

P S B 9 0 ] L . P l a n k / E . S t e i n / D . B i s c h o , A c -

c u r a c y a n d A d a p t i v i t y i n t h e N u m e r -

i c a l A n a l y s i s o f T h i n W a l l e d S t r u c -

t u r e s , C o m p . M e t h . i n A p p l . M e c h .

E n g . , 8 2 ( 1 9 9 0 )

R e i 9 5 ] U . R e i f , A u n i e d a p p r o a c h t o s u b d i -

v i s i o n a l g o r i t h m s n e a r e x t r a o r d i n a r y

v e r t i c e s , C A G D 1 2 ( 1 9 9 5 ) , p p . 1 5 3 {

1 7 4

R i c 9 6 ] J . R i c c i u s / K . S c h w e i z e r h o f / M .

B a u m a n n , C o m b i n a t i o n o f a d a p t i v -

i t y a n d m e s h s m o o t h i n g f o r t h e n i t e

e l e m e n t a n a l y s i s o f s h e l l s w i t h i n t e r -

s e c t i o n s , t o a p p e a r i n I n t . j . n u m e r .

m e t h o d s e n g .

S F 8 9 a ] J . S i m o / D . F o x , O n a S t r e s s

R e s u l t a n t G e o m e t r i c a l l y E x a c t S h e l l

M o d e l . P a r t I : F o r m u l a t i o n a n d

O p t i m a l P a r a m e t r i z a t i o n , C o m p .

M e t h . A p p l . M e c h . E n g . 7 2 ( 1 9 8 9 )

S F 8 9 b ] J . S i m o / D . F o x , O n a S t r e s s

R e s u l t a n t G e o m e t r i c a l l y E x a c t S h e l l

M o d e l . P a r t I I : T h e L i n e a r T h e -

o r y : C o m p u t a t i o n a l A s p e c t s , C o m p .

M e t h . A p p l . M e c h . E n g . 7 3 ( 1 9 8 9 )

S c h 9 3 ] K . S c h w e i z e r h o f / J . R i c c i u s /

M . B a u m a n n , V e r b e s s e r u n g v o n F i -

n i t e E l e m e n t B e r e c h n u n g e n d u r c h

A d a p t i v i t a t u n d N e t z g l a t t u n g a m

B e i s p i e l e b e n e r u n d g e k r u m m t e r

F l " a c h e n t r a g w e r k e , T e c h n i c a l r e -

p o r t , U n i v e r s i t a t K a r l s r u h e , I n s t i -

t u t f u r W i s s e n s c h a f t l i c h e s R e c h n e n

u n d M a t h e m a t i s c h e M o d e l l b i l d u n g ,

1 9 9 3

S c h 9 5 ] R . S c h r e i b e r , K o m b i n a t o r i s c h o p -

t i m i e r t e K o n s t r u k t i o n f a s t - r e g u l r e r

P o l y g o n n e t z e T h e s i s , U n i v e r s i t a t

S t u t t g a r t , 1 9 9 5

S c o 6 9 ] A . S c o r d e l i s / K . L o , C o m p u t e r a n a l -

y s i s o f c y l i n d r i c a l s h e l l s , J . o f A m e r -

i c a n C o n c r e t e I n s t i t u t e , 6 1 ( 1 9 9 6 ) ,

5 3 9 { 5 6 1

S R 9 1 ] E . S t e i n / W . R u s t , M e s h a d a p -

t a t i o n s f o r l i n e a r 2 D n i t e - e l e m e n t

d i s c r e t i z a t i o n s i n s t r u c t u r a l m e c h a n -

i c s , e s p e c i a l l y i n t h i n s h e l l a n a l y s i s ,

J . C o m p u t . A p p l . M a t h . , 3 6 ( 1 9 9 1 ) ,

1 0 7 { 1 2 9

Z S S 9 6 ] D . Z o r i n / P . S c h r o e d e r / W .

S w e l d e n s , I n t e r p o l a t i n g S u b d i v i s i o n

f o r M e s h e s o f A r b i t r a r y T o p o l o g y t o

a p p e a r i n S I G G R A P H P r o c e e d i n g s

1 9 9 6

1 0

10.1.1.61.7862

Documents