×
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

Checking the FEA quality
5

Checking the FEA quality

Checking the FEA quality

(OP)
Hi every body,

I would like to know how do you check the quality of your FEA for static linear stress analysis?

I know that an error estimation method was proposed by Zienkiewichz on 1987 and many other publications came later for this purpose. I am not sure if there is any adaptive mesh refinement strategy available based on this norm (in any commercial FEA code)! I use MSC.Nastran, but I am not sure if this error estimation method is available in Nastran too.

Usually, there are many small geometrical features in my models, so I do some geometrical clean up and then I mesh my model with HyperMesh. The mesh quality is very good (about 1 million degrees of freedoms for my quadratic tet meshes).  My boss recently asked me, how we can be sure about the stress quality in the stress concentration zones? I thought that the best way is checking with error estimation norms. But, I would like to ask you first, may be some one have a better idea or easier way for this purpose. By the way, I found this forum very helpful.  

Regards,

Ali
http://www.geocities.com/fea_tek/asd/

RE: Checking the FEA quality

FEA is a tool and give as an approximation of the stress in the analysed parts.
If you will refine your mesh in the hot spot sectors you will have a different configuration of stress: much small area but higher stress, this thing is not necessarily more truth or much closer from reality.
In HM, your are lucky, you can use QI features. Try to use the Quality Index feature and ask to HM to optimize your mesh with this criterion.
You have some steps to do for a quality check:
*merge nodes with a convenient tolerance (pay attention to the value of tolerance, is very model dependent, start with very small value),
*use “edge” to detect disconnected elements, nodes, or duplicated,
*use “shrink” to detect duplicated (shell) elements,
*warping (shell), Jacobian, aspect ratio and others criteria if you want (I don’t known very well the criteria for the solid element, I work more with shell-plate and beam elements)
Some model imperfections can generate stress concentrations (unconnected nodes).
Hope this help,
gelu

RE: Checking the FEA quality

feajob,

One way to check model quality is to run convergence studies.  This is very time consuming, and as gelu mentions, it does not necessarily substanciate the results, just the quality of the mesh.  gelu provides many good checks to do on the elements themselves before you even run the model, or after you perform some of the checks below.

Another way is to plot the unaveraged results using element fill options.  In this case, you should get a fringe plot pattern that looks very similar to your averaged/continuous plot.  I.E. smooth transition between elements.  If it looks like patchwork, then your model could use more refinement.  Also, check the difference plot (instead of average) and see how much variation there is between adjacent elements.  The smaller the value, the better.  However, it is relative to the ultimate stress in that region.

Finally, you should check a net section against theoretical results just as a sanity check that loads and boundary conditions have been applied correctly.  If you have the ability to strain gauge an actual product at some point, verify strains against your model.  These are the only ways to truely verify the model accuracy.

Hope this is of use.

jetmaker

RE: Checking the FEA quality

feajob is concerned about the validity of the results from a model, and their accuracy. Generally you can check your results by plotting the unaveraged element nodal stresses and comparing the contours between adjoining elements. If they contours are close then it's reasonable to assume that your results will be fairly accurate. If you want to be precise about it then re-run your analysis with a greater mesh density. Preferably do this with sub-models of the region as a million degrees of freedom in a model seems to be a little large. If you can mesh using brick elements, and the analysis run-time is relatively short, then it's easy to run the model with linear elements and then change them to quadratic elements, thereby keeping the same number of elements but increasing the accuracy for your comparison.
For stress concentrations though it's generally better to model the whole part with a coarse mesh and then sub-model the region of interest, but include the stress concentration, and let rip with your automatic mesh generator, while checking for compatibility at the sub model boundaries with the coarse model.
Nafems might give better guidance if you wish to check their site or literature.

corus

RE: Checking the FEA quality

(OP)
Thank you gelu for your response. I am agree with you HyperMesh has some good tools for checking the element quality. I use them, how ever, as you mentioned FEA give us just an approximation of the real solution. In the stress concentration area, if I refine and refine my mesh then I will get closer to the real solution.

