Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
CQUAD4 in bending test
#1
Hi, 

I am getting some unexpected results running the attached cantilevered composite plate model in Mystran 11.0. :
  • The model is a flat plate in the X-Y plane, with dimensions of 1.5m x 0.15m and fully fixed at one of the short ends. 
  • A Z-direction shear force is applied through an RBE3 at the unsupported end (Grid 737).
  • I would expect a Z-deflection at Grid 737 close to T3 = 0.29625 (simple CLT calculation and in agreement with other FEM solvers)
  • But I am getting T3 = 0.03613559 at Grid 737 instead.
  • AUTO SPC forces are zero (constraining Z-rotations / drilling DOF on the CQUADs only)
Any ideas what could trigger this behaviour?

Do you have a simple test case for shell (CQUAD4) element bending that I can run to test the installation? 
Thanks.


Attached Files
.txt   shell_bending.txt (Size: 122.08 KB / Downloads: 31)
Reply
#2
Can you check the reaction forces to make sure you are applying the correct load? There is a 0.1 scaling factor in the FORCE card so you may want to double check that. Have you tried an isotropic material to see if the results are correct for that first? I will let Bill know as well.
Reply
#3
I'll take a look at it and get back to you
Reply
#4
(08-27-2020, 06:26 PM)Admin Wrote: Can you check the reaction forces to make sure you are applying the correct load? There is a 0.1 scaling factor in the FORCE card so you may want to double check that. Have you tried an isotropic material to see if the results are correct for that first? I will let Bill know as well.

Reaction forces are correct - factor of 0.1 is intended. 
The element node numbering is such that the vectors from grids 1 to 2 point in the x-direction for all elements - this is matched in the comparison models - so I don't think this is an input error? Besides, the results are not any closer if the plies are set to 90deg in the PCOMP cards. 

As you suggested, I reran the model using an isotropic material - results are much closer (<2% error).
Comparing Zmax in the T3 deflection at grid 737:

Analytical (Bernoulli beam):
I = (1.25e-04*12)**3*0.15/12 = 4.21875e-11
F = 0.1
L = 1.5
E = 70.0E9
Zmax  = F*L**3/(3*E*I) = 0.038095

Mystran (see file attached)
Zmax = 3.762719E-02


Attached Files
.txt   isotropic_shell_bending_1.txt (Size: 122.1 KB / Downloads: 14)
Reply
#5
I just did a quick calculation treating the plate as a beam (should be good assumption since the aspect ratio is 20 to 1) and get a displacement at the end that is within 7% of the MYSTRAN plate problem result. I did run the problem in MYSTRAN and get exactly what you got, the displacement of 0.03613558. See the attached jpg which shows the beam displacement within 7% of the plate solution. Can you tell me what other sources gave you such a different answer?


Attached Files Thumbnail(s)
   
Reply
#6
O_Stodieck, can you show your hand calculations for the original problem as well? The composite version.
Reply
#7
Sorry, I used the wrong width (.075 rather than .15). But where did the 70E+09 for modulus come from?
Reply
#8
(08-28-2020, 12:05 AM)drbillc Wrote: Sorry, I used the wrong width (.075 rather than .15). But where did the 70E+09 for modulus come from?

Bill, there are two decks in the different posts and it was changed between them. Its a bit confusing, but maybe it was an attempt to try to make the composite more "isotropic-like". Not sure.
Reply
#9
O_Stodieck, I was looking for a good QUAD4 example to send you and, in the process, uncovered an error in MYSTRAN that crept in between versions 10.1 and 11. It is a plate problem with a pressure load that has a theoretical solution so that you can see how well MYSTRAN gets that answer. I will fix the error and then send you the bdf, but in the meantime MYSTRAN may crash with a shell pressure load.
Reply
#10
(08-28-2020, 12:02 AM)Admin Wrote: O_Stodieck, can you show your hand calculations for the original problem as well? The composite version.

Sure - I have posted a python implementation of the CLT approach here https://gist.github.com/ostodieck/d0411f...2db4ccf512
I hope this is self-explanatory enough - let me know if something is unclear. 
The laminate and materials are the same as in the first post, and the tip deflection is  0.29625 m as expected. 
Thanks for checking this!
Reply


Forum Jump:


Users browsing this thread: 1 Guest(s)