champgeoquad "ChampGeo" "hm16-10000.geoquad"
scalaire "AlphaMin" 0.
scalaire "AlphaMax" 1.
scalaire "MaxIt" 26
//
scalaire "Epais" 0.0001
scalaire "Normal" f(Epais)=-Epais*(z*z-y*y-x*x)/((x*x+y*y+z*z)*sqrt(x*x+y*y+z*z))
scalaire "p_0x" f(Normal)=Normal*x
scalaire "p_0y" f(Normal)=Normal*y
scalaire "p_0z" f(Normal)=Normal*z
scalaire "zoro" 0
empilementv3d "p_0" [p_0x,p_0y,p_0z]
vectoriel3d "PointEvaluation" [0.,0.,1.]
// si L/h = 100
//scalaire "UZTheo" -7.73688e-3
// si L/h = 1000
//scalaire "UZTheo" -8.08931e-3
// si L/h = 10000
scalaire "UZTheo" -8.19934e-3
scalaire "zero1d" 0
vectoriel3d "zero3d" [0,0,0]
scalaire "E1" 2.e5
scalaire "Nu1" 0.3
v3dcrourav "U" [0,0,0]
v3dcrourav "UNc" [0,0,0]
scalaire "DU" f(UX,UNcX)=UX-UNcX
scalaire "DV" f(UY,UNcY)=UY-UNcY
scalaire "DW" f(UZ,UNcZ)=UZ-UNcZ
empilementv3d "DiffUUNc" [DU,DV,DW]
scalaire "G1" f(Nu1,E1)=E1/(2.0*(1+Nu1))
to4elasticite "E" [E1,E1,E1,Nu1,Nu1,Nu1,G1,G1,G1]