Design Optimization of Beam Elements

Given: The bucket of the front-end loader consists of two identical mechanisms. Each mechanism carries a load of 12 kN (24 kN total for the front-end loader) as shown in the figure below. Include the live load and gravity in the analysis.

Find: A preliminary design with:

  1. the optimum cross sections for the three links AB, BCE, and DEFG, based on minimizing the weight of the structure and limiting the stress to 100 MPa
  2. the size required for cylinders CD and FH.

Figure 1

 

Since we are after a preliminary design and only considering in-plane effects, the shape of the cross-section is not overly important. Thus, beam elements with a rectangular cross section will be used. (Once the solution is found, a different shape with similar area and section modulus can be substituted, if needed.)

The optimization procedure is performed in the following basic steps:

You can create the geometry and input in either of two ways:

Part One – Build the Model:

The model will be created with beam elements for the arms (since they need to resist bending) and truss elements for the cylinder (since they will only transmit an axial force). The initial sections are as follows:

Link Part Number Surface Number Element Type Initial Cross Section Material
       

width out-of-plane (b) x height in-plane (h)

 
DEFG 1 2 Beam 24 x 250 mm AISI 1005 Steel
Bucket 2 3 Beam

rigid but small area

AISI 1005 Steel
AB 3 2 Beam 14 x 150 mm AISI 1005 Steel
BCE 4 3 Beam 20 x 150 mm AISI 1005 Steel
CD and FH 5 any Truss

1900 mm 2

AISI 1020 Steel, cold rolled
Table 1: FEA Properties

                                                     

Different part numbers will be used as a convenience for displaying or hiding each link in the Results environment. To simulate the pinned connections between the beams, end releases will be used on the bucket (at points A and D) and link BCE (at points B and E).

The model is planned as in Figure 2. The orientation of the model can be in any direction since beam elements are three dimensional, but let's use the XZ orientation of the model shown in Figure 2. (This makes the Front view natural.) The different part numbers are shown (P) based on the above table.

For beam elements, the orientation of axis 2 needs to be decided to enter the cross sectional properties. The desired orientation scheme is as follows:

  • Axis 2 (weak axis of bending) lies in the plane of the model
  • Axis 3 (strong axis of bending) points out of the model plane (parallel to the global Y-axis).

Since the orientation of beam elements is controlled by the surface number of the lines, this dictates the surface numbers (S), shown in Table 1 and in Figure 2. (See the paragraph Beam Element Orientation on the page Setting Up and Performing the Analysis: Analysis-Specific Information: Linear: Element Types and Parameters: Beam Elements for details on orienting beam elements.)

Figure 2: The FEA Model

 Part numbers (P), axis 2 orientation decided by the user, and surface numbers (S) dictated by the orientation, are added to the figure.

