reset
#{ R = 0.5 }
#{ t = 0.02 }
#{ H = 1.0 }

# create overall domain
create surface rectangle width { R + t } height { H } zplane 

# move such that we have left side at rotation axis
move Surface 1  x { (R + t)/2 }

# cut off the pipe
webcut body 1  with plane xplane offset { R }

# merge for conforming mesh
merge curve all

# regions
block 1 add surface 3 
block 1 name "S_a"
block 2 add surface 2 
block 2 name "S_s"

# interfaces
nodeset 1 add curve 5  
nodeset 1 name "L_as"
nodeset 2 add curve 6 7  
nodeset 2 name "L_s_end"

# mesh
# mech domain
curve 6 7  interval 4
curve 6 7  scheme equal
curve 5 4  size 0.015
curve 5 4  scheme equal
# acou domain
curve 10  tangent opposite
curve 9 10  scheme bias fine size 0.02 coarse size 0.075 
mesh curve all
surface 3  scheme pave
mesh surface 3 
mesh surface all
smooth surface 3 


# save cdb
export ansys cdb "pipe.cdb"  overwrite