Pure Julia Implementation of Geometrically Exact Beam Theory
Author byuflowlab
15 Stars
Updated Last
2 Years Ago
Started In
June 2020


Pure Julia Implementation of Geometrically Exact Beam Theory

Author: Taylor McDonnell

GXBeam is a pure Julia implementation of Geometrically Exact Beam Theory, based on the open source code GEBT and its associated papers[1][2].

As a sample of one of the many things this package can do, here's a time domain simulation of the dynamic response of a joined wing subjected to a simulated gust, scaled up in order to visualize the deflections:

Package Features

  • Performs multiple types of analyses including:
    • Linear/Nonlinear static analyses
    • Linear/Nonlinear steady-state dynamic analyses
    • Linear/Nonlinear eigenvalue analyses (by linearizing about a steady state condition)
    • Linear/Nonlinear time-marching dynamic analyses
  • Accurately models arbitrary systems of interconnected highly flexible composite beams.
    • Captures all geometric nonlinearities due to large deflections and rotations
    • Capable of using the full 6x6 Timoshenko beam stiffness matrix
    • Singularity-free rotational deflections of any magnitude using only 3 rotational parameters
  • Models arbitrary time-varying distributed forces/moments on beam elements using:
    • Dead forces/moments (which do not rotate as the beam element rotates)
    • Follower forces/moments (which rotate as the beam element rotates)
  • Models arbitrary time-varying prescribed forces/moments and/or displacements/rotations at the connection points between beam elements using:
    • Dead forces/moments (which do not rotate as the point rotates)
    • Follower forces/moments (which rotate as the point rotates)
  • Capable of using arbitrary units (as long as they are compatible)
  • Simple result visualization using WriteVTK
  • Built-in DifferentialEquations interface for time domain simulations.
  • Extensively validated against published analytical and computational results. See the examples in the documentation.


Enter the package manager by typing ] and then run the following:

pkg> add GXBeam


This code has been optimized to be highly performant, primarily by maintaining type stability and minimizing allocations. As a result the performance of this package rivals (or even outperforms) that of the Fortran implementation in GEBT. At this point, differences in performance between the two codes can be primarily attributed to the performance of the sparse linear system solver in each.


See the documentation

Note that while the theoretical basis for this code is identical to Wenbin Yu's code, some of the implementation details vary.


[1] Yu, W., & Blair, M. (2012). GEBT: A general-purpose nonlinear analysis tool for composite beams. Composite Structures, 94(9), 2677-2689.

[2] Wang, Q., & Yu, W. (2017). Geometrically nonlinear analysis of composite beams using Wiener-Milenković parameters. Journal of Renewable and Sustainable Energy, 9(3), 033306.