Draw the Geometry:

  1. Start a new FEA model ( New) and
    1. Set the analysis type to Linear: Static Stress with Linear Material Models.
    2. Click the Override Default Units button.

      Since the dimensions are given in millimeters and the load in Newtons, lets choose units of Newtons, millimeters, and seconds for the model...

      • Set the Unit System drop-down to Custom.
      • Set the corresponding units of Force, Length, and Time to N, mm, and s, respectively.
      • Click the OK button to set the units.
    3. Click the New button and enter a model name.
  2. Start drawing the mechanism at point G. Add a line using Draw Draw Line.
    1. Remove the check from the Use as Construction check box.
    2. Set the Part number and Surface number to 1 and 2, respectively, per Figure 2.
    3. Start the line at (0,0,0), and then activate the Use Relative check box. Per Figure 1, put the next points at the following relative positions. (To enter a coordinate, put the cursor in the X, Y, or Z box (or DX, DY, DZ), type the value, and press the Enter key.)
      Point  Relative Distance From Previous Point
      F (-480,0,120)
      E (-420,0,240)
      D (-300,0,140)
      Table 2

      Keep the Define Geometry dialog box open. We will be creating additional lines.

  3. The lines may be off of the screen at this point. Use View Navigate Orientation Front View to re-orient the model view and to fit it to the display area.
  4. On the add line dialog, set the Part number and Surface number to 2 and 3, respectively, per Figure 2.
    • Enter the next point A at a relative coordinate of (0,0,320).
  5. Set the Part number and Surface number to 3 and 2, respectively, per Figure 2.
    • Enter the next point B at a relative coordinate of (400,0,0).
  6. Set the Part number and Surface number to 4 and 3, respectively, per Figure 2.
    1. Enter the next point C at a relative coordinate of (0,0,-160).
    2. The final segment of link BCE can be snapped to point E using the mouse. (The cursor will change to a padlock when the mouse is near point E, and clicking will snap the line segment to the point.)
    3. Click the X button to close the Define Geometry dialog box.
  7. Use View Navigate Enclose to show the entire model.
    Tip:
    • If you used the wrong part or surface number when drawing the lines, simply select the line segment (Selection Select Lines), right-click, and choose Edit Attributes. Then set the correct part or surface number.
    • If you enter a wrong dimension, you can undo the lines (on the Quick Access toolbar, click Undo), or select the wrong lines, delete them, and redraw the lines.
  8. Beam elements should be divided into multiple segments to get a stress distribution along the length of the member. Therefore, select all of the lines (Selection Select Lines, then Selection Select All).
    1. Use Draw Modify Divide to split the selected lines into 4 lines each.
    2. Click the OK button to perform the divide operation.
  9. Draw the lines representing the two cylinders using Draw Draw Line.
    1. Change the Part number to 5.
    2. Activate the Single Line option, since the two cylinders are not connected together.
    3. Use the mouse to draw a line from D to C.
    4. Use the mouse to start the next line at point F.
    5. Enter the final point at a relative distance of (600,0,-600).
    6. This completes the geometry, so click the X button to close the Define Geometry dialog box.
  10. Use View Navigate Enclose, if necessary, to show the entire model.

Specify the FEA Parameters:

  1. Enter a name for each of the parts. Select the Part <Unnamed> entry in the browser (tree view), right-click, and choose Rename. Type in the name of each link, per Table 1.
  2. Set the element type for each part per Table 1. Select the Element Type <Unknown> entry in the browser, right-click, and choose the appropriate element type. The element type for parts 1 through 4 can be set simultaneously by holding the Ctrl key and selecting each entry, then right-clicking on one of the selected entries.
  3. Set the element definition for each part per Table 1. Select the Element Definition entry in the browser, right-click, and choose Edit Element Definition. For the beam elements, highlight a cell in the Sectional Properties spreadsheet; this will enable the Cross-Section Libraries button. When the button is clicked, the Cross-Section Libraries dialog will appear. Set the Section database pull-down on the left to User Defined, and set the pull-down on the right to Rectangular. Based on the dimensions in Table 1 and the orientation of the beam elements, enter the height dimension for h and the width dimension for b. After entering the dimensions, click OK and OK to set the element definition. For the bucket (part 2), make the part rigid but with a small cross-sectional area. (The reason for a small area will be apparent when we do the design optimization.) In the analysis, the bucket will experience compression due to the 12000 N load and bending due to the moment. It will not experience any torsion. So, simply enter an area of 10 and moments of inertia I2 and I3 values of 100E6 in the Sectional Properties spreadsheet. Set the section modulus S2 and S3 to 0 so that no bending stresses will be calculated for the bucket. (The area is estimated from the load, length, and modulus of the bucket part to keep a reasonable deflection, and the inertia is a few times larger than the other cross sections.)
  4. Set the material properties for each part per Table 1. Select the Material <Unnamed> entry in the browser, right-click, and choose Edit Material. Parts 1 through 4 can be set simultaneously by holding the Ctrl key and selecting each entry, then right-click one of the entries. (Since the cross-sectional area of the bucket is small compared to the other sections, the bucket will be virtually weightless even though the material includes a mass density. If you truly wanted to make the weight 0, you can edit the Material for part 2 and set the mass density to 0.)
  5. Apply the constraints. Select points G and H (Selection Select Vertices, click point G, hold the Shift key, and click point H).
    1. Right-click in the display area and choose Add Nodal General Constraint.
    2. Click the No Translation button.
    3. Click the OK button.
  6. The above constraints do not create a statically stable 3D model; the model is free to rotate about the Z axis. Add one more constraint by selecting the vertex where the load will be applied, right-click, and choose Add Nodal General Constraint.
    • Activate the Ty check box and click OK. This condition will prevent the out-of-plane rotation.
  7. With the same vertex selected, right-click and choose Add Nodal Force.
    1. Enter a Magnitude of -12000.
    2. Set the direction to Z.
    3. Click the OK button to apply the force.
  8. With the same vertex selected, right-click and choose Add Nodal Moment.
    1. Enter -12000*300 in the Magnitude field. The product will be calculated and
    2. Set the direction to Y. Note that the equation in the Magnitude field is evaluated as soon as you click the direction radio button and 3600000 appears in this field.
    3. Click the OK button to apply the moment.
  9. Apply gravity to the model. Click Setup Loads Gravity.
    1. Click the Set for standard gravity button and confirm that the direction is correct (Z multiplier is set to -1).
    2. Go to the Gravity/Acceleration tab and set the Accel/Gravity multiplier to 1.
  10. By default, the connection between beam elements is a rigid connection. Therefore, end releases need to be added at points A, B, D, and E to simulate the pinned condition. Based on the chosen orientation, we know that the moment about axis 3 (R3) needs to be released, but we do not know a priori whether to release the moment on the I end or the J end of the element. This can be determined by displaying the beam orientation marks (View Visibility Object Visibility Element Axis 1 and Element Axis 2). This shows us that R3 needs to be released on the I end of the beams at points D and B, and the J end of the beams at points A and E. (See Figure 3.) So,
    1. Select the beam element in the bucket at point D and in link BCE at point B (Selection Select Lines, click the first line, hold the Shift key and click the second line),
    2. Right-click and choose Add Beam End Releases.
    3. Activate the R3 check box for the I Node.
    4. Click the OK button to apply the end release.
    5. Repeat the process for the elements at points A and E and activate the R3 check box for the J Node.
    6. When done properly, the end release symbol will appear on the proper end of each beam element as shown in Figure 3.

      Figure 3: View Beam Orientations to set the Proper End Releases

       The orientation marks at the center of each element show which direction is axis 1, and this determines which node is I and which is J. (Axis 1 runs from node I to node J.)
  11. The beam orientation marks can be turned off (View Visibility Object Visibility Element Axis 1 and Element Axis 2).

