Pushover Analysis with 1D Beam Element

Modified on Thu, 11 Jul at 2:37 PM

What is in this Article?

1. Introduction

2. Problem Description

3. Before We Start

4. Implementation

5. Results


     Terminologies

Collection

The folder containing data of an asset/model on Akselos Cloud or users’ computer

Components

Components created by the componentization process to use with Akselos Integra 

Node port

A mesh node or beam end used to connect 1D beam elements to other component types 

Nodeset

A mesh node assigned with a certain ID number

Model ribbon

An Akselos Modeler ribbon containing tools for model assembling and management

Ribbons

The top-sided toolbars of Akselos Modeler

Property Tree

A panel at the left bottom of Akselos Modeler, where shows properties of user selection


1. Introduction

Pushover analysis for beam models is now supported by the Akselos solver. Pushover capabilities within Akselos solver can provide insight into the structural aspects by enabling the users to conduct a full elastoplastic pushover analysis of a structure.

In this tutorial, we will do the pushover analysis for the frame model with the 1D beam element. It will guide you on how to create the beam frame model via the Akselos Modeler tool. You will learn how to create, apply loads, set up time steps, boundary conditions on the model, and solve it.

 Figure 1.1. Pushover analysis with 1D beam element workflow


2. Problem Description

In this tutorial, we use the 1D beam element to build the beam frame model which is constrained at 4 legs, under self-weight, and forces applied on 4 points on top. FV is the constant load. FH is the pushover load which is increased by step until any degree-of-freedom displacement first equals or exceeds 0.25m then the solver will terminate. This time also means the structure is considered collapsed. The model schematic and beam cross-section are shown in the figure below.

 Figure 2.1. The beam frame problem

The beam geometry, beam properties, and load cases used in this example model are as follows:

  • Beam cross-section (m): Tubular with below dimensions:
    • Cross-section 1          
      • Diameter: 1.016 m
      • Thickness: 0.035 m
    • Cross-section 2          :
      • Diameter: 1.0 m
      • Thickness: 0.04 m
    • Cross-section 3          :
      • Diameter: 2.0 m
      • Thickness: 0.06 m
  • Material: Elasticity with the below properties:
    • Mass density: 8000 kg/m3
    • Poisson ratio: 0.3
    • Young Modulus: 210 GPa
    • Tangent Modulus: 0.42 GPa
    • Yield Stress: 330 MPa
  • Applied loads:
    • Self-weight
    • Constant Load: FV = 30 MN (Point Load)
    • Pushover Load: FH = 30 MN (Point Load)

Figure 2.2. The beam Cross-section


3. Before We Start

To follow the instructions below, please notice the point below:

  

Note:
▶ If you have not installed Akselos Modeler yet, please see Akselos Modeler - Installing, Updating and Managing the Akselos Simulation Software to download and install the software.
▶ If you have trouble accessing to any of those prerequisites. Please contact us at: support.akselos.com

  


4. Implementation

In the Implementation section, we will show you a step-by-step in detail guide on how to analyze the pushover problem in Akselos Modeler.

Figure 4.1. Pushover analysis with 1D beam element workflow

Step 1: Create Nodes

Figure 4.2. Create Model – Create Nodes

Akselos Modeler provides an Edit Beams tool that helps users create nodes and beams. The tool also can modify beams such as split/merge or change the beam type.

  • On the left panel, select the Edit Beam tool.

Figure 4.3. Opening the Edit Beams tool

Figure 4.4. The Edit Beams tool interface

To build this beam model, we need to create 13 nodes at the needed positions. There are two options to perform this action:

  • Option 1: Create nodes from table, click on the icon to insert the row  Insert until there are 13 rows on the table  Fill X, Y, Z values into them as in the figure below  Click on the Create Nodes button below the table.

Figure 4.5. Creating Nodes

  • Option 2: Create nodes from csv file.
    1. Firstly, you need to prepare a csv file with the correct format as the table in Option 1 or download here - Model Nodes.

Figure 4.6. Data format in csv file

  1. In the Edit Beams interface, click on  icon  find and open csv file.

Figure 4.7. Data format in csv file

After finishing one of these options, 13 nodes will appear on the Graphic Window.

Figure 4.8. Thirteen nodes on the Graphic Window after creating

Step 2: Create Beams Model From Nodes

Figure 4.9. Create Model – Create Beams Model From Nodes

From those nodes, you can create the beam frame.

  • Click on the Create Beam From Nodes button  Select nodes on the Graphic Window (every two nodes selected will create one beam). Do the same until finish the model with 28 beams.

 

