XIAN-FEM-2026June/三维matlab代码/matlab 3D二阶矢量本征问题/BF_CurlE.m

43 lines
861 B
Matlab

function out=BF_CurlE(u,v,w,n)
if n==1
out=[0,0,-2+9*u];
elseif n==2
out=[0,0,-2+9*v];
elseif n==3
out=[0,2-9*u,0];
elseif n==4
out=[0,2-9*w,0];
elseif n==5
out=[0,-2+9*u,2-9*u];
elseif n==6
out=[0,7-9*u-9*v-9*w,-7+9*u+9*v+9*w];
elseif n==7
out=[-2+9*v,0,0];
elseif n==8
out=[-2+9*w,0,0];
elseif n==9
out=[2-9*v,0,-2+9*v];
elseif n==10
out=[-7+9*u+9*v+9*w,0,7-9*u-9*v-9*w];
elseif n==11
out=[-2+9*w,2-9*w,0];
elseif n==12
out=[7-9*u-9*v-9*w,-7+9*u+9*v+9*w,0];
elseif n==13
out=[-3*u,-3*v,6*w];
elseif n==14
out=[3*u,-6*v,3*w];
elseif n==15
out=[3*u,3*v,6-9*u-9*v-6*w];
elseif n==16
out=[-3*u,6*v,3-9*v-3*w];
elseif n==17
out=[-3*u,-6+9*u+6*v+9*w,-3*w];
elseif n==18
out=[3*u,-3+3*v+9*w,-6*w];
elseif n==19
out=[6-6*u-9*v-9*w,3*v,3*w];
elseif n==20
out=[3-3*u-9*w,-3*v,6*w];
end
end