Part Two – Run the Analysis:

In most cases, the input for design optimization is based on the results of the model as setup. In any event, it is always a good idea to perform the analysis first to make sure it runs properly before doing a more lengthy design optimization analysis. (This model is small enough that a Check Model is not really needed. Check the setup after performing the analysis.)

  1. Analysis Analysis Run Simulation.
  2. When it completes and the model appears in the Results environment, the first thing to check is the beam orientation to confirm that the weak axis of bending (axis 2 in our setup) and strong axis are oriented as intended. Select parts 1, 3, and 4 in the browser, right-click, and choose 3D Visualization. The rectangular cross section is shown for each element. Use View Navigate Orbit to rotate the model around. If the beams are oriented properly, the height of the rectangle should be in the XZ plane and the narrower width in the Y direction. Turn the visualization off when finished.
  3. Use Results Contours Displacement Displacement Magnitude to check the deflections. The maximum deflection of 3 mm is acceptable.
  4. Use Results Contours Stress Beam and Truss Worst to check the stress. The artificial stress in the bucket produces a large compressive stress which is not of interest. Right-click part 2 in the browser and clear Visibility. The stress in the remaining parts range from 51 N/mm 2 (51 MPa) compression to 60 N/mm 2 (60 MPa) tension. With the limit of 100 MPa, there is considerable room to reduce the cross section.
  5. One last thing to confirm: that the end releases are on the proper end of the elements to create the pinned connection at points A, B, D, and E. View the moment about the strong axis of bending (Results Contours Other Results Element Forces Moment About Axis 3). To check the moment on the fly, use Results Inquire Probes Probe. As you hold the mouse near a node, the results at that node are displayed. Confirm that the results near points A, B, D, and E are 0. (At node E, there should be a 0 moment in link BCE, and a large bending moment in link DEFG.)

Part Three – Specify the Design Variables:

The parameters that are adjusted by the design optimization are known as the design variables, and these are specified in the Element Definition dialog. Let's optimize the height of the cross sections for the links indicated in the problem statement.

  1. Return to the FEA Editor (Tools Environments FEA Editor).
  2. Right-click the Element Definition heading of part 1 in the browser and choose Edit Element Definition. Highlight a cell in the Sectional Properties spreadsheet, then click the Cross-Section Libraries button. Right-click in the h field and choose Set as Design Variable. The DV symbol appears to show that the height of the rectangular cross section will be a design variable. Click OK and OK to save the input.
  3. Repeat the above step for parts 3 and 4.

