5.4 BEAM in PLANE WITH BEAMS NO.13

Copy the example file B4_X.DXF to Z88X.DXF.

B4_X.DXF ---> Z88X.DXF input file for CAD converter Z88X

CAD:

Import Z88X.DXF into your CAD program and look at it. Usually you would have designed this example in a CAD system (makes not much sense because this example is extremely simple) and then exported as Z88X.DXF.

Z88: (in reduced form, more detailed instructions cf. examples 5.1, 5.2 and 5.3)

Z88X, conversion, "from Z88X.DXF to Z88I*.TXT"

Z88P, looking at structure, structure file Z88I1.TXT

Z88F calculates deflections

Z88D calculates stresses

Z88E calculates nodal forces

Z88P, plot FE structure, now also deflected (FUX, FUY, FUZ per 10.)

This example deals with a beam, fixed on both sides, and loaded with 1,648 N in the middle in downward direction. This mechanical problem is covered in every mechanical and civil engineering handbook. Geometry: Length 1,000 mm, cross-cut 50 x 10 mm. Thus: A = 500 mm**2, Izz = 4,167 mm**4, ezz = 5 mm.

The deflection curve has inflection points, we therefore take 4 beams No.13. Nodes 1 and 5 will be fixed and node 3 is loaded.

You would calculate analytically:

f in the middle: F*L**3/(192*E*I) = 10 mm

f in the inflection points: fw = f/2 = 5 mm

The bending moments on the left, middle, on the right: F*L/8 = 206,000 Nmm

The slope angle in the inflection points: Phi = atan (3*f/L) = 0.029991 rad

When interpreting the results of Z88O2.TXT (deflections) and Z88O4.TXT (nodal forces and moments) refer to the sign definition of chapter 3.13. Especially Z88O4.TXT, node 3: The force F(2) = force in Y direction is the sum of the forces of elements 2 and 3, due to extrinsic force. The force F(3) = bending moment is not a summary of elements 2 and 3, because it is an intrinsic moment, not an extrinsic load ! Also the signs of the load F(3) at node 1 and F(3) at node 5 are correct, refer to chapter 4.13. Keep in mind that the classical mechanical science sometimes uses different conventions.

5.4.1 Input:

This example shows that a FEA basically needs nodes in all locations where you want to get results. As the beam is fixed left and right, the maximum of displacements appears in the middle for x = L/2, but the bending curve features two inflection points for x = L/4 and x = 3L/4. To calculate results for this locations, the structure must be subdivided with nodes in x = 0, x = L/4, x = L/2 and x = 3L/4.

Only the file input is shown here because CAD use is not worth here.

Z88I1.TXT so becomes:

2 5 4 15 1 0 1 0 (2-D,5 nodes,4 ele, 5 DOF,1 mat info, KFLAG 0, IBFLAG , IPFLAG 01)

1 3 0 0 (1.node, 3 DOF, X and Y coordinate)

2 3 250 0

3 3 500 0

4 3 750 0

5 3 1000 0

1 13 (1. element, type beam in plane No.13)

1 2 (coincidence for 1. element)

2 13

2 3

3 13

3 4

4 13

4 5

1 4 206000 0.3 1 500 0 0 4167 5 0 0 (mat info for ele 1 to 4, Young's, Poisson's, INTORD (any), QPARA = area, Ixx=0, exx=0, Izz, ezz, It=0, Wt=0)

The node 1 is fixed in all degrees of freedom at the boundary conditions. It is important to fix especially the DOF 1 = displacement in X direction so that the structure cannot move. Node 5 is fixed in DOF 2 = displacement in Y direction and DOF 3 = rotation around Z axis. You could also fix DOF 1 for node 5, if you wish. But in reality one of the bearings or supports will allow for thermal expansion. This was taken into account in Z88I2.TXT.

Here Z88I2.TXT:

6 (6 Boundary conditions)

1 1 2 0 (Node 1, DOF 1 gets a displacment of 0 = DOF 1 fixed)

1 2 2 0 (Node 1, DOF 2 fixed)

1 3 2 0 (Node 1, DOF 3 fixed (restraining moment)

3 2 1 -1648 (Node 3, DOF 2 gets load of -1,648 N)

5 2 2 0

5 3 2 0

The parameter file for the stress processor Z88I3.TXT can have any content (cf. sections 3.5 and 4.13), because Gauss points, radial and tangential stresses as well as calculation of the von Mises stresses has no significance for Beams No.13.

5.4.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 Beams No.13.

The nodal force processor Z88E internally uses the calculated deflections of Z88F and stores Z88O4.TXT, the computed nodal forces.

The following picture of the plot program shows the deflected structure for FUX, FUY and FUZ = 10 each (magnifications of the deflections).

Attention to the results of the nodal force calculation: Node 3: The force F(2) = force in Y direction is the sum of the forces of elements 2 and 3, due to extrinsic force. The force F(3) = bending moment is not a summary of elements 2 and 3, because it is an intrinsic moment, not an extrinsic load ! Also the signs of the load F(3) at node 1 and F(3) at node 5 are correct, refer to chapter 4.13. Keep in mind that the classical mechanical science sometimes uses different conventions.

Additional remark: Such simple examples are well suitable to become aware of the sign omen definitions. Experiment with this example and calculate other bend cases from good handbooks. Frameworks with Beams No.2 are calculated accordingly. However, a real spatial structure then must be available: At least one Z coordinate must not equal 0.

View of undeflected and deflected structure

Take into account: The plot program Z88P connects the nodes with straight lines, although the deflection curve represents a cubic parable in the case of a Beam No.13 or No.2. This means: Z88P shows the deformations correctly for the node, but straight lines are between the nodes. Therefore, no deflection curve is shown. If you want to plot a real nice deflection curve with Z88P, then use basically more nodes, e.g. 15 to 20 nodes for this example (the cubic bending curve is then featured by a couple of straight lines).