1st order versus 2nd order tetrahedral elements
1st order versus 2nd order tetrahedral elements
(OP)
Hi,
I usually deal with 3d complex geometries and I use tet10 (2nd order solid elements) for static analysis. My models are very big (talking about 800,000 nodes). We have good hard-wares and I get my results in les than 1 hour. But, now, I have to do some contact analysis. Since, 2nd order elements are not suggested for contact analysis and my models are big. I have converted my 2nd order elements to 1st order elements. In this way, for the same number of elements, I have 6 times less DOFs in my model.
I do contact analysis in an efficient way, but, as you know, tet4 is stiffer than tet10. So, I am worry about the quality of my results. I found this paper:
http://www.mscsoftware.com/support/library/conf/au...
It confirms my anxiety. At this point, I don’t know, what should be sacrificed? Accuracy or efficiency? I am thinking about a mix model. Only 1st order elements in contact area and 2nd order elements elsewhere, but, it won’t be efficient (solving iteratively big models!). Please let me know your suggestions.
Regards,
AAY
I usually deal with 3d complex geometries and I use tet10 (2nd order solid elements) for static analysis. My models are very big (talking about 800,000 nodes). We have good hard-wares and I get my results in les than 1 hour. But, now, I have to do some contact analysis. Since, 2nd order elements are not suggested for contact analysis and my models are big. I have converted my 2nd order elements to 1st order elements. In this way, for the same number of elements, I have 6 times less DOFs in my model.
I do contact analysis in an efficient way, but, as you know, tet4 is stiffer than tet10. So, I am worry about the quality of my results. I found this paper:
http://www.mscsoftware.com/support/library/conf/au...
It confirms my anxiety. At this point, I don’t know, what should be sacrificed? Accuracy or efficiency? I am thinking about a mix model. Only 1st order elements in contact area and 2nd order elements elsewhere, but, it won’t be efficient (solving iteratively big models!). Please let me know your suggestions.
Regards,
AAY





