Next: Using macros to model
Up: Modelling the geometry
Previous: Modelling the geometry
  Contents
Both plungers are topologically the same. They consist of a tube where
inside of the tube a circular cylinder with a rounded cap sits in.
Figure 5.2:
Details of the technical drawing, showing the plungers.
|
Figure 5.3:
An outline of the body of revolution that shall decribe the body of a plunger.
|
A plunger is a body of revolution. gd1 can model this directly.
We edit our inputfile so that it contains:
#
# a plunger
#
define(PlungerRadius0, 110e-3/2 )
define(PlungerInnerRadius, 100e-3/2 )
define(PlungerCurvature, 16e-3)
-gbor
material= 10
originprime= (0,0,0)
zprimedirection= (0,0,1)
rprimedirection= (1,0,0)
range= (0,360)
clear
# point= (z,r)
point= ( 0,0 )
point= ( 0, PlungerInnerRadius-PlungerCurvature )
arc, radius= PlungerCurvature, size= small, type= counterclockwise
point= ( -PlungerCurvature, PlungerInnerRadius )
point= ( -170e-3, PlungerInnerRadius )
point= ( -170e-3, 0 )
show= now,
doit
The figure 5.3 shows an outline of the body of revolution that
this decribes.
This plunger has its axis direction in z-direction.
Our plungers shall have their axis lying in the x-y-plane,
with an angle of -90+22.5 and 17 degrees.
In order to have the axis of our plunger direct in the proper direction, we
change the values of zprimedirection, rprimedirection
.
These two vectors define the local z', r'
coordinate-system, in which
the body of revolution is described.
We edit our inputfile:
define(PlungerAngle, (-90+22.5)*@pi/180 )
originprime= (0,0,0)
zprimedirection= ( -cos(PlungerAngle) ,\
-sin(PlungerAngle) ,\
0 )
rprimedirection= ( 0, 0, 1 )
The resulting outline is shown in figure 5.4.
Figure 5.4:
The plunger now has its axis showing in the right direction.
|
We are not done yet: The plunger is not yet at the right position.
The origin of the plunger shall not be at (x,y,z)=(0,0,0), but at
(x,y,z)
=(
), with
degrees.
We change our inputfile:
define(PlungerRadius0, OuterRadius-50e-3 )
define(PlungerAngle, (-90+22.5)*@pi/180 )
originprime= ( cos(PlungerAngle)*PlungerRadius0 ,\
sin(PlungerAngle)*PlungerRadius0 ,\
0 )
zprimedirection= ( -cos(PlungerAngle) ,\
-sin(PlungerAngle) ,\
0 )
rprimedirection= ( 0, 0, 1 )
When we feed gd1 with this inputfile, we do not see the plunger in the
plot of the material-discretisation.
The reason is: The plunger is outside of the specified computational volume.
Since the geometry with the plunger does no longer have all three
symmetry-planes, we have to compute in a much larger volume.
The only symmetry plane left is the plane z=0.
So we change the specifications for the boundaries of the computational
volume to:
###
### We define the borders of the computational volume,
### we define the default mesh-spacing,
### and we define the conditions at the borders:
###
-mesh
spacing= InnerRadius/15
pxlow= -1.1*OuterRadius
pylow= -1.1*OuterRadius
pzlow = -(GapLength/2+TaperLength+9e-2)
pxhigh= 1.1*OuterRadius
pyhigh= 1.1*OuterRadius
pzhigh= 0
#
# The conditions to use at the borders of the computational volume:
#
cxlow= electric, cxhigh= electric
cylow= electric, cyhigh= electric
czlow= electric, czhigh= electric
The so edited inputfile can be found as
"/usr/local/gd1/Tutorial-SRRC/wPlunger00.gdf".
When we feed gd1 with this inputfile (gd1 wPlunger00.gdf), we see
a picture similiar as the one shown in figure 5.5
Figure 5.5:
Screenshot of the desktop when the inputfile wPlunger00.gdf has been fed
into gd1.
|
Subsections
Next: Using macros to model
Up: Modelling the geometry
Previous: Modelling the geometry
  Contents