fcFEM - FEA from start to finish
Moderator: bernd
Forum rules
and Helpful information for the FEM forum
and Helpful information for the FEM forum
Re: fcFEM - FEA from start to finish
... and how about uniform pressure on a curved surface. E=1000MPa, p=1MPa, D=10mm. Theoretical and calculated increase in diameter 0.01mm
Re: fcFEM - FEA from start to finish
To be precise: calculated change in diameter: 0.0100001941635mm, so an error of 0.001941635%. Not bad considering the 10-node terahedron has 6-node parabolic (and not spherical) faces.
Re: fcFEM - FEA from start to finish
That will be absolutely amazing! I'm currently involved in scripting compilation of CalculiX and seeing things like fortran '77, spooles/arpack code from the 90's hacked together with some shady perl scripts and shell macros is like digging in an old basement. There might be dead bodies I know that CalculiX won't be replaced any time soon, but having internal solver would be so much better!
Re: fcFEM - FEA from start to finish
... and much easier to extend with special features.
Re: fcFEM - FEA from start to finish
+1 and On one side CalculiX has a good community and lot of people use it, but the development still is a one man show of guido (even after over 20 years of development (AFAIK)). This is not future save IMHO.PrzemoF wrote: ↑Fri Feb 08, 2019 7:44 am That will be absolutely amazing! I'm currently involved in scripting compilation of CalculiX and seeing things like fortran '77, spooles/arpack code from the 90's hacked together with some shady perl scripts and shell macros is like digging in an old basement. There might be dead bodies I know that CalculiX won't be replaced any time soon, but having internal solver would be so much better!
Re: fcFEM - FEA from start to finish
I think it is open, but writing code in fortran 77 is some serious mental exercise. A random code sample from CalculiX:
For me it's just one step away from https://en.wikipedia.org/wiki/Brainfuck
Code: Select all
if(nelem.eq.nelem1) then
!
ang1s=(1.41d0-0.00594*alpha1)*alpha1*pi/180
!
cang1s=dcos(ang1s)
!
if(alpha1.le.22.5) then
lam11=0.0712*alpha1**0.7041+0.37
lam12=0.0592*alpha1**0.7029+0.37
else
lam11=1.d0
lam12=0.9d0
endif
zeta=lam11+(2.d0*lam12-lam11)*(V1V0*a0a1)**2
& -2d0*lam12*V1V0*a0a1*cang1s
zeta=zeta*(W0W1)**2
!
elseif(nelem.eq.nelem2) then
Re: fcFEM - FEA from start to finish
I am not planning to replicate Calculix. It’s just that I have many wishes on my todo list that would be impossible to implement without full access to and understanding of 3rd party code. I am now implementing buckling and stress stiffening (the opposite of bucking) and interface elements. This would take me months Calculix and days in something as simple as fcFEM.
Re: fcFEM - FEA from start to finish
i think Brainfuck (i had never heard about it ) is worsePrzemoF wrote: ↑Fri Feb 08, 2019 11:31 am For me it's just one step away from https://en.wikipedia.org/wiki/Brainfuck
would be great if Calculix comes at github. would speedup the development i guess