×
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

Peak stress in solid FE models
4

Peak stress in solid FE models

Peak stress in solid FE models

(OP)
Hi all,

Does anyone have any guidelines or philosophies with respect to high peak stresses?

There are several guidelines to for design stress and maximum allowable stress for classic "manual" calculation methods, but I can't find anything for FEM/FEA maximum. How to deal with peak stresses?

Both as an engineer and as a checker I regularly find statements as

* "The high peak stresses should be ignored"
* "The maximum stress is local therefore allowable"
* "The hot-spots are negligible"
* "Due to the local nature of the high stress etc..."
* "Next to the red region there is a green region" (!!!)

As an engineer, this is not very satisfactory. Anybody else finds these situations? Is there a recognised guideline to handle peak stresses? How to evaluate them?

Hoping for a healthy discussion!

Regards,

RE: Peak stress in solid FE models

for myself, I'm ok with "ignoring" localised stress peaks ... tho' i don't think i'd use such wording in a report.  i think i'd say that the peaks are acceptable as they are highly localised, and a results tof the linear elastic assumptions of the model.  In reality such stress peaks would not occur due to localised yielding.  these days i should really use NL codes but ...

RE: Peak stress in solid FE models

There are various industry guidelines for handling highly localized stress.  A few of them are:

>10% stress gradient across a single element (may need to refine local mesh and rerun analysis)