Figure 4.10. Creating beams from nodes

Step 3: Split Beams

Figure 4.11. Create Model – Split Beams

You can skip this step but to increase the accuracy of the result, splitting the beam into many elements is also important.

  • In the Edit Beams interface, change to the Split/Merge tab  select all beams on the Graphic Window  Input 8 into the Number of Segment box → Leave the default for other options and click on the Split button.

Figure 4.12. Splitting beams

Figure 4.13. The model after splitting

Step 4: Add Cross-sections

Figure 4.14. Assign Properties – Add Cross-sections

In this step, we will guide you to create and define three cross-sections and then apply them to the model.

  • In the Model ribbon on the left panel, right-click on Model  Insert Property  Cross-Sections. The group of cross-sections will be created with one default item. You can add more Cross-sections by right-clicking on the group → Insert Property  Crosssectiongroupitem or delete them by clicking on then pressing the Delete button. Each item has an ID number on the right, it starts from 1 with the default item. The ID number is used to assign the cross-section to the model.

 

Figure 4.15. Creating CrossSectionGroup

  • Click on Cross-Section 1 then see the Property Tree  Define the dimensions.

Figure 4.16. Defining Cross-section properties

  • Insert two cross-sections and define them.

Figure 4.17. Defining Cross-sections 2 & 3

To assign the recently created Cross-sections to the model, change to Mesh Elements Mode:

  • Select the top beams as in the figure below  In the Property Tree, expand Fields and input the ID of CrossSection 2 into the cross_section_id box.

Figure 4.18. Assigning CrossSection 2 to the top beams

  • Select the beams of four legs as in the figure below  In the Property Tree, expand Fields and input the ID of CrossSection 3 into the cross_section_id box.

Figure 4.19. Assigning CrossSection 3 to the top beams

Step 5: Add Materials

Figure 4.20. Assign Properties – Add Materials

  • In the Model ribbon on the left panel, right-click on Model  Insert Property  Materials. The group of materials will be created with one default item. You can add more Materials by right-clicking on the group → Insert Property  Materialgroupitem or delete them by clicking on then pressing the Delete button. Each item has an ID number on the right, it starts from 1 with the default item. The ID number is used to assign the material to the model.

 

Figure 4.21. Creating MaterialGroup

  • Click on Material 1 then see the Property Tree  Define the properties of elasticity.

Figure 4.22. Defining elasticity properties of Material 1

  • In the Property Tree of Material 1, right-click on value → Insert New… → plasticity  Expand hardening_model and input values to tangent_modulus and yield_stress.

Figure 4.23. Inserting and defining plasticity properties of Material 1

  • To assign the newly created material to the model, change to Subdomains Selection Mode  Select model on the Graphic Window  See the Property Tree and fill the ID of the created material to the material_id box. The default is 1.

Figure 4.24. Assigning material to model

After defining and assigning material and cross-sections, we can visualize it on the Graphic Window by turning on Show Beam Cross-section in Configure view settings. Note that we need to assign both for the Show Beam Cross-section tool to work properly. You can turn off it to easily work in the next step.

 

Figure 4.25. Showing beam model with cross-section on the Graphic Window

Step 6: Constrain Model

Figure 4.26. Model Set-up – Constrain model

In this tutorial, four legs of the frame are constrained. We will constrain 4 nodes at the bottom of the model. Before doing it, we need to create Node Port at these nodes and then fix them.

  • Firstly, we need to turn on the Unconnected Ports and Boundary Conditions options in Configure view settings. The ports and boundaries of the model are visualized on the Graphic Window.

Figure 4.27. Turning on Unconnected Ports and Boundary Conditions options

  • Select the node at the bottom Right-click on the Graphic Window → Create Node Port. Do the same for the remaining three nodes. Note that we have to do it one by one because the tool does not support when selecting multiple nodes.

Figure 4.28. Creating Node Port

Figure 4.29. Four node ports are created

  • To constrain them, select four node ports created above see the Property Tree and select (x, y, z, theta_x, theta_y, theta_z) in Constraint Type. The software will automatically do the same for other ports. After that, the black points will appear at four ports, it performs that four ports are constrained successfully.

 

Figure 4.30. Constraining node ports

Figure 4.31. The model with four constrained nodes

Step 7: Set Time Array & Time Steps

Figure 4.32. Model Set-up – Set Time Array & Time Steps

