Artículos

URI permanente para esta colección

https://hdl.handle.net/10784/1349

Examinar

Mostrando 1 - 2 de 2

MultiZ: A Library for Computation of High-order Derivatives Using Multicomplex or Multidual Numbers
(Association for Computing Machinery (ACM), 2020-01-01) Aguirre-Mesa A.M.; Garcia M.J.; Millwater H.; Mecánica Aplicada
Multicomplex and multidual numbers are two generalizations of complex numbers with multiple imaginary axes, useful for numerical computation of derivatives with machine precision. The similarities between multicomplex and multidual algebras allowed us to create a unified library to use either one for sensitivity analysis. This library can be used to compute arbitrary order derivates of functions of a single variable or multiple variables. The storage of matrix representations of multicomplex and multidual numbers is avoided using a combination of one-dimensional resizable arrays and an indexation method based on binary bitwise operations. To provide high computational efficiency and low memory usage, the multiplication of hypercomplex numbers up to sixth order is carried out using a hard-coded algorithm. For higher hypercomplex orders, the library uses by default a multiplication method based on binary bitwise operations. The computation of algebraic and transcendental functions is achieved using a Taylor series approximation. Fortran and Python versions were developed, and extensions to other languages are self-evident. © 2020 ACM.
A stiffness derivative local hypercomplex-variable finite element method for computing the energy release rate
(Elsevier BV, 2019-01-01) Aguirre-Mesa A.M.; Ramirez-Tamayo D.; Garcia M.J.; Montoya A.; Millwater H.; Mecánica Aplicada
A “local” hypercomplex-variable finite element method, L-ZFEM, is proposed for the computation of the energy release rate (ERR) using the stiffness derivative equation. This approach is analogous to the stiffness derivative method proposed by Parks and Hellen but has superior numerical accuracy. In addition, this method is significantly more efficient than the previously published “global” hypercomplex-variable finite element method, ZFEM, in that the global hypercomplex system of FE equations is not assembled nor solved. Instead, the displacement field is computed using a traditional, real-valued finite element method, and the numerical derivative of the stiffness matrix at the element level is only computed for a group of local, surrounding elements to the crack tip by using a Taylor series expansion based on complex numbers or dual numbers. The ERR is then determined as a sum of the element contributions. Derivatives of the ERR with respect to an arbitrary model parameter such as a crack extension, material property, or geometric feature are also available using a combination of the global and local methods, GL-ZFEM. Both L-ZFEM and GL-ZFEM were implemented into the commercial finite element software Abaqus through user defined element subroutines. Numerical results show that the ERR obtained by L-ZFEM has the same accuracy as that estimated through the global ZFEM or the J-integral methods but exhibits superior computational efficiency. © 2019 Elsevier Ltd

Examinar

Examinando Artículos por Autor "Aguirre-Mesa A.M."

Resultados por página

Opciones de ordenación