Open joneumann opened 5 years ago
Hi Johannes,
sorry for the late reply. We are currently completely reworking oFEM and I didn't really check the Git recently.
If you are still interested send a reply and I will try to help you.
Best regards, Michael
Hi oFEM-Team, is the project still alive?
Is the gong.inp from your publication available somewhere?
I tried to run testelastic.m, but failed miserably. It seems that you need at least to change the code to the following lines:
% oFEM elem = ofem.finiteelement.P1; gaussRule = ofem.gaussianquadrature(mesh, elem); solver = ofem.development.elastic(mesh,elem,gaussRule);
However, then the code fails in obj.mesh.jacobiandata (running myplate2.inp) because there are all-zero lines in obj.el. The lowest element number is 4961, hence there are 4960 all-zero lines.
Cheers, Johannes