Unexpected results from ANSYS workbench

11 replies · 4K views · last reply
Akarsh
Messages
6
Reaction score
0
Hi,
I have problem of finding the deformation of skin using a probe. I have identified the different layers of skin and made the model accordingly along with a probe in contact with the top layer of skin. The model is basically 15x15mm model with extrudes of different thickness for each layer.
The layers are stratum corneum, epidermis, dermis, hypodermis and muscle respectively. Each layer has different mechanical properties and is given isotropic elasticity condition.
The force of 0.1N is given to the probe and the bottom and the sides of the skin model are given fixed support condition.
The deformation in the first four layers are acceptable but the muscle layer acts like a cube of immense stiffness in a sense that no stress, no deformation, no strain occurs in the muscle layer. I cannot seem to identify the cause of this problem and i need some HELP figuring this out!

PS: I have attached the picture of the model, the sectional view of the results and the material properties of each layer
The model was created on SOLIDWORKS 2015
 
Attachments
  • ANSYS 16.0 mesh view of a 2D skin model: a cylindrical probe sits on a rectangular block of stacked meshed layers, with a 0 to 10 mm scale bar and axis triad showing Y up and Z to the right.
    Entire model with 5 layers and 1 probe.PNG
    17.9 KB · Views: 868
  • ANSYS 16.0 view of a finite element skin model with a red spherical probe touching the top layer, showing a scale bar from 0.000 to 0.700 mm with layers of differing colour below.
    Layers of skin.PNG
    12.6 KB · Views: 753
  • ANSYS Workbench contour plot of normal elastic strain (Y axis) in a layered skin model under a probe, with a color scale from -0.087835 min to 0.013158 max; strain is concentrated in the upper layers, and the bottom layer shows no strain.
    Normal strain (Y-axis).PNG
    40.6 KB · Views: 962
  • ANSYS total deformation contour plot in mm of a skin-layer cross-section under a rectangular probe, colored from blue (0 mm) to red (0.12938 mm max); deformation is largest directly beneath the probe and fades toward the blue, nearly undeformed lower region.
    Sectional view of deformation.PNG
    27.3 KB · Views: 1,404
  • Table of user defined materials with modulus of elasticity in MPa and Poisson's ratio: Stratum corneum 500, 0.49; Viable Dermis 1.5, 0.49; Dermis 0.02, 0.49; Hypodermis 0.002, 0.49; Muscle 0.25, 0.46.
    Material properties.PNG
    10.5 KB · Views: 831
I notice that your top layer has a modulus that is two orders of magnitude or more above that of any of the lower layers. So I am wondering if the composite restraining force at the tensile stress of primarily that layer plus the two underlying layers with all edges constrained is sufficient to resist any more deflection of the probe from its penetration point without requiring any support from the muscle layer.

I would suggest that you do a simple diaphragm analysis without the bottom restraint on the muscle and a similar diaphragm analysis without the muscle layer to see the results of those two cases.
 
With regard to the above, increasing the area of your model beyond the current 15x15 mm size will reduce the tensile restraining effect of your peripheral edge restraint on the relatively high modulus top layer.

Also a test you could perform on your model would be to temporarily remove the top three layers and place a test point load on only that remaining muscle layer to see how it responds in order to determine if there is some basic issue with its modeling or meshing.
 
  • Like
Likes   Reactions: Akarsh
@JBA i tried doing the same analysis without the muscle layer and keeping the bottom of the hypodermis layer fixed. I got the similar results except that there is no muscle layer in the results.
PS in this result i did not add the friction between the skin and the probe. If i did add it then i would i get the same result as the one's i initially posted
 
Attachments
  • ANSYS Static Structural total deformation contour plot in mm, showing a red rectangular probe pressed into a layered blue-to-red deformed block, with a colour legend from 0 Min to 0.11303 Max and a scale bar of 0 to 8.000 mm.
    Deformation.PNG
    40.6 KB · Views: 703
