# test en reference a l'exo 51 de MSD mmk2

#MeshFile Sphere_Creuse-fin.msh
MeshFile Sphere_Creuse-grossier.msh

# Definition de la sphere creuse
CutMeshSphere 0 0 0 1.0
CutMeshSphere 0 0 0 2.0
CutElasticDomain 52 1e9 0
# Attention : l'effet poisson peut générer des interpénétrations suivant le déplacement imposé

#-------------------------------------------------------
# rho = sqrt(x^2+y^2+z^2)
# phi = acos(z/sqrt(x^2+y^2+z^2))
# theta = acos(x/sqrt(x^2+y^2))               si y>=0
#         2*3.14159265-acos(x/sqrt(x^2+y^2))  si y<0
#       = (1+sign(y))/2*acos(x/sqrt(x^2+y^2))-(sign(y)-1)/2*(2*3.14159265-acos(x/sqrt(x^2+y^2)))
#-------------------------------------------------------

# Ur = Ar + B/r^2 = U
#    = sqrt(x^2+y^2+z^2)+1/(x^2+y^2+z^2)
# U  = Ur*sin(phi)*cos(theta) ex
#    + Ur*sin(phi)*sin(theta) ey
#    + Ur*cos(phi)            ez

StableLagMultFieldFunction 2 0 (sqrt(x^2+y^2+z^2)+1/(x^2+y^2+z^2))*sin(acos(z/sqrt(x^2+y^2+z^2)))*cos((1+sign(y))/2*acos(x/sqrt(x^2+y^2))-(sign(y)-1)/2*(2*3.14159265-acos(x/sqrt(x^2+y^2)))) x y z
StableLagMultFieldFunction 2 1 (sqrt(x^2+y^2+z^2)+1/(x^2+y^2+z^2))*sin(acos(z/sqrt(x^2+y^2+z^2)))*sin((1+sign(y))/2*acos(x/sqrt(x^2+y^2))-(sign(y)-1)/2*(2*3.14159265-acos(x/sqrt(x^2+y^2)))) x y z
StableLagMultFieldFunction 2 2 (sqrt(x^2+y^2+z^2)+1/(x^2+y^2+z^2))*cos(acos(z/sqrt(x^2+y^2+z^2))) x y z

Solver 3