I couldn't spot what's the malfunctioning is but the submitted code is configured to work with LU decomposition. P.S.: There is also the preconditioned conjugate gradient embedded and implemented in the source code, but for some reasons preconditioner is filled with zeros and so the solver throws exception. the coordinates, modulus of elasticity and force units. I triple checked the code for unit consistency, e.g. Xiong Zhang, Tsinghua University, FORTRAN 90). The thing that needs to be noted from results is : At the node 4, vertical -345.53 force has been applied an expected vertical displacement should be approx -3.08 but I get -0.4403 from the solution. OpenSTAP is a Finite Element Method (FEM) Solver based on STAP (KJ Bathe, FORTRAN IV), and modified version STAP90 (Prof. The program mainly uses the "Archibald Archambaud" second link for obtaining the stiffness matrix of triangular element, composes the global assembly matrix and uses the LU decomposition to solve the linear algebraic equations. There is also the input file (input_file.dat) which feeds the input for 893 noded 1672 triangular element plate. I've managed to purge the code, at the attachements you can find the source code written in fortran.
0 Comments
Leave a Reply. |