A block forward substitution method for solving the hypercomplex finite element system of equations

dc.contributor.authorAguirre-Mesa, Andres M.
dc.contributor.authorGarcia, Manuel J.
dc.contributor.authorAristizabal, Mauricio
dc.contributor.authorWagner, David
dc.contributor.authorRamirez-Tamayo, Daniel
dc.contributor.authorMontoya, Arturo
dc.contributor.authorMillwater, Harry
dc.date.accessioned2023-05-23T17:13:55Z
dc.date.available2023-05-23T17:13:55Z
dc.date.issued2021-12-15
dc.description.abstractThe hypercomplex finite element method, ZFEM, allows the analyst to compute highly-accurate arbitrary-order shape, material property, and loading derivatives by augmenting the traditional finite element method with multiple imaginary degrees of freedom. In ZFEM, the real variables are converted to hypercomplex variables such as multicomplex, multidual, or quaternions. By uplifting the real variables to hypercomplex, derivatives are computed in an automated fashion using a standard finite element formulation. The use of multicomplex or multidual numbers provides higher-order derivatives. The drawback of ZFEM is that it increases the number of degrees of freedom of the real variable system by a factor 2^n, where n is the order of the required derivative. In consequence, ZFEM increases the memory consumption and the solution time of the system of equations compared to the real variable system. The block forward substitution method (BFS), proposed in this work, addresses the memory and runtime issues. This new method solves the original real-valued FEM system once. Then, the derivatives are computed using pseudo-loads with the original system of equations. In contrast with the conventional solution method of ZFEM, BFS computes the hypercomplex contributions to the stiffness matrix element-wise, and it never assembles nor solves the full hypercomplex system of equations. In effect, the BFS method generalizes the first-order semi-analytical complex variable method to any order derivative. The BFS method provides the capability to allow a combination of real-variable and hypercomplex-variable elements within the same model. The numerical results indicate that a first-order derivative can be obtained with 1% to 8% additional computational time of the real-variable analysis. This allows the computation of multiple first-order derivatives by post-processing of a single FEM analysis. Additionally, it was shown that fourth-order shape sensitivities can be computed in less than 5% additional runtime of the real-variable FEM analysis.en_US
dc.description.departmentMechanical Engineeringen_US
dc.description.departmentCivil and Environmental Engineering, and Construction Management
dc.description.sponsorshipUnited States Department of Defense; Universidad EAFIT, Colombia and COLCIENCIAS’ Scholarship Programsen_US
dc.identifier.citationAguirre-Mesa, A. M., Garcia, M. J., Aristizabal, M., Wagner, D., Ramirez-Tamayo, D., Montoya, A., & Millwater, H. (2021). A block forward substitution method for solving the hypercomplex finite element system of equations. Computer Methods in Applied Mechanics and Engineering, 387, 114195. doi:https://doi.org/10.1016/j.cma.2021.114195en_US
dc.identifier.issn1879-2138
dc.identifier.otherhttps://doi.org/10.1016/j.cma.2021.114195
dc.identifier.urihttps://hdl.handle.net/20.500.12588/1851
dc.language.isoen_USen_US
dc.publisherElsevieren_US
dc.relation.ispartofseriesComputer Methods in Applied Mechanics and Engineering;387
dc.rightsAttribution 3.0 United States*
dc.rights.urihttp://creativecommons.org/licenses/by/3.0/us/*
dc.subjectsemi-analytical complex-variable methoden_US
dc.subjecthypercomplex numbersen_US
dc.subjectdual numbersen_US
dc.subjecthyperdual numbersen_US
dc.subjectcomplex-variable finite element methoden_US
dc.subjectcomplex Taylor series expansionen_US
dc.titleA block forward substitution method for solving the hypercomplex finite element system of equationsen_US
dc.typeArticleen_US

Files

Original bundle

Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
Aguirre-Mesa 2021 - A block forward substitution method for solving the hypercomplex finite element system of equations.pdf
Size:
1.18 MB
Format:
Adobe Portable Document Format
Description:

License bundle

Now showing 1 - 1 of 1
No Thumbnail Available
Name:
license.txt
Size:
1.86 KB
Format:
Item-specific license agreed upon to submission
Description: