Example of 3D skeleton analysis by Python
The latest version of the program is posted below, so please refer to it.
Introduction
I will introduce an example calculated by a three-dimensional skeleton analysis program using Python. Here, I tried to reproduce the grid girder analysis result and the plane frame analysis result. I was wondering if the 3D coordinate conversion was working, but as a result, it seems to be working.
Until now, plane skeleton analysis and grid girder analysis were used properly according to the structure to be designed, but if the three-dimensional skeleton analysis program can be used freely, this one can play both roles. It can be done, so it will be convenient.
But in reality, the grid girder analysis will be done by 3D skeleton analysis, but the plane skeleton analysis will use a plane program. This is because creating input data is a little troublesome for 3D objects.
The program used was the one introduced last time. http://qiita.com/damyarou/items/2b83b9b20181c1b8e35d
The program, input data, analysis procedure and section force drawing script by GMT are uploaded to the following Gist. https://gist.github.com/damyarou/8b99e5e57576db5e71bf630fa1b396a2
About the program
The three-dimensional frame structure analysis program by the finite element method is Python3, but the section force diagram is created by GMT (Generic Mapping Tools). The reason is that the 3D plot of Python-matplotlib is not cool. matplotlib is good for 2D plots, but I feel that it suddenly becomes cartoonish when it comes to 3D. In that respect, GMT is troublesome to handle, but it is possible to draw in a unified manner in both 2D and 3D, and the product feels good. Both are a matter of taste. In order to plot in GMT, a Python program that reads the output data of the analysis program and creates a data file for plotting in GMT is created separately.
Analysis of in-plane grid girders in the global coordinate system XY
Azimuth looking at the 3D plot: 150 degrees from north, 30 degrees elevation Code angle: All 0
Analysis of multi-layer ramen in the whole coordinate system Y-Z plane
Azimuth looking at the 3D plot: 120 degrees from north, 30 degrees elevation Code angle: All 0
Analysis of multi-layer ramen in the whole coordinate system X-Z plane
The dimensions, cross-sectional characteristics, and loads of the structural system are exactly the same as the above-mentioned model in the Y-Z plane of the overall coordinate system, but when making a three-dimensional model, the structure is placed in the X-Z plane. In this case, the same cross-sectional force diagram as the above model is obtained by setting the cord angle to -90 degrees for the column member.
Azimuth looking at the 3D plot: 150 degrees from north, 30 degrees elevation Code angle: 0 for beam members, -90 degrees for column members
It is necessary to consider the relationship between the column member and the beam member when setting the code angle. The following documents will be helpful.
http://www.archi.hiro.kindai.ac.jp/laboratory/SAL/dfujii/Report/fem/fem_11.pdf
that's all
Recommended Posts