This is not really my cup of tea but I would have a question: Do Ansys and Code Aster have knowledge of the system beyond its discrete finite element representation? I.e. do they know from the CAD geometry that the normal vector at point of load is exactly (1,0,0)? If you have just access to the discrete mesh this information is not trivially available. If you have access to the CAD model you can always retrieve the exact normal and you're not really comparing the same thing!
Also the ShellSolver of Elmer may use exact normal directionality data. This is unfortunately not something that is easily available since it is usually lost in the process of meshing. You will find numerous consistency tests where "mesh.director" data is given to enable more exact solution:
Code: Select all
Shell_And_Solid2/shell_and_solid/mesh.director
Shell_And_Solid/shell_and_solid/mesh.director
Shell_BenchmarkCase1_Quad/p1mesh3/mesh.director
Shell_BenchmarkCase1_Tria/tria3/mesh.director
Shell_BenchmarkCase2_Quad/p1mesh3/mesh.director
Shell_BenchmarkCase2_Tria/tria3/mesh.director
Shell_Cantilever/strip/mesh.director
Shell_Eigenanalysis_Cylinder/fan_blade/mesh.director
Shell_Eigenanalysis_Spherical/piece_of_sphere/mesh.director
Shell_OpenHemisphere/hemisphere/mesh.director
Shell_PinchedCylinder/Mesh/mesh.director
Shell_PlateCase_Quad/plates/mesh.director
Shell_with_Beam_Sections/strip_with_beams/mesh.director
Shell_with_Solid_BenchmarkCase1b/cylinder/mesh.director
Shell_with_Solid_BenchmarkCase1/cylinder/mesh.director
Shell_with_Solid_BenchmarkCase2/cylinder/mesh.director
Shell_with_Solid_Eigenanalysis/cylinder/mesh.director