Small issue - for a 3 mm disc it should be 3e-3 not 3e-4 in lines:
import numpy as np
my_model.mesh = F.Mesh1D(vertices=np.linspace(0, 3e-4, num=1001))
it should be:
import numpy as np
my_model.mesh = F.Mesh1D(vertices=np.linspace(0, 3e-3, num=1001))
Small issue - for a 3 mm disc it should be 3e-3 not 3e-4 in lines:
it should be: