Closed form integration of element stiffness matrices using a computer algebra system

1995 ◽  
Vol 56 (4) ◽  
pp. 529-539 ◽  
Author(s):  
C.K. Yew ◽  
J.T. Boyle ◽  
D. MacKenzie
Author(s):  
Sara E. McCaslin ◽  
Panos S. Shiakolas ◽  
Brian H. Dennis ◽  
Kent L. Lawrence

Closed-form solutions for straight-sided tetrahedral element stiffness matrices used in finite element analysis have been proven more efficient than numerically integrated solutions. These closed-form solutions are symbolically integrated using computer algebra systems such as Mathematica or Maple. However, even with memory and processing speed available on desktop computers today, major hindrances exist when attempting to symbolically evaluate the stiffness matrices for high order elements. This research proposes a new approach to obtaining closed-form solutions. Results are presented that demonstrate the feasibility of obtaining the stiffness matrices for high order tetrahedral elements through p-level 9 by use of parallel processing tools in Mathematica 7. Comparisons are made between serial and parallel approaches based on memory required to generate a solution. The serial approach requires more memory and can only generate closed-form solutions up to 7th order. The parallel processing approach presented requires less memory and can generate solutions up to 9th order.


Author(s):  
Panos S. Shiakolas ◽  
David C. Wilhite ◽  
Sara E. McCaslin

Computer algebra systems (CAS) have been advantageously employed to generate closed form expressions for finite elements. The advantages relate to the time improvements or savings realized by employing closed form generated expressions as compared to numerical integration. However as the element order increases, the size of the closed form generated expressions become unmanageable causing the source code files to possibly become unusable due to their size. One approach to reducing the size of the source files is to take advantage of the utilities found in CAS to identify common expressions or sub-expressions. In this manuscript we present on-going research by comparing two widely used CAS, Mathematica and Maple, as they relate to identifying common expressions in low order tetrahedral finite element stiffness matrices generated in symbolic form, associated time savings and possible issues. The results indicate that the use of CAS could be advantageously employed to identify common subexpressions through pattern matching to further reduce the size of the generated source files and realize time improvements during execution of the source codes. In addition, the developed procedures could be easily applied to higher order elements with much larger number of entries of closed form expressions where even more savings could be realized.


1998 ◽  
Vol 37 (03) ◽  
pp. 235-238 ◽  
Author(s):  
M. El-Taha ◽  
D. E. Clark

AbstractA Logistic-Normal random variable (Y) is obtained from a Normal random variable (X) by the relation Y = (ex)/(1 + ex). In Monte-Carlo analysis of decision trees, Logistic-Normal random variates may be used to model the branching probabilities. In some cases, the probabilities to be modeled may not be independent, and a method for generating correlated Logistic-Normal random variates would be useful. A technique for generating correlated Normal random variates has been previously described. Using Taylor Series approximations and the algebraic definitions of variance and covariance, we describe methods for estimating the means, variances, and covariances of Normal random variates which, after translation using the above formula, will result in Logistic-Normal random variates having approximately the desired means, variances, and covariances. Multiple simulations of the method using the Mathematica computer algebra system show satisfactory agreement with the theoretical results.


1995 ◽  
Vol 10 (3) ◽  
pp. 329-337 ◽  
Author(s):  
John Hutton ◽  
James Hutton

2020 ◽  
Vol 9 (1) ◽  
pp. 53
Author(s):  
Kamhar Ngado ◽  
Rosnawati Rosnawati ◽  
Heri Retnawati ◽  
Sri Andayani

2021 ◽  
Vol 254 ◽  
pp. 02006
Author(s):  
Liubov Feshchenko ◽  
Gleb Vodinchar

The paper describes a technology for the automated compilation of equations for shell models of turbulence in the computer algebra system Maple. A general form of equations for the coefficients of nonlinear interactions is given, which will ensure that the required combination of quadratic invariants and power-law solutions is fulfilled in the model. Described the codes for the Maple system allowing to generate and solve systems of equations for the coefficients. The proposed technology allows you to quickly and accurately generate classes of shell models with the desired properties.


Sign in / Sign up

Export Citation Format

Share Document