CC BY
34
Наука та прогрес транспорту. Вкник Дншропетровськош нацюнального ушверситету залiзничного транспорту, 2016, № 2 (62)
UDC 624.04
M. SYBIS1*, A. SMOCZKIEWICZ-WOJCIECHOWSKA2*, A. SZYMCZAK-GRACZYK3*
1 Institute of Civil and Geoingeenering, Poznan University of Life Sciences, ul. Pi^tkowska 94, 60-649 Poznan, Poland, email msybis@up.poznan.pl, ORCID 0000-0002-0032-6313
2*Institute of Civil and Geoingeenering, Poznan University of Life Sciences, ul. Pi^tkowska 94, 60-649 Poznan, Poland, email asw72@up.poznan.pl, ORCID 0000-0002-7202-6678
3*Institute of Civil and Geoingeenering, Poznan University of Life Sciences, ul. Pi^tkowska 94, 60-649 Poznan, Poland, email agraczyk@up.poznan.pl, ORCID 0000-0002-1187-9087
IMPACT OF MATRIX INVERSION ON THE COMPLEXITY OF THE FINITE ELEMENT METHOD
Purpose. The development of a wide construction market and a desire to design innovative architectural building constructions has resulted in the need to create complex numerical models of objects having increasingly higher computational complexity. The purpose of this work is to show that choosing a proper method for solving the set of equations can improve the calculation time (reduce the complexity) by a few levels of magnitude. Methodology. The article presents an analysis of the impact of matrix inversion algorithm on the deflection calculation in the beam, using the finite element method (FEM). Based on the literature analysis, common methods of calculating set of equations were determined. From the found solutions the Gaussian elimination, LU and Cholesky decomposition methods have been implemented to determine the effect of the matrix inversion algorithm used for solving the equations set on the number of computational operations performed. In addition, each of the implemented method has been further optimized thereby reducing the number of necessary arithmetic operations. Findings. These optimizations have been performed on the use of certain properties of the matrix, such as symmetry or significant number of zero elements in the matrix. The results of the analysis are presented for the division of the beam to 5, 50, 100 and 200 nodes, for which the deflection has been calculated. Originality. The main achievement of this work is that it shows the impact of the used methodology on the complexity of solving the problem (or equivalently, time needed to obtain results). Practical value. The difference between the best (the less complex) and the worst (the most complex) is in the row of few orders of magnitude. This result shows that choosing wrong methodology may enlarge time needed to perform calculation significantly.
Keywords: finite element method; FEM; LU; Cholesky; Gaussian elimination; decomposition methods; optimizations
Introduction
The development of a wide construction market and a desire to design innovative architectural building constructions has resulted in the need to create complex numerical models of objects having increasingly higher computational complexity.
There are a lot of numerical methods designed to provide an approximate solutions to the considered problem including finite element method (FEM), boundary element method (BEM) or finite difference method (FDM).
The article is organized as follows. The first section describes the main characteristics of Finite Element Methods for the beam element. The second point concerns the computational complexity of the FEM. In the next point, an estimate of the complexity of the individual stages of the FEM is presented. In the fourth section the
basic methods used to solve set of linear equations are presented. The next point is dedicated to the presentation and discussion on the complexity of the implemented algorithms, depending on the implemented matrix inversion method and the number of nodes used in calculations. Summary of the article is included in section six.
Purpose
The purpose of this work is to show that choosing a proper method for solving the set of equations can improve the calculation time (reduce the complexity) by a few levels of magnitude.
Methodology
FEM for the beam element: Finite element method is one of the most popular tools used in engineering calculations. One of the basic
Наука та прогрес транспорту. Вкник Дншропетровськош нацюнального ушверситету залiзничного транспорту, 2016, № 2 (62)
assumptions of this method is the abandonment of the analytical division of the area for the division into a finite number of elements [12]. The calculation takes place only at specific points, called nodes, and the results for the rest of the area are approximated using, so-called, shape functions on the basis of the results obtained for each node [10], [2], [4].
The basic steps in the FEM for the beam element can be described as follows [7], [8]:
- In the first step, we divide the beam in the given number of finite elements and select nodes. We accept both the local and the global coordinate system,
- Number of components of the displacement vector qe need to be set,
- Shape functions N need to be assumed:
N = [[, N2, N3, N4 ]
(1)
Shape functions for the beam construction can be assumed as a third degree polynomial written in the classical form as:
N1 = -3-(2x3 - 3lx2 +13 )
N2 =
-3 (( - 2l2x2 + xl3 )
N3 = -3 (-2 x3 + 3lx2 )
N4 = 1 (( -12x2 )
or in the matrix form:
N = -3Г 2x3 - 3lx2 +13, Ix3 - 2l2x2 + xl3 l3 L
(2)
(3)
(4)
(5)
-2x3 + 3lx2, lx3 -12x2 ],
(6)
The stiffness matrices Ke are determined in each node
Ke =
EL T
12 6l -12 6l
6l 4l2 -6l 2l2
-12 -6l 12 -6l
6l 2l2 -6l 4l2
(7)
Calculation of substitute nodal loads matrix
pb =}NTbydx = 1-[6,1,6, -l] (8) 0 12
Calculation of the B matrix that describes the deflection at each point of the element with the following equation:
B = LN (9)
where L is the differentiation operator described as follows:
L = -y-r dx
(10)
Substituting equations (6) and (10) to (9) gives:
B=LN= -
[12 x - 6l ,6lx - 4l2
-12 x + 6l,6lx - 2l2 ], (11)
The stiffness and nodal load matrices are aggregated to the global system with application of the boundary conditions.
We solve the equation for the balance of the structure (12) identifying unknown deflections u:
Ku = Pb
(12)
Computational complexity: A key issue regarding the performance of all algorithms, including FEM [7], is the computational complexity or, in other words, computing requirements. We call computational requirements all the necessary arithmetic operations that need to be performed during calculations.
The concept of computational complexity is also associated with the issue of the resources availability. Resources can be defined in the form of time or memory [8], [9].
- Time complexity - the measurement of the time complexity is performed by the calculation of the number of basic operations that depend on the size of the input data. Measuring the actual clock time can vary depending on the implementation of the algorithm, computing machine or a compiler used.
- Memory complexity - the memory comp-lexity expresses the amount of memory used, expressed as the number of memory cells in a function of the size of input data or expressed in bits or bytes, to perform specific computing.
2
Наука та прогрес транспорту. Вкник Дшпропетровського нацюнального ушверситету залiзничного транспорту, 2016, № 2 (62)
The creators of complexity theory are Juris Hartmanis and Richard Stearns [3]. It can be said that the larger the scale of the problem and thus the larger size of the input data, the more resources it will need to perform computing operations [12].
A common way to compare the complexity of various algorithms is to determine their asymptotic growth rate. It is a measure of the time requirements with increasing size of the input data. It can be said that it describes how quickly time requirements are changing.
To describe the asymptotic growth rate so-called O notation is used. It can be said that a function f has a complexity of O(g(n)) if there is a positive constant value c and non-negative value N satisfying the condition:
f (n)< cg (n) for n > N.
(13)
In the remainder of this article the above notation will be used to determine the computational complexity of different stages of FEM algorithm.
Estimating the complexity of the FEM algorithm: In this section we will investigate the complexity of the different stages of the FEM algorithm and we will assess its impact on the complexity of the entire algorithm.
a) Dividing the Q area into finite elements.
In the case of beam element, this step has a very low complexity, which can be estimated as O(n).
b) Adoption of shape functions vector Ne and degrees of freedom qe.
The complexity of this stage is constant and is equal to O(1). It is due to the fact that regardless the size of the problem (e.g. number of nodes), this step is performed only once.
a) Calculation of elementary stiffness matrices Ke and nodal loads vectors Qe.
This step, similar like the previous one, has the complexity of the order of O(1). This is because the number of selected nodes does not affect the complexity of calculating the elementary matrices Ke and Qe .
b) Aggregation of the stiffness matrix K and the loads vector Q for the entire object.
The complexity of this stage is linear and thus can be noted as O(n). This is related to the need to add the Ke matrix to the aggregated matrix K n-times.
b) Boundary conditions.
Because taking into account boundary conditions takes place only at the ends of the beam it not depend on the number of nodes. It is therefore assumed that this step has a complexity of the order of O(1).
a) Solving Kq = Q equation On the basis of the literature it can be stated that the complexity of this stage is in the order of O(n3).
Methods of solving equations set: As shown in section 4 the most complex part of each of the algorithms is a step designated to solve the set of equations. This problem is even more complex in the case of finite element methods because this method has to generally deal with matrices of considerable sizes, which entails the need of conducting significant number of mathematical operations. Below, collection of the most popular methods for solving linear sets of equations (and matrix inversion), which can be found in the literature is presented.
The method of algebraic complements: Application of this method requires the appointment of an inverse matrix which subsequently need to be multiplied by the vector of the right side values. The algorithm for inverse matrix calculation can be presented in a following steps [3]:
- Calculation of the determinant of the matrix
A,
- Calculation of the algebra complements of all elements of a matrix,
- Transposition of the matrix containing algebraic complements,
- Determination of the inverse matrix.
These steps can be represented by the following formula:
A"1 =-
det A
D Di
D21 D22
Dn1 Dn2
D d2
1n
Dn
(14)
where detA is the determinant of a matrix A and the
D„
D
12
Dnn represent the algebraic
complement values of consecutive elements of matrix A [12].
Gauss-Jordan method: This method is another method whose immediate goal is to determine the inverse matrix (further step in the form of
1
Наука та прогрес транспорту. Вкник Дншропетровськош нацюнального ушверситету залiзничного транспорту, 2016, № 2 (62)
multiplication of inverse matrix and right values vector need to be performed). Gaussian elimination is also called matrix method of the attached identity matrix.
To find the inverse of a matrix A the following set of equations need to be solved:
AB = I
A = L • L
(17)
where: L - the lower-triangular matrix with positive values on the main diagonal, referred to as the «square root» of positive definite matrix A.
In equation (17), the unknowns are the elements of the matrix L:
L =
I
l21 l22
□ □ □ □ □
' . □
ln1 ln 2 ' " ln
(18)
a11 AT 1 21 " l11 0 1 "lu Il 1 21
_ A21 A22 _ _ L21 L22 _ 0 LT 22 _
12
lii
l11I21
it T It + T IT
11 21 21 21 22 22
thus
(15) and
where: A - output matrix, B - inverse matrix, I -identity matrix. It is necessary for both submatrices [A 11 ] to multiply them by matrix B, to get the matrix
[AB | IB] (16)
Taking into account that B = Athe [A 11] -— [11 A-1 J is obtained. In order to get the
inverse of matrix A the sub-matrix A in the [ A11 ]
need to be converted to unitary sub-matrix with the use of elementary operations on the rows.
Gauss elimination: This method was developed by a famous German mathematician Carl Friedrich Gauss [6]. Its purpose is to bring the matrix to form a stepwise matrix (Lower or Upper triangular matrix) using elementary operations on the rows. It should be also checked the existence of solution by using the Kronecker-Capelli theorem.
Cholesky decomposition: Every positively definite matrix A can be written as [11, 5, 13]
l11 = van
L = — A
21 i 21
l11
It can be also noticed that:
T T1 = A - T T1 = A - -
22 22 22 21 21 22
- A A1
21 21
(19)
(20) (21)
(22)
Solving of equations Ax = b is performed in two stages:
Ax = b LLTx = b Ly = b LTx = y Solving the Ly=b equation
У1 =
(23)
(24)
(25)
(26)
(27)
У,- =
bj -Z к=11кУк
L
dla j = 2,3,...,n (28)
Solving the L x = y equation
~ = Уп
ln
(29)
x. =
yi -Z
к=i+1likxk
L
i = n -1,n - 2,...,1 (30)
Using the formulas for calculating the matrix product it is possible to describe the relations that allow to calculate individual elements of the matrix L. As a result:
Expression (22) is so-called the Schur complement. Due to the positive value of the expression for positive definite matrix, decomposition algorithm can be applied to the matrix diminished by first row and first column.
LU decomposition: The method of LU decomposition is to produce a matrix A as a product of two triangular matrices: lower L and upper U with the addition of zero elements, above and below, respectively, the main diagonal of the
l
Наука та прогрес транспорту. Вкник Дшпропетровського нацюнального ушверситету залiзничного транспорту, 2016, № 2 (62)
matrix [1]. This formally can be written as:
A = LU
where:
L =
l
0
l21 l22
ln1 ln 2
(31)
(32)
and
U =
0 u
22
2n
0 0
(33)
Analyzed set of equations can be described as: LUx = y (34)
Solving the equation (34) is achieved by solving two sets of equations:
Lz = y (35)
Ux = z (36)
In this way, we obtain a searched vector x.
Findings
This section provides an analysis of the computational requirements of different matrix inversion algorithms for finite element method.
The calculations were performed for the fixed-end beam of the length equal to lx = 6.5 m, loaded with a uniformly distributed load q equal to 15 kN / m, which diagram is shown in Fig. 1.
example, it has been taken into account that the matrices are symmetrical and contain a large number of coefficients equal to 0, which meant that some multiplications and additions can be omitted when performing calculations. The obtained results have been shown in the latter part of this section.
Fig. 1. Scheme of the investigated beam
For each method, the calculations of deflection of the beam, which is presented in Fig. 2, divided into 5, 50, 100 and 200 nodes have been performed. The process of calculation of the equation (12) has been carried out using the Gauss method, Cholesky decomposition and LU decomposition. In addition, each method has been further optimized. This optimization consisted in the use of some knowledge of matrices. For
Fig. 2. Deflection of the beam for 5 (top) and 50 (bottom) nodes calculated with the use of FEM and accurate methods
Table 1 presents the number of mathematical operations performed during calculation of the deflection of the beam divided into 5 nodes, depending on the implemented method for matrix inversion (solving equation set).
Second column shows the assignment operations necessary to be performed. Assignment operations are the operations of saving the calculated values in the memory of the computer (after the value is calculated it needs to be stored in the memory). The third column shows the number of addition operations and in subsequent columns, number of multiplications, divisions and square root calculations are presented. Last column presents the total number of required operations (sum of the operation of the preceding columns).
Table 1
Number of mathematical operations for different implemented algorithms, FEM methodology, division into 5 nodes
Method Gauss Opt. Gauss LU Opt. LU Cholesky Opt. Cholesky
Assignments 106 74 72 60 98 84
Additions 85 56 120 56 91 79
Multi-plications 85 56 120 56 65 54
Divisions 21 18 21 18 32 31
Square roots 0 0 0 0 6 6
Total 297 204 333 190 292 254
Hereafter the name basic methods has been used regarding the Gauss, LU decomposition and Cholesky decomposition methods while the term optimized methods has been used regarding the optimized versions of these algorithms.
Comparing the Gauss method to the optimized Gauss method it can be noticed that the number of arithmetical operations performed is reduced by about 31%. For the LU decomposition method and its optimized version complexity reduction is at the level of 43%. In the case of optimized Cholesky decomposition method compared to its non optimized version the profits is equal to approximately 13%.
While taking into account the non optimized versions of the algorithms the number of operations performed by a Cholesky method is the smallest. However, while taking the optimized versions into account more preferably is to use LU and Gauss method.
For 50 nodes, we can observe a significant reduction in the number of mathematical operations required to perform for optimized algorithms compared to basic methods. The obtained results are shown in Table 2. Reduction of the number of arithmetic operations in the Gauss method and LU decomposition compared to their non optimized methods is about 99.5%. The difference between these two optimized methods is however low and is equal to about 2.5% in favor of the LU method. Basic Cholesky method gives results more than 17% better (in terms of required complexity) compared to other basic methods. Optimization of this method gives a gain of about 96%. However, the situation is changing if optimized methods are
Наука та прогрес транспорту. Вкник Дншропетровськош нацюнального ушверситету з^зничного транспорту, 2016, № 2 (62)
compared because other methods are more than 80% less complex.
Table 2
Number of mathematical operations for different implemented algorithms, FEM methodology, division into 50 nodes
Method Gauss Opt. Gauss LU Opt. LU Cholesky Opt. Cholesky
Assignments 304096 1784 18432 1320 161408 2244
Additions 299440 1406 446880 1406 299536 14884
Multiplications 299440 1406 446880 1406 156560 1674
Divisions 4656 378 5646 378 143072 13306
Square roots 0 0 0 0 96 96
Total 907632 4974 916848 4510 760672 32204
For 100 nodes (see Table 3) application of the optimized methods result in a reduction in the number of arithmetic operations performed in relation to the basic methods, even by three orders of magnitude, giving a gain of 99.9%. It is also worth noting that optimized LU decomposition method requires 10 % less operations than the optimized Gauss method and 92% fewer than the Cholesky method.
Table 3
Number of mathematical operations for different implemented algorithms, FEM methodology, division into 100 nodes
Method Gauss Opt. Gauss LU Opt. LU Cholesky Opt. Cholesky
Assignments 2548196 3684 76832 2720 1312808 4644
Additions 2528890 2906 3783780 2906 2529086 59834
Multiplications 2528890 2906 3783780 2906 1293110 3474
Divisions 19306 778 19306 778 1236172 56556
Square roots 0 0 0 0 196 196
Total 7625282 10274 7663698 9310 6371372 124704
Table 4
Number of mathematical operations for different implemented algorithms, FEM methodology, division into 200 nodes
Method Gauss Opt. Gauss LU
Assign-ments 20856396 7484 313632
Additions 20777790 5906 31127580
Multiplications 20777790 5906 31127580
Divisions 78606 1578 78606
Square roots 0 0 0
Total 62490582 20874 62647398
End of table 4
Method Opt. LU Cholesky Opt. Cholesky
Assign-ments 5520 10585608 9444
Additions 5906 20778186 239734
Multiplications 5906 10506210 7074
Divisions 1578 10272372 233056
Square roots 0 396 396
Total 18910 52142772 489704
By analyzing the amount of necessary operations to be performed for 200 nodes (Table 4) it can be noticed that the difference between the optimized versions and the basic algorithms is in the order of four orders of magnitude (profit is more than 99.99%) for the Gauss method and LU decomposition. Profit for the Cholesky decomposition between the basic version and its optimized version reaches the value of 99%.
Similar as for 100 nodes the difference between the optimal versions of LU and Gauss is still around 10% (in favor of the LU decomposition method). Optimized Cholesky method offers the highest complexity compared with other optimized algorithms where number of required operations is more than 95% lower.
The basic version of Cholesky method is about 16% less complex compared to other methods considered.
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету з^зничного транспорту, 2016, № 2 (62)
Originality and practical value
The main achievement of this work is that it shows the impact of the used methodology on the complexity of solving the problem (or equivalently, time needed to obtain results). The difference between the best (the less complex) and the worst (the most complex) is in the row of few orders of magnitude. This result shows that choosing wrong methodology may enlarge time needed to perform calculation significantly.
Conclusions
In the article the impact of the method used for calculating the equations set on the number of arithmetic operations performed is investigated. Deflection of the beam when divided into 5, 50, 100 and 200 nodes using the finite element method has been calculated. To solve the set of equations the Gauss method, LU and Cholesky decomposition methods were implemented. In addition, for each method an optimized version of the algorithm was also investigated, what was intended to further reduce the number of computational operations performed. It was found that in order to maximally reduce the number of necessary to conduct operations at first the advantage of all known properties of the resulting matrix should be taken. Afterwards, the selection of a method that aims at the highest possible reduction of complexity in solving the system of equations should be made.
In the case of optimized methods number of operations can be reduced, by several orders of magnitude. It is also important to choose a proper method because in the basic versions it is preferably to choose Cholesky method while the worst choice is the LU decomposition. However, the situation changes for the optimized methods. The best results are observed for LU decomposition method giving a gain of 3-10% compared to the Gaussian method, and more than 90% compared to the for the Cholesky method.
It was also found that the highest differences between the complexity of the methods are observed if the beam is divided into a larger number of nodes.
Наука та прогрес транспорту. Вкник Дншропетровського нащонального ушверситету залiзничного транспорту, 2016, № 2 (62)
LIST OF REFERNCE LINKS
1. A parallel domain decomposition method for 3D unsteady incompressible flows at high Reynolds number / R. Chen, Y. Wu, Z. Yan [et al.] // J. of Scientific Computing. - 2013. - Vol. 58. - Iss. 2.
- P. 275-289. doi: 10.1007/s10915-013-9732-x.
2. Balonek, K. Wprowadzenie do metody elementow skonczonych / K. Balonek, S. Gozdur ; Akademia Gorniczo Hutnicza. - Krakow, 2009.
3. Drwal, M. Algorytmy algebry numerycznej. Rozwi^zywanie ukladu rownan liniowych / M. Drwal. - Wroclaw : Politechnika Wroclawska, 2012.
4. Dzierzankowski, G. Samouczek Metody Elementow Skonczonych / G. Dzierzankowski, M. Sitek.
- Warszawa : Oficyna Wydawnicza Politechniki Warszawskiej, 2012. - 106 p.
5. Geijn, R. A. van de. Notes on Cholesky Factorization / R. A. van de Geijn // Report TX 78712 University of Texas at Austin, Department of Computer Science. - Austin, 2011. - P. 1-16.
6. Introduction to Algorithms / T. H. Cormen, C. E. Leiserson, C. Stein, R. L. Rivest. -Cambridge : The MIT Press, 1990. - 1292 p.
7. Khawaja, H. Application of a 2-D approximation technique for solving stress analyses problem in FEM / H. Khawaja // Intern. J. of Multiphysics. -2015.- Vol. 9. - Iss. 4. - P. 317-324. doi: 10.1260/1750-9548.9.4.317.
8. Metody komputerowe w liniowej mechanice konstrukcji / C. Cichon, W. Cecot, J. Krok, P. Plucinski. - Krakow : Politechnika Krakowska, 2009. - 428 p.
9. Papadimitriu, Ch. H. Zlozonosc obliczeniowa / Ch. H. Papadimitriu. - Warszawa : Wydaw-nictwa Naukowo-Techniczne, 2002. - 540 p.
10. Rakowski, G. Metoda Elementow Skonczonych w mechanic konstrukcji / G. Rakowski, Z. Kacprzyk. - Warszawa : Oficyna Wydawnicza Politechniki Warszawskiej, 2005. - 434 p.
11. Vandenberghe, L. Applied Numerical Computing (lecture) [Electronic resource] / L. Vandenberghe. - Available at: http://www.ee.ucla.edu/vandenbe/-ee103.html. - Title from the screen. - Accessed 10.04.2016.
12. Sybis, M. Analiza zlozonosci metody roznic skonczonych (MRS) oraz metody elementow skonczonych (MES) na przykladzie elementu belkowego / M. Sybis, A. Smoczkiewicz-Woj-ciechowska, A. Szymczak-Graczyk // Projek-towanie, eksploatacja, diagnostyka i naprawy wybranych obiektow budownictwa ogolnego i hydrotechnicznego : Monografia / Instytut Budownictwa i Geoinzynierii Uniwersytet Przyrodniczy. - Poznan, 2015. - P. 135-151.
13. Younis, G. Practical method to solve large least squares problems using Cholesky decomposition / G. Younis // Geodesy and Cartography. - 2015. -Vol. 41. - Iss. 3. - P. 113-118. doi: 10.3846/20296991.2015.1086118.
М. СЫБИС1*, А. СМОЧКЕВИЧ-ВОЙЦЕХОВСКА2*, А. ШИМЧАК-ГРАЧИК3
1 Каф. «Гражданская и инженерная геология», Познаньский университет естественных наук, ул. Пиатковска, 94, Познань, Польша, 60-649, эл. почта msybis@up.poznan.pl, ОЯСГО 0000-0002-0032-6313
2*Каф. «Гражданская и инженерная геология», Познаньский университет естественных наук, ул. Пиатковска, 94, Познань, Польша, 60-649, эл. почта asw72@up.poznan.pl, ОЯСГО 0000-0002-7202-6678
3*Каф. «Гражданская и инженерная геология», Познаньский университет естественных наук, ул. Пиатковска, 94, Познань, Польша, 60-649, эл. почта agraczyk@up.poznan.pl, ОЯСГО 0000-0002-1187-9087
ВЛИЯНИЕ ОБРАЩЕНИЯ МАТРИЦЫ НА СЛОЖНОСТЬ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ
Цель. Развитие обширного строительного рынка и стремление проектировать инновационные архитектурные строительные конструкции привели к необходимости создания сложных численных моделей объектов, у которых возрастает сложность вычислений. Цель данного исследования - показать, что выбор правильного метода для решения системы уравнений может улучшить время расчета (уменьшить сложность) на несколько уровней величины (амплитуды). Методика. В статье представлен анализ влияния алгоритма обращения матрицы на расчет прогиба в балке с использованием метода конечных элементов (МКЭ). На основе анализа литературы были определены общие методы расчетных систем уравнений. Используя найденные решения, были применены методы исключения переменных Гаусса, разложения ЬИ матрицы и разложения Холецкого с целью определения влияния алгоритма обращения матрицы, используемой для решения комплекса уравнений, на количество выполненных вычислительных операций.
Наука та прогрес транспорту. Вкник Дншропетровського нацюнального ушверситету залiзничного транспорту, 2016, № 2 (62)
Кроме того, каждый из реализованных методов был в дальнейшем оптимизирован, тем самым уменьшая количество необходимых арифметических операций. Результаты. Эти оптимизации были выполнены с использованием определенных свойств матрицы, таких как симметрия или значительное число нулевых элементов в матрице. Результаты анализа представлены для разделения балки на 5, 50, 100 и 200 узлов, для которых был рассчитан прогиб. Научная новизна. Главным достижением данной работы является то, что она показывает влияние используемой методики на сложность решения задачи (или, что одно и то же, на время, необходимое для получения результатов). Практическая значимость. Разница между лучшим (менее сложным) решением задачи и худшим (более сложным) находится в ряде нескольких порядков величины. Результаты показывают, что при выборе неправильного метода может значительно увеличиться время, необходимое для выполнения расчета.
Ключевые слова: метод конечных элементов; разложение LU матрицы; Холецкий; исключение методом Гаусса; методы разложения; оптимизации
М. СИБ1С1*, А. СМОЧКЕВИЧ-ВОЙЦЕХОВСЬКА2*, А. ШИМЧАК-ГРАЧ1К3
1 Каф. «Громадянська та шженерна геологш», Познаньський ушверситет природничих наук, вул. Пиатковська, 94, Познань, Польща, 60-649, ел. пошта msybis@up.poznan.pl, ORCID 0000-0002-0032-6313
2*Каф. «Громадянська та шженерна геологш», Познаньський ушверситет природничих наук, вул. Пиатковська, 94, Познань, Польща, 60-649, ел. пошта asw72@up.poznan.pl, ORCID 0000-0002-7202-6678
3*Каф. «Громадянська та шженерна геологш», Познаньський ушверситет природничих наук, вул. Пиатковська, 94, Познань, Польща, 60-649, ел. пошта agraczyk@up.poznan.pl, ORCID 0000-0002-1187-9087
ВПЛИВ ЗВЕРНЕННЯ МАТРИЦ НА СКЛАДН1СТЬ МЕТОДУ К1НЦЕВИХ ЕЛЕМЕНТ1В
Мета. Розвиток великого буд1вельного ринку та прагнення проектувати шноващйш архггектурш буд1вельш конструкцп призвели до необхщносп створення складних чисельних моделей об'екпв, у яких зростае складшсть обчислень. Мета даного дослщження - показати, що виб1р правильного методу для виршення системи р1внянь може полшшити час розрахунку (зменшити складшсть) на кшька порядив величини (амплиуди). Методика. У статп представлений анал1з впливу алгоритму звернення матрищ на розрахунок прогину в балщ з використанням методу ск1нченних елеменпв (МСЕ). На основ1 анал1зу лггератури були визначеш загальш методи розрахункових систем р1внянь. Використовуючи знайдеш ршення, були застосоваш методи виключення змшних Гаусса, розкладання LU матрищ i розкладання Холецкого з метою визначення впливу алгоритму звернення матрищ, використовувано! для виршення комплексу рiвнянь, на шльшсть виконаних обчислювальних операцiй. Крiм того, кожен iз реалiзованих методiв був надалi оптимiзований, тим самим зменшуючи кшьшсть необхвдних арифметичних операцiй. Результати. Ц оптимшцд були виконанi з використанням певних властивостей матрищ, таких як симетрiя або значне число нульових елеменпв у матрицi. Результати аналiзу представленi для подiлу балки на 5, 50, 100 i 200 вузлiв, для яких був розрахований прогин. Наукова новизна. Головним досягненням дано! роботи е те, що вона показуе вплив використовувано! методики на складшсть виршення задачi (або, що одне й те ж, на час, необхвдний для отримання результата). Практична значимiсть. Рiзниця мiж кращим (менш складним) рiшенням задачi i гiршим (бiльш складним) знаходиться в рядi дек1лькох порядков величини. Результати показують, що при виборi неправильного методу може значно збшьшитися час, необхвдний для виконання розрахунку.
Ключовi слова: метод шнцевих елементiв; розкладання LU матрищ; Холецький; виключення методом Гаусса; методи розкладання; оптимшцд
REFERNCES
1. Chen R., Wu Y., Yan Z., Zhao Y., Cai X.-C. A parallel domain decomposition method for 3D unsteady incompressible flows at high Reynolds number. Journal of Scientific Computing, 2014, vol. 58, no. 2, pp. 275-289. doi: 10.1007/s10915-013-9732-x.
2. Balonek K., Gozdur S. Wprowadzenie do metody elementow skonczonych. Akademia Gorniczo Hutnicza Publ., Krakow, 2009.
Наука та прогрес транспорту. Вкник Дншропетровського нащонального ушверситету залiзничного транспорту, 2016, № 2 (62)
3. Drwal M. Algorytmy algebry numerycznej. Rozwi^zywanie ukladu rownan liniowych. Wroclaw, Politechnika Wroclawska Publ., 2012.
4. Dzierzankowski G., Sitek M. Samouczek Metody Elementow Skonczonych. Warszawa, Oficyna Wydawnicza Politechniki Warszawskiej Publ., 2012. 106 p.
5. Geijn R.A. van de. Notes on Cholesky Factorization. Report TX 78712 University of Texas at Austin, Department of Computer Science Publ., Austin, 2011. P. 1-16.
6. Cormen T.H., Leiserson C. E., Stein C., Rivest R. L. Introduction to Algorithms. Cambridge, The MIT Press Publ., 1990. 1292 p.
7. Khawaja H. Application of a 2-D approximation technique for solving stress analyses problem in FEM. Intern. Journal of Multiphysics, 2015, vol. 9, issue 4, pp. 317-324.
8. Cichon C., Cecot W., Krok J., Plucinski P. Metody komputerowe w liniowej mechanice konstrukcji. Krakow, Politechnika Krakowska Publ., 2009. 428 p. doi: 10.1260/1750-9548.9.4.317.
9. Papadimitriu Ch.H. Zlozonosc obliczeniowa. Warszawa, Wydawnictwa Naukowo-Techniczne Publ., 2002. 540 p.
10. Rakowski G., Kacprzyk Z. Metoda Elementow Skonczonych w mechanic konstrukcji. Warszawa, Oficyna Wydawnicza Politechniki Warszawskiej Publ., 2005. 434 p.
11. Vandenberghe, L. Applied Numerical Computing (lecture). Available at: http://www.ee.ucla.edu/vandenbe/ee103.html (Accessed 10 April 2016).
12. Sybis M., Smoczkiewicz-Wojciechowska A., Szymczak-Graczyk A. Analiza zlozonosci metody roznic skonczonych (MRS) oraz metody elementow skonczonych (MES) na przykladzie elementu belkowego. Projektowanie, eksploatacja, diagnostyka i naprawy wybranych obiektow budownictwa ogolnego
i hydrotechnicznego: Monografia. Instytut Budownictwa i Geoinzynierii Uniwersytet Przyrodniczy Publ., Poznan, 2015. P. 135-151.
13. Younis G. Practical method to solve large least squares problems using Cholesky decomposition. Geodesy and Cartography, 2015, vol. 41, issue 3, pp. 113-118. doi: 10.3846/20296991.2015.1086118.
Prof. D. O. Bannikov, Sc. Tech. (Ukraine) recommended this article to be published
Received: 29 Jan., 2016
Accepted: 31 March, 2016
