reset
# Geometrie of the bar
#{ length = 1}
#{ hight = 0.01}
#{ width = 0.02}

# Number of nodes for each direction
#{ NodesX = 500}
#{ NodesY = 5}
#{ NodesZ = 10}

# Geometrie data of the pressure load on top of the bar
#{pressureMiddle = 0.7}
#{pressureWidth = 0.1}

# Length of each volume in x-driection
#{ length1 = pressureMiddle-pressureWidth/2}
#{ length2 = pressureWidth}
#{ length3 = length-pressureMiddle-pressureWidth/2}

brick x {length} y {hight} z {width}
move Volume 1  x 0.5 y 0 z 0 include_merged 
webcut volume 1  with plane xplane offset {length1}
webcut volume 1  with plane xplane offset {length1 + length2}
imprint all
merge all

curve 25 26 28 27  interval {NodesX*length1/length}
curve 25 26 28 27  scheme equal
mesh curve 25 26 28 27 

curve 34 36 35 33  interval {NodesX*length3/length}
curve 34 36 35 33  scheme equal
mesh curve 34 36 35 33  

curve 43 44 41 42  interval {NodesX*length2/length}
curve 43 44 41 42  scheme equal
mesh curve 43 44 41 42 

curve 15 29 3 1 5 31 13 7  interval {NodesY}
curve 15 29 3 1 5 31 13 7  scheme equal
mesh curve 15 29 3 1 5 31 13 7 

curve 12 30 16 11 10 32 14 9 interval {NodesZ}
curve 12 30 16 11 10 32 14 9 scheme equal
mesh curve 12 30 16 11 10 32 14 9 

mesh volume 1 2 3

block 1 add volume 1 2 3
block 1 name "V_beam"

nodeset 1 add surface 4 6
nodeset 1 name "S_fix"

nodeset 2 add surface 26
nodeset 2 name "S_pressureY"

nodeset 3 add surface 24
nodeset 3 name "S_pressureZ"

export ansys cdb "bendingBeam.cdb" overwrite