Pyramid elements C3D5 & C3D13

21

under the GNU profiler to recognize in which routines the C3D15 have been used, further I also did a complete search for all places where C3D15 had been implemented in

Clever way! That would have given you a lot more confidence than trying to walk through it all by hand.

I don’t think it’s too big of a deal if there are some bugs. They’ll get discovered soon enough and can be fixed one by one. If it makes it into official CCX, I’ll be making my own set of test cases too.

Thanks for spending the time on such a useful (from my POV) piece of work!

22

@vicmw
During the implementing of the pyramid element I got by incident implement an error in the 2’order wedge element which I just have discovered with a test of the pre-tension card.

Wedge_error
Wedge_error707×600 85.6 KB

I have updated the code so the result now should be correct.
Wedge_fixed
Wedge_fixed733×624 87.5 KB

I actual don’t who are using my pyramid code, but if you should need the updated source code and the windows exe you can get it at the link in the beginning in this thread.

23

I may want to test them later this year… thank you @fgr!

24

I’m not using it but I have code to generate / read solution with these new elements for CCX and hope Guido includes it.

I just ran my big set of test cases on it and found these failures. I only had a quick look so let me know if you need more details.

Thermal doesn’t work but I’m not sure if you intended that or not.

CCX S30 density.inp CCX_S30_density.inp - Pastebin.com C3D5 Small difference with with gravity compared to collapsed hex and Mecway’s internal pyramids. Not sure if it’s worse or better.

CCX S30 Pressure pyramids.inp CCX_S30_Pressure_pyramids.inp - Pastebin.com Bad stress and displacement

CCX S30 solid element force.inp CCX_S30_solid_element_force.inp - Pastebin.com Slight increase in error of reaction force over collapsed hex pyramids.

CCX L30 compression only support cylinder rotating.inp CCX_L30_compression_only_support_cylinder_rotating.inp - Pastebin.com C3D13 won’t converge whereas it would for collapsed hex pyramids.

1 Like
25

When I recently e-mail Guido, I asked him to take a look at this code. He said that he would do it.

3 Likes
26

@vicmw
I really appreciate your test files. It helps a lot to have to alternative input in order to make verification test. :+1:

I’ll return when I have looked into your test files, which I really have had time to look into yet, but in advance I also would like to refer from “Advanced Finite Element Methods” by Carlos Felippa, that some discrepancy should be taken into account when running tests with few elements.

1 Like