Graphique en log (graduations logarithmiques)

1. Charger les données à mettre en X et Y

rayon = c(1e-5,1e-4,1e-3,1e-1,0.025,0.05,0.1,0.125,0.15,0.175,0.2,0.3,0.4,0.5,0.6,0.7,

0.8,0.9,1,1.5,2,2.5,3,3.5,4,5,6,7,8,9,10,11,12,13,14,15,20,100,1000,10000,1e5,1e6,1e7,1e8,1e9,1e10)

aire = rayon^2*4*pi

volume = rayon^3*4/3*pi

volume

ratio = aire/volume

ratio

2. Tracer le graphique en log

y_lab = expression(Ratio~~Surface/Volume~~(en~~nm^-1)) # ~~ pour les espaces

xlab=expression(Volume~~(en~~nm^3))

par(mar=c(5,6,3,4))

plot(volume,ratio,log="xy", xlim=c(0.00001, 1e11),ylim=c(0.001, 1000), pch=16,col="red",type="l",lwd=3,xlab=xlab,ylab=y_lab,cex.axis=1.3,cex.lab =1.5)

# le paramètre log="xy" implique un affichage en loge des 2 axes.

3. Tracer le graphique

panel.first = grid(lty = 2, lwd = 2) 

panel.first = grid(3, lty = 3, lwd = 2) 

lines(smooth.spline(volume,ratio),col="red",lwd=2)

abline(h=1)

4. Enregistrer les points à ajouter

#########################################################

# NANOTUBES et autres

#########################################################

# Nanotube : 1 à 10nm de diamètre, jusqu'à 1µm de long.

r = 5/2

H = 1e3

surface_nano = 2*pi*r*H+2*pi*r^2

volume_nano = pi*r^2*H

ratio_nano = surface_nano/volume_nano

#########################################################

# AMIANTE

#########################################################

# Amiante : 3 µm de diamètre, 100 à 200 µm de long : peuvent entrer dans les alvéoles du rat (source : INRS)

# Chez l'homme 1,5µm de diamètre pour max 360 µm Prenons diamètre 1 µm et H de 110 µm

r = 1e3/2

H = 110e3

surface_amiante = 2*pi*r*H+2*pi*r^2

volume_amiante = pi*r^2*H

ratio_amiante = surface_amiante/volume_amiante

#########################################################

# PLOMB

#########################################################

# 180 pm de diamètre

r = 90/2/1000

surface_pb = r^2*4*pi

volume_pb = r^3*4/3*pi

ratio_pb = surface_pb/volume_pb

#########################################################

# TiO2 nano

#########################################################

# 21 nm de diamètre

r = 21/2

surface_tinano = r^2*4*pi

volume_tinano = r^3*4/3*pi

ratio_tinano = surface_tinano/volume_tinano

#########################################################

# TiO2

#########################################################

# 250 nm de diamètre

r = 250/2

surface_ti = r^2*4*pi

volume_ti = r^3*4/3*pi

ratio_ti = surface_ti/volume_ti

5. Ajouter les points

points(volume_nano,ratio_nano,pch=8,col="brown",cex=3,lwd=3)

points(volume_amiante,ratio_amiante,pch=18,col="#40411E",cex=3,lwd=3)

points(volume_pb,ratio_pb,pch=16,col="#333333",cex=2,lwd=2)

points(volume_ti,ratio_ti,pch=15,col="#C03000",cex=4,lwd=4)

points(volume_tinano,ratio_tinano,pch=15,col="#C03000",cex=2,lwd=4)