5.3
TRANSMISSION CAM WITH CAM ELEMENTS NO.5
Copy the example file
B3_X.DXF to Z88X.DXF.
B3_X.DXF ---> Z88X.DXF
input file for CAD converter Z88X
CAD:
You should only
look within this example at the CAD FE structure without producing it.
This comes
with later examples. Import Z88X.DXF into your CAD program and view it.
Usually
you would draw or model the structure in your CAD system. Do not change
anything and leave your CAD program without saving, converting etc. If
you do
not have any suitable CAD system, then drop this step.
Z88X, conversion from Z88X.DXF to
Z88I1.TXT, Z88I2.TXT and Z88I3.TXT.
Windows: In the Z88 commander
launch Z88X, press Button DXF --> Z88I* (default), Run
Button.
LINUX/UNIX:
In
the Z88 commander press button DXF --> Z88I* under CAD converters.
Z88O, look at finite element structure. Proceed as follows:
Windows: In the Z88
commander launch Z88O, Run Button. Use
the Wireframe Mode. Switch on Mouse
(4th button from left). You may also
LINUX/UNIX: In the Z88
commander launch Z88O, Run Button. Use
the Wireframe Mode. Switch on Mouse (7th button from
top). You may
also
Windows: In the Z88 commander
press button Z88F, CD Button (is
default), Run Button.
LINUX/UNIX:
In
the Z88 commander press button Z88F -C.
Windows: In the Z88 commander
press button Z88D, Run Button.
LINUX/UNIX:
In
the Z88 commander press button Z88D.
Windows: In the Z88 commander
press button Z88E, Run Button.
LINUX/UNIX:
In
the Z88 commander press button Z88E.
Windows, LINUX/UNIX: Choose Deflected.
The
deflections are magnified per default by the factor 100 which is not
enough for
this example. Enter
1,000 each for FUX, FUY
and FUZ in Menu > Factors >
Deflections.
Basically, the
calculation
and displaying of von Mises stresses is not provided in Z88 for cams No.5, because newer literal sources
state correctly that reduced stresses for cams and other machinery
parts under
dynamic loads do not only depend on the normal and direct stresses
(which are
computed by Z88), but also on stress concentration factors (impossible
to
calculate in Z88 and other FEA systems) and other factors.
Task: A transmission cam is designed as
follows:
* Cam section, D= 30 mm, L = 30 mm , fixed bearing at the left end
* Gear wheel 1, reference circle D = 45 mm, L = 20 mm
* Cam section, D = 35 mm, L = 60 mm, moveable bearing in the middle
* Gear wheel 2, reference circle D = 60 mm, L = 15 mm
* Cam section, D = 40 mm, L = 60 mm, moveable bearing at the right end
For the loads we picture
the cam with the following coordinate system: If we look onto the cam
as the main
view, then the origin should be at the left end in the middle of the
cam. X
runs along the cam, Z runs to the upper direction, Y runs in the rear.
Gear wheel 1 gets the
following loads in the (physical) point X1 = 40, Y1 = -22.5, Z1 = 0:
Fx1 =
-10,801 N, Fy1 = 6,809 N, Fz1 = 18,708 N. Fx1 results in a bending
moment M1
around the Z axis of -243,023 Nmm.
Gear wheel 2 gets the
following loads in the (physical) point X2 = 117.5, Y2 = 0, Z2 = 30:
Fx2 =
8,101 N, Fy2 = -14,031 N, Fz2 = -5,107 N. Fx2 results in a bending
moment M2
around the Y axis of -243,030 Nmm.
This results in loads in XY
and XZ plane. The "physical" points do not exist in the FE
calculation, of course, because a cam element is formed analytically
only of
two points along an axis. The Y and Z coordinates are always 0.
The cam is subdivided into eight cam elements No.5 = 9 nodes. The bearings are assumed in the nodes 1, 5 and 9. Very important: Node 1 is fixed in addition in the degree of freedom 4 (the torsion degree of freedom) in order to compute the torsion angle between the two gears. Otherwise, the structure is statically underdefined !
5.3.1
Input
This example can
almost be entered easier by editor into a file than with CAD. The CAD
use has
real advantages for the examples 1, 2, 5 and 6. Both ways are shown
below:
With
CAD program:
Proceed
according to the description of chapter 2.7.
Do not forget to write the element
information on the layer Z88EIO by TEXT function:
FE 1 5 (1st finite element
type 5)
FE 2 5 (2nd finite element type 5)
FE 3 5 (3rd finite element type 5)
FE 4 5 (4th finite element type 5)
FE 5 5 (5th finite element type 5)
FE 6 5 (6th finite element type 5)
FE 7 5 (7th finite element type 5)
FE 8 5 (8th finite element type 5)
Write the general
information and material information on the layer Z88GEN :
Z88I1.TXT 3 9 8 54 3 0 0 0 0
(3-Dim, 9 nodes, 8 elements, 54 DOF, 3 mat infos, flags 0 )
MAT 1 1 3 206000 0.3 1 30 (1st mat info for ele 1 to ele 3,
Young's,Poisson's,QPARA)
MAT 2 4 6 206000 0.3 1 35 (2nd mat info for ele 4 to ele 6,
Young's,Poisson's,QPARA)
MAT 3 7 7 206000 0.3 1 40 (3rd mat info for ele 7 to ele 7,
Young's,Poisson's,QPARA)
(INTORD is set here to
1, has no influence)
As cam elements No.5 are
structure elements (thus not subdividable like finite elements), the
mesh
generator cannot be used. You can immediately write the boundary
conditions
with the TEXT function on the layer Z88RBD:
Z88I2.TXT 18 (18
Boundary conditions altogether)
RBD 1 1 1 2 0 (1.BC: Node 1, DOF 1 (=X) fixed)
RBD 2 1 2 2 0 (2.BC: Node 1, DOF 2 (=Y) fixed)
RBD 3 1 3 2 0 (3.BC: Node 1, DOF 3 (=Z) fixed)
RBD 4 1 4 2 0 (4.BC: Node 1, DOF 4 (=torsion) fixed)
RBD 5 3 1 1 -10801 (5.BC: Node 3, DOF 1 (=X), load -10,801 N)
RBD 6 3 2 1 +6809 (6.BC: Node 3, DOF 2 (=Y), load 6,809 N)
RBD 7 3 3 1 +18708 (7.BC: Node 3, DOF 3 (=Z), load 18,708 N)
RBD 8 3 4 1 -420930 (8.BC: Node 3, DOF 4 (torsion) -420,930 Nmm)
RBD 9 3 6 1 -243023 (9.BC: Node 3, DOF 6 (bend. moment around
Z),-243,023Nmm)
RBD 10 5 2 2 0
RBD 11 5 3 2 0
RBD 12 7 1 1 +8101
RBD 13 7 2 1 -14031
RBD 14 7 3 1 -5107
RBD 15 7 4 1 +420930
RBD 16 7 5 1 -243030
RBD 17 9 2 2 0
RBD 18 9 3 2 0
... And write on the layer
Z88GEN onto any free place of your drawing the stress parameters for
the stress
calculation:
Z88I3.TXT 0 0 0 (any
stress parameters for Trusses No.4)
Export the drawing as DXF
file with the name Z88X.DXF and then launch the CAD converter Z88X with
the
option "from Z88X.DXF to Z88I*.TXT" (DXF -> I*). The CAD converter
will produce the input files Z88I1.TXT, Z88I2.TXT, Z88I3.TXT.
With
an editor:
Enter the structure data into Z88I1.TXT by
editor (cf. section 3.2):
3 9 8 54 3
0 0 0 0 (3D,
9 Node, 8 Ele, 54 DOF, 3 E-Gesetze, Flags 0)
1 6 0
0 0 (Node 1, 6 DOF, X-, Y- und Z-Koordinate)
2 6 30 0
0 (Node 2, 6 DOF, X-, Y- und
Z-Koordinate)
3 6 40 0 0
4 6 50 0 0
5 6 80 0 0
6 6 110 0 0
7 6 117.5 0 0
8 6 125 0
0
9 6 185 0
0
1 5 (Element 1, cam No.5)
1 2 (Coincidence Ele 1)
2 5 (Element 2, type 5)
2 3 (coincidence Ele 2)
..... (Elemente 3 to 7 dropped here)
8 5
8 9
1 3 206000
0.3 1
30 (mat info from Ele 1 to
3,Young's,Poisson's, QPARA= 30)
4 6 206000
0.3 1
35 (mat info from Ele 4 to
6,Young's,Poisson's, QPARA= 35)
7 7 206000
0.3 1
40 (mat info from Ele 7 to
7,Young's,Poisson's, QPARA= 40)
INTORD is set here to 1,
has no influence.
The boundary conditions Z88I2.TXT:
18 (18 Boundary
conditions)
1 1
2 0
(1.BC: Node 1, DOF 1 (=X) fixed)
1 2
2 0
(2.BC: Node 1, DOF 2 (=Y) fixed)
1 3 2
0 (3.BC:
Node 1, DOF 3 (=Z) fixed)
1 4 2
0 (4.BC:
Node 1, DOF 4 (=torsion) fixed)
3 1 1 -10801
(5.BC: Node 3, DOF 1 (=X), load
-10,801 N)
3 2 1
+6809
(6.BC: Node 3, DOF 2 (=Y), load 6,809 N)
3 3 1
+18708
(7.BC: Node 3, DOF 3 (=Z), load 18,708 N)
3 4 1
-420930
(8.BC: Node 3, DOF 4 (torsion) -420,930 Nmm)
3 6 1
-243023
(9.BC: Node 3, DOF 6 (bend. moment around Z),-243,023Nmm)
5 2
2
0
5 3 2
0
7 1 1
+8101
7 2 1
-14031
7 3 1
-5107
7 4 1
+420930
7 5 1
-243030
9 2 2
0
9 3
2 0
The parameter file for the
stress processor Z88I3.TXT
can have any content (cf. sections 3.5 and 4.4), because Gauss points,
radial
and tangential stresses as well as calculation of the von Mises
stresses has no
significance for Cam Elements No.5.
CAD
and editor:
Because now the structure data Z88I1.TXT, the boundary conditions
Z88I2.TXT and
the header file for the stress processor Z88I3.TXT (with any content)
do exist,
you can launch
>Z88F Cholesky
solver for
computing the deflections
>Z88D
stress processor
>Z88E
nodal force processor
5.3.2
Results
The Cholesky solver Z88F
provides the following output files:
Z88O0.TXT stores the processed structure
data. For documentation purposes.
Z88O1.TXT stores the processed boundary conditions: For
documentation
purposes.
Z88O2.TXT, the displacements, the main task and solution of the
FEA
problem.
The stress processor Z88D internally uses the calculated
displacements
from Z88F and stores Z88O3.TXT, the calculated stresses. The
results in
Z88O3.TXT do not depend on the header parameters in Z88I3.TXT for Cam
Elements
No.5.
The nodal force processor Z88E internally uses the calculated
deflections of Z88F and stores Z88O4.TXT, the computed nodal
forces. Keep
in mind, that the "forces" of the DOF 4, 5 and 6 are really moments,
because the DOF 4, 5 and 6 are rotations.
The following pictures of
the plot program show the deflected structure for FUX, FUY and FUZ =
1,000 each
(magnifications of the deflections):
View of XY pane, deflected with nodal numbers
View of XZ pane, deflected with coarse plot of the Z deflections (Button Deflections Z). Compare this to the exact values in Z88O2.TXT.