Question

How to plot the magnetic flux density for the concentric hollow cylinders using matlab

Answer 1

% Magnetic Flux distribution using Finite Element Method % of a Solenoid vertically placed on XY plane." % Field distribution is calculated on XZ plane @ origin (Y-0).. % We will use Biot-Savart Law for this % Refer - http://en.wikipedia.org/wiki/Biot-Savart-law,.. clear all clc d-2; % Radius of the Loop." % Direction of current through the Loop 1(ANTI-CLOCK), -1(CLOCK) % Number sections/elements in the conductor % Number of Turns in the Solenoid." % Gap between Turns.. % Length of each element.. N 30 ; 8 0.75 dl-2*pi d/N; % XYZ Coordinates/Location of each element from the origin(0,0,0), Center of the Loop is taken as origin.. dtht 360/N; tht (0-dtht/2): dtht : (360-dtht/ 2); % Angle of each element w.r.to origin.. % Loop elements divided w.r.to degrees." XC= d."cosd (tht)."ones(1,N); yC-d. sind(tht).*ones (1,N); % X coordinate % Y coordinate

% zc will change for each turn, -z to +z, and will be varied during iteration 2C ones(1,h); h (g/2) (T-1) : g : (g/2) (T-1); % Length (Projection) & Direction of each current element in Vector form % Length of each element on X axis.. % Length of each element on Y axis.. % Length of each element is zero on z axis.. Lz -zeros (1,N); % Points/Locations in space (here XZ plane) where B is to be computed NP125; % Detector points.. % Dimensions of detector space.., arbitrary.. zPmax 1.5* (T*g); XP linspace(-xPmax,XPmax, NP); zP - linspace(-zPmax, zPmax, NP); [xXP zzP] meshgrid(xP,2P); % Divide space with NP points." % Creating the Mesh.. % Initialize B.. Zeros (NP,NP); By-zeros (NP, NP); B2 zeros (NP,NP);

% Computation of Magnetic Field (B) using superposition principle." % Compute B at each detector points due to each small cond elements & integrate them for p 1: T .. for q 1:N % Displacement Vector along X direction from Conductor.. % No detector points on Y direction.. % zc vertical location changes per Turn. r sqrt(rx."2+ry."2+rz."2); % Displacement Magnitude for an element on the conductor. rA3 rx ry rz I % m - direction of current element end

end B -B/max(max(B)); % Plotting % Magnitude of B % Normalizing % figure(1); % pcolor(XXP,zzPS); % colormap(jet); % shading interp; axis equal; % title( "Magnetic Field Distibution"); % colorbar; figure(2); surf(xxP,z2P,B, 'Facecolor', 'interp',... EdgeColor, 'none',... Facelighting', 'phong'); daspect([1 1 1]); axis tight; view(-10,30); camlight right;

colormap (jet); grid off axis off; colorbar; title('Magnetic Field Distibution) % figure (3); % quiver(xx?zzP,Bx, Bz); % colormap (lines); % axis([-2*d 2*d -T"g T"g]); % title('Magnetic Field Distibution'); % xlabel('<--x -->'');ylabel('<-- z -->'); % zoom on;