How to build a FEA program in the shortest time?

More Fan Shen's questions See All

What's the different between Geometrically Exact and Absolute Nodal Coordinate Formulation ?

In my view, GE is in the purpose of raise an analytical model for highly flexible problem, however, there's some inaccuracy in it. ANCF is a FEM method completely, which can obtain an accurate...

01 February 2016 7,419 2 View

Is there any study in high-dimensional matrix?

3D, or 4D matrix, what form their multiply operation or inverse matrix will be? Is there any paper study them or use the studies about them?

01 February 2015 4,816 12 View

Is there any open-source FEA software that is good for non-linear plates or membranes?

I want to do a simulation about a plate or membrane-based model. Because of the non-linearity, I tried Abaqus - that's good. In many papers, NASTRAN was used, and the answer was similar to Abaqus....

03 April 2014 9,594 2 View

How to choose a Free FEA program to follow?

There are many free FEM programs under GPL, such as DUNE, FreeCASE, z88, Hermes, etc. Which one should we choose to follow? Have there been any articles which compare them?

11 December 2013 5,626 5 View

error in running a model in Abaqus

Excessive rotation of nodes in node set ErrNodeExcessRotation-Step1

17 August 2021 0 0 View

error in modelling model in Abaqus

Excessive rotation of nodes in node set ErrNodeExcessRotation-Step1

17 August 2021 0 0 View

MGL tools

20 June 2021 0 0 View

what is the simplest definition of ion exchange capacity

What this unit identifies for the ion exchange capacity ca. 0.92 meq. g-1. I am a newcomer in the membrane area. I just started. I highly appreciate your help

08 May 2021 0 0 View

what is the simplest definition of ion exchange capacity

What this unit identifies for the ion exchange capacity ca. 0.92 meq. g-1. I am a newcomer in the membrane area. I just started. I highly appreciate your help

08 May 2021 0 0 View

i have some problem with chtMultiRegionSimpleFoam

--> FOAM FATAL ERROR: Maximum number of iterations exceeded From function Foam::scalar Foam::species::thermo::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo = Foam::species::thermo >, Foam::sensibleEnthalpy>] in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::heRhoThermo >, Foam::sensibleEnthalpy> > > >::calculate() at ??:? #3 Foam::heRhoThermo >, Foam::sensibleEnthalpy> > > >::correct() at ??:? #4 ? at ??:? #5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #6 ? at ??:?

24 March 2021 0 0 View

How to compare two groups with only two measurements?

Hiiiii everyone! I have an enquiry on statistical analysis. I was looking for many forum and it's still cannot solve my problem. I want to compare means of two groups of data but only with two...

03 March 2021 8,796 3 View

Can we classify HV and LV lines based on line length?

I have dataset which shows the length of power lines. I need to classify the lines based on the line length. Is there a rule to classify the High voltage (HV) and low voltage (LV) lines based on...

03 March 2021 4,116 4 View

Enhanced Yellow Fluorescent Protein (Aequorea victoria) DNA sequence?

I am on the lookout for the Enhanced Yellow Fluorescent Protein (Aequorea victoria) DNA sequence. Does anyone know where I can find it? Thank you in advance

03 March 2021 3,568 1 View

What kind of fluid could I use for a pitfall trap that also does not invalidate molecular testing?

Hi, I want to start testing pitfall trap to obtain ants samples, but I need to conduct molecular analysis on those insects. So, what kind of fluid can I use? Ethanol expires too early and I need...

03 March 2021 5,978 5 View

Kari Tammi Popular answer

I'd be too lazy to code own, since there are so many available. How about trying Elmer (http://www.csc.fi/english/pages/elmer )? It has full capability for multiphysics problems and even parallelisation.

Just my few cents, WBR, Kari

Alberto Ortona

Check this:

http://www.springer.com/engineering/computational+intelligence+and+complexity/book/978-1-4020-9199-5

Vinícius Gabriel Segala Simionatto

Hi, Fan Shen,

I don't believe that writing your FEA program in Matlab is a bad idea. Matlab provides some handy solutions for quick software development. Also, I, myself, have already written a FEA software in Matlab.

However, there are two drawbacks. The first is that if your system is large you'll need more performance than Matlab can provide you. Then, for high performance applications, Fortran would be the best choice, although it is not a user friendly language at all.

The second drawback is that Matlab supports object oriented programming, but it is not as well implemented as the simple m-code, so if you want to develop a object oriented program, things would be a little more confusing, and in this very case I'd tell you to try using Java or Python, or evn C++. However, in these three languages you'd have to start from the very beginning, by implementing mathematical data structures and so on.

Best wishes,

Vinícius

Kari Tammi

Troy Tamas

I've also built a FEM solver by hand, but wouldn't recommend it.

There's a lot of packages out there that will give you a good head start. I thought this one looked pretty good to get going quick in your choice of language:

http://download.gna.org/getfem/html/homepage/index.html

Also, "Trilinos" is a more comprehensive set of packages by Sandia Labs that should get you going through pretty much anything, with probably a bit higher of a learning curve:

http://trilinos.sandia.gov/packages/

However, I would look seriously at the Sundance package. It looks like it could get you going quick, if you're on a Linux machine. (I tried once to build once on Windows, but it's a pain in the ass.)

Get your linear solver and (maybe) mesh generator in a high-performance language like C++, and the speed of all other operations should be insignificant.

Best,

Troy

Fan Shen

I decided to write my FEM with MATLAB, and have a try with the Object-Oriented style, thank Vinicius and Kambiz a lot. And I am also glad to receive wonderful ideas from Kari and Troy, my work had to run on Windows, but they will still give me great help. Also thanks to Alberto, I would consider that book. Thanks again!