Peak stress in area of little or no interest (engineering judgement tells you that this isn't the load path, so poor FEA geometry is causing the problem).

peak stress across <10% of the local cross-section

Check with your particular industry.  There should be a reason why the statements that you mentioned were made...

RE: Peak stress in solid FE models

I usually find these occur because of poor details or poor modelling, and generally try to amend details to reduce these "hotpots".
If they persist a plastic analysis is usually very quick to run, and values of plastic strain will often be found to be very low.

eelco71...I agree with you and dont like the statements described, which are often made to avoid proper consideration of results.

RE: Peak stress in solid FE models

The statements:

* "The high peak stresses should be ignored"
* "The maximum stress is local therefore allowable"
* "The hot-spots are negligible"
* "Due to the local nature of the high stress etc..."
* "Next to the red region there is a green region" (!!!)

without a detailed explanation show that the FE analyst has no clue what they are doing.  And yes, these are all too common.

Evaluation of peak stresses depends on:
- material (not all exhibit plasticity)
- structure
- stress state at peak (tension, compression, combined, etc)
- loading type (static, cyclic, thermal, etc)
- application
- applicable codes and regulations

So there is no universal guide.  Please provide more detail of the specific situation.

RE: Peak stress in solid FE models

Peak stresses aren't ignored, unless it's a particular aspect of your modelling that you're aware of such as a concentrated force or restraint in an area that you're not interested in.
Usually you consider the peak stresses for fatigue damage but if it's purely a static loading with no cyclic loading then you can ignore it. Refer to design standards such as those given in pressure vessel design codes, usually in the Design by Analysis sections, which give details of stress classification and their appropriate limits.   

corus

RE: Peak stress in solid FE models

eelco-

If you're looking for a code or standard which provides some acceptance criteria for dealing with peak stresses in metal structures, I'd suggest you meander on over to your local pressure vessel engineering outfit (or tech library). Take a look at ASME Boiler and Pressure Vessel Code Section VIII Div. 2 Appendix 4 (pre 2007 ed). The 2007 ed (part 5)  and its companion, API 579-1 / ASME FFS-1 (Annex B1) deal with nonlinear analyses better but are first edition with the customary errors associated withfirst editions. For what its worth, the 1968-2006 VIII-2 is more of an ASD philosophy while the 2007 version has jumped to a LRFD type approach.

jt

RE: Peak stress in solid FE models

(OP)
Few reactions to the healthy discussion above and a wrap-up:

The statements mentioned are just as meaningful as just shouting that the analyst has no clue of what he is doing. It would be good if there could be a proper objective unambiguous way to justify these stresses.

In my experience there are several types of peak stresses:

* due to modelling in regions that you are truly not interested in. Those stresses can be ignored.
* due to modelling in regions where you want to know the stress. To get accurate stresses you should change your model by changing element type (eg quadratic volume elements) or changing mesh size
* truly occurring stresses, where changing the mesh size does not significantly change the result. Element results and average results are more or less similar.

I hope to constrain this discussion to the latter cases.

Classic allowable stress according to most codes is of order magnitude 60%-75% yield. Although I understand how to input plasticity and how to interpret the results, it does not solve the problem as plasticity only occurs around yield stress.

In case it wasn’t clear, material in my case is usually steel, could be S355 or A572 gr 50. Plastic properties are obtained using Ramberg Osgood. Application is not a pressure vessel.

To me it is clear that in some cases “high peak” stresses in “small local” regions cause no problem because of redistribution of stress due to plasticity. At the same time I realise that due to the undefined nature of “high peak” and “small local” this becomes very very slippery.

Anyone has any idea on how to define maximum values of “high peak” and “small local”?

Thanks!

RE: Peak stress in solid FE models

Hi eelco71,

I usually use Neuber's Rule to transform the peak stresses of a Linear FEA. In other words, I return those points back to the real strain-stress curve by Neuber's Rule.

You can get more details by checking the following thread:
Linear Finite Element Analysis & Neuber’s Rule
thread727-210076: Linear Finite Element Analysis & Neuber’s Rule

But, if I have enough time then I will go with Nonlinear FEA.

A.A.Y.

RE: Peak stress in solid FE models

If we have expereince with the product and know the failure modes, some of the stresses could be neglected, with proper explanation. However if I have any doubt I go for non-linear analysis using abaqus. Run time is longer, but the if the model size is not large it is not a problem.

Before doing the non-linear analysis, I would debug the linear model.

Gurmeet

RE: Peak stress in solid FE models

(OP)
feajob, gurmeet,

I have mastered the non-linear plasticity, it takes a bit longer to run, but that's OK.

If the peak stress in linear calculation was higher than yield, it will now be smaller but still higher than yield. Not surprising.

Now what? Is higher than yield now OK? Are the hot-spots now negligible because next to the red region there is a green region?

The statement by GBor is intresting peak stress less than 10% of the local cross section is a good definition of "local". Is there any reference?

What we need now is a good definition of "peak stress"...

RE: Peak stress in solid FE models

2
Hi,
eelco71, imho you should properly understand a pair of points:
1- why you run a plastic analysis in cases you find higher-than-elastic-yield stresses
2- what is the real meaning of a "local", a "secondary" and a "peak" stress.

Let's begin from the simplest of the two points:
2- even if you don't deal with pressure vessels, ASME B.P.V.C. VIII-2 App.4 / 5 (edition 2005 - I don't have the 2007 edition which, as it can be seen in another thread, has considerably changed) give very good indications / explanations on what and where you can consider local, secondary, peak. By the way: for an elastic linear analysis, ASME criteria allow for limits which are considerably higher than yield for "local primary + secondary stresses" (limit "Sps"), but in this case of course you are obliged to perform a fatigue analysis any time it is needed.

1- As for the plastic analysis... except if you already have exact knowledge of a case-figure exactly similar to the one you are analyzing, seeing a load redistribution will give you very little information "per se". First of all, if you really want to go in the "inelastic route", then congratulations, it's certainly the very best way to go, but it will be VERY time-consuming. In fact, you will have to check not only for plastic collapse (using a purely-plastic constitutive law), i.e. check that the redistributed stress can be held by the net section without "frankly" collapsing, but also for PROGRESSIVE collapse, i.e. if you have a variable load you must check that successive load / unload cycles won't make the plasticized zone grow indefinitely until the section will no longer be able to sustain the load itself. The check is now "in control of strain" instead of "in control of stress", since you will have to check that no point of the part has a total strain higher than a given limit (it depends on the Norm you use).
These things are well defined in the first code which explicitely allowed for the elasto-plastic "direct" analysis without stress categorization, i.e. EN-13445. As far as I know, the 2007 edition of ASME updates its methods to almost the same level of completeness.
Warning: depending on the field you work in, you may have VERY BIG problems to have your customers accept the "inelastic route"...

