Discussion Forum

Dear All,

i am working on a MRI coil design . I want to extract the edge numbers of MRI Coils . I want to integrate (integration formula is different for each coil) magnetic due to each coil(there are around 56 coils).
I want to extract edge number of each coil , and enter it in the next optimization iteration.

i m using
clear edg
edg.ind = [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23, ...
24,25,26,26,27,25,28,24,29,23,30,22,31,21,32,20,33,19,34,18,35,17,36, ...
16,37,15,38,14,39,13,40,12,41,11,42,10,43,9,44,8,45,7,46,6,47,5,48,4, ...
49,3,50,2,51,1,52,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43, ...
44,45,46,47,48,49,50,51,52];

% Edge expressions
edg.expr = {'ut1',{'a*r26','a*r25','a*r24','a*r23','a*r22','a*r21','a*r20', ...
'a*r19','a*r18','a*r17','a*r16','a*r15','a*r14','a*r13','a*r12', ...
'a*r11','a*r10','a*r9','a*r8','a*r7','a*r6','a*r5','a*r4','a*r3', ...
'a*r2','a*r1','-a*r1','-a*r2','-a*r3','-a*r4','-a*r5','-a*r6','-a*r7','-a*r8', ...
'-a*r9','-a*r10','-a*r11','-a*r12','-a*r13','-a*r14','-a*r15','-a*r16','-a*r17', ...
'-a*r18','-a*r19','-a*r20','-a*r21','-a*r22','-a*r23','-a*r24','-a*r25','-a*r26'}};
fem.edg = edg;

% Scalar expressions
fem.expr = {'a','((mu*I)/(4*pi))*((t1x*(y-y0)-t1y*(x-x0))/(((x-x0)^2+(y-y0)^2+(z-z0)^2)^1.5))'};



commands to specify the expression on respective coil. and finally i m using "postint"
command to integrate

How could i extract edg.ind of the respective coil.

is there any other way to integrate magnetic field ?

thanks in advance