5.1
FORK WRENCH WITH PLANE STRESS ELE. NO.7
Copy the example files B1_*
into Z88 entry files Z88* (has been already carried out on the Z88 CDs or
Internet packages for your immediate start):
B1_X.DXF ---> Z88X.DXF
input file for CAD converter Z88X
B1_2.TXT ---> Z88I2.TXT boundary conditions for Cholesky solver Z88F
B1_3.TXT ---> Z88I3.TXT header parameters for stress processor Z88D
Simply
proceed with the following steps to get familiar with Z88:
CAD:
As for this
first example, you should only look at the CAD super 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 super 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.
Z88:
Z88X,
conversion from Z88X.DXF to Z88NI.TXT. Windows: Compute > Z88X
> Type Conversion > 6 from Z88X.DXF to Z88NI.TXT > Compute > Go,
UNIX: pushbutton DXF <-> Z88 with radiobutton DXF ->
NI (Z88-Commander) or z88x -nifx (console or X-term).
Z88P, looking at the super structure.
The error message should not disturb you, because Z88P still has no .STO file
and therefore expects Z88I1.TXT to load (which is not yet present). However,
you want to load Z88NI.TXT: Windows: Plot > Z88P > File
> Structure File > Z88NI.TXT, UNIX: with the Z88-Commander
pushbutton Plot feature and radiobutton Z88P or enter from an
X-term z88p, enter z88ni.txt in textfield Struc., press
Return. Now try Z88O.
Z88N, mesh generator, reads
Z88NI.TXT and produces Z88I1.TXT. Windows: Compute > Z88N >
Compute > Go. UNIX: pushbutton Z88N (Z88-Commander) or z88n
(console or X-term).
Z88P, look at the finite element
structure. Z88P now has a Z88P.STO file from the former Run and expects
Z88NI.TXT to load. Now you want to load Z88I1.TXT. Windows: Plot >
Z88P > File > Structure File > Z88I1.TXT, UNIX:, with the
Z88-Commander pushbutton Plot feature and radiobutton Z88P or
enter from an X-term z88p, enter z88i1.txt in textfield Struc.,
press Return. This could have been avoided if you would have deleted the file
Z88P.STO before start of Z88P, because Z88P.STO stores all parameters of the
last Z88P run. If Z88P.STO does not exist, Z88P works with default values, for
example Z88I1.TXT as structure file.
Z88F, calculates displacements. Use the
Compute mode: Windows: Compute > Z88F > Mode > Compute mode,
> Compute > Go, UNIX: pushbutton SOLVER with radiobutton Z88F -C
(Z88-Commander) or z88f -c (console or X-term).
Z88D, calculates stresses. Windows:
Compute > Z88D > Compute > Go, UNIX: pushbutton Z88D
(Z88-Commander) or z88d (console or X-term).
Z88P, look at the deflected finite
element structure. Z88P has a .STO file and thus features correct parameters.
The displacements are multiplied per default by the factor 100, which is a bit
too large for this example. Windows: Plot > Z88P > Factors >
Deflections > enter 10 for FUX and FUY each,
> Structure > Deflected UNIX: Z88-Commander: pushbutton Plot
feature and radiobutton Z88P or enter from an X-term z88p,
enter 10 into both textfields FUX and FUY, press either a
Return for each textfield or pushbutton Regen. Click radiobutton Deflected.
Because Z88D has been already run, you can display the von Mises stresses right
now. Go to Undeflected structure. Windows: > Mises Stresses
> Show von Mises stresses. UNIX: Togglebutton Mises .
Furthermore you could produce a plotter file, perhaps undeflected and without
von Mises stresses. Windows: > Output > Plotter. UNIX:
Pushbutton Plot. The resulting plot file Z88O6.TXT contains HP-GL
commands.
Z88E, nodal forces calculation. Windows:
Compute > Z88E > Compute > Go, UNIX: pushbutton Z88E
(Z88-Commander) or enter z88e from a console or X-term.
Your
task:
A fork wrench should be loaded with the screw's tightness torque. A couple of
forces are applied in the wrench's mouth according to the torque and the fixed
points are assumed to be at the locations where the mechanic's hand grips the
wrench. In fact, these clever boundary conditions are doing the same task as
(in reality !) the fixed points in the mouth and the forces applied to the
grip, but are much easier to handle.
The fork wrench should be
modelled by 7 super elements Plane Stress No.7. The mesh generator should produce 66 finite
elements from the super elements. The element thickness is 10 mm each. Mesh
generation: Local and global axises are not the same direction in this example:
Local x direction at super element 1 defines by the local nodes 1 and 2 which
correspond to the global nodes 1 and 3. The local y direction of SE 1 is determined
by local nodes 1 and 4 which correspond to the global nodes 1 and 7. Further
take into account: Super elements which have a joint side must have an
absolutely identical subdivision at this side. Thus, SE 1 and SE 2 share the
line 3-4-5: The subdivisions in y direction must be exactly the same. Here 3
subdivisions, respectively.
Now calculate this example
as indicated above. After that, one can experiment: Subdivide the SE 7 in Z88NI.TXT as a meaningful variation as
follows:
7 7 ("Super element 7 is of type 7, i.e.
Plane Stress Element No.7")
6 L 3 E (
"Subdivide SE 7 into finite elements Plane Stress No.7 and subdivide into
x direction 6 times geometrically ascending and in y direction 3 times
equidistant") Of course, the SE 1 to SE 5 as well could each be
condensed in direction of the screw:
1 7
3 L 3 E
2 7
3 L 3 E
.... continue ....
Note: As it is obvious for the input
files, you can add comments after all required data are entered in every line.
Separate the last date from the comment by at least one blank. You can do this
just the same in your own files. A maximum of altogether 250 characters per
line is permitted.
5.1.1
Input
This example works
with a super structure, i.e. a very rough FE mesh. The mesh generator should
generate a FE structure from the super structure. Thus, the first task is to
design the mesh generator input file Z88NI.TXT. Chapter 2.6 outlines the procedure
if working with CAD. If you work without a CAD system, you design the file
Z88NI.TXT by editor or word processing program. The super structure shall look
as follows:
With
CAD program:
Follow the description of chapter 2.7. Do not forget to write the super
element information on the layer Z88EIO by TEXT function. Thus
SE 1 7 7 3 E 3 E ( 1st
SE, SE type 7, FE type 7, subdivide into x 3 times equid., into y 3 times
equid. )
SE 2 7 7 3 E 3 E ( 2nd SE, SE type 7, FE type 7, subdivide into X 3 times
equid., into y 3 times equid. )
SE 3 7 7 3 E 3 E
SE 4 7 7 3 E 3 E
SE 5 7 7 3 E 3 E
SE 6 7 7 1 E 3 E
SE 7 7 7 6 E 3 E
...and write the general
information and material information onto the layer Z88GEN :
Z88NI.TXT 2 38 7 76 1 0 0 0
0 0 ( 2-DIM, 38 nodes, 7 superelements, 76 DOF, 1 mat info, flags 0 )
MAT 1 1 7 206000 0.3 3 10 ( 1.mat info from SE 1 to SE 7: Young's modulus,
Poisson's ratio, INTORD, thickness, )
Export the drawing as DXF
file with the name Z88X.DXF and start the CAD converter Z88X with the option "from Z88X.DXF
to Z88NI.TXT" (DXF -> NI). Z88X will produce the mesh generator input
file Z88NI.TXT. You should have a look at it with Z88P.
With
editor:
Write the mesh
generator input file Z88NI.TXT (cf. chapter 3.3) with an editor:
2 38 7 76 1 0 0 0 0 0 (2-DIM, 38 nodes, 7 superelements,
76 DOF, 1 mat info line, flags 0)
1 2 22.040 32.175
(Node 1, 2 DOF, X and Y coordinates)
2 2 31.913 28.798
(Node 2, 2 DOF, X and Y coordinates)
3 2 43.781 24.826
4 2 43.880 32.373
5 2 43.980 39.424
...... (Coordinates for nodes 6... 36 not represented)
37 2 202.847 27.507
38 2 144.905 42.403
1 7 (SE 1 of the type Plane Stress
No.7)
1 3 5 7 2 4 6 8 (Coincidence
for 1st SE)
2 7 (SE 2 of the type Plane
Stress No. 7)
3 10 12 5 9 11 13 4 (Coincidence
for 2nd SE)
..... (Coincidence for elements 3 .. 6 droped here)
7 7
30 35 37 32 34 36 38 31
1 7 206000 0.3
3 10 (mat info from SE 1 to SE
7:Young,Poisson,INTORD,thickness)
1 7 (Subdivide 1st SE into FE type 7
and
3 E 3 E subdivide
into x 3 times equidistant + into y 3 times equidistant)
2 7 (Subdivide 2nd SE into FE type 7
and
3 E 3 E subdivide
into x 3 times equidistant + into y 3 times equidistant)
3 7
3 E 3 E
4 7
3 E 3 E
5 7
3 E 3 E
6 7
1 E 3 E
7 7
6 E 3 E
With
CAD program and editor:
Start the mesh generator
Z88N for producing
the final Z88 structure file Z88I1.TXT. Look at it either
Enlarge the wrench's mouth
by zooming for defining the two nodes which will get the load representing the
torque (to simplify matters it is assumed, that the screw gets only selectively
a couple of forces as torque at the corners and that the screw itself and not
the wrench revolves):
We find the nodes 11 and
143. The pictures printed here were produced directly by Z88P.
In the same way both the nodes for fixing the wrench are determined and the
boundary conditions are entered in the plot program or CAD system:
In
the CAD program:
Switch to the
layer Z88RBD and write with the TEXT function into any free place:
Z88I2.TXT 16 (16
Boundary conditions altogether)
RBD 1 11 2 1 -7143 (1st BC: Node 11, DOF 2, Force -7,143 N assumed)
RBD 2 143 2 1 7143 (2nd BC: Node 143, DOF 2, Force 7,143 N assumed)
RBD 3 216 1 2 0 (3rd BC: Node 216, DOF 1, Displacement 0 (= fixed) assumed)
RBD 4 216 2 2 0
RBD 5 220 1 2 0
RBD 6 220 2 2 0
RBD 7 227 1 2 0
RBD 8 227 2 2 0
RBD 9 231 1 2 0
RBD 10 231 2 2 0
RBD 11 238 1 2 0
RBD 12 238 2 2 0
RBD 13 242 1 2 0
RBD 14 242 2 2 0
RBD 15 249 1 2 0
RBD 16 249 2 2 0
with
an Editor:
Design the
boundary condition file Z88I2.TXT by editing:
16 (16 Boundary
conditions altogether)
11 2 1 -7143 (1st BC: Node 11, DOF 2, Force -7,143
N assumed)
143 2 1 7143
(2nd BC: Node 143, DOF 2, Force 7,143 N assumed)
216 1 2 0 (3rd BC: Node 216, DOF 1, Displacement 0
(= fixed) assumed)
216 2 2 0
220 1 2 0
220 2 2 0
227 1 2 0
227 2 2 0
231 1 2 0
231 2 2 0
238 1 2 0
238 2 2 0
242 1 2 0
242 2 2 0
249 1 2 0
249 2 2 0
Input for stress
calculation:
In
the CAD program:
Switch to the
layer Z88GEN and write with the TEXT function into any free place:
Z88I3.TXT 3 0 1 ( 3x3 Gauss points for stresses, KFLAG 0,
von Mises stresses)
Export the drawing as DXF
file with the name Z88X.DXF, then start the CAD converter Z88X with the option
"from Z88X.DXF to Z88I*.TXT" (DXF -> I*). The CAD converter
produces the three Z88 input files Z88I1.TXT, Z88I2.TXT, Z88I3.TXT.
With
an editor:
Enter in the
parameter file for the stress processor Z88I3.TXT (cf. Chapter 3.5):
3 0 1 ( 3x3 Gauss points for stresses, KFLAG 0,
von Mises stresses)
Now launch the Cholesky
solver Z88F and
then the stress processor Z88D. You will see during the run of Z88F, that 14.848 memory places (8
bytes each) are needed in the total stiffness matrix GS. NKOI, i.e. memory places in the coincidence
vector KOI, is
printed as 540 (4 Bytes each). Well, this also matches Z88.DYN. Where does the
number 540 come from? 66 finite elements of the type Plane Stress No.7 with 8
nodes each, makes 66*8 = 528. The number 540 results because Z88F always calculates
20 nodes for security reasons for the last finite element. Thus, NKOI becomes
here: 65*8 + 20 = 540.
You calculate the nodal
forces with Z88E.
5.1.2
Results
The Cholesky solver
Z88F provides the following output files:
Z88O0.TXT stores the processed structure data. It is mainly intended for
documentation purposes, but also shows if your input file Z88NI.TXT for the
mesh generator did what you meant it to do.
Z88O1.TXT stores the processed boundary conditions: For documentation
purposes. And: Was your boundary conditions input in Z88I2.TXT correctly
interpreted ?
Z88O2.TXT, the displacements, the main task and solution of the FEA
problem.
The stress processor Z88D uses internally the calculated displacements
from Z88F and stores Z88O3.TXT, the calculated stresses. The results in
Z88O3.TXT depend on the header parameters in Z88I3.TXT.
The following picture of
the plot program shows the deflected structure for FUX and FUY = 30 each
(magnifications of the deflections):
The nodal force processor Z88E
uses internally the calculated deflections of Z88F and stores Z88O4.TXT,
the computed nodal forces.
The OpenGL plot program Z88O shows the average von Mises stresses per node. These stresses are usually somewhat
lower than the stresses computed in the Gauss points shown by the other plot program
Z88P.