DELTA 46 2842 607 SVNKRTRRfunction [vn,fn] = drndt(rn,appliedstress,MU,NU,mob,rann) %force on dislocations fn = nodal_force(rn,appliedstress,MU,NU,rann); %mobility function vn = feval('mobrelax',rn,fn,mob); ENDREP DELTA SVNfunction vn = mobrelax(rn,fn,mob) N=length(rn(:,1)); b=zeros(N,3); b(1:N,1:3)=rn(1:N,3:5); %based on force compute dislocation velocity (linear mobility) vn = fn * mob; %project velocity onto glide plane for i=1:N, b_norm=sqrt(b(i,1)^2+b(i,2)^2); % b_norm=norm(b(i,1:2)); if b_norm~=0 b_unit=b(i,1:2)/b_norm; vn(i,:)=dot(vn(i,:),b_unit)*b_unit; end endENDREP DELTA 46 4130 259 SVNeYg~gx^Sr +B/KN6#R`+ɥ\xov͠ւTmZdފF!wbw'|m-ϲ4hCG4%^Aۉ R)PQ ۳rv%]xVaƝq@uv=$1@5u@8mu#hI2Z|:H:ǥՋҏ!_߾ ~`JxSWyuҳ{ݲ5V44ct΍FAGDCOWkőϪ (vxrfC_kNKcbưxY8-''Oeĕu`U*'pвxc1 zk U|Pe3S<_ƫsz>7l L8ze W