Download ZINC 3.2 (Finite element code) User Manual

ZINC 3.2 user manual
NPL Report
part 2
region 2
type sphere
fab 0.3
surface region 1
coat 1 6
end
part 3
region 3
type boundyup
end
1
1
1
6
xxxxxxxxxxxxxxxxxxxxxxxxxxx
6
6
xxxxxxxxxxxxxxxxxxxxxxxxxxx
1 xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx6
xxxxxxxxxxxxxxxxxxxxxxxxxxx
1
1
1
1
1
1
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
6
2
6
2
6
1
6
1
2
6
6
1
6
6
1
3
1
1
1
1
1
1
1
1
1
1
1
4
(a)
3
part 4
region 4
type boundydn
end
3
2
3
3
3
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxx
2
2
1
1
4
4
4
4
2
2
2
2
2
2
2
2
1
2
1
2
2
2
4
1
(b)
Figure 3.4: (a) Use of the coat command to change the node region numbers at an interface;
(b) use of open parts (boundyup and boundydn) to change the node region numbers
to clamped displacements at the top and bottom of the simulation, and so on.
Figure 3.5 shows the corresponding file.zin for Figure 3.4 set up as an electrostatic problem. This file will be described in more detail in Section 5.2 but for now note how region 1 and
region 2 elements are given permittivity values of ǫ0 (permittivity of free space, 8.854 × 10−12
F/m) and 2ǫ0 . Also, region 3 and region 4 nodes are set to voltage values 1 V and 0 V respectively. Region 1 and 2 nodes are not set at all, so these nodes (the interior nodes in the problem)
are allowed to vary in order to solve the electrostatic equations.
By the end of the part fitting process, all elements and nodes should be given an region
number. If this is not the case, Zmesh will report an error. It is common practice for the first
part to be a box part which covers the whole simulation area. Other parts are then inserted
into this part.
We have so far described most of the part commands: type, region, fab, surface,
coat. The other two commands are shift and rotate, allowing the part to be shifted and rotated
into the required position. Parts are conceptually constructed at the origin (0, 0, 0) with major
axes generally along x, y, z. Thus, the “box” part is constructed with principle axes along the
x,y,z directions. It is centred at the origin and its extents are given by the 3 numbers following
the fab (fabricate) command. The meaning of the fab parameters depends on the part and is
given in Table 3.1.
The rotation command has the form
rotate rx ry rz [string]
If string is absent then a rotation is performed about the (global) x-axis, then the y-axis then
the z-axis by the angle given, in degrees. If a different order is required, this can be given using
the string. Eg using “yzx” would cause rotation to occur about y, then z then x. Since these are
Euler angles, rotation order is important. Furthermore, rotation is performed before the object
is shifted irrespective of the order of commands in the part block.
The shift command has the form
Page 15 of 63