May be I didn’t explain well my question. I am looking for a way to do this job more efficiently. As jetmaker mentioned this job may be very time consuming and our boss won’t be happy. Jetmaker solution seems more practical, because as we know FEA provides discontinuities for stress fields across the elements (no matter which code you use). The stress are initially calculated  at all integration points then they are extrapolated into all nodes (using a smoothing + averaging method). Jetmaker, at this stage of project, we don’t have opportunity to measure stains with stain gauges. Comparing the FEA solution with theoretical solution is a good approach, but this solution it is not always  available for a complex 3D model.

Corus, sub modelling is a good approach, but still, it may be time consuming. I use Patran and HyperMesh, but I think that checking the compatibilities  between coarse mesh and fine mesh is a time consuming job.

John Barlow is a guy in Rolls-Royce (stress department) did a good job many years ago. He reported in his paper the Optimal Stress Locations in FEA. He demonstrated that FEA provides higher accuracy on some points of an element. Later, Zienkiewicz and Zhu developed a global and local smoothing methods based on this optimal points. In fact, they built estimated stress field and compared this field with FEA stress field. As a result, they can estimate the element error for whole mesh. MSC.Marc is using this norm for refining mesh of nonlinear analysis.

I though that may be some one is using this norm for evaluation the quality of FEA solution in linear FEA.

Thank you for your responses.
Ali
http://www.geocities.com/fea_tek/asd/

RE: Checking the FEA quality

feajob,

Sorry, not familiar with the approach you describe.  My only understanding is that displacements, reported at the Gaussian points are the most accurate results presented by FEA.  Other results (i.e stress, strain) are derived from those displacements.  Another check has something to do with strain energy.  Something along the lines of checking the model strain energy and seeing if any element is picking up a greater proportion of that energy.  Can't remember the specific details, but I'll get back to you if I remember.

jetmaker

RE: Checking the FEA quality

(OP)
jetmaker,

FEA solves at the first step the following equation:

[k]{u}={f}

where
[k] denotes stiffness matrix of our model
{u} nodal displacements
{f} nodal forces

Evidently, the distributed forces can be transformed into the nodal concentrated force.

So, it finds nodal displacement first. Then, at the second step, it calculate the strain field by derivation of the displacement field at integration points.