The incremental steps in pushover analysis represent the progression of loading from one limit state to another as the applied lateral forces increase. The magnitude of these increments depends on factors like the required accuracy, the complexity of the structure, and the specific objectives of the analysis. Smaller increments generally provide more precise results but lead to longer computation. In this step, we will guide you on how to set up time steps in Akselos Modeler for the load applied in the next step. 

  • Time Array: The model has both constant loads and pushover loads, time_array should be an array with two entries: e.g. "time_array": "step_array([0.0, 1.5])*s". The solver initializes constant loads in the first step (at t = 0), and the pushover loads are added to the model in the second step.
    • To create a time array property, right-click on Model on the left pannel Insert NewGeneric Property Property.

Figure 4.33. Inserting Model Property

    • Rename the recently created property to time_arrayand enter into the box as in the figure below.

Figure 4.34. Time array property

  • Time Steps:
    • There are three child items in this property tree: load_substep_type, number_of_substeps, and time. The load_substep_type can be changed from a drop-down list of options, but the number_of_substeps must be given by the users. It should be noted that the time_array established previously is used to enter time.
    • There are three child items in this property tree: load_substep, number_of_substeps, and time.
    • The default value of max_arclen_substep is 100, which is sufficient in most cases. However, for some test cases, this value may need to be increased for the solver to trace a complete load-displacement curve.
    • The max_displacement_size stopping control method is optional, it uses the user input value as a criterion to stop solving. This input value represents the maximum allowable displacement (of any node in the structure) before the structure is considered collapsed.
    • To change the time steps settings, select default_time_stepping on the left panel.

Figure 4.35. Default_time_stepping

    • In the Property Tree, right-click on default_time_step  Insert new → number_of_substeps (InPropertyArray).

Figure 4.36. Inserting number_of_substeps

    • Change load_substep_type to ArcLengthContinuationProperties.

Figure 4.37. Setting load_substep_type

    • Right-click on load_substep_type  Insert new → path_predictor_type (PathPredictorType).

Figure 4.38. Inserting path_predictor_type

    • Right-click on auto_step_properties  Insert New  max_displacement_size (m).

Figure 4.39. Inserting max_displacement_size

    • Right-click auto_step_properties  Insert New  max_substeps (PropertyArray).

Figure 4.40. Inserting max_substeps

    • Change all values as in the figure below.

Figure 4.41. Setting default_time_stepping properties

Step 8: Apply Loads

Figure 4.42. Model Set-Up – Apply Loads

In this tutorial, there are three load cases applied to the mode (Self-Weight, Point-Load-Constant, Point-Load-Pushover). Self-Weight is applied on a subdomain selection and two Point-Load are applied on four nodes. To create a load applying on the node, we need to create Node Set and store them in Stored Selection.

  • Create Node Set:
    • Change to the Collection ribbon → Show component tree → Right-click on Boundary Sets → Create Boundary Set → Node.

Figure 4.43. Creating Node Set

    • Click on New Node  Select four nodes on top as the figure below.

Figure 4.44. Selecting four nodes for Node Set

    • We need to turn on the PointLoad operator for this Node Set. Click on New Node  see the Property Tree and turn on PointLoad in Operators.

Figure 4.45. Turning on the PointLoad operator

    •  After all, we need to save the change of component. Right-click on the component  Save Changes

Figure 4.46. Saving component

  • Stored Selection for Self-Weight Load Case:
    • Change to the Model ribbon  On the left panel, right-click on Stored Selections  Create Subdomains Selection.

Figure 4.47. Creating Subdomains Selection

    • Click on Subdomain Selection 1  Select the model on the Graphic Window.

 

 

Figure 4.48. Storing Subdomain for Self-Weight load case

  • Stored Selection for Point-Load Load Case:
    • Right-click on Stored Selections → Create Boundary Sets Selection.