Regards

RE: Peak stress in solid FE models

(OP)
Hi CBRN, Thank you for your helpfull post. You hit the core of the problem. I studied the ASME code you described, there it mentioned the acceptance of "peak stress" in "local" regions due to model discontinuities.

I could not find a numerical definition of "peak" and "local". The plastic calculation I perform is just to try to get a better grip on what is happening, I can switch it on or off in a blink. Wether the calculation is elastic or plastic, the evaluation still comes to the same thing. If higher than allowable stresses are acceptable over a small region of the model, then where is the limit?

The last (plastic) model I had had a maximum stress of 101% yield, with approximately 4% of the load bearing area above the allowable 75% yield. So yes, it is probably OK, but now I need to find a renown document to prove it...

At this point I cannot even convince my colleagues, let alone clients and - in case something happens - insurance.

Note that most of our items are for one temporary use only, no fatigue needs to be concidered.

Thanks,

RE: Peak stress in solid FE models

eelco71,

Your questions exactly echo my concerns about FEA.  i.e. When is a hot-spot OK?

Maybe the answer is to do some hand calcs for certain regions, to convince yourself or your clients that things are probably OK.

I've heard some very experienced FE people dismiss many peak stresses when they were near a load or other stiff boundary condition.  Usually, these people have extensive testing experience to corroborate this.

One thing I have also seen, was (very qualified) people extracting nodal forces from the FE then performing a hand calc (such as My/I) to determine local stresses.

More often than not, in my industry, the controlling failure mechanism is buckling, which is not even captured by linear static FEA, thereby requiring still more hand calcs.

tg

RE: Peak stress in solid FE models

Hi,
ASME code has a way to guide you in order to reckon if a stress can be treated as "local": the main dimension of the area where the stress is 10% higher than the "general" limit should not extend more than SQRT(R*t) and two successive "local" regions can not be closer than 2.5*SQRT(R*t). It's in ASME VIII-2 App.4, 4.112 (i).

Regards

RE: Peak stress in solid FE models

eelco71

In my opinion here is how the non-linear analysis helps. If you see high local stresses more than yield. After the non-linar analysis one has a more realistic picture of stresses and also the extent of high stress zone. In your question you mentioned that you are talking about static loads and fatigue is not involved.

Once we rule out fatigue, the failure mode is yielding through the section. ASME code is a very good document if you are involved with pressure vessel design. However if you are designing products other than pressure vessels, which also look very different, I do not think it is necesary to go back to ASME code after the non-linear analysis.

Stress concentrations are sudden changes in geometry, such as corners, fillets, holes or steps etc. If the gross stresses are much below yield, the high stress zone would be confinded to a small neighbourhood near the stress concentration. The word 'small' is compared to the geometry of the part; thickness, diamter etc. (a 10%) figure sounds pretty reasonable.

If there is still some doubt remaining, one could calculate the limit load for your part. I usually consider that load as the limit load, which causes yielding through the thickness.

I would avoid putting one stress concentration on the top of the other; a hole on fillet.

Gurmeet  

RE: Peak stress in solid FE models

