Problem with thin plate calculation

I apologize for any translation errors.

I created a simple plate model with holes and a uniform load. I can’t obtain nonlinear analysis results for the thin plate.

B x L = 1.2 x 0.8 m, t = 5 mm, P = 1 kPa

The linear analysis of the plate gives good results.

Converting the model to a solid gives incorrect results in both the nonlinear and linear ranges.

I tried to create a better mesh, but I can’t.

Reference results from without holes:

Results from PrePoMax for 2D:

Results from PrePoMax for 3D:

Szklane 3D.pmx (6.3 MB)

Szklane 2D.pmx (3.2 MB)

Rigid body constraints shouldn’t be used with shells in nonlinear analyses, this leads to non-convergence and even an immediate error in CalculiX 2.23.

Regarding the solid model, you have C3D6 and C3D8 elements - they can overstiffen the model in bending. Try at least C3D8R or C3D8I, but I would just switch to second-order elements altogether. Some global refinement and better meshing algorithm (Gmsh with hex-dominated mesh generation via extrusion or sweep with quasi-structured quad algorithm) may also help. In your model, it seems that you just used the Thicken Shell Mesh tool on the ready shell mesh, right ? This is ok if you generate a good shell mesh using Shell Gmsh and surface partitions to make use of the transfinite 3- and 4-sided faces.

What supports do you want to model here ? Why do you leave the rotational DOFs of the rigid body rerefence nodes unconstrained and why don’t you apply the BCs directly on the edges or faces of the holes ? Is there a description of the model from another FEM software available for reference anywhere ? And are the reference results from linear or nonlinear analysis (what nonlinearities ?) ?

You could also utilize symmetry here and model just a quarter of this plate. Just be careful with symmetry BCs on shells - fixing the drilling DOF leads to overconstraint and even non-convergence.

1 Like

The reference result is from a linear analysis and is consistent with many results I’ve obtained using various programs.

The support is a rotary ball bearing. Nonlinear analysis should yield smaller deformations in this case.

There might be differences in how supports are defined. Here, you are using a rigid body constraint, so making the edge/surface of the hole completely rigid and following a single reference node whose translations are fixed (but rotations aren’t). Some software may use flexible spiders instead (CalculiX has them too; they are called distributing couplings, but they have various limitations and aren’t available in PrePoMax yet).

I assume that material properties are the same in all cases (you used a material from PrePoMax’s library, that’s why I mention it).

Then, for nonlinear cases, you are only considering geometric nonlinearity ?

You could try with this model (refined 2nd order mesh, using symmetry):

Plate with hole symm - solid.pmx (7.6 MB)

Plate with hole symm - shell.pmx (679.1 KB)

Now it’s most likely just a matter of using the same supports as in other software (keeping in mind the differences between solid and shell elements).

1 Like

I’m aware of the RIGID BODY limitation. I replaced the hole with a rigid plate and TIE connection. The deformation results are acceptable. The model cannot use symmetry because it’s training for asymmetric tasks.

I’m using version 2.2. You’re using 2.4 of PrePoMax.

The material is chosen at random. Ultimately, it will be glass.

The edges will be on rigid rotules and those further on on a quasi-rigid ball joint.

Do you mean some area around the hole used for rigid body constraint ? Then It could be done with a continuous mesh too:

They may improve if you refine the mesh and use second-order elements or at least ones with reduced integration.

Right, you may not be able to open those files in older versions.

Yeah, I’ve noticed the name of the file “szklane”. I’m Polish too - if you want to talk in this language, we can move to private chat (on the public forum we should keep using English).

1 Like

I downloaded version 2.4. You’ve fixed the hole, so there should be free support in the rotula. I know how to do this and I use support, but that’s not the point. Interesting. Are you using a “Quasi-structured quad” without Rercombine and have squares in the slab?

Thanks for your time :-). I have to remember that RIGID BODY is for solids, because it mixes plates with shells differently,

I fixed the edge/surface of the hole and, as you said, thought you could try BCs more appropriate for your desired setup.

Yeah, that’s the only 2D algo not requiring recombine. Others need it. Also, it provides the best structured quad meshes in most cases, but it can’t be used with extruded meshes (you would need a sweep mesh, but it has some additional limitations) to generate hexahedral solid elements.

Unfortunately, it’s a significant limitation of CalculiX, and it’s enforced (even for linear analyses) in its new 2.23 version: Rigid body cannot be used with shells at all in ccx 2.23 · Issue #134 · Dhondtguido/CalculiX · GitHub

1 Like

Also, for verification of linear cases with analytical results, you can take the portion of the plate between the holes and treat it as a plate fixed at corners: Plate Supported On All Corners With Uniform Load

This yields the deflection of around 7.2 mm in your case.

I use the model to learn how to model real-world structures. The model must be accurate, and I must compare it with commercial engineering software.