Part Four – Specify the Design Optimization Parameters:

  1. Use Analysis Analysis Optimization to access the Design Optimization dialog. The Design Variable tab shows the design variables set in the model. The Current Value column is filled in with the dimension entered in the Element Definition. For each of the design variables, enter a Lower Limit of 50; we do not want the cross section to be unreasonably small. Since the stress is below the allowable, there is no need to allow a larger upper limit, so enter an Upper Limit equal to the current value.
  2. The Performance tab is where the objective of the analysis and any constraints are entered. In this example, we want to minimize the weight (which generally minimizes the cost) and limit the stress. Click the Add Row button. By default, this adds an objective of minimizing the volume, which is equivalent to minimizing the weight. (Objectives are the quantities being minimized or maximized.)
    1. Since only one objective row can be specified, the Part to be minimized must be 0 or All.
      Tip: Design optimization adjusts the dimensions that are specified as design variables. Decreasing any part dimension reduces the volume of that part, which in turn reduces the total volume for the assembly.
    2. To get the input for the Current Value, switch back to the Autodesk Simulation Mechanical or Multiphysics interface and run the Analysis Analysis Weight and Center-of-Gravity command. This gives the total volume as 1.256e7.
      • Enter 1.256e7 in the Current Value column of the Design Optimization dialog box.
      Tip: Previously, the specified cross-sectional area of the bucket (part 2) was kept small, even though we wanted the bucket to be rigid. The purpose of using a small area is so that its volume would be small compared to the volume being optimized. If the volume of the bucket were to be very large relative to the other parts of the model, then a small change in volume of parts 1, 3, or 4 might represent an insignificant change in the total volume. We want to avoid this potentially misleading result.
    3. The Limit Value of the volume can be estimated from the Weight and Center of Gravity dialog by scaling the volume of each part by the ratio of the design variable's Lower Limit/Current Value and summing the result. This yields a total minimum volume of the model of 4.7616E6. Enter this value into the Limit Value column.
  3. Click Close on the Weight and Center of Gravity dialog.
  4. Back in the Design Optimization dialog, Performance tab, click the Add Row button three times to add three new constraints. We want to set an upper limit on the stress for each part that is being optimized.
    1. Change the Objective/Constraint to Max Stress.
    2. Set the Type to Upper Limit.
    3. Change the Part numbers to 1, 3, and 4 respectively.
    4. The Current Value for each part can be found from the Results environment in the Autodesk Simulation Mechanical or Multiphysics interface. (Click the Results tab.) Display each of the three parts individually to determine the current stress for each (select the parts to hide in the browser, right-click, and clear Visibility).
      • With the Results Contours Stress Beam and Truss Worst result shown for one part, the current value of stress for the current cross-section can be obtained. (Be sure to use the largest absolute value. The maximum stress can be negative or positive.)
      • Clear the Visibility of this part and select Visibility of the next part. Repeat until you have determined the stress for parts 1, 3, and 4.
    5. The Limit Value of 100 (MPa) maximum stress is per the problem statement, and the Current Value of the part stresses were determined in the preceding step. Enter these values, each in the appropriate row and column. The correct values are listed in Table 3 below.
      Objective/ Constraint Loadcase Type Part Current Value Limit Value
      Volume All Minimize All 1.256E7 4.7616E6
      Max Stress All Upper Limit 1 59.86 100
      Max Stress All Upper Limit 3 5.42 100
      Max Stress All Upper Limit 4 27.38 100
      Table 3: Performance Tab Input
  5. The default input on the Parameters tab will be used, and the results checked to confirm that they were acceptable.

