Creating the finite element method (FEM) with python ~ vba → python translation ~

Series table of contents

--Creating the finite element method (FEM) with python ~ vba → python translation ~ -Creating the finite element method (FEM) with python ~ damyarou's messing around ~ --Creating a finite element method (FEM) with python-matrix solution- (planned)

Overview of vba → python translation

If you explain the concept roughly Let's create the finite element method (FEM) in python

First of all, let's say ~ vba → python translation ~ [this book](https://rpx.a8.net/svt/ejp?a8mat=2TINC2+D4Q7VE+2HOM+BW8O1&rakuten=y&a8ejpredirect=http%3A%2F%2Fhb.afl .rakuten.co.jp% 2Fhgc% 2F0eac8dc2.9a477d4e.0eac8dc3.0aa56a48% 2Fa08100810103_2TINC2_D4Q7VE_2HOM_BW8O1% 3Fpc% 3Dhttp% 253A% 252F% 252Fbooks.rakuten.co.jp% 252Fbooks.rakuten.co.jp%252Fbooks. .rakuten.co.jp% 252Frb% 252F1700147% 252F) Excel sample 5-1 Three-dimensional skeleton analysis .xls (hereinafter sample) will be rewritten in python.

[Introduction to Excel Finite Element Method (Excel Civil Engineering Course)](https://rpx.a8.net/svt/ejp?a8mat=2TINC2+D4Q7VE+2HOM+BW8O1&rakuten=y&a8ejpredirect=http%3A%2F%2Fhb .afl.rakuten.co.jp% 2Fhgc% 2F0eac8dc2.9a477d4e.0eac8dc3.0aa56a48% 2Fa08100810103_2TINC2_D4Q7VE_2HOM_BW8O1% 3Fpc% 3Dhttp% 253A% 252F% 252Fbooks.rakuten.co.jp % 252F books.rakuten.co.jp% 252Frb% 252F1700147% 252F)

The translation is as follows [Sample](https://rpx.a8.net/svt/ejp?a8mat=2TINC2+D4Q7VE+2HOM+BW8O1&rakuten=y&a8ejpredirect=http%3A%2F%2Fhb.afl.rakuten.co. jp% 2Fhgc% 2F0eac8dc2.9a477d4e.0eac8dc3.0aa56a48%2Fa08100810103_2TINC2_D4Q7VE_2HOM_BW8O1% 3Fpc% 3Dhttp% 253A% 252F%252Fbooks.rakuten.co.jp%252Frb%252F252F2 jp% 252Frb% 252F1700147% 252F) → python will be translated.

5-1 3D skeleton analysis.xls

Private Sub cmd calculation_Click()

    Dim I As Long
    Dim NEQ As Integer

    For I =1 To array upper limit
        AJCB(I) = 0
    Next

Call data entry
Call SKY Matrix(NEQ)
    For I = 1 To NEQ
        FORCE(I) = 0
    Next I
Call distributed load distribution
Call external force add
    Call decomp(NEQ)
    Call redbak(NEQ)
Call displacement calculation
Call force and moment calculation

Call result output

    MsgBox "Analysis is finished"

End Sub

framecalc.py

from Calcrate import Calcrate

if __name__ == "__main__":

    cal = Calcrate()
    NEQ = 0

    cal.Data entry()
    cal.SKY matrix(NEQ)
    cal.Distributed load distribution()
    cal.External force add()
    cal.decomp(NEQ)
    cal.redbak(NEQ)
    cal.Displacement calculation()
    cal.Force and moment calculation()

    cal.Result output()

    print("Analysis is finished")

Environmental condition

My environment is as follows.

OS

• windows 7 Home premium Sp1

Development environment

• python 3.5.2
• Anaconda 4.2.0(64-bit)
• pandas 0.18.1
• numpy 1.12.0
• tensorflow 0.12.0re1

Let's move it first

First, download the source code.

[Sample](https://rpx.a8.net/svt/ejp?a8mat=2TINC2+D4Q7VE+2HOM+BW8O1&rakuten=y&a8ejpredirect=http%3A%2F%2Fhb.afl.rakuten.co.jp%2Fhgc%2F0eac2 9a477d4e.0eac8dc3.0aa56a48%2Fa08100810103_2TINC2_D4Q7VE_2HOM_BW8O1% 3Fpc% 3Dhttp% 253A% 252F% 252Fbooks.rakuten.co.jp%252Frb%252F1700147%252F%26m%3Dhttp% 253A%252F%26m%3Dhttp% 252F) Please download 5-1 3D skeleton analysis.xls according to the guidance in the book. We will translate this Excel VBA into python.

The source translated into python is here. $ git clone -b "#1" https://github.com/sasaco/FramePython.git

When the environment is ready, move to the source code directory and hit framecalc.py to start the analysis. If the message "Analysis is complete" is displayed, the analysis was performed correctly.

python

$ cd FramePython

$ python framecalc.py

$Analysis is finished

python program specifications

About data input / output

[Sample](https://rpx.a8.net/svt/ejp?a8mat=2TINC2+D4Q7VE+2HOM+BW8O1&rakuten=y&a8ejpredirect=http%3A%2F%2Fhb.afl.rakuten.co.jp%2Fhgc%2F0eac8dc2 9a477d4e.0eac8dc3.0aa56a48% 2Fa08100810103_2TINC2_D4Q7VE_2HOM_BW8O1% 3Fpc% 3Dhttp% 253A% 252F% 252Fbooks.rakuten.co.jp% 252Frb% 252F1700147% 252F% 26m% 3Dhttp% 253A%252F%26m%3Dhttp% 252F) Excel uses sheet ** [input 1] [input 2] ** as input data. The analysis result is stored in the sheet ** [Calculation / Output] **.

The translated python program does not work with the excel sheet, so [Sample](https://rpx.a8.net/svt/ejp?a8mat=2TINC2+D4Q7VE+2HOM+BW8O1&rakuten=y&a8ejpredirect=http%3A%2F%2Fhb.afl.rakuten.co.jp%2Fhgc%2F0eac8dc2 9a477d4e.0eac8dc3.0aa56a48% 2Fa08100810103_2TINC2_D4Q7VE_2HOM_BW8O1% 3Fpc% 3Dhttp% 253A% 252F% 252Fbooks.rakuten.co.jp%252Frb%252F1700147%252F%26m%3Dhttp% 253A%252F%26m%3Dhttp% 252F) The file obtained by converting the Excel sheet ** [Input 1] [Input 2] [Calculation / Output] ** to csv is used for input / output.

-Input 1.csv -Input 2.csv -[Calculation / Output.csv](https://github.com/sasaco/FramePython/blob/%231/%E8%A8%88%E7%AE%97%E3%83%BB%E5%87%BA % E5% 8A% 9B.csv)

About array subscripts

The sequence of vba is 1 ~

5-1 3D skeleton analysis.xls

Range(1,1)

The python array is 0 ~

Calcrate.py

Range[0][0]

So I had to be careful when translating.

About drawing function

[Sample](https://rpx.a8.net/svt/ejp?a8mat=2TINC2+D4Q7VE+2HOM+BW8O1&rakuten=y&a8ejpredirect=http%3A%2F%2Fhb.afl.rakuten.co.jp%2Fhgc%2F0eac8dc2 9a477d4e.0eac8dc3.0aa56a48% 2Fa08100810103_2TINC2_D4Q7VE_2HOM_BW8O1% 3Fpc% 3Dhttp% 253A% 252F% 252Fbooks.rakuten.co.jp% 252Frb% 252F1700147% 252F% 26m% 3Dhttp% 253A%252F%26m%3Dhttp% 252F) Excel drawing function is not ported to python

Comparison of analysis results

For the time being, the translation work of vba → python has been completed, so I compared the analysis results of both. A small value of almost 0 has an error, but isn't it a fair reproducibility?

Sample Excel analysis results (excerpt)

Analysis result of python program (excerpt)

next time

For the time being, I was able to get similar results.

next time, --Separate the analysis part class ――The matrix solution is like python

I am thinking.

This article is an ongoing project, so Not completed yet.

The current program state is in the master branch. https://github.com/sasaco/FramePython/tree/master You may have started fixing some of the above

Please contact us if you are interested.