View source: R/find_densification_limit.R
find_densification_limit | R Documentation |
For porous (foam-like) materials, estimates the end of the plateau segment where the compression force starts to raise dramatically. In foams, this marks the end of pore space compression in favour of increasing wall material loading (Gibson&Ashby, 1997).
find_densification_limit(theData,touch_point,elastic_limit,gel_thickness=NULL)
theData |
Smoothened texture analyzer data, as produced by the smooth_texture_analyzer_data function; should contain only a one-way movement (down or up) for the chuck |
touch_point |
Point where the chuck touches the gel. This can be found from the Force - Distance diagram. where the force curve in a Force - Distance diagram raises abruptly |
elastic_limit |
Output of find_elastic_limit, which should be run before this function. |
gel_thickness |
If an estimation of the overall gel thickness is known, this is used to generate a better initial guess of the densification limit. |
The function uses both the slope and the actual stress (force per area, pressure) to find the initial touch point, hence the theData
argument must have at least the columns Distance, pressure and pressureSlope. The estimation is based on the idea to find the point where the measured pressure starts to rise substantially above the expected plateau pressure. Quantitatively, this is judged by comparing measured pressure to extrapolated plateau pressure (linear extrapolation based on minimal slope), with the requirement that the excess pressure exceeds the actual plateau pressure value.
Numerical value, indicating the position of the densification limit
Thomas Braschler
Gibson, L. J. & Ashby, M. F. Cellular Solids: Structure and properties. (Cambridge University Press, 1997).
smooth_texture_analyzer_data
data(sampleGelSmooth) touch_point=find_initial_touch_point(sampleGelSmooth[sampleGelSmooth$direction=="down",]) elastic_limit = find_elastic_limit(sampleGelSmooth,touch_point) densification_limit=find_densification_limit(sampleGelSmooth[sampleGelSmooth$direction=="down",],touch_point,elastic_limit) plot(Force ~ Distance, sampleGelSmooth[sampleGelSmooth$direction=="down",]) lines(rep(touch_point,2),c(0,11),col="black") lines(rep(elastic_limit,2),c(0,11),col="black") lines(rep(densification_limit,2),c(0,11),col="red") text(touch_point,11.5,labels="Touch point",adj=0.5) text(elastic_limit,11.5,labels="Elastic limit",adj=0.5) text(densification_limit,11.5,labels="Densification limit",adj=0.5,col="red")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.