Finally, at the third step it calculates the stress field at integration point using law of material (Young's modulus and poisson ratio).

That's why displacement field is a continuous field in FEA, while stress and strain fields are discontinuous.

Ali







RE: Checking the FEA quality

HI all,

i´m generally (like all of you) after stresses at stress concentration zones to estimate fatigue behavior.

My experience is that in fatigue there are so many unknown parameters that the error in linear elastic stress prediction is much smaller than errors from boundary conditions and material data and history of specimen.

My criterias are very simple and maybe not good enough for all fea-users.

-I just mesh my model very coarse (stiffnes should be represented good enough with an coarse mesh) and refine it at locations of interest.
-Then i take a look at the averaged stress distribution at the stress concentratoin zones (should be smooth and the differnce in stress along a element side should be smaller then 10%)

-sometimes i use adaptive meshing to see the convergence (mesh independent result), but i always model the notches in the model with a small radius. Therefore i get convergent results

regards

mfinke

RE: Checking the FEA quality

(OP)
Hi mfinke,

It is an interesting and practical criteria. I guess you apply 10% criteria for 2D problems? Are you doing 3D too?

Thanks,
Ali

RE: Checking the FEA quality

Ali,

The results of your analysis may differ from the real answer for many reasons such as:
accuracy of loads
accuracy of material properties( tested vs. literature)
is the mesh fine enough?

It appears that you are looking for information on the last point mentioned above.

One simple way to check that the mesh is fine enough in area of interest is to either halve or double the mesh density in that area. If your results do not change, it indicates that your mesh is fine enough. However if your results change in that area, it indicates that the mesh needs to be refined further. Please note that when you are focussing on an area of stress concentration you do not have to increase the mesh density in the whole model but only in the area of interest.

Gurmeet

RE: Checking the FEA quality

Simple problem.
Simple solution.
Automatic solution feedback using the P-method.
Mechanica perhaps?
see Forum828

RE: Checking the FEA quality

Hi FeaJob

In my experience on Fatigue studies, when commenting on mesh quality, you are referring to the independence of the results to the mesh size, formulation etc.  

I tend to refer to mesh size, angle, warp etc as 'mathematical checks', most pre-processors tend to perform these checks for free during the mesh creation process, and although its important to get these within certain limits for the solver (nastran can get upset with large aspect ratio limits for example) as loong as these are within sensible limits you will find that variation in these parameters don't affect your results.

Element formulation and density are for me the two more important things to concentrate on.  For fatigue analysis, a fine mesh size <=5mm is a good starting point for panel areas - with concentric circles around holes, of finer density around the hole, then broadening out as you move away from the hole.  Another important point is to always model quads in the area under fatigue study - some solvers can get really upset with tris and highlight them as problems, when common sense can eliminate them.  Increasing element order is obviously another good thing - but with fine mesh and high order elements, you need the computing resource to get results in a reasonable time.

Personally, I tend to prefer to mesh the component in consideration to its loading, although I do use the Auto mesh as often as possible - I do cut up geometry to alter its alignment and direction so it looks 'neat'.

Once you've got a fine quad mesh of reasonable quality in the affected area - and run the job, I prefer to look at the results in conjunction with engineering judgement - when I say engineering judgement I mean :-

Firstly I look at the entire model as a whole. I'm looking for a difference in the contour plot across a group 3 shell elements of less than 10% of the total contour value (this shows me that the distribution is smooth and not affected by the way the component has been modelled.  If its greater than this - I tend to remodel it in the affected area. Re-modelling in this context means altering mesh density and alignment.

Then I start looking at areas of concentration which I can't explain with the above.  In these cases I will do some rudimentary calculations and see if the results are in the same ballpark.  These calculations derive from the k-u.F formula.  If they look out - then I look at re-meshing the affected area - changing the mesh alignment to the loading direction, not necessarily the density or order of the mesh.

Anything that remains after these checks, is usually a concern for fatigue.  I then approach texts from the authors already mentioned or my peers for detailed advice.  Another author of note is R.A. Higgins, has some good points about the fatigue failure mechanism.

Depending on the component under study - you may have to sub-model it as corus suggests.  This can be appropriate if you are modelling in solids - or across rapid changes in section for example.  If you are modelling in bricks (solids) be careful to get at least three bricks across the thickness though - if you don't you're results will be a little suspect.

Hope this helps

sean

RE: Checking the FEA quality

(OP)
Hi,

Sean, Thank you for your reply. Your points are interesting. I can use some of them, but not all of them. Because, I do 3D FEA. Even after many geometrical clean up, my geometry is still very complex, so, I don't any other choice than 3D tetrahedral meshes (quadratic). I usually end up with more than 1,200,000 DOFs. I have good hardware and usually I can solve my linear problems in 45 minutes.

gurmeet2003, I am agree with you, the best way to check the quality of a FE solution is refining mesh in stress concentration area. However, it may not be the best cost effective solution, when you are dealing with large and complex models.

d3ets, I am agree with you too. P-elements are helpful, because you are involve directly with mesh refinement procedure. In fact, promechanica increase the p degree of elements and it increase accuracy without refining your coarse mesh, but, as far as, I remember, there are some limitation. For example, you can not solve complex contact problem with p-elements? Is it correct?

Thank you for your comments,
Ali
http://www.geocities.com/fea_tek/asd/

RE: Checking the FEA quality

You can solve contact problems in Mechanica, however you do need a fine mesh in the contact area. Mechanica has an automatic method for doing this, or you can manually seed the contact area.

RE: Checking the FEA quality

feajob,

back to your original question:

I agree with gelu (check for coincident nodes, boundaries: free edges, duplicate elements, elements distortion: taper, skew, warp …,  degenerated elements, check association of elements and nodes, element normals, connectivity, shade FEM) and with jetmaker (plotting the unaveraged stresses and look for a smooth gradient).

You can also check the difference. The difference is determined as the delta between the maximum and minimum stress contributor to each node. The fringe plot of this max difference indicates the quality of the mesh and the location where this mesh needs to be refined by comparing its values with the actual values of the results.

If you use Nastran open the *.f06 file and check:

- OLOADS and SPCFORCES have to equal out
- model mass and CG (verify against known values), PARAM, GRDPNT,0 - not just for dynamic analysis, I find it good to verify that I modeled the structure realistically
-  EPSILON must be < 10E-9
- MAXRATIO > 10E5 may indicate trouble

The EPSILON is a measure of the numeric conditioning of the model. An acceptable value of EPSILON depends on the model complexity and the machine that it runs on. An EPSILON value of less than 10E-9 is generally considered acceptable, for a clean and converged model EPSILON is around 10E-15. A large value of EPSILON is evidence of numerical ill-conditioning (e.g. poor constraints or mechanism).

MAXRATIO is the maximum ratio of the matrix diagonal to the factor diagonal. A large MAXRATIO is an indication of a potential modeling problem.

RE: Checking the FEA quality

(OP)
AnnHer,

In fact, I use Nastran as a solver. Please let me know how you can get MAXRATIO? Because, I see only EPSILON in my *.f06 file.

BTW, can you explain more about "the delta between the maximum and minimum stress contributor to each node". How  can you make a fringe plot of this delta in Patran?


Thank you for your comments,
A.A.Y
http://www.geocities.com/fea_tek/asd/

RE: Checking the FEA quality

Regarding the issue of convergence, I recommend taking a look at the FAQ P-Elements (FAQ828-810) to, hopefully, clear up some confusion surrounding that topic. The subject of actual results quality is still, unfortunately, a separate issue. Just because a mesh is of good quality, or convergence on strain energies (or other measure) has been achieved, it does not guarantee that the results are accurate. The model could also be prepared poorly (having cracks, or sharp corners) that can represent a singularity; no convergence check will produce accurate answers in this case.
Other sources of error can be in the modeling loads and boundary conditions. One can also have errors in the type of element used (membranes rather than shells, plane strain when plane stress was correct, linear triangles or tetrahedrons instead of parabolic elements, etc.).

As has been mentioned before, a Finite Element Model is just that, a model. There is always potential for errors with any model.

Best regards,

Matthew Ian Loew
"I don't grow up. In me is the small child of my early days" -- M.C. Escher

Please see FAQ731-376 for tips on how to make the best use of Eng-Tips Fora.

RE: Checking the FEA quality

Yes, as we know the output of a finite element analysis is always as good as the input and the input is a lot of guesstimate especially when it comes to boundary conditions, loads, connections or lovely homogeneous material. Results are always to be seen relative. And regarding the issue of convergence, it depends on your needs. One does not need a perfectly converged model for a concept design study. One does also not need a perfectly converged model if e.g. the max. stresses reach  10% of the yield strength.


feajob,

PARAM, MAXRATIO, 1.E7     specified in the Case Control Section

MAXRATIO Default=1.E7 The ratios of terms on the diagonal of the stiffness matrix to the corresponding terms on the diagonal of the triangular factor are computed. If, for any row, this ratio is greater than MAXRATIO, the matrix will be considered to be nearly singular (having mechanisms). If any diagonal terms of the factor are negative, the stiffness matrix is considered implausible (nonpositive definite). The ratios greater than MAXRATIO and less than zero and their associated external grid identities will be printed out. The program will then take appropriate action as directed by the parameter BAILOUT.

By default
PARAM, BAILOUT, 0 and it should stay this way

If you experience FATAL MESSAGE 9050 and you don’t know what went wrong, you may set PARAM, BAILOUT, -1 to force the run through. But this you do just to check your model! Don’t believe your results, since thy might not be accurate.

About the Difference in Patran:

Go to Results: Create/Fringe and select your results. Click on Plot Option and go to Averaging Definition. Domain: All entities, Method: Difference

Check it out.
A

RE: Checking the FEA quality

AnnHer,

I am not sure I would trust comparing results from models that had unknown or different convergences. I think that the issue of convergence is routinely ignored by many analysts. I use MECHANICA when I do structural analysis, so I have this capability built in.

I have found some articles that may be of interest:
  • P. Kurowski Easily made errors mar FEA results, Machine Design, 13 September 2001
  • P. Kurowski More errors that mar FEA results, Machine Design 21 March 2002
  • V. Apanovitch, P. Kurowski Discretization Error Control in Precise Solids Method, part 1, BENCHmark, April 2001
  • V. Apanovitch, P. Kurowski Discretization Error Control in Precise Solids Method, part 2, BENCHmark, July 2001
  • P. Kurowski,  B.Szabo  How to find Errors in the Finite Element Models, Machine Design, 25 September 1997



  • I will (eventually) create a FAQ on this topic so we can maintain this important information.

    Best regards,

    Matthew Ian Loew
    "I don't grow up. In me is the small child of my early days" -- M.C. Escher

    Please see FAQ731-376 for tips on how to make the best use of Eng-Tips Fora.

    RE: Checking the FEA quality

    MLoew,

    How does this capability in MECHANICA work?

    I fully agree, I would not trust comparing results from models that had unknown or different convergences either. But we are still talking linear and convergence in respect to mesh discretization – regardless h- or p-elements and assuming a good mesh quality. Perfectly converged for me means 0% difference in results due to mesh discretization. It is the duty of every analyst to check convergence, but there is often a limit, since we unfortunately don’t always get paid for this.

    Assuming a concept design: An analyst’s model setup is well done, he made good assumptions and did some convergence studies. He ends up with an error of 10% in the results due to convergence, but stresses just reach 10% of the yield strength. Fine.

    Assuming the same, but final design and resulting stresses reaching much closer to the yield strength. Definitely get your result difference due to convergence as close to 0 as possible.

    This again looks different for an analyst doing e.g. a material or geometry parameter study. The analyst has to decide and QUOTE for the convergence requirements.

    Good links.

    A

    RE: Checking the FEA quality

    AnnHer,

    Please take a look at the FAQ P-Elements (FAQ828-810) and a couple of the links above. Briefly, in MECHANICA, the polynomial level of the elements is increased until convergence of a predetermined level (typically 10% on strain energies) is reached for a multi-pass adaptive analysis (MPA). There is also a single pass adaptive (SPA) method that presumes the increase in p-order to achieve convergence.

    Best regards,

    Matthew Ian Loew
    "I don't grow up. In me is the small child of my early days" -- M.C. Escher

    Please see FAQ731-376 for tips on how to make the best use of Eng-Tips Fora.

    RE: Checking the FEA quality

    A star for MLoew - I will be circulating the links to these papers to a few colleagues who will benefit greatly!

    RE: Checking the FEA quality

    JulianHardy,

    Thanks for the ! I included these links in the FAQ Errors in FEM/FEA (FAQ727-980) along with some other links I found.

    Best regards,

    Matthew Ian Loew
    "I don't grow up. In me is the small child of my early days" -- M.C. Escher

    Please see FAQ731-376 for tips on how to make the best use of Eng-Tips Fora.

    RE: Checking the FEA quality

    When I check FEA models I mainly look at the following:

    1) Are the loads realistic? Shock load, fatique etc.
    2) Is the model over/under contrained? In other words does the displacement plots make sense.
    3) Is the CAD geometry good? Check for cracks, overlapping surfaces.
    4) Correct material property assumptions?
    5) Correct failure mode assumptions? Brittle or ductile etc.
    6) Benchmark against something similar that you know is OK.
    7) Has the solution converged? Adequate mesh density.

    If your model has sayisfied 1-6 then you should be in the right ball park. I hardly ever do mesh convergence studies unless I am sure I have satisfied all the above criteria first.

    RE: Checking the FEA quality

    Dear feajob:

    I am with a company that manufactures measuring systems for optical full-field strain and stress analysis. This is a  3-D non-contact measurement.
    They enable fast acquisition of the strains and stress in the observed surface section with regard to magnitude and direction.

    In order to present the acquired data of all different regions simultaneously, 3-D measuring results are mapped e.g. onto a 3-D CAD model of the tested component. Using known location point coordinates, both in the Field of Measurement as well as CAD model coordinate systems, the image re-sectioning process is accomplished. This 3-D data mapping capability provides a unique, direct comparison between the experimentally collected strain fields and the analytically calculated strain data by CAE.


    You might have a look at www.dantec-ettemeyer.com and follow the link to the strain and stress product lines.

    That might provide help with your callenging task.

    Regards
    Jim

     

    Btw.: Here, this approach might be superior to the strain gauge because the latter can only determine the strain at the point where it is mounted, averaged over its length and in the specified measurement direction. So, easy full-field hot spot detection becomes reality.


     

    RE: Checking the FEA quality

    (OP)
    Hi Jim,

    Thank you for your reply. It is an interesting approach, I wlll send your link to our test engineer.

    AAY

    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