A paper by R.Seshadri (http://cat.inist.fr/?aModele=afficheN&cpsidt=16873436) gives a method for determining the acceptability of hot spots by caclulating remaining strength factors. Although the method is for hot 'hot spots', ie. those caused by temperature, the method assesses not only the stress but also the extent of the 'local' hot spot and its acceptibility. This could be extended to hot spots as defined by geometric stress concentrations. This would then give you an assessment based upon not only the value of stress, but also the size of the region. It's a 'Level 2' assessment, ie. a hand calculation method, but it could probably be extended to results from a linear FE model where the full non-linear elastic-plastic 'Level 3' capabilities weren't available. There's probably a thesis there for anyone wishing to take it up.

corus

RE: Peak stress in solid FE models

(OP)
CBRN, thanks again,

if you do the math with max dimension of stressed region of sqrt(Rt) (and assume it is circular) and the cross sectional area of pipes it turns out that the region could be maximum 1/8 = 12.5% of total area. Intresting, now this could be applied to any area.

Few more issues with this definition:

1:
the code states that the region is a peak stress only when it is 110% of allowable stress. So how about stresses between 100% and 110%?? This stress is higher than allowable, but it is not classified as peak stress.

2:
Is there any limit to the height of the peak stress? even if the region is smaller than the sqrt(Rt)=12.5% I would expect that there would also be a maximum peak stress.

3:
How about when the area is loaded with bending moments as well as axial forces? Peak stresses (due to geometrical discontinuities) will occur at the outer fibres (My/I, but then including stress concentrations)

Calculating the limit load as Gurmeet suggested could be intresting. Just investigating when the entire section is yielding.


I know, this is getting complicated, that is why I feel FEM together with correct interpretation should provide a solid answer.

Thanks all,

RE: Peak stress in solid FE models

Hi,
1- there is a little confusion. In the chapter I mentioned, ASME code is talking about LOCAL, not peak stresses. It means that, when you see a region where the stress intensity is higher than the "primary general stress" limit, you can compare it to a higher (1.5 times higher if no secondary effects are involved) limit, provided that you verify that this "high-stress" zone effectively has the characteristics to be considered "local", i.e. the "SQRT(R*t) for over-110%" criterion. If it doesn't respect that, i.e. if it extends more, then for safety's sake you must consider it "general" stress (not local any more) and so consider the verification failed. Consider the hot-spot stress you have, and "draw" (physically or ideally) two markers at 0.5*SQRT(R*t) from the hot-spot. Now "draw" two other markers where the stress begins to be higher than 110% of the "primary general" limit. Are these two last markers inside the first ones? If yes, the region IS a local stress; if not, bad luck, you have to consider it failed for the primary general stress criterion. IF it is LOCAL, THEN the hot-spot stress must be lower than the "primary local" stress limit, i.e. 1.5 times the primary general.
If you are in a zone where secondary effects arise, then the limit is higher again. Secondary stresses are self-limiting, i.e. they will tend to be "confined" by what happens around. In other terms, they will never be the direct cause of a frank failure of the section where they happen. In order to know that, ASME gives some examples. Other ways to reckon a secondary stress zone is... experience. Another way is the plastic analysis, especially in your case where, as far as I understand, the load will be applied only ONCE and then never removed-and-reapplied again.
So, in YOUR SPECIFIC case, most probably the plastic check against direct collapse is a sufficient verification. You will have to input a "pure-plastic" constitutive law if you want to be safe, as stated by EN-13445.
2- No, not with purely-linear-elastic constitutive laws. You will have to ensure that you can not initiate a damage (crack propagation theory) if you have a traction stress. However, if the peak stress gets too high, most probably you will see a "more-basic" criterion fail BEFORE (primary general, or primary local. If the hot-spot is extremely high (even higher than ultimate stress) but ALL other criteria are respected in the same time, then 99% probably the hot-spot is not "true", i.e. it is a numerical artifact.
Due to the fact you can do it with little expense, in this case a pure-plastic analysis is a good verification.
3- I don't see a conceptual difference in this case.

I'd add: if you run a plastic analysis, the determination of the limit load is easy: it's the lowest load for which the solution doesn't converge any more (be careful to be somewhat restrictive in the convergence criteria, or it will take very long !!!)

Regards

RE: Peak stress in solid FE models

Hello,

I have experienced a similar problem recently.  A static stress assessment of a steel structure calculated a stress concentration where several members are joined together.  I originally modelled the structure with shells, but had to build a solid submodel to find the results in more detail. The results are shown in the attachment below.

As can be seen, there is a stress concentration at the fillet radius.  The yield is 355MPa, and the results indicate the maximum stress concentration is ~550MPa.  I do not believe that the entire section will yield, due to the support offered by the surrounding structure.  The problem is, how much support is offered by the surrounding structure, and therefore when does the stress concentration turn unacceptable?



I also have past experience with ASME III DBA, and have performed non-linear plastic analysis.  A limit load assessment using linear-perfectly plastic properties can be used to gradually increase the load until convergence does not occur.

A full non-linear plastic shakedown analysis of a component requires that the number of loadings above 1.5Sm are known.  The model is subject to repetitive cycles of the intended loadings, and the strain is calcualated.  This is then compared against the ASME limit below:

According to ASME III, Subsection NB-3228.4:

The limits on thermal stress ratchet in shell (NB3222.5) and progressive distortion of non-integral connections (NB-3227.3) need not be satisfied at a specific location, if, at the location, the procedures of (a) through (c) below are used.

(a)    In evaluating stresses for comparison with the remaining stress limits, the stresses shall be calculated on an elastic basis.
(b)    In lieu of satisfying the specific requirements of NB-3221.2, NB-3222.2, NB 3222.5 and NB-3227.3 at a specific location, the structural action shall be calculated on a plastic basis, and the design shall be considered to be acceptable if shakedown occurs (as opposed to continuing deformation).  However, this shakedown requirement need not be satisfied for materials having a minimum specified yield strength to specified minimum ultimate strength ratio of less than 0.70 provided the maximum accumulated local strain at any point, as a result of cyclic operation to which plastic analysis is applied, does not exceed 5.0%.  In all cases, the deformation which occurs shall not exceed the specified limits.
(c)    In evaluating stresses for comparison with fatigue allowables, the numerically maximum principal total strain range shall be multiplied by one-half the modulus of elasticity of the material (Section II, Part D, Subpart 2, Tables TM) at the mean value of the temperature of the cycle.


I have long sought where the 5% strain value comes from, and also as to what type of strain it refers.  I believe it to be a 'membrane' strain though the section, as this section guards against progressive plastic collapse.




Incidentally, you can also perform a fatigue assessment using the results of the non-linear assessment, to remove conservatism.


I don't want to perform a full non-linear shakedown analysis for this static asssessment since it takes a long time/would have difficulty defining the inteded loads over the lifetime.  I'm also not sure if I should mix rules from a pressure vessel code with a strutural design code.

What would be useful would be some guideline for when a stress concentration in a load bearing memeber becomes a mode of failure.  Any suggestions?

RE: Peak stress in solid FE models

I just want to correct something that chaboche wrote in his previous post.  He wrote that you can perform the fatigue analysis using the results of the non-linear analysis.  If chaboche is referring the using the strains from said analysis, and comparing them to the ASME fatigue charts, then that would be WRONG!!!

In TF-EPFEA and SG-DA, we have been slowly working our way through how to do fatigue with non-linear analysis, and what I can tell you is that that is NOT the way to do things.  The fatigue charts are based on elastically calculated stresses (hence the use of a K_e factor to adjust for plasticity), and comparing them with a non-linear analysis would simply not correlate.

RE: Peak stress in solid FE models

What is surprising in the last picture is that the results from the sub model don't appear to relate to the global model, in particular of the location of the maximum stress. I'd expect the results at the limits of the sub model to be the same (or close to) the results of the global model. I'd check that the displacements from the global model have been correctly transformed to the sub model.

You ask when a stress concentration becomes a mode of failure. I presume you mean by regard to some yield criteria and not by fatigue damage which the peak stresses at a concentration are assessed against. By defining the stress as that occuring at a concentration (or a feature) then you have already classified the stress and hence its limits. In terms of stresses away from the maximum stress at the concentration then you can look at the stress distribution up to and including the peak stress and attempt to fit a straight line through the stresses that will remove the non-linear contribution to the curve from the concentration. It is largely based upon your own judgement and so it would be better to err on the conservative side in your assessment. The resulting stress would be classed as secondary at the structural discontinuity. In the shell model this stress concentration of the fillet will be excluded and hence it would be better to assess the stresses as either secondary or primary based on their location. If in doubt then be conservative in your assessment and choose the lower limit.

corus

RE: Peak stress in solid FE models

I am interested to learn about the work TSG4 has done on non-linear fatigue.

Traditionally, the stress range during a cycle of a thermal transient can be calculated and compared against S-N data to give an allowable number of cycles.  If pasticity effects are included in the FE model, why cannot the strain range be compared against the S-N data instead?  This would remove the need for a conservative plasticity correction factor.


In response to corus, thanks for your comments.  The reason for the difference in location of maximum stress between the global (shells) and submodel (solids) is due to the cut boundary interpolation.  The maximum stress has shifted to the submodel boundary, where several of the shells are joined together.  This has created a highly localised concentration.  I have checked and the overall displacements between the models are in agreeement.  Another contributing factor is that the shells weren't defined exactly through the centre of the solids i.e. some shells are located at the top/bottom of the solids.  This was unfortunately decided before it was evident that a submodel would be required.

Regarding this location, you are correct that I wish to determine the strength of the section against a yield criterion, as opposed to the fatigue damage.  I can see why the stress would be defined as a secondary stress according to a pressure vessel design code.  Most structural steelwork codes however don't seem to classify stresses in the same way though.  AISC or DNV, for instance just give limits on the calculated stresses without classifying them as primary or secondary (or maybe I'm looking in the wrong sections?).  Is it therefore still valid to draw a line through the results to remove the peak?

A good way to show the stresses through the section is to plot those which are acceptable and those which are not.  This is shown in the figure at the link below.  It is then to be judged as to whether the concentration extends significantly enough to cause yielding of the entire section.







RE: Peak stress in solid FE models

The problem with structural codes, as I see it, is that they're based upon hand calculations of simple direct and bending stresses and don't consider the extent to which finite element methods can predict stresses. In fact I don't think they consider thermal stresses from what I remember. The pressure vessel codes do give more thorough ways of assessing stresses which the structural codes I've seen don't seem to possess.

Another aspect of peak stresses in transient thermal analyses to consider is the non-linear stress distribution through the thickness, which your 3D solid model will pick up. Your stress distribution might therefore be composed of this thermal stress component as well as the stress concentration from a geometric feature. You may need to calculate the equivalent linear bending stresses through the thickness of each section, and thus remove that peak component to leave only the secondary stress (range) and peak stress from the geometric feature. Effectively you'd have to extrapolate your results in two directions to remove any peak stress components.

In your case it's clear that the stresses at the fillet are due to the stress concentration of the geometry (plus perhaps a thermal peak stress). Stresses at 1.0t from the weld toe should be considered as nominal stresses and assessed accordingly.

corus

RE: Peak stress in solid FE models

I would like to make one more comment on calculation of limit loads. This may be obvious to some. One knows the target factor of safety. In non-linear zone the factor of safety should be based on load and not stress. For example let us say that the factor of safety is 3. One could run another analysis with 3 times the load and look at the results. If the part has not yielded through the section, one may consider it safe. On the other hand if the run does not converge, it may indicate the part does not have the adequate factor of safety.

Gurmeet

RE: Peak stress in solid FE models

Hi,
yes, and in addition to what Gurmeet says, I'd add (just because I've seen it asked sometimes in Eng-tips threads...): when you perform a plastic analysis against direct collapse, you can not consider stresses in the plasticized zone: of course, they will be clipped to the yield value of the material (slightly more due to numerical reasons). You will have to:
- either compare the strains to related limits
- either perform a "limit-load" verification as suggested.

The important thing is that, in most cases, the direct collapse verification is subjected to strain limits by codes like EN-13445, not to limit-load safety factors (though these can be far more "understandable").

Regards

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