Part Four – Run the Design Optimization Analysis:

  1. Go to Analysis Perform Design Optimization.
    • When the dialog appears, click the Analyze button to begin the optimization process.
  2. When the analysis completes, check the last lines of information in the display window (the optimization log file). The critical items are described in Table 4.
    01 ** The change of constraints are less than limit!
    02 ** Finished Optimization after 5 Iterations
    03 ** Optimum Design Variables **
    04 ===========================================================================
    05 Design Current Value Lower Limit Start Value Upper Limit
    06 ---------------------------------------------------------------------------
    07 1 1.8922e+002 5.0000e+001 2.5000e+002 2.5000e+002
    08 2 5.0000e+001 5.0000e+001 1.5000e+002 1.5000e+002
    09 3 7.6111e+001 5.0000e+001 1.5000e+002 1.5000e+002
    10 ===========================================================================
    11 ** Optimum Objective/Constraint Value **
    12 ===========================================================================
    13 Obj/Const Current Value Limit Value Change/Violation
    14 ---------------------------------------------------------------------------
    15 Objective 9.38600e+006 1.25600e+007 25.270701 % Decreased
    16 1st Const 9.96355e+001 1.00000e+002 -0.364452 % Not Violated
    17 2nd Const 1.62569e+001 1.00000e+002 -83.743147 % Not Violated
    18 3rd Const 9.99214e+001 1.00000e+002 -0.078580 % Not Violated
    19 ===========================================================================
    20 ** Finished Optimization !! **
    Table 4: The Last Iteration of Design Optimization Analysis
    Line numbers were added to help with the following description.
    The analysis has converged (line 01) because the change in the user-entered constraints from one iteration to the next was less than the user-entered tolerance, and this was reached in fewer iterations than the entered maximum number of iterations (line 02). Line 15 shows the objective (volume, being minimized). The Change/Violation column indicates the optimum design is 25.3% less than the user-entered current value. The first constraint (line 16, stress) had an upper limit of 100 MPa specified. The final optimum value was 99.6 MPa. The second constraint (line 17, stress) was greatly under upper limit of 100 MPa. The stress is not higher because the design variable (the height of the cross section, line 08) is at the lower limit of 50 mm; thus, the cross section cannot be any smaller, so the stress cannot be any higher than 16.3 MPa. The third constraint (line 18, stress) is also at the limit of 100 MPa. (Text shown in bold is highlighted to correspond with this description.)
  3. Click the Done button to close the analysis window. Since the results are acceptable, no changes are required. If the outcome were not acceptable, such as one of the constraints had been violated, then the design optimization could be re-run with different upper and lower limits on the constraints (Design Variable tab) or with different optimization parameters (Parameters tab).

Part Six – Review the Optimum Results

  1. From the Design Optimization Results menu, choose Design Optimization History. Then the Results frame contains each objective, constraint, and design variable in the drop-down menu. Select each item to see a graph of how the value changed with each iteration.
  2. Use the Results menu and choose Optimum Result. This will start a new instance of Autodesk Simulation and load a copy of the original model but with the design variables set to the optimum value calculated by the design optimization.
  3. In the Autodesk Simulation interface, click the Results tab to view the optimum model results. Following the sequence of steps performed on the original model, we find that the displacement magnitude is slightly higher (4.5 mm) than before, but still acceptable. Naturally, the stress results are acceptable. Thus, the answers to item 1 of the problem description are summarized in the following table.
    Link Part Number Optimum Cross Section
       

    width out-of-plane ("b")

    x height in-plane ("h")

    DEFG 1 24 x 189 mm
    AB 3 14 x 50 mm
    BCE 4 20 x 76 mm
    Table 5: Optimum Cross Section
  4. To get the cylinder sizes (item 2 of the problem description), hide all parts except for part 5. (Select parts 1 through 4 in the browser, right-click one of them, and clear Visibility.)
  5. Set the results type to Results Contours Other Results Element Forces Axial Force.
  6. The forces in the cylinders are as follows, determined either from the legend or by inquiring on the elements (Results Inquire Probes Probe).
    Cylinder Axial Force
    CD 21400 N
    FH 72450 N
    Table 6: Cylinder Forces
    (to 4 Significant Figures)

    The size of the cylinders can be calculated based on an assumed fluid pressure.

Two archives of the model and results are located in the Models subfolder of the Autodesk Simulation Mechanical or Multiphysics installation directory:

  • front-end loader input results.ach is the original model with the design variables defined
  • front-end loader_OPT.ach is the model with geometry based on the optimum results.

This example is based on example 6.149 from Beer, Ferdinand P. and Johnston, E. Russell, Vector Mechanics for Engineers, Statics and Dynamics, McGraw Hill, 7th Edition, page 341.