When and how to "ignore" a singularity
When and how to "ignore" a singularity
(OP)
I am modelling a shipping container wall with rectangular openings. I'm using risa 3d, which, from my research on this site is a "lower" powered FEM program, compared to what yall mechanical guys use. I'm using shell elements for the wall and beam elements at the edges in the model. I have singularities at the corners of the openings which do show an increased stress as I decrease the mesh size in these areas. There may be singularities at the corner supports as well, but the stresses are so much lower than the opeing corners, that they don't show up. So, how do I deal with these infinite stress locations and get usefull results from this model. I'm trying to determine deflection, a realistic max stress in the corrugated wall, and perfomance of the boundary beam elements. Unfortunately I am way ahead of the corresponding real world load tests which we will be doing in the future.
thanks
an ps I have learned so much from this forum and it is indispensible to a one person office, thanks
thanks
an ps I have learned so much from this forum and it is indispensible to a one person office, thanks





RE: When and how to "ignore" a singularity
RISA3D is great for certain applications...just not sure this is one of them.
Garland E. Borowski, PE
Star Aviation
RE: When and how to "ignore" a singularity
RE: When and how to "ignore" a singularity
RE: When and how to "ignore" a singularity
some tests ...
1) fill in the hole (i don't think this'll fix the problem)
2) replace the beams with rods (just for giggles, to see what'll happen)
3) add full plate properties to the panel elements (maybe wack up their stiffness to see what happens)
4) add a beam (maybe several) across the panel (joining the cut-out frame to the edge of the panel) ... constrain the bending freedoms
are you only applying in-plane loads (there was a thread about the compression strength of corrugated container walls recently) ? ... even if you have no out-of-plane loads (and so nothing driving stresses into these out-of plane stiffnesses) FE codes don't like zero stiffness
RE: When and how to "ignore" a singularity
Crisb:
I'm modellilng a 40' long wall so a fillet or chamfer at the hole corner seems like a problem of another scale. I'm pretty much stuck with 4node shell elements and Risa3d does not do any non-linear analysis. I think this is ok, because we don't want to load the wall beyond any signifigant deformation (something that would crack applied finishes).
rb1957
1. I started with a no holes in the wall model.
2. I'm not sure of the difference between a beam and a rod. I have a choice of beam or plate elements. I can apply many different properties to either.
3. Little unclear on the full plate properties. My plate elements are as wide as the distance between vertical bends in the metal and about the same height.
4. I have messed around with this but it seems like, and I'm not entirely sure, that nothing is really stiff enough to erase the singularity which sorta makes sense.
I am only loading in the plane of the plate elements. No wind/lateral loads on the current model. Intuition says the model is acting correctly, except for the singularity.
RE: When and how to "ignore" a singularity
corus
RE: When and how to "ignore" a singularity
That is what I thought and we are not doing a fatique analysis. Can you guide me to some literature that would back this up? I am beginning to understand why it exists and why I can "ignore" it, but I need to have more of an understanding in order to explain it to my colleagues.
thanks
RE: When and how to "ignore" a singularity
From BS EN 13445 (the later version of BS 5500):
C.2.5
peak stress
that part of stress which is additive to the respective primary and secondary stresses, to form the total
stress
NOTE 1 Peak stresses do not cause any noticeable distortion and are only important to fatigue and brittle
fracture in conjunction with primary and secondary stresses.
NOTE 2 Peak stresses also comprise deviations from nominal stresses at hole edges within tube-hole fields
due to pressure and temperature, in which case the nominal stresses are derived from equilibrium of forces
considerations.
corus
RE: When and how to "ignore" a singularity
The problem with this type of analysis is that we start with the assumption that the stress tensor is symmetric. Most FEM programs make this assumption. The shear stress is assumed to be symmetric, Txy = Tyx. Resissner (QAM ~ 1947) demonstrated how to handle the edge stresses of a plate that had a zero shear stress along one corner and a finite shear on the other face.
Grioli (Acad. Press 1962) provides a detailed elastic analysis when it is assumed that the stress tensor is nonsymmetric. This author Citerley, (J. Franklin Inst 1966) demonstrated how to impose this non-symmetry condition using Lagrangian multipliers.
The engineering approach would be to establish the basic state of stress adjacent to the re-entrant corner using any verified FEM program and then apply well established stress concentration factors obtained from Peterson's Stress concentration book.
RE: When and how to "ignore" a singularity
Is Grioli saying that isotropic linear elastic materials have non symmetric stress tensors? That seems to contradict the fundamental assumption, the isotropy. It's as if you make the assumption that the process is 'isentropic' then claim the entropy did increase.
Wish I had the book so I could see for myself just what Grioli is talking about. Is this the Grioli book?
Mathematical theory of elastic equilibrium: (recent results) (Ergebnisse der angewandten Mathematik) (Unknown Binding)
It's Academic Press, 1962.
RE: When and how to "ignore" a singularity
The symmetry of Cauchy's (i.e. true) stress comes from the conservation of the angular momentum (assuming no couple tractions) and has nothing to do with any material formulation.
Some stress measures are not symmetric, for example the first Piola-Kirchhoff stress tensor (sometimes referred to either as the nominal stress or as the transpose of the nominal stress) is not symmetric.
Some of the symmetries of the elasticity tensor [E_ijkl] are shown assuming the symmetry of the stress and strain tensors.
RE: When and how to "ignore" a singularity
This method was employed to the end grain problem of a cylinder supported at its outer surface. The cylinder is cooled and thermal stresses result. At the wall interface, assuming a symmetric stress condition, infinite stresses are predicted. The same is true for a bi-metallic strip. The Picone method suggests otherwise.
Xfer: Thank you for a more complete explanation regarding symmetry.
RE: When and how to "ignore" a singularity
I have some reading to do now. I will post my current analysis soon and perhaps we can discuss some specifics. This is a new area for myself and I really appreciate everyones thoughts. I am analyzing thin corrugated metal and I have discovered another stress concentration where the corrugation transisions to flat. I haven't found any direct reference for dealing with this type of discontinuity yet although I am confident that the metal will deform and we aren't dealing with any cyclic loading or crack propogation, yet. Again, thanks to all for the time and I hope this discussion continues.
RE: When and how to "ignore" a singularity
Not to try to cover what the excellent posts above discuss...but one aspect of RISA's graphic output (the stress contours) is that you can identify only certain ranges of stress on the plot window to get a better look at stress gradations in the low stress regions and ignore, temporarily, the hot spots at the corners.
On the plot window, under the selection of the type of stress/moment/etc. it lets you set a range. So this alters the contours to read more detail in just the range areas you are concerned with.
RE: When and how to "ignore" a singularity
If modeled correctly, any FEM program should be able to give the correct results.