A light-service, cargo crane is shown in Figure 6–11. You have been asked to determine the static deflections of the crane when it carries a load of 10 kN. You should also identify the critical members and joints in the structure: i.e., those with the highest stresses and loads. Because this is a static analysis you will analyze the cargo crane using Abaqus/Standard.
The crane consists of two truss structures joined together by cross bracing. The two main members in each truss structure are steel box beams (box cross-sections). Each truss structure is stiffened by internal bracing, which is welded to the main members. The cross bracing connecting the two truss structures is bolted to the truss structures. These connections can transmit little, if any, moment and, therefore, are treated as pinned joints. Both the internal bracing and cross bracing use steel box beams with smaller cross-sections than the main members of the truss structures. The two truss structures are connected at their ends (at point E) in such a way that allows independent movement in the 3-direction and all of the rotations, while constraining the displacements in the 1- and 2-directions to be the same. The crane is welded firmly to a massive structure at points A, B, C, and D. The dimensions of the crane are shown in Figure 6–12. In the following figures, truss A is the structure consisting of members AE, BE, and their internal bracing; and truss B consists of members CE, DE, and their internal bracing.
The ratio of the typical cross-section dimension to global axial length in the main members of the crane is much less than 1/15. The ratio is approximately 1/15 in the shortest member used for internal bracing. Therefore, it is valid to use beam elements to model the crane.
You should use the default global rectangular Cartesian coordinate system shown in Figure 6–11 and Figure 6–12. Locate the origin of the coordinate system midway between points A and D. If you build your model with a different origin or orientation of the coordinate system, ensure that the input data in your model reflect your coordinate system and not the one shown here.
The cargo crane will be modeled with three-dimensional, slender, cubic beam elements (B33). The cubic interpolation in these elements allows us to use a single element for each member and still obtain accurate results under the applied bending load. The mesh that you use in the simulation is shown in Figure 6–13.
The welded joints in the crane provide complete continuity of the translations and rotations from one element to the next. You, therefore, need only a single node at each welded joint in the model. The bolted joints, which connect the cross bracing to the truss structures, and the connection at the tip of the truss structures are different. Since these joints do not provide complete continuity for all nodal degrees of freedom, separate nodes are needed for each element at the connection. Appropriate constraints between these separate nodes must then be given by using the *MPC, *BOUNDARY, or *EQUATION options. The *MPC and *EQUATION options are discussed in more detail later.
The node numbers for the various members of the cargo crane model are shown in Figure 6–14.
These are the node numbers from the input file given in “Cargo crane,” Section A.4. Separate nodes have been defined on the cross-bracing elements and the truss structures that they connect. Separate nodes are also needed at the end of each truss structure, point E in Figure 6–11. The node numbers in your model may be different from those shown here.The element numbers for the various members of the cargo crane model are shown in Figure 6–15.
These are the element numbers from the input file given in “Cargo crane,” Section A.4. The element numbers in your model may be different from those shown here.The full input file for this example is crane.inp, and it is available in “Cargo crane,” Section A.4. The Abaqus input options used to create the nodes and elements shown on the preceding pages can be found in “Cargo crane,” Section A.4. If you wish to create the entire model using Abaqus/CAE, refer to “Example: cargo crane,” Section 6.4 of Getting Started with Abaqus: Interactive Edition.
This section describes how the model data are described in the input file for this example. These data include the descriptions for the input file heading, its nodes and elements, beam sections and orientations, constraints, and boundary conditions.
Heading
The heading used in this example provides a short description of the model and the units used:
*HEADING 3-D model of light-service cargo crane S.I. Units (m, kg, N, sec)
Nodal coordinates and element connectivity
Define the nodal coordinates in a *NODE option block. If you decide to do this with an editor, you may want to use the mesh generation commands found in “Cargo crane,” Section A.4. In this example, a node set called ATTACH is created; this node set contains the nodes at points A, B, C, and D, the points at which the crane is attached to the parent structure.
*NSET, NSET=ATTACH 100, 107, 200, 207
Create elements in your model that correspond to the elements shown in Figure 6–15, but remember that your numbering may be different. (Having the same numbering will make defining some modeling features easier, however.) There are several element sets that you will need in this simulation. As the elements are defined, they are grouped into following element sets:
OUTA 100, 101, 102, 103, 104, 105, 106 BRACEA 110, 111, 112, 113, 114 OUTB 200, 201, 202, 203, 204, 205, 206 BRACEB 210, 211, 212, 213, 214 CROSSEL 300, 301, 302, 303, 304, 305, 306, 307where these element numbers refer to those elements shown in Figure 6–15. The element sets OUTA and OUTB contain the main outer members for the two truss structures. Element sets BRACEA and BRACEB contain the elements modeling the internal bracing within each truss structure. Element set CROSSEL contains the cross bracing that connects the two truss structures.
Beam element properties
Since the material behavior in this simulation is assumed to be linear elastic, use the *BEAM GENERAL SECTION option to define the section properties. All of the beams in this structure have a box-shaped cross-section.
A box-section is specified using the parameter SECTION=BOX. The first data line contains the section dimensions, which are the dimensions a, b, ,
,
, and
shown in Figure 6–16 for a box section. The dimensions shown in Figure 6–16 are for the main members of the two trusses in the crane.
The beam section axes for the main members should be oriented such that the beam 1-axis is orthogonal to the plane of the truss structures shown in the elevation view (Figure 6–12) and the beam 2-axis is orthogonal to the elements in that plane. Specify this orientation by giving the approximate direction of the beam 1-axis (the -vector) on the second data line of the *BEAM GENERAL SECTION option. To get the correct normal,
, in this case, you need to provide a very accurate
. It is somewhat easier to provide an approximate direction, which would be the negative 3-direction. However, given the logic that Abaqus uses to determine
given
and
, the normal
is rotated slightly from its proper orientation if we use this approximate
. You can specify the same
-direction for all elements in each of the two truss structures. The third data line contains the elastic and shear moduli, assuming a mild strength steel with
= 200.0 GPa,
= 0.25, and
= 80.0 GPa. These modeling data are included in the input file in the following option blocks:
*BEAM GENERAL SECTION, SECTION=BOX, ELSET=OUTA 0.10,0.05,0.005,0.005,0.005,0.005 -0.1118, 0.0, -0.9936 200.E9,80.E9 *BEAM GENERAL SECTION, SECTION=BOX, ELSET=OUTB 0.10,0.05,0.005,0.005,0.005,0.005 -0.1118, 0.0, 0.9936 200.E9,80.E9
The dimensions of the beam sections for the bracing members are shown in Figure 6–17. Both the cross bracing and the bracing within each truss structure have the same beam section geometry, but they do not share the same orientation of the beam section axes. Therefore, separate *BEAM GENERAL SECTION options must be used. The bracing is made from the same steel as the main members.
The approximate -vector for the internal truss bracing is the same as for the main members of the respective truss structures. The following input defines the element properties of this bracing:
*BEAM GENERAL SECTION, SECTION=BOX, ELSET=BRACEA 0.03,0.03,0.003,0.003,0.003,0.003 -0.1118, 0.0, -0.9936 200.E9,80.E9 *BEAM GENERAL SECTION, SECTION=BOX, ELSET=BRACEB 0.03,0.03,0.003,0.003,0.003,0.003 -0.1118, 0.0, 0.9936 200.E9,80.E9
We make some assumptions so that the orientation of the cross bracing is somewhat easier to specify. All of the beam normals (-vectors) should lie approximately in the plane of the plan view of the cargo crane (see Figure 6–12). This plane is skewed slightly from the global 1–3 plane. Again, a simple method for defining such an orientation is to provide an approximate
-vector that is orthogonal to this plane on the element property options. The vector should be nearly parallel to the global 2-direction. Since the angle between the normals from one element to the next is always greater than 20°, the normals are not averaged at the nodes.
Depending on the exact orientation of the members in the cross bracing, it is possible that we would have to define the normals individually for each of the cross-bracing elements. Such an exercise would be very similar to what you have already done by defining the normals for the two truss structures. Since the square cross-bracing members are subjected to primarily axial loading, their deformation is not sensitive to cross-section orientation. Therefore, we accept the default normals that Abaqus calculates to be correct. The approximate -vector for the cross bracing is aligned with the y-axis. The following option block specifies the cross-bracing:
*BEAM GENERAL SECTION, SECTION=BOX, ELSET=CROSSEL 0.03,0.03,0.003,0.003,0.003,0.003 0.0,1.0,0.0 200.E9,80.E9
Beam section orientation
In this model you will have a modeling error if you provide data that only define the orientation of the approximate -vector. The averaging of beam normals at the nodes (see “Beam element curvature,” Section 6.1.3) causes Abaqus to use incorrect geometry for the cargo crane model. To see this, you can use Abaqus/Viewer to display the beam section axes and beam tangent vectors (see “Postprocessing,” Section 6.4.7). However, the normals in the crane model appear to be correct in Abaqus/Viewer; yet, they are, in fact, slightly incorrect. You can also find such modeling mistakes by examining the averaged nodal normals that are printed in the data (.dat) file. Some of the normals in the incorrect model of the cargo crane are shown in the following output:
The problem is the normal vector for node 102, which does not match those at the other nodes defining the lower, main member in truss A (see Figure 6–14). Four elements (101, 102, 112, and 113) contain node 102. When the averaging of beam normals at nodes produces multiple independent normals, the additional normals at the node are also printed in the data file (see “Beam element curvature,” Section 6.1.3, for details). The correct geometry for the crane model requires three independent beam normals at node 102: one each for the bracing elements 112 and 113 and a single normal for elements 101 and 102. The normal shown above for node 102 is not the normal needed for elements 101 and 102. If it were, it would match the normals shown for nodes 100, 101, or 103. Nor is it the correct normal for either of the other two elements, whose normals are printed in the data file, as shown below.
N O R M A L D E F I N I T I O N S ELEMENT NODE NORMAL ELEMENT NODE NORMAL 100 101 -0.1820 0.9831 2.0481E-02 101 101 -0.1820 0.9831 2.0481E-02 101 102 -0.1820 0.9831 2.0481E-02 102 102 -0.1820 0.9831 2.0481E-02 102 103 -0.1820 0.9831 2.0481E-02 103 103 -0.1820 0.9831 2.0481E-02 103 104 -0.1820 0.9831 2.0504E-02 104 105 6.1600E-02 -0.9981 -6.9312E-03In this table element 113 shows no normal for either of its nodes, node 102 and node 105. Thus, the normal shown for node 102 above in the NODE DEFINITIONS table was the average of the normals from elements 101, 102, and 113. Using the Abaqus logic for averaging normals, we could have predicted that the normals at the nodes of element 113 would be averaged with the normals for the adjacent elements. For this problem the important part of the averaging logic is that normals that subtend an angle less than 20° with the reference normal are averaged with the reference normal to define a new reference normal. In the case of the normals at node 102, the original reference normal is the normal for elements 101 and 102. Since the normal for element 113 at node 102 subtends an angle less than 20° with the original reference normal, it is averaged with the normals for elements 101 and 102 at node 102 to define the new reference normal at that node. On the other hand, since the normal for element 112 subtends an angle of approximately 30° with the original reference normal, it has an independent normal at node 102, as shown in the data file.
This incorrect average normal means that elements 101, 102, and 113 have a section geometry that exhibits curvature about the beam axis from one end of the element to the other, which is not the intended geometry. You must use the *NORMAL option to define the normals at node 102 for element 113 explicitly. Explicitly specifying the normal directions prevents Abaqus from applying its averaging algorithm. You must also use *NORMAL for the corresponding elements on the opposite side of the crane: element 213, node 202 of truss B.
There is also a problem with the normals at nodes 104 and 204 at the tip of each truss structure, again because the angle between elements 103 and 104 is less than 20°. Since we are modeling straight beam elements, the normals are constant at both nodes in each element. Thus, the *NORMAL option block that you should put in your input file has six data lines. If you used the numbering scheme shown in Figure 6–14 and Figure 6–15, the following option block should be added to your input file:
Multi-point constraints
The cross bracing, unlike the internal truss bracing, is bolted to the truss members. You can assume that these bolted connections are unable to transmit rotations or torsion. The duplicate nodes that were defined at these locations are needed to define this constraint. In Abaqus such constraints between nodes can be defined using multi-point constraints (MPCs) or constraint equations.
Multi-point constraints allow constraints to be imposed between different degrees of freedom of the model. A large library of MPCs is available in Abaqus. (See “Linear constraint equations,” Section 35.2.1 of the Abaqus Analysis User's Guide, for a complete list and a description of each one.) The format of the *MPC option is
*MPC <type of MPC>,<node 1 or node set 1>,<node 2 or node set 2>, ......
You can define multiple constraints of the same type with just a single data line by using node sets. The MPC type needed to model the bolted connection is PIN. The pinned joint created by this MPC constrains the displacements at two nodes to be equal, but the rotations, if they exist at the nodes, remain independent.
There are many bolted joints in the crane model. The following is the complete *MPC option block for the model from “Cargo crane,” Section A.4:
*MPC PIN,101,301 PIN,102,302 PIN,103,303 PIN,105,305 PIN,106,306 PIN,201,401 PIN,202,402 PIN,203,403 PIN,205,405 PIN,206,406
Add a similar option block to your model, changing the node numbers to correspond to those in your model. If all of the nodes on the two truss structures had been grouped into a node set called TRUSNODE and all of the nodes on the cross bracing had been grouped into a node set called CROSNODE, the option block could have been shortened to the following:
*NSET, NSET=TRUSNODE, UNSORTED 101,102,103,105,106,201,202,203,205,206 *NSET, NSET=CROSNODE, UNSORTED 301,302,303,305,306,401,402,403,405,406 *MPC PIN, TRUSNODE, CROSNODE
If a node set is provided as the first item after the MPC type, the second item can be either another node set or a single node. When the data line of an *MPC option contains a node set and then a single node, as shown below, Abaqus creates an MPC constraint between each node in the set and the individual node specified. For example, the following option block would create a pinned joint between node 301 and each node in node set TRUSNODE.
*MPC PIN, TRUSNODE, 301
Constraint equations
Constraints between nodal degrees of freedom can also be specified with linear equations by using the *EQUATION option. The form of each equation is
<Exactly four terms must be given on each data line except the final one, which can have fewer terms.>,<
>,<
>,<
>,<
>,<
>,
,<
>,<
>,<
>
In the crane model the tips of the two trusses are connected together such that degrees of freedom 1 and 2 (the translations in the 1- and 2-directions) of each tip node are equal, while the other degrees of freedom (3–6) at the nodes are independent. We need two linear constraints, one equating degree of freedom 1 at node 104 to degree of freedom 1 at node 204:
You may have to change the node numbers if you created this model using a preprocessor. The following option block defines the appropriate constraints at point E in the crane model (see Figure 6–11):
*EQUATION 2 104,1,1.0, 204,1,-1.0 2 104,2,1.0, 204,2,-1.0
The degrees of freedom at the first node defined in an *MPC or *EQUATION are eliminated from the stiffness matrix. Therefore, these nodes should not appear in other MPCs or constraint equations. Nor should boundary conditions be applied to the eliminated degrees of freedom.
Boundary conditions
The crane is attached firmly to the parent structure. The following *BOUNDARY option block constrains all of the nodes at the attachment points, which should have been grouped into node set ATTACH:
*BOUNDARY ATTACH, ENCASTRE
The following options specify a static, linear perturbation simulation:
*STEP, PERTURBATION Static tip load on crane *STATIC
Loading
A concentrated load of 10 kN is applied in the negative y-direction to node 104. Since there is a constraint equation connecting the y-displacement of nodes 104 and 204, the load is carried equally by both nodes. The following *CLOAD option block provides an equal load on both nodes:
*CLOAD 104,2,-1.0E4
Output requests
Write the displacements (U) and reaction forces and moments (RF) at the nodes, and the section forces and moments (SF) in the elements to the output database for postprocessing with Abaqus/Viewer, as shown in the following option blocks:
*OUTPUT, FIELD *NODE OUTPUT U, RF *ELEMENT OUTPUT SF, *END STEP
Store the input in a file called crane.inp. Run the analysis using the command
abaqus job=crane
Start Abaqus/Viewer by typing the following command at the operating system prompt:
abaqus viewer odb=craneAbaqus/Viewer plots the undeformed shape of the crane model.
Plotting the deformed model shape
To begin this exercise, plot the deformed model shape with the undeformed model shape superimposed on it. Specify a nondefault view using (0, 0, 1) as the X-, Y-, and Z-coordinates of the viewpoint vector and (0, 1, 0) as the X-, Y-, and Z-coordinates of the up vector.
Tip:
You can also display the model using this view by clicking from the Views toolbar.
The undeformed shape of the crane superimposed upon the deformed shape is shown in Figure 6–18.
Using display groups to plot element and node sets
You can use display groups to plot existing node and element sets; you can also create display groups by selecting nodes or elements directly from the viewport. You will create a display group containing only the elements associated with the main members in truss structure A.
To create and plot a display group:
In the Results Tree, expand the Sections container underneath the output database file named crane.odb.
To facilitate your selection, change the view back to the default isometric view using the tool in the Views toolbar.
Tip:
If the Views toolbar is not visible, select ViewToolbars
Views from the main menu bar.
In succession, click the items in the container until the elements associated with the main members in truss A are highlighted in the viewport. Click mouse button 3 on this item and select Replace from the menu that appears.
Abaqus/Viewer now displays only this group of elements.
To save this group, double-click Display Groups in the Results Tree; or use the tool in the Display Group toolbar.
The Create Display Group dialog box appears.
In the Create Display Group dialog box, click Save As and enter MainA as the name for your display group.
Click Dismiss to close the Create Display Group dialog box.
This display group now appears underneath the Display Groups container in the Results Tree.
Beam cross-section orientation
You will now plot the section axes and beam tangents on the undeformed model shape.
To plot the beam section axes:
From the main menu bar, select PlotUndeformed Shape; or use the
tool in the toolbox to display only the undeformed model shape.
From the main menu bar, select OptionsCommon; then, click the Normals tab in the dialog box that appears.
Toggle on Show normals, and accept the default setting of On elements.
In the Style area at the bottom of the Normals page, specify the Length to be Long.
Click OK.
The section axes and beam tangents are displayed on the undeformed shape.
Rendering beam profiles
You will now display an idealized representation of the beam profile and contour the stress results.
To render beam profiles:
From the main menu bar, select ViewODB Display Options.
The ODB Display Options dialog box appears.
In the General tabbed page, toggle on Render beam profiles and accept the default scale factor of 1.
Click OK.
Abaqus/Viewer displays beam profiles with the appropriate dimensions and in the correct orientations. Figure 6–20 shows the beam profiles on the whole model. Your changes are saved for the duration of the session.
Click to contour the Mises stress on the rendered profile.
Creating a hard copy
You can save the current image to a file for hard-copy output.
To create a PostScript file of the current image:
From the main menu bar, select FilePrint.
The Print dialog box appears.
From the Settings area in the Print dialog box, select Black&White as the Rendition type; and toggle on File as the Destination.
Select PS as the Format, and enter beam.ps as the File name.
Click .
The PostScript Options dialog box appears.
From the PostScript Options dialog box, select 600 dpi as the Resolution; and toggle off Print date.
Click OK to apply your selections and to close the dialog box.
In the Print dialog box, click OK.
Abaqus/Viewer creates a PostScript file of the current image and saves it in your working directory as beam.ps. You can print this file using your system's command for printing PostScript files.
Displacement summary
Write a summary of the displacements of all nodes in display group MainA to a file named crane.rpt. The peak displacement at the tip of the crane in the 2-direction is 0.0188 m.
Section forces and moments
Abaqus can provide output for structural elements in terms of forces and moments acting on the cross-section at a given point. These section forces and moments are defined in the local beam coordinate system. Toggle off the rendering of beam profiles, then contour the section moment about the beam 1-axis in the elements in display group MainA. For clarity, reset the view so that the elements are displayed in the 1–2 plane.
To create a “bending moment”-type contour plot:
From the list of variable types on the left side of the Field Output toolbar, select Primary.
From the list of output variables in the center of the toolbar, select SM.
Abaqus/Viewer automatically selects SM1, the first component name in the list on the right side of the Field Output toolbar, and displays a contour plot of the bending moment about the beam 1-axis on the deformed model shape. The deformation scale factor is chosen automatically since geometric nonlinearity was not considered in the analysis.
Open the Common Plot Options dialog box, and select a Uniform deformation scale factor of 1.0.
Color contour plots of this type typically are not very useful for one-dimensional elements such as beams. A more useful plot is a “bending moment”-type plot, which you can produce using the contour options.
From the main menu bar, select OptionsContour; or use the Contour Options
tool in the toolbox.
The Contour Plot Options dialog box appears; by default, the Basic tab is selected.
In the Contour Type field, toggle on Show tick marks for line elements.
Click OK.
The plot shown in Figure 6–21 appears. The magnitude of the variable at each node is now indicated by the position at which the contour curve intersects a “tick mark” drawn perpendicular to the element. This “bending moment”-type plot can be used for any variable (not just bending moments) for any one-dimensional element, including trusses and axisymmetric shells as well as beams.