RE: 1st order versus 2nd order tetrahedral elements
Firstly, Abaqus have a specially formulated tet10 element that allows contact (C3D10M)
Secondly, experienced Abaqus users won't touch this element with a barge pole! It's stress recovery is very inaccurate (Abaqus at first denied there was a problem to me, until the very next slide on their ppt display said they were working on the problem!).
Thirdly, experienced Abaqus users (at least three seperate consultants/companies) use mixed meshing, 8 noded bricks for the contact and 10 noded bricks elsewhere. They connect 4 bricks to 2 tet10 elements, thus ensuring continuity of displacement across the join, and living with the dis-continuity in stress.
Your model using 2nd order elements does strike me as too large to be solved iteratively. Why don't you produce and solve a coarser mesh, which will still provide accurate displacements and hence loading paths, then use the results from this "coarse" model to be fed into fine meshed sub-models as boundary conditions?
RE: 1st order versus 2nd order tetrahedral elements
corus
RE: 1st order versus 2nd order tetrahedral elements
Corus, using brick elements for a landing gear structure is not an easy task. We can do that, but, it takes too much time and my boss won't be happy.
Johnhors, Corus, I am not using Abaqus. I am using MSC products (Nastran, Marc). Well, I never did submodeling in these softwares, I can try this idea. I think that it would be equivalent of finding the interface nodal forces(on a small portion of my model, where contact happens) with Nastran in a linear analysis and applying them on the interface nodes of my small model for non-linear analysis with Marc. Please correct me, if I am wrong.
AAY
http://www.geocities.com/fea_tek/asd/
RE: 1st order versus 2nd order tetrahedral elements
Yes you are quite right on the necessity to use tet elements for meshing of the complex parts of a landing gear.
My procedure for sub-modelling is the reverse of what said!
I run a coarse mesh of the assembly with full contact, which I know will not produce good stress results, it will however produce good displacements which I can use in a sub-model.
When I sub-model, I first sub-divide a complex part into several or many connected volumes (sharing common boundary faces). I make planar cuts through the volume around areas of concern. I then have nice convenient flat boundaries available for sub-modelling these regions. From the analysis of the complete coarse mesh which can be a non-linear contact problem, I obtain the nodal displacements on these previously defined planar cut faces and apply these as enforced displacements on the boundaries of the sub-model.
The technique is very easy to do using any FE system. Checking is a breeze, just run the sub-model without enhancing the mesh and you should get identical results to the whole model, you can then safely refine the mesh of the sub-model to obtain a converged stress.
RE: 1st order versus 2nd order tetrahedral elements
With reference to brick meshing, I'd tell the manager that if he wants the job doing properly then he'll have to wait. If he wants a pretty picture of no real value then just use a coarse tet mesh for your analysis, and look for another place of work where they do things properly.
corus
RE: 1st order versus 2nd order tetrahedral elements
On the other hand, our models are very complex. Using brick elements, is a very time consuming task. I don't know how we can justify that? We can refine more our tet10 mesh , or use brick elements at some places with tet10 mesh (as johnhors suggested) and get the good solution.
RE: 1st order versus 2nd order tetrahedral elements
I may have the wrong end of the stick here but be careful about generalising your assumptions on tet elements. First order tets are useless in anything but pure thermal analyses, however second order tets are very good for modelling complicated geometry and for use in structural problems AS LONG AS YOU UNDERSTAND THEIR LIMITATIONS. It sounds as though your geometry wil probably have to use these according to your description - this is not a problem. If you need to mesh a structure in which the stresses/response is dominated by bending consider that second order tets are quite stiff (1st order tets are useless in this respect). You will need to use about two/three tets through thickness for adequate resolution of stress results, more if your analysis includes non-linear material props. For displacement resolution only you may require less through thickness. If your problem is dominated by planar loads (membrane stresses as in in-plane bending for example) you could get away with less than three. Whatever, don't be afraid to use tet elements for contact problems. Most CAD geometries we import can only be tet meshed in most cases. We could get mapped or hex-dominant meshes, but because of the time taken to produce this offsets the benefit of a hex mesh. Once meshed wth tets, these components are used to carry out complex contact problems. The key is knowing how to mesh your component(s), knowing the limitations of your elements and knowing exactly how to tweek the contact to obtain best results.
Cheers,
-- drej --
RE: 1st order versus 2nd order tetrahedral elements
Again your comments are sound. With regard to brick or tet10 element usage, neither will converge to a "correct" solution, just that bricks will converge to solution that is more "correct" than tets and with far fewer elements on a "simplistic" test model. I said "simplistic" because there is no meshing software in this world that could mesh a fully detailed landing gear leg with reasonably shaped brick elements in "reasonable" numbers.
The last time I used brick elements on a landing gear leg (in days before tet meshers) was using a dumb Tektronix terminal on a VAX/VMS mainframe. I built geometry from scratch (no CAD import available, in fact no 3D solid existed). I ignored all small details (chamfers, fillets, grease holes... etcetera). The volume had to consist of thousands of simple shaped map-meshable regions i.e. bricks or wedges. A drawing refresh when a lot of lines were being displayed would give you long enough to take a coffee break (OpenGL was a distant dream). How long did I take to complete the model? Answer, a whole year!!!!
I'm sure with the improved software and hardware, the time would now be considerably less, but probably at least 2 months, BUT that is still ignoring all the small details! Such a model would only be good enough for ultimate strength analysis and completely useless for fatigue analysis.
Using tet elements (in large numbers) is the only practicable solution for meshing very complex structures.
RE: 1st order versus 2nd order tetrahedral elements
Johnhors, thank you for your comments. It is an interesting approach. Since, the convergence rate of displacement fields are higher than those of stress fields. You assume that displacement fields obtained from a coarse mesh is good enough to be applied on a small portion of your mesh (sub-model) with a very fine mesh. I guess that the number nodes on the common boundary faces must be the same? So, you refine your sub-model without touching the boundary nodes on common faces. Now, I would like to have a rough idea about the following ratio for your sub-model:
(DOFs of your final fine mesh) / (DOFs of your initial coarse mesh) = ?
Thank you,
AAY
RE: 1st order versus 2nd order tetrahedral elements
On a second point, quadrialteral elements have a problem with the mid-side nodes in contact which can cause problems in convergence and produce odd results at the mid-point positions. As tet elements can only be used in their quadrilateral form I find it is better to use linear elements where ever possible in contact problems. Even though it means spending a little time on reviewing the details needed to be modelled, on the mesh quality, and a little time on deciding on what is to be obtained from the model.
corus
RE: 1st order versus 2nd order tetrahedral elements
Yes you are right again, the mesh density on the common boundary faces remains the same (I mesh these with a density greater than that used in the rest of the coarse model - a kind of half way house between coarse and fine meshing).
As for DOF's, this is determined by hardware resources, on the machines I use, any model below 400,000 nodes will run with ease, the largest I have run is 440,000 nodes using either Lusas or Abaqus. So the coarse mesh of the whole component is usually about at the max capability i.e. 1.2 to 1.5 million DOF's. Sub-models after refinement approach the same size. So the ratio that you want is the ratio of the volume of the whole model to the sub-model, this can be anything between 10 to 500.
RE: 1st order versus 2nd order tetrahedral elements
I'm afraid that omitting any geometric features is politically unacceptable! If I do remove a grease hole I would have to justify my action. How do you do that? Answer - run the complete geometry model intact! Sub-modelling is the only feasible way to produce models good enough for fatigue analysis that include all possible stress raisers.
RE: 1st order versus 2nd order tetrahedral elements
corus
RE: 1st order versus 2nd order tetrahedral elements
Where I work the Catia model is sacrosanct. The model that I use for analysis is the same one used by manufacturing. There is only one current model allowed per component for auditing purposes. I fought and lost the battle to make allowances for engineering judgement some years back. This is now the common (enforced) practice in aerospace as a whole. Using engineering judgement to omit geometric features for an analysis would be viewed upon in a very dim light by litigation lawyers, if the occasion arose (and it does!). You simply have no choice but to mesh Catia solid models as-is.
RE: 1st order versus 2nd order tetrahedral elements
It's probably fairer to assume that lawyers are unable to reason and ensure that you have plenty of space on your machine for the massive models that must occur. It makes you wonder how the Wright brothers managed without computers.
corus
RE: 1st order versus 2nd order tetrahedral elements
AAY
http://www.geocities.com/fea_tek/asd/
RE: 1st order versus 2nd order tetrahedral elements
I know this might be a little late, but in MSC/PATRAN Fields you can easily convert (interpolate) displacement fileds from a global to a local model
See:
http://support.mscsoftware.com/kb/results_kb.cfm?S...
Good luck
RE: 1st order versus 2nd order tetrahedral elements
AAY