From what I see comparing the results, it would appear that the diaphragm tension in the top (and to some degree, the second layer) is the controlling factor due to its high modulus. The modulus in the third and fourth layers is so low that they are almost fluid in their behavior.

To test if this is true, rerun the problem with a reduced modulus (possibly 300) for the top layer, that will transfer more of the loading to the supporting underlying layers and to the muscle layer. If it is true, you should see a deeper probe penetration and some deflection of the muscle layer.
 
  • Like
Likes   Reactions: Akarsh
@JBA i removed the top layers and placed the probe on the muscle layer and saw whether it would deform or not.
It did deform and caused some normal stress
I know the mesh is not refined properly as i just wanted to see whether i could get any deformation for the muscle layer
 
Attachments
  • ANSYS Static Structural screenshot of normal stress in the Y axis in MPa, showing a coarse triangular mesh on a rectangular block with a green cylindrical part on top, with a colour legend from 0.25754 Max to -0.35299 Min and Max and Min markers, scale bar in mm.
    normal stress.PNG
    33.9 KB · Views: 751
  • ANSYS Workbench total deformation contour plot in mm, showing a red rigid probe pressed into a blue block with a cyan-to-red deformation zone beneath its tip; legend scale from 0 to 0.073854 mm.
    Deformation1.PNG
    38.5 KB · Views: 703
@JBA i changed the young's modulus of the top layer to 30MPa and ran the simulation,
The results were that there was more deflection in the model but the muscle layer was still not deformed at all

I tried changing the young's modulus of the muscle layer to a value lower than the hypodermis (the layer prior to the muscle layer). The results from this trial shows deformation in the muscle layer
 
Looking at your initial strain diagram you can see that as the loading progresses down through the layers the strain area is also expanding and therefore as this area of pressure is expanded it reaches the point that the area is so large and evenly distributed that it is deforming the entire surface only a minuscule amount.

If you look at that initial strain display values chart you will see -00248... is the minimum negative strain value given and .00193... is the minimum positive strain value given and there is no distinguishable color variation in the color band between those two values. As result, if the strain across the surface of the muscle layer lies in the zone between those values, it is not discernible even if it is not actually 0.000.

Considering all of the above, it would appear that is the most probable answer to what you have perceived to be a problem, that does not actually exist, in your analysis.
 
  • Like
Likes   Reactions: Akarsh
So when the force is applied through a series of other layers the strain would be lesser?
 
As a test of my theory I ran the below simple strain calculation for your .1 N force distributed over the 15x15 mm surface area of the muscle material (think of it as a the end of a square bar of your muscle material) and the results compared with the below copy of your original strain diagram with a solid color strain range from + .0019 to - .0092 tend to confirm that a similar distributed loading on you muscle area might fall between those two boundaries. (please forgive the units conversions, I am simply more confident using USA units)
Spreadsheet calculation table listing force, muscle surface area in mm² and in², stress in psi, modulus in MPa and psi, and strain in in/in and mm/mm, with values such as -0.100 N and strain -0.0045156 mm/mm.

ANSYS static structural plot of normal elastic strain (Y axis) in a block under a loaded punch, with a colour legend from 0.013158 max to -0.087835 min and a 0 to 7 mm scale bar.
 
Last edited:
In response to your last question, It is the extreme softness of the third and forth layers that causes them to act like a fluid trapped between the top two layers and the bottom muscle layer and spread the probe load across almost the entire surface of the muscle layer that results in that layer's very small, but not zero, deflection.

If you look carefully at the above strain plot you will see that the horizontal grid line between the fourth layer and the muscle layer is not actually flat, it has one negative offset just inside the second grid line from the outside end on both ends; and second negative offset just outside of each vertical grid line next to the centerline that indicate a small but present strain in the top surface of the muscle layer.
 
  • Like
Likes   Reactions: Akarsh