¶metre aa=15.0d-6 !rayon du fil en metre. xLfil=3.0d-2 !longueur du fil en metre. xMs=6.6d5 !aimantation à saturation A/m psi=1.0d0 !angle d,anisotropie en degré Hk=4.0d1 !champ d'anisotropie intrinsèque en A/m sigma=8.0d5 !conductivité en (Ohm*m)-1 Hsurface=0.0d0 !champ d'anisotropie supplémentaire à lasurface (rapporté à Hk). C champ est appliqué sur un nombre de points = npsurf à partir de la surface Htorsion=0.0d0 !champ de tortion à la surface rapporté à Hk. Ce champ est linéaire à r/a. xnu=1.7565d11 !rapport gyromagnetique en Hz/Tesla alpha=1.0d-2 !constante d'amortissement de Gilbert Aex=1.0d-11 !constante d'échange en J/m &end &donnee courant=1.0d-3 ! amplitude du courant sinusoidal traversant le fil en A(rms) nformeI=1 !fonction du courant nformeI=1-->sinusoidale, nformeI=2-->carrée nformeI=3-->pulse frequence=1.0d8 !frequence du courant en Hertz instants=500000 !nombre d'intervalle de temps dans une periode npoint=100 !nombre de point considérés sur la direction radiale. le centre du fil et la surface exclus npsurf=1 !nbre de points qui appartriennent à la surface (points pour lesquels le champ Hsurface est appliqué). npsur >= 1 (si pas de champ de surace alors il faut mettre Hsurface=0) ngrad=0 !voir ngradfonc ngradfonc=1 !fonction distribution de l'anisotropie de surface sous la surface dure. si ngradfonc=1 --> polynome de degré 3 s'etalant sur ngrad intervalles entre points ! si ngradfonc=2 --> fonction exp f(r)=Hsurface*exp(ngrad(r-1)) rapamort=1 !rapport du coefficient d'armotissement de la surface sur celui de l'interieur du fil. spin libre --> rapamort=1 nharm=4 !nombre d'harmoniques considérées ll=1 !si ll=1 --> équation de Landau-Lifshitz, si ll#1 --> équation de Gilbert xprecision=1.0d-2 !precision sur la périodicité du mouvement de l'aimantation Mperiode=100 !nombre Max de periode pour chaque valeur de Hz nconveMax=100 !nombre de periode de récolte de données si la convergence n'est pas atteinte après ''Mperiode'' periodes nRK=1 !nombre de pas de runge-Kutta pour une meme valeur de champ effectif. HzHys1=30.0d0 !champ statique pour lequel sont tracées les composantes de M et H en fonction du temps et de la position radiale npHys=31 !voir nhys. npHys peu prendre les valeurs de 1 (centre du fil) jusqu'à npoint+2 (surface) nHzdyn=1 ! si nHzdyn=1 la composante sur z du champ dynamique est prise en compte ndomaine=1 !si ndomaine=1 un seul domaine sera calculé. si ndomaine#1 deux domaine seront calculés nbreHys=50 !nombre de points de l'axe du temps dans les fichiers de sortie (doit être inférieur ou égal à instants) codemag=0 !champ de démagnétisation axial codemag=Nd*Ms en A/m (ne pas mettre le signe -) &end &champ Hzi=50.0d0 !Hz initial rapporté à Hk Hzf=0.0d0 !Hz final // pasHz=-2.0d0 ! incrémentation de balayage de Hz rapporté à Hk si pasHz=0-->il n'y a pas de balayage et Hz considéré est Hzi Hzizone2=10.0d0 ! début de la zone 2 de Hz (à plus haute résolution) Hzfzone2=0.0d0 ! fin de la zone 2 paszone2=-0.2d0 ! pas de la zone 2 iretour=0 !si iretour=0 ==> balayage du champ dans une seule direction. Si iretour#0 ==> le champ est balayé en aller-retour pascourant=0.0d0 !pas du courant en Arms. si pascourant#0.0 --> courant est incrémenté jusqu'à courantf courantf=5.0d-4 nfrequence=0 !si nfrequence=0-->pas de balayage de la frequence, si =1 ou 2-->balayage linéaire (F(i+1)=F(i)+pasfreq) ou exponentiel (F(i+1)=F(i)*pasfreq) respectivement Ffinal=1.0d10 pasfreq=1.2d0 &end &bobine ! les paramètres suivants ne sont pas pris en compte si nHzdyn différent de 1. Ceux-ci concernent le champ magnétique dynamique axial appliqué à la surface latérale du microfil Hz0=0.0d0 !amplitude rapportée à Hk du champ dynamique axial Hz créé à la surface du conducteur par une bobine longitudinal. si pas de champ mettre Hz0=0.0d0 nfHz=1 !frequence de Hz dynamoque rapporté à la fréquence du courant dephas=0.0d0 !déphasage en degrés de Hz dynamique par rapport au courant. &end end