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.
+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 amThat 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!
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