I am currently working on the following scenario - a cantilever beam (rectangular section) meshed with solid elements subjected to a concentrated force at the free end. Bending occurs in the elasto-plastic regime (elastic perfectly-plastic material model). My goal is to obtain the deflection as close to the analytical one as possible using solid elements. Here are the current outcomes:
- analytical solution: 23,625 mm
- beam elements in Abaqus (40 x B22): 23,69 mm
- solid elements in Abaqus (7000 x C3D8R): 23,36 mm
- beam elements in CalculiX (100 x B32): 68,3 mm !!!
- beam elements in CalculiX (100 x B32R): 22,6 mm
- solid elements in CalculiX (7000 x C3D8R): 23,27 mm
- solid elements in CalculiX (26451 x C3D10): 23 mm
As you can see, the result obtained with B32R elements in CalculiX is not so good and with B32 elements (full integration) I get completely wrong value. That’s probably because of shear locking.
Solid elements provide significantly better results (especially C3D8R) but they could be a bit more accurate.
Anyway, what I would like to ask you is whether there are some additional element types that could be better for this type of application and how far I should go with mesh refinement.
Also, does it matter if I use rigid body constraint or coupling to transfer the concentrated force to the front surface of the beam ?
Here’s the input data:
- dimensions: 200x140x2000 mm
- material - steel: Young’s modulus 210 GPa, yield stress 235 MPa
- load: 155 kN