WANtaroHP: Python FEM Index for English

  Last update: 22 May 2017 at 11:42:52 (since 08 December 2016) pic


Thank you for your visit to WANtaroHP

Nickname of webmaster of this site is WANtaro. Generally in Japan, the website is called Home Page. The page title WANtaroHP means the website (Home Page) of WANtaro.

The purpose of this site is to store the FEM programs and auxiliary programs by Python made by WANtaro. From this site, you may be able to find sample programs of FEM. Sample programs are provided as text files of source code by Python3. If you want to execute these programs, you need to install Python3, numpy, scipy and matplotlib by yourself. Furthermore, GMT 5 (Generic Mapping Tools) and gnuplot are used for drawings, and ImageMagick and TeX are used for making materials for this website.


Programming environment
ItemDescription
Machine MacBook Pro (Retina, 13-inch, Mid 2014)
Processor 2.6GHz Intel Core 5
Memory 8GB 1600MHz
Storage SSD 128GB, SD card 128GB
OS macOS Sierra Version 10.12.3
Software Python 3.6.0 (numpy, scipy, matplotlib),

GMT 5, gnuplot 5, ImageMagick,

TeX, Geany, CotEditor

Published FEM programs are shown below.


Published FEM programs
FEM programType of elementSolver of equation
2D Frame Analysis 2 nodes / element, 3 DOF / node, beam element x=numpy.linalg.solve(A,b)
2D Truss Analysis 2 nodes / element, 2 DOF / node, truss element x=numpy.linalg.solve(A,b)
Grid Girder Analysis 2 nodes / element, 3 DOF / node, beam element x=numpy.linalg.solve(A,b)
2D Stress Analysis 4 nodes - 4 Gauss points / element, 2 DOF / nodex=numpy.linalg.solve(A,b)
Axisymmetric Stress Analysis 4 nodes - 4 Gauss points / element, 2 DOF / nodex=numpy.linalg.solve(A,b)
2D Saturated-Unsaturated Seepage Flow Analysis4 nodes - 4 Gauss points / element, 1 DOF / nodex=numpy.linalg.solve(A,b)
2D Thermal Conductivity Analysis 4 nodes - 4 Gauss points / element, 1 DOF / nodeA-1=numpy.linalg.inv(A)
2D Frame Response Time History Analysis 2 nodes / element, 3 DOF / node, beam element A-1=numpy.linalg.inv(A)
2D Frame Vibration Eigenvalue Analysis 2 nodes / element, 3 DOF / node, beam element w,vr=scipy.linalg.eig(K,M)
2D Frame Geometrically Nonlinear Analysis 2 nodes / element, 3 DOF / node, beam element x=numpy.linalg.solve(A,b)
2D Frame Buckling Analysis 2 nodes / element, 3 DOF / node, beam element w,vr=scipy.linalg.eig(K,G)
1D Thermal Conductivity Analysis 2 nodes / element, 1 DOF / node, linear elementA-1=numpy.linalg.inv(A)



About

I am a civil engineer and work for a civil engineering consultant company. My main job is the design and construction supervision of civil structures in hydro-power plants such as dams, pressure tunnels, steel penstocks and power stations.

Just now, I am staying in Malaysia and designing a power station building in a hydro power plant. Although my mother country is Japan and I like Japan of course, I like Malaysia also and I am enjoing Malaysian life.

By the way, in this website, some programs made by me are shown. These programs have sometimes used for my job. Generally, versatile programs that were developed by the company engineers have been used for structural analysis in the jobs. Because those have been maintained by company engineers and user can ask them about his troubles if one has any troubles about those programs.

On the other hand, if we have small own programs, it is very convenient in case of doing small model analysis. Therefore, I want to have own programs for my job and my study. In addition, I think that many people use a spreadsheet software like Exel to solve complicated problems in their jobs these days. In my opinion, to use a spreadsheet software is almost the same as to make program using program languages such as Python or Fortran or C . And it is easy for me to make programs and to solve the problems using program language comparison with solving the problems using a spreadsheet software. These are the reasons why I make programs by myself.

Although I used Fortran for programing previously, I am using Python mainly now. Because the numerical calculation libraries of Python is very useful and convenient for numerical calculation.

Please enjoy to see my website.

Although this site is for my private purpose, if you can get something that is helpful to you from this site, I will feel happiness.


Thank you.



Notice

  • This website is written by HTML5, and JavaScript is called in order to enable IE to understand HTML5 elements.
  • MathJax is used for indication of mathematical equations in this website. Waiting time for indication varies depending on the PC environment.
  • Information contained on this site may involve inaccurate descriptions and typos.
  • Webmaster may change the articles of this website without notice.
  • Webmaster of this site assumes no responsibility whatever for any direct, indirect, special, incidental, consequential damages and any other damages resulting from the use of this site.
  • You can link and unlink to this site freely.


Contact

e-mail adress of WANtaro
wantaro@jcom.home.ne.jp
damyarou@yahoo.co.jp


inserted by FC2 system