Figure 4.49. Creating Boundary Set Selection

    • Click on Subdomain Selection 1  Select node-set (New Node created in the previous stepon the Graphic Window. Because the position of this node set coincides with another one, to select the right node-set, move the mouse to one of its nodes, use the Tab button to change to the right name shown at the bottom of the Graphic Window then click once.

Figure 4.50. Storing node sets for Point-Load load case

  • Create Self-Weight Load Case:
    • Right-click on Load Cases → Create Load Case

Figure 4.51. Creating Load Case

    • Right-click on Load Case 1 → Create Load

Figure 4.52. Creating load in load case

    • Click on Load 1  see the Property Tree and set them as in the figure below.

Figure 4.53. Defining Self-Weight load

Figure 4.54. The load vector of self-weight appeared after applying

  • Create Point-Load-Constant Load Case:
    • Right-click on Load Cases → Create Load Case.

Figure 4.55. Creating Load Case

    • Right-click on Load Case 2 → Create Load

Figure 4.56. Creating load in load case

    • Click on Load 1 of Load Case 2  see the Property Tree and set them as in the figure below.

Figure 4.57. Defining Point-Load

 

Figure 4.58. The load vectors of the point load (constant load) appeared after applying

  • Create Point-Load-Pushover Load Case: We can copy from Load Case 2 and modify some values to use.
    • Right-click on Load Case 2 → Clone.

Figure 4.59. Creating load case by cloning from another one

    • Expand Load Case 3  click on Load 1 → Modify the values in the Property Tree as in the figure below.

Figure 4.60. Creating load in load case

To easily manage the load cases, we should rename them following the load (Right-click on it  Rename this Load Case or click on it  press F2 button)

  • Load Case 1 = Self-Weight
  • Load Case 2 = Point-Load-Constant
  • Load Case 3 =Point-Load-Pushover

Figure 4.61. Renaming three load cases

Step 9: Create Solver Option

Figure 4.62. Solve – Create Solver Option

For the 1D beam element problem, we can only solve it with the FEA solver.

  • Under Solver Options on the left panel, there is RB-FEA Options, it is the default and can be used with some changes. Click on RB-FEA Options  See the Property Tree and change Solver Strategy to FEA. After changing, the name of the default solver will automatically change to FEA Options.

Figure 4.64. Setting FEA Solver

  • In the Solver Option, you can choose the result fields that you want to visualize. The default results of the solution are displacements in x, y, z-directions. In this example, we choose some additional fields which are Axial Force, Shear Forces, Bending Moments, and Plastic Indicator. Turning on three options as in the figure below.

Figure 4.65. Adding some fields on Solver Option

Step 10: Create Scenario 

Figure 4.66. Solve – Create Scenario

  • Under Scenario on the left panel, there is Default Scenario and we can use it. Click on Default Scenario  See the Property Tree and set the coefficients of Point-Load and Self-Weight to 1.0.

Figure 4.67. Setting Scenario

  • Set Solve List: Since there is only one scenario, select the default in Solve List to check that the default scenario is turned on.

Figure 4.68. Setting Solve list

Step 11: Sync Collection and Solve

Figure 4.69. Solve – Sync collection & Solve

Now you can proceed to solve the model. First, you have to save aks file to save all setting up the model  sync the collection to Akselos Dashboard  Run solving.

  • Save aks File:
    • Click on the File tab Save.

Figure 4.70. Saving aks file

    • The Select destination to save the file window appears, enter the name into the File name box then click the Save button. You will receive a success message when the file is successfully saved.

Figure 4.71. Naming and saving aks file

  • Sync Collection:
    • Change to the Collection ribbon Click on Repository...  Sync.

Figure 4.72. Syncing collection

    • Click the Commit button when this window appears.

Figure 4.73. Committing for synchronization

  • Run Solving: Change to the Solutions ribbon  Click on the Solve button. The solving request will be submitted to Akselos Cloud. You will see the status of this request on the Solution tree. You can see the results once the solving is done, which is usually within seconds with Akselos solver. 

Figure 4.74. Solving model

Figure 4.75. Solution after solving.


5. Results

Figure 5.1. Validate – Compare with theoretical results.

After receiving the solution, you can study its solution fields. There are some solution fields (visualization & results). You can review it to compare it with yours.

Figure 5.2. Axial Force result

Figure 5.3. Displacement x-direction result

Figure 5.4. Displacement y-direction result

Figure 5.5. Displacement z-direction result

Figure 5.6. Plastic Indicator result

Figure 5.7. Shear Force results

Figure 5.8. Bending Moment results

  • Validation:

The horizontal response calculated by the Akselos solver for the current jacket structure is compared to reference results [1] in the figure below. In terms of ultimate load capacity, the results appear to be in good agreement.

 

Figure 5.9. Load-displacement curves for the Four-legged jacket structure

[1] Engseth, A. Finite Element Collapse Analysis of Tubular Steel Offshore Structures, Div. of Marine Structures, The University of Trondheim, The Norwegian Institute of Technology, Trondheim, 1984.


Related Articles & Similar Case Studies:

Was this article helpful?

That’s Great!

Thank you for your feedback

Sorry! We couldn't be helpful

Thank you for your feedback

Let us know how can we improve this article!

Select at least one of the reasons
CAPTCHA verification is required.

Feedback sent

We appreciate your effort and will try to fix the article