Hexahedral Mesh Refinement in Salome

     In this video I show some ways to refine the hexahedral mesh of the beam we did before. This procedure is meant to be more clever than simple making the "number of segments" parameter bigger and bigger.

Thanks again!

Cantilever beam finite element analysis using open source software

     In this video I'm modeling, meshing and computing a cantilever beam using open source software. For the model and mesh, I've used Salome7, and the finite element computations were carried out with Calculix. This is meant to be an practical introductory video to free fea software. Theres an insane amount of theoretical books about this subject out there, take a look on them before designing your own cranes (LoL). 

The code for the calculix input deck:

Cheers!

Plotting Beam Diagrams with Python

     Today I'm sharing how to plot Beam Diagrams using python. More specifically the code is about a cantilever beam subjected to gravity load, a special case of uniformly distributed load. There's a bunch of cool python features in this post, like the very convenient use of "x>" to evaluate a function, or part of it, only after a specified value in the array.

For example, the code: (my blog font shows this symbol =, as equals, this is a minus -)

import numpy as np
x=np.linspace(1,10,10)
y=2*x*(x>5)
print x
print y

returns these two arrays.

[  1.   2.   3.   4.   5.   6.   7.     8.   9.    10.]
[  0.   0.   0.   0.   0.  12.  14.  16.  18.    20.]

     Values with "false" boolean condition returned 0, and the values with "true" condition were evaluated. This is very useful when computing singularities functions, also known as Macaulay method, due to the British mathematician William Herrick Macaulay.  Oddly this method was first proposed by the German mathematician Alfred Clebsch. 

    Other features that worth take a look is how to make very neat plots using matplotlib subplots method.

The code:



The plot:

See ya!