Convergence analysis for plate with a hole
Convergence analysis for plate with a hole
(OP)
Hi all,
I have a plate of 300*300 mm and a hole of 20 mm. I have applied a tensile pressure load of 200 MPa. I want to find the stress concentration factor. I have developed a quarter-symmetric model. I have used first-order plate elements. And this, is what I find when I do the convergence analysis.
No. of elements Stress Concentration factor
29 1.236666667
116 1.474666667
464 1.544666667
1856 2.188666667
7424 2.599333333
29696 3.733333333
118784 3.509333333
The SCF for a plate with a hole for the dimensions I have should be 2.81 (from Strength of Materials text book). But, I am unable to converge it to that value. What is happening here? Any sort of stress stiffenning?
Thanks
Kumar.
I have a plate of 300*300 mm and a hole of 20 mm. I have applied a tensile pressure load of 200 MPa. I want to find the stress concentration factor. I have developed a quarter-symmetric model. I have used first-order plate elements. And this, is what I find when I do the convergence analysis.
No. of elements Stress Concentration factor
29 1.236666667
116 1.474666667
464 1.544666667
1856 2.188666667
7424 2.599333333
29696 3.733333333
118784 3.509333333
The SCF for a plate with a hole for the dimensions I have should be 2.81 (from Strength of Materials text book). But, I am unable to converge it to that value. What is happening here? Any sort of stress stiffenning?
Thanks
Kumar.





RE: Convergence analysis for plate with a hole
Second, are you sure you are applying the theory correctly? When using the S.C. are you subtracting out the hole diameter from the area when calculating the nominal stress? The nominal stress in this case is actually 200*300/280 MPa. After that you are doing something wrong, this should be a simple problem to get right.
RE: Convergence analysis for plate with a hole
RE: Convergence analysis for plate with a hole
This is stress singualrity phenomena and not stress stiffening effect.
Just do stress singularity search and find threads related to this topic.
Vineet
vineet
RE: Convergence analysis for plate with a hole
Regarding stress stiffening effect, it could be the case, because thin walled elements will undergo this phenomenon. But what i probably don't understand, that you model your thin structure under pressure with plate elements, it should be modeled with shell elements, since the structure will have curvature under pressure load. Please see the theory about the difference between "plate" and "shell", any suggestion?
regards
RE: Convergence analysis for plate with a hole
I have refined my elements using h-method. As you can see from the data, each element has been split in the mid-side to obtain 4 elements. So, 29 elements in the first run has become 29*4=116 in the second run.
I have applied a uniform tensile load of 200 MPa on one edge of the plate. It is not a bending load. Hence, I think the use of plate element is justified. Any arguments and suggestions are most welcome.
RE: Convergence analysis for plate with a hole
Other aspect is that the stress after 6. step (elements = 29696), the stress is 3.733333333, whereby the theoretical stress should be 2.81, i.e. above the theoretical solution (the model is flexible then the reality). Normally FEM Model is stiffener then the reality/theory, so it should have small stress and the stress would increase with the inreasing of elements number. The step number 1 to 5 seems to be right, but 6. and 7. step is strange. I lose the theoretical explanation for those steps. May be you better do FeaGuru's suggestion and compare the solution to each other, Furthermore you have better control over the quality of the mesh.
regards
RE: Convergence analysis for plate with a hole
Apart from that, if I understand correctly your refinement procedure, in the last step you should have elements with a side of about 0.7 mm: this could be still a bit high to get a good approximation for the local stress. Moreover it is not clear how you account, in refining, for the hole curved edge: you could have locally elements with bad shapes.
What about trying a further step? However you should refine the mesh only around the hole, using there only nearly square elements in a polar, rather than cartesian, pattern.
prex
motori@xcalcsREMOVE.com
http://www.xcalcs.com
Online tools for structural design
RE: Convergence analysis for plate with a hole
Have you applied right boundary conditions for your quarter symmertric model?
I feel stress singularity due to boundary conditions may be behind this phenomena.
Regards,
elogesh
RE: Convergence analysis for plate with a hole
I have applied the right boundary conditions, since the results of the quarter-symmetric model matched with the full scale model.
I will go ahead and refine only the area around the hole and try to get a converging solution. Any more suggestions and tips are highly welcome.
Thanks
Kumar.
RE: Convergence analysis for plate with a hole
Now you are a little bit closer to your goal, but still a difference exists.
If you want more help may you describe how are your boundary conditions applied?
prex
motori@xcalcsREMOVE.com
http://www.xcalcs.com
Online tools for structural design
RE: Convergence analysis for plate with a hole
So the solution should converge way before 100000 elements.
Maybe in your model, when splitting the elements, the boundary conditions were constantly applied only to the original corner nodes, not the the new ones (additional stress peaks)??
RE: Convergence analysis for plate with a hole
Can you elaborate about the model and the boundary conditions used for the analysis?
Whether it is full model or quarter sysmmetry model?
Regards,
amar
RE: Convergence analysis for plate with a hole
Denzil, could you tell me what was your maximum stress and the nominal stress and the way you modelled.
y
|
|
* * *
* *
* *
* *
* *
* * *--->x
Thanks
Kumar.
RE: Convergence analysis for plate with a hole
I modelled a quarter of the plate (150x150 mm, R10mm quarter hole at the corner) using MSC/Patran automatic paver mesh. Element size along the hole edge was 0.5 mm and 6 mm at the distance from the hole.
Boundary conditions were: SymmetryX (1,5,6) along left edge, SymmetryY (2,4,6) along lower edge, and Z (3) at every node.
Upper edge of the model was subjected to constant pressure=1. Thus, the average stress in the weak section is 150mm/140mm*1MPa = 1.0714 MPa.
Maximum stress from FE analysis was 3.04 MPa i.e. peak factor 2.837.
Regars Denzil
RE: Convergence analysis for plate with a hole
reason for avaoidance of plate elements in the present problem is that...even if u use a plate element at lower refinement u should be intrested in getting stresses from higher order derivatives and at higher refinements u should be intrested in getting stresses from lower order derivatives....the end result being for the present problem its better if u go in for a plane stress element so that u reach a constant strain state earlier and that your results converge faster and not llike as we see in ur solution.
Refer to paper by barlow on stress extraction and refinement
may be i shall let u now the journal and vol in next post
hope it helps
regds
rAJ
Raj
RE: Convergence analysis for plate with a hole
to move. Call Node 1 the location of the highest
stress. From Kumar's drawing, that should be the
top of the hole, on the x=0 axis. If one fixed
the y displacement to be zero at Node 1, that
would not allow Node 1 to move down (Poisson
effect) when you pulled to the right. Does
the model deform as expected? What do the stress
contours look like? Refinement seems to be
the key answer, nevertheless a sparse mesh
does not explain why the stresses keep going up.
A bad boundary condition like a point constraint
does explain why the stresses keep going up.
RE: Convergence analysis for plate with a hole
Here's what I did, using plane stress quads in p-version code StressCheck, took about 10 minutes to build mesh and get result. Solution converges quickly as there are no singularities. Kt=3.038 based on far field, 2.835 based
on net section stress. Only 1344 D.O.F. needed for less than 1% error. Same numbers for plane strain of course.
Now fix the node at the top of the hole in y-direction.
=================================
| DOF | Max. Sx
================================
| 27 | 1.964885E+000
| 95 | 2.526961E+000
| 203 | 2.793709E+000
| 351 | 3.009912E+000
| 539 | 3.115877E+000
| 767 | 3.250458E+000
| 1035 | 3.200530E+000
| 1343 | 3.320644E+000
==================================
Ugly, isn't it?
s