×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

(OP)
Dear all,

I have a peculiar case in which the exact opposite happens from what is expected. I am modelling crack propagation using cohesive elements. In the beginning I am just modelling a mixed-mode bending test in plane strain to calibrate my model parameters. The model is bending dominated, much like a cantilever with tip-load. Trying different mesh sizes and elements for the bulk material I found the following:

1. Quadratic plane strain elements (CPE8) deliver the most accurate solution. (Closest to analytic solution)
2. Linear elements (CPE4) are softer than quadratic elements
3. With increasing mesh density the solution with linear elements approaches the quadratic solution, thus the model becomes stiffer.

This is exactly contrary to what I know, that linear elements are softer in bending than quadratic ones and that the response stiffens with mesh refinement.

Can anyone explain this?

RE: Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

There is nothing wrong in what you've found; your prior knowledge, it turns out, is inaccurate.

When a mathematical model (in this case, an FE model) is determined to be "stiff", it means there are not enough degrees of freedom available for the variation in the field variable (such as stress). So, when you increase mesh density or increase the order of the elements, you provide the model with more degrees of freedom making it "less stiff". This "stiffness" is different from the typical definition of stiffness being equal to delta_load/delta_displacement in the linear region of the experimental data from, say, a uniaxial tensile test.

http://www.eng-tips.com/faqs.cfm?fid=376
http://www.eng-tips.com/faqs.cfm?fid=1083

RE: Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

(OP)
Dear IceBreaker,

Its the other way around: More elements --> Stiffer

RE: Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

If you are disagreeing with IceBreakerSours point then I'm sorry but I have to agree with IceBreakerSours. Below is a quote from an article in MachineDesign (http://machinedesign.com/article/a-closer-look-at-...). I don't have any of my texts to copy a page that describes, more or less, the same thing.

If you are just reiterating what you are seeing in your model then I guess just ignore this response.

QUOTE: "Consider a cantilever beam meshed with three different densities to see how the maximum deflection changes. To make this comparison fair, even the first mesh (the rough one) must be able to model bending stiffness.

The beam deflects further with finer meshing. This proves that the finite-element model becomes "softer" with a rise in element density. The effect arises because artificial constraints become less imposing with mesh refinement. ..."

Dan

Han primo incensus

RE: Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

(OP)
Dear DanStro,

I do not disagree with IceBreaker, neither do I disagree with you. What both of you wrote is also what I have experienced so far.

But you both have misread my post, or I have not clearly expressed myself. My model behaves the other way around from what you describe:

- Common knowledge: Finer mesh --> softer response

- My model: Finer mesh --> stiffer response

I do not understand why. That is my question.

RE: Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

(OP)
@IceBreaker,

I am looking at plots of applied force vs tip displacement. The more elements I use through the thickness of the specimen, the more load is measured for the same applied displacement, thus the response is stiffer (steeper plot).

Thank you btw.

RE: Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

(OP)
I have attached a plot to show what I mean. t stands for number of elements through the thickness, all else being equal.

http://files.engineering.com/getfile.aspx?folder=f...

RE: Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

a) Firstly, check the .msg and .dat files when t=1. My guess is that there may be some warning messages. Even if they are not, output at t=1 is incorrect because one linear element can only allow constant strain through it. In a problem such as yours, there is a variation in strain going through the thickness. Bottomline: t=1 is not even a correct solution. This is a classic "stiff" problem but the solution (from the curve) is, if I use your terminology, compliant/soft.

b) Secondly, as I said before, this may be just a matter of terminology used by mathematicians (FEA developers) vs. engineers. Stiffness of the structure and a mathematical problem/model being "stiff" are not the same thing. You are converging to the analytical solution as you increase mesh density, which is what is expected. So, if one is observing your load-displacement curves, then yes, the structure is becoming stiffer as you increase mesh density but that is not the same thing as the problem being "stiff" or not.

http://www.eng-tips.com/faqs.cfm?fid=376
http://www.eng-tips.com/faqs.cfm?fid=1083

RE: Linear elements stiffer in bending than quadratic ones (not the other way around as it should be)

(OP)
dear IceBreaker, thx for your answers

then I am confused as to the terminologies of stiff in FEA.

From the document that DanStro posted I quote: "The beam deflects further with finer meshing"

In my weird model the opposite holds true: My model deflects less with finer meshing. Definition of stiffness be as it may, but that fact remains that my model behaves the contrary to common believe.

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources