I solved a plate problem by Maple (similar to Matlab) program using Classical Plate Theory (CPT) and First Shear Deformable Theory (SDT). I attached the picture.
After I get the results from my Maple codes, I want to compare the results with FE results so I set up FE model in FEMAP and run the analysis with NX Nastran.
The FE results are somewhat close to my Maple results, but they are not identical. For example, for the deflection of a node, Maple gives me 0.68'' whereas the FE result gives 0.715''. The main question that I want to ask is why is there such a different in result assuming my Maple program is good (I double check with my professor). If anyone has compared the FE result with the theoretical results before, I would like to hear from your perspective of why you think there is a difference in my case.
I did some research and know that NX Nastran by default uses CQUAD4 element for plate which is based on Mindline-Reissner theory with transverse flexibility. And this simulation is linear static (SOL 101). I have tried the same problem with different load configurations, but there is always a difference.
I would appreciate if someone can comment on this.
Solved! Go to Solution.
The Finite Element Method is an approximate method, the accuracy of the solution depends of various factors: element type, mesh density, loadings, boundary conditions, analysis type, etc.. Here you are a few suggestions to improve accuracy:
In summary, using Finite Element Method you can have MANY solutions, all are solutions, yes, the valid one will depending of what you are looking for: stress, displacements, reaction force, buckling load factor, etc.. You can "play" with many factors, and all have a "price" to pay (solution time, model size, model creation time, etc..), enjoy!.
Thank you blas for your suggestions. I will try some of them out.
Right now, I am gonna use Patran/ Nastran and possibly other softwares to compare the results.
There is another question I want to ask. CQUAD4 element is based on Mindlin-Reissner plate theory with transverse shear flexibility which means that the users can turn on or off the transverse shear. In Femap, the transverse shear can be turned off in plate property. I want to know how does "NX Nastran" get rid of transverse shear for CQUAD4 elements. If the transverse shear is taken out, the Mindlin-Reissner shouldn't be valid anymore. Then what would CQUAD4 element become?
I'm trying to look into NX Nastran theoretical manual without much luck on CQUAD4 derivation. Let me know what you think.
The CQUAD4 element of NX NASTRAN is my favorite between CTRIA6, CQUAD8, CQUAD4R, etc.., is an isoparametric quadrilateral element with optional coupling of bending and membrane stiffnesses. The formulation of the CQUAD4 element is based on the Mindlin-Reissner shell theory. The element do not provide direct elastic stiffness for the rotational degrees-of-freedom which are normal to the surface of the element.
Depending on the 2-D element, you define the properties for plate and shell elements with either PSHELL, PCOMP, PCOMPG, or PLPLANE bulk entries. Transverse shear flexibility may be included for all bending elements on an optional basis. Differential stiffness matrices are generated for all shell elements except CQUADR and CTRIAR.
The PSHELL entry allows to define the material ID for the membrane properties, the bending properties, the transverse shear properties, the bending-membrane coupling properties, and the bending and transverse shear parameters. By choosing the appropriate materials and parameters, virtually any plate configuration may be obtained.
The CQUAD4 element can model in-plane, bending, and transverse shear behavior. The element’s behavior is controlled by the presence or absence of a material ID number in the appropriate field(s) on the PSHELL entry. To add transverse shear flexibility to bending, fill in MID2 and MID3.
Adding transverse shear flexibility means that using MID3 adds a shear term in the element’s stiffness formulation. Therefore, a plate element with an MID3 entry will deflect more (if transverse shear is present) than an element without an MID3 entry. For very thin plates, this shear term adds very little to the deflection result. For thicker plates, the contribution of transverse shear to deflection becomes more pronounced, just as it does with short, deep beams.
You should use the CQUAD4 element when the surfaces you are meshing are reasonably flat and the geometry is nearly rectangular. For these conditions, the quadrilateral elements give more accurate results for the same mesh size than any other element type, OK?.