% COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) flclear fem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.4'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 248; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2007/10/10 16:07:51 $'; fem.version = vrsn; % Geometry g1=circ2(0.4,'base','center','pos',[-0.6,0]); g2=circ2('250e-6','base','center','pos',{'0','0'},'rot','0'); [g3]=geomcopy({g2}); [g4]=geomcopy({g3}); g4=move(g4,[0.0010,0]); [g5]=geomcopy({g3}); g5=move(g5,[-0.0010,0]); [g6,g7,g8]=geomcopy({g2,g4,g5}); [g9,g10,g11]=geomcopy({g6,g7,g8}); g9=move(g9,[0,-0.0010]); g10=move(g10,[0,-0.0010]); g11=move(g11,[0,-0.0010]); [g12,g13,g14,g15,g16,g17]=geomcopy({g2,g4,g5,g9,g10,g11}); [g18,g19,g20,g21,g22,g23]=geomcopy({g12,g13,g14,g15,g16,g17}); g18=move(g18,[0,0.0020]); g19=move(g19,[0,0.0020]); g20=move(g20,[0,0.0020]); g21=move(g21,[0,0.0020]); g22=move(g22,[0,0.0020]); g23=move(g23,[0,0.0020]); [g24,g25,g26,g27,g28,g29,g30,g31,g32,g33,g34,g35]=geomcopy({g2,g4,g5,g9,g10,g11,g18,g19,g20,g21,g22,g23}); [g36,g37,g38,g39,g40,g41,g42,g43,g44,g45,g46,g47]=geomcopy({g24,g25,g26,g27,g28,g29,g30,g31,g32,g33,g34,g35}); g36=move(g36,[0.0030,0]); g37=move(g37,[0.0030,0]); g38=move(g38,[0.0030,0]); g39=move(g39,[0.0030,0]); g40=move(g40,[0.0030,0]); g41=move(g41,[0.0030,0]); g42=move(g42,[0.0030,0]); g43=move(g43,[0.0030,0]); g44=move(g44,[0.0030,0]); g45=move(g45,[0.0030,0]); g46=move(g46,[0.0030,0]); g47=move(g47,[0.0030,0]); [g48,g49,g50,g51,g52,g53,g54,g55,g56,g57,g58,g59,g60,g61,g62,g63,g64,g65,g66,g67,g68,g69,g70,g71]=geomcopy({g2,g4,g5,g9,g10,g11,g18,g19,g20,g21,g22,g23,g36,g37,g38,g39,g40,g41,g42,g43,g44,g45,g46,g47}); [g72,g73,g74,g75,g76,g77,g78,g79,g80,g81,g82,g83,g84,g85,g86,g87,g88,g89,g90,g91,g92,g93,g94,g95]=geomcopy({g48,g49,g50,g51,g52,g53,g54,g55,g56,g57,g58,g59,g60,g61,g62,g63,g64,g65,g66,g67,g68,g69,g70,g71}); g72=move(g72,[0,-0.0040]); g73=move(g73,[0,-0.0040]); g74=move(g74,[0,-0.0040]); g75=move(g75,[0,-0.0040]); g76=move(g76,[0,-0.0040]); g77=move(g77,[0,-0.0040]); g78=move(g78,[0,-0.0040]); g79=move(g79,[0,-0.0040]); g80=move(g80,[0,-0.0040]); g81=move(g81,[0,-0.0040]); g82=move(g82,[0,-0.0040]); g83=move(g83,[0,-0.0040]); g84=move(g84,[0,-0.0040]); g85=move(g85,[0,-0.0040]); g86=move(g86,[0,-0.0040]); g87=move(g87,[0,-0.0040]); g88=move(g88,[0,-0.0040]); g89=move(g89,[0,-0.0040]); g90=move(g90,[0,-0.0040]); g91=move(g91,[0,-0.0040]); g92=move(g92,[0,-0.0040]); g93=move(g93,[0,-0.0040]); g94=move(g94,[0,-0.0040]); g95=move(g95,[0,-0.0040]); [g96,g97,g98,g99,g100,g101,g102,g103,g104,g105,g106,g107,g108,g109,g110,g111,g112,g113,g114,g115,g116,g117,g118,g119,g120,g121,g122,g123,g124,g125,g126,g127]=geomcopy({g4,g10,g19,g22,g36,g37,g38,g39,g40,g41,g42,g43,g44,g45,g46,g47,g73,g76,g79,g82,g84,g85,g86,g87,g88,g89,g90,g91,g92,g93,g94,g95}); [g128,g129,g130,g131,g132,g133,g134,g135,g136,g137,g138,g139,g140,g141,g142,g143,g144,g145,g146,g147,g148,g149,g150,g151,g152,g153,g154,g155,g156,g157,g158,g159]=geomcopy({g96,g97,g98,g99,g100,g101,g102,g103,g104,g105,g106,g107,g108,g109,g110,g111,g112,g113,g114,g115,g116,g117,g118,g119,g120,g121,g122,g123,g124,g125,g126,g127}); g128=move(g128,[0.0040,0]); g129=move(g129,[0.0040,0]); g130=move(g130,[0.0040,0]); g131=move(g131,[0.0040,0]); g132=move(g132,[0.0040,0]); g133=move(g133,[0.0040,0]); g134=move(g134,[0.0040,0]); g135=move(g135,[0.0040,0]); g136=move(g136,[0.0040,0]); g137=move(g137,[0.0040,0]); g138=move(g138,[0.0040,0]); g139=move(g139,[0.0040,0]); g140=move(g140,[0.0040,0]); g141=move(g141,[0.0040,0]); g142=move(g142,[0.0040,0]); g143=move(g143,[0.0040,0]); g144=move(g144,[0.0040,0]); g145=move(g145,[0.0040,0]); g146=move(g146,[0.0040,0]); g147=move(g147,[0.0040,0]); g148=move(g148,[0.0040,0]); g149=move(g149,[0.0040,0]); g150=move(g150,[0.0040,0]); g151=move(g151,[0.0040,0]); g152=move(g152,[0.0040,0]); g153=move(g153,[0.0040,0]); g154=move(g154,[0.0040,0]); g155=move(g155,[0.0040,0]); g156=move(g156,[0.0040,0]); g157=move(g157,[0.0040,0]); g158=move(g158,[0.0040,0]); g159=move(g159,[0.0040,0]); [g160,g161,g162,g163,g164,g165,g166,g167,g168,g169,g170,g171,g172,g173,g174,g175,g176,g177,g178,g179]=geomcopy({g72,g73,g74,g75,g76,g77,g84,g85,g86,g87,g88,g89,g144,g145,g148,g149,g150,g151,g152,g153}); [g180,g181,g182,g183,g184,g185,g186,g187,g188,g189,g190,g191,g192,g193,g194,g195,g196,g197,g198,g199]=geomcopy({g160,g161,g162,g163,g164,g165,g166,g167,g168,g169,g170,g171,g172,g173,g174,g175,g176,g177,g178,g179}); g180=move(g180,[0,-0.0020]); g181=move(g181,[0,-0.0020]); g182=move(g182,[0,-0.0020]); g183=move(g183,[0,-0.0020]); g184=move(g184,[0,-0.0020]); g185=move(g185,[0,-0.0020]); g186=move(g186,[0,-0.0020]); g187=move(g187,[0,-0.0020]); g188=move(g188,[0,-0.0020]); g189=move(g189,[0,-0.0020]); g190=move(g190,[0,-0.0020]); g191=move(g191,[0,-0.0020]); g192=move(g192,[0,-0.0020]); g193=move(g193,[0,-0.0020]); g194=move(g194,[0,-0.0020]); g195=move(g195,[0,-0.0020]); g196=move(g196,[0,-0.0020]); g197=move(g197,[0,-0.0020]); g198=move(g198,[0,-0.0020]); g199=move(g199,[0,-0.0020]); g2=move(g2,[-0.0050,0]); g4=move(g4,[-0.0050,0]); g5=move(g5,[-0.0050,0]); g9=move(g9,[-0.0050,0]); g10=move(g10,[-0.0050,0]); g11=move(g11,[-0.0050,0]); g18=move(g18,[-0.0050,0]); g19=move(g19,[-0.0050,0]); g20=move(g20,[-0.0050,0]); g21=move(g21,[-0.0050,0]); g22=move(g22,[-0.0050,0]); g23=move(g23,[-0.0050,0]); g36=move(g36,[-0.0050,0]); g37=move(g37,[-0.0050,0]); g38=move(g38,[-0.0050,0]); g39=move(g39,[-0.0050,0]); g40=move(g40,[-0.0050,0]); g41=move(g41,[-0.0050,0]); g42=move(g42,[-0.0050,0]); g43=move(g43,[-0.0050,0]); g44=move(g44,[-0.0050,0]); g45=move(g45,[-0.0050,0]); g46=move(g46,[-0.0050,0]); g47=move(g47,[-0.0050,0]); g72=move(g72,[-0.0050,0]); g73=move(g73,[-0.0050,0]); g74=move(g74,[-0.0050,0]); g75=move(g75,[-0.0050,0]); g76=move(g76,[-0.0050,0]); g77=move(g77,[-0.0050,0]); g78=move(g78,[-0.0050,0]); g79=move(g79,[-0.0050,0]); g80=move(g80,[-0.0050,0]); g81=move(g81,[-0.0050,0]); g82=move(g82,[-0.0050,0]); g83=move(g83,[-0.0050,0]); g84=move(g84,[-0.0050,0]); g85=move(g85,[-0.0050,0]); g86=move(g86,[-0.0050,0]); g87=move(g87,[-0.0050,0]); g88=move(g88,[-0.0050,0]); g89=move(g89,[-0.0050,0]); g90=move(g90,[-0.0050,0]); g91=move(g91,[-0.0050,0]); g92=move(g92,[-0.0050,0]); g93=move(g93,[-0.0050,0]); g94=move(g94,[-0.0050,0]); g95=move(g95,[-0.0050,0]); g128=move(g128,[-0.0050,0]); g129=move(g129,[-0.0050,0]); g130=move(g130,[-0.0050,0]); g131=move(g131,[-0.0050,0]); g132=move(g132,[-0.0050,0]); g133=move(g133,[-0.0050,0]); g134=move(g134,[-0.0050,0]); g135=move(g135,[-0.0050,0]); g136=move(g136,[-0.0050,0]); g137=move(g137,[-0.0050,0]); g138=move(g138,[-0.0050,0]); g139=move(g139,[-0.0050,0]); g140=move(g140,[-0.0050,0]); g141=move(g141,[-0.0050,0]); g142=move(g142,[-0.0050,0]); g143=move(g143,[-0.0050,0]); g144=move(g144,[-0.0050,0]); g145=move(g145,[-0.0050,0]); g146=move(g146,[-0.0050,0]); g147=move(g147,[-0.0050,0]); g148=move(g148,[-0.0050,0]); g149=move(g149,[-0.0050,0]); g150=move(g150,[-0.0050,0]); g151=move(g151,[-0.0050,0]); g152=move(g152,[-0.0050,0]); g153=move(g153,[-0.0050,0]); g154=move(g154,[-0.0050,0]); g155=move(g155,[-0.0050,0]); g156=move(g156,[-0.0050,0]); g157=move(g157,[-0.0050,0]); g158=move(g158,[-0.0050,0]); g159=move(g159,[-0.0050,0]); g180=move(g180,[-0.0050,0]); g181=move(g181,[-0.0050,0]); g182=move(g182,[-0.0050,0]); g183=move(g183,[-0.0050,0]); g184=move(g184,[-0.0050,0]); g185=move(g185,[-0.0050,0]); g186=move(g186,[-0.0050,0]); g187=move(g187,[-0.0050,0]); g188=move(g188,[-0.0050,0]); g189=move(g189,[-0.0050,0]); g190=move(g190,[-0.0050,0]); g191=move(g191,[-0.0050,0]); g192=move(g192,[-0.0050,0]); g193=move(g193,[-0.0050,0]); g194=move(g194,[-0.0050,0]); g195=move(g195,[-0.0050,0]); g196=move(g196,[-0.0050,0]); g197=move(g197,[-0.0050,0]); g198=move(g198,[-0.0050,0]); g199=move(g199,[-0.0050,0]); g2=move(g2,[0.0030,0.0020]); g4=move(g4,[0.0030,0.0020]); g5=move(g5,[0.0030,0.0020]); g9=move(g9,[0.0030,0.0020]); g10=move(g10,[0.0030,0.0020]); g11=move(g11,[0.0030,0.0020]); g18=move(g18,[0.0030,0.0020]); g19=move(g19,[0.0030,0.0020]); g20=move(g20,[0.0030,0.0020]); g21=move(g21,[0.0030,0.0020]); g22=move(g22,[0.0030,0.0020]); g23=move(g23,[0.0030,0.0020]); g36=move(g36,[0.0030,0.0020]); g37=move(g37,[0.0030,0.0020]); g38=move(g38,[0.0030,0.0020]); g39=move(g39,[0.0030,0.0020]); g40=move(g40,[0.0030,0.0020]); g41=move(g41,[0.0030,0.0020]); g42=move(g42,[0.0030,0.0020]); g43=move(g43,[0.0030,0.0020]); g44=move(g44,[0.0030,0.0020]); g45=move(g45,[0.0030,0.0020]); g46=move(g46,[0.0030,0.0020]); g47=move(g47,[0.0030,0.0020]); g72=move(g72,[0.0030,0.0020]); g73=move(g73,[0.0030,0.0020]); g74=move(g74,[0.0030,0.0020]); g75=move(g75,[0.0030,0.0020]); g76=move(g76,[0.0030,0.0020]); g77=move(g77,[0.0030,0.0020]); g78=move(g78,[0.0030,0.0020]); g79=move(g79,[0.0030,0.0020]); g80=move(g80,[0.0030,0.0020]); g81=move(g81,[0.0030,0.0020]); g82=move(g82,[0.0030,0.0020]); g83=move(g83,[0.0030,0.0020]); g84=move(g84,[0.0030,0.0020]); g85=move(g85,[0.0030,0.0020]); g86=move(g86,[0.0030,0.0020]); g87=move(g87,[0.0030,0.0020]); g88=move(g88,[0.0030,0.0020]); g89=move(g89,[0.0030,0.0020]); g90=move(g90,[0.0030,0.0020]); g91=move(g91,[0.0030,0.0020]); g92=move(g92,[0.0030,0.0020]); g93=move(g93,[0.0030,0.0020]); g94=move(g94,[0.0030,0.0020]); g95=move(g95,[0.0030,0.0020]); g128=move(g128,[0.0030,0.0020]); g129=move(g129,[0.0030,0.0020]); g130=move(g130,[0.0030,0.0020]); g131=move(g131,[0.0030,0.0020]); g132=move(g132,[0.0030,0.0020]); g133=move(g133,[0.0030,0.0020]); g134=move(g134,[0.0030,0.0020]); g135=move(g135,[0.0030,0.0020]); g136=move(g136,[0.0030,0.0020]); g137=move(g137,[0.0030,0.0020]); g138=move(g138,[0.0030,0.0020]); g139=move(g139,[0.0030,0.0020]); g140=move(g140,[0.0030,0.0020]); g141=move(g141,[0.0030,0.0020]); g142=move(g142,[0.0030,0.0020]); g143=move(g143,[0.0030,0.0020]); g144=move(g144,[0.0030,0.0020]); g145=move(g145,[0.0030,0.0020]); g146=move(g146,[0.0030,0.0020]); g147=move(g147,[0.0030,0.0020]); g148=move(g148,[0.0030,0.0020]); g149=move(g149,[0.0030,0.0020]); g150=move(g150,[0.0030,0.0020]); g151=move(g151,[0.0030,0.0020]); g152=move(g152,[0.0030,0.0020]); g153=move(g153,[0.0030,0.0020]); g154=move(g154,[0.0030,0.0020]); g155=move(g155,[0.0030,0.0020]); g156=move(g156,[0.0030,0.0020]); g157=move(g157,[0.0030,0.0020]); g158=move(g158,[0.0030,0.0020]); g159=move(g159,[0.0030,0.0020]); g180=move(g180,[0.0030,0.0020]); g181=move(g181,[0.0030,0.0020]); g182=move(g182,[0.0030,0.0020]); g183=move(g183,[0.0030,0.0020]); g184=move(g184,[0.0030,0.0020]); g185=move(g185,[0.0030,0.0020]); g186=move(g186,[0.0030,0.0020]); g187=move(g187,[0.0030,0.0020]); g188=move(g188,[0.0030,0.0020]); g189=move(g189,[0.0030,0.0020]); g190=move(g190,[0.0030,0.0020]); g191=move(g191,[0.0030,0.0020]); g192=move(g192,[0.0030,0.0020]); g193=move(g193,[0.0030,0.0020]); g194=move(g194,[0.0030,0.0020]); g195=move(g195,[0.0030,0.0020]); g196=move(g196,[0.0030,0.0020]); g197=move(g197,[0.0030,0.0020]); g198=move(g198,[0.0030,0.0020]); g199=move(g199,[0.0030,0.0020]); g2=move(g2,[-0.0010,-0.0010]); g4=move(g4,[-0.0010,-0.0010]); g5=move(g5,[-0.0010,-0.0010]); g9=move(g9,[-0.0010,-0.0010]); g10=move(g10,[-0.0010,-0.0010]); g11=move(g11,[-0.0010,-0.0010]); g18=move(g18,[-0.0010,-0.0010]); g19=move(g19,[-0.0010,-0.0010]); g20=move(g20,[-0.0010,-0.0010]); g21=move(g21,[-0.0010,-0.0010]); g22=move(g22,[-0.0010,-0.0010]); g23=move(g23,[-0.0010,-0.0010]); g36=move(g36,[-0.0010,-0.0010]); g37=move(g37,[-0.0010,-0.0010]); g38=move(g38,[-0.0010,-0.0010]); g39=move(g39,[-0.0010,-0.0010]); g40=move(g40,[-0.0010,-0.0010]); g41=move(g41,[-0.0010,-0.0010]); g42=move(g42,[-0.0010,-0.0010]); g43=move(g43,[-0.0010,-0.0010]); g44=move(g44,[-0.0010,-0.0010]); g45=move(g45,[-0.0010,-0.0010]); g46=move(g46,[-0.0010,-0.0010]); g47=move(g47,[-0.0010,-0.0010]); g72=move(g72,[-0.0010,-0.0010]); g73=move(g73,[-0.0010,-0.0010]); g74=move(g74,[-0.0010,-0.0010]); g75=move(g75,[-0.0010,-0.0010]); g76=move(g76,[-0.0010,-0.0010]); g77=move(g77,[-0.0010,-0.0010]); g78=move(g78,[-0.0010,-0.0010]); g79=move(g79,[-0.0010,-0.0010]); g80=move(g80,[-0.0010,-0.0010]); g81=move(g81,[-0.0010,-0.0010]); g82=move(g82,[-0.0010,-0.0010]); g83=move(g83,[-0.0010,-0.0010]); g84=move(g84,[-0.0010,-0.0010]); g85=move(g85,[-0.0010,-0.0010]); g86=move(g86,[-0.0010,-0.0010]); g87=move(g87,[-0.0010,-0.0010]); g88=move(g88,[-0.0010,-0.0010]); g89=move(g89,[-0.0010,-0.0010]); g90=move(g90,[-0.0010,-0.0010]); g91=move(g91,[-0.0010,-0.0010]); g92=move(g92,[-0.0010,-0.0010]); g93=move(g93,[-0.0010,-0.0010]); g94=move(g94,[-0.0010,-0.0010]); g95=move(g95,[-0.0010,-0.0010]); g128=move(g128,[-0.0010,-0.0010]); g129=move(g129,[-0.0010,-0.0010]); g130=move(g130,[-0.0010,-0.0010]); g131=move(g131,[-0.0010,-0.0010]); g132=move(g132,[-0.0010,-0.0010]); g133=move(g133,[-0.0010,-0.0010]); g134=move(g134,[-0.0010,-0.0010]); g135=move(g135,[-0.0010,-0.0010]); g136=move(g136,[-0.0010,-0.0010]); g137=move(g137,[-0.0010,-0.0010]); g138=move(g138,[-0.0010,-0.0010]); g139=move(g139,[-0.0010,-0.0010]); g140=move(g140,[-0.0010,-0.0010]); g141=move(g141,[-0.0010,-0.0010]); g142=move(g142,[-0.0010,-0.0010]); g143=move(g143,[-0.0010,-0.0010]); g144=move(g144,[-0.0010,-0.0010]); g145=move(g145,[-0.0010,-0.0010]); g146=move(g146,[-0.0010,-0.0010]); g147=move(g147,[-0.0010,-0.0010]); g148=move(g148,[-0.0010,-0.0010]); g149=move(g149,[-0.0010,-0.0010]); g150=move(g150,[-0.0010,-0.0010]); g151=move(g151,[-0.0010,-0.0010]); g152=move(g152,[-0.0010,-0.0010]); g153=move(g153,[-0.0010,-0.0010]); g154=move(g154,[-0.0010,-0.0010]); g155=move(g155,[-0.0010,-0.0010]); g156=move(g156,[-0.0010,-0.0010]); g157=move(g157,[-0.0010,-0.0010]); g158=move(g158,[-0.0010,-0.0010]); g159=move(g159,[-0.0010,-0.0010]); g180=move(g180,[-0.0010,-0.0010]); g181=move(g181,[-0.0010,-0.0010]); g182=move(g182,[-0.0010,-0.0010]); g183=move(g183,[-0.0010,-0.0010]); g184=move(g184,[-0.0010,-0.0010]); g185=move(g185,[-0.0010,-0.0010]); g186=move(g186,[-0.0010,-0.0010]); g187=move(g187,[-0.0010,-0.0010]); g188=move(g188,[-0.0010,-0.0010]); g189=move(g189,[-0.0010,-0.0010]); g190=move(g190,[-0.0010,-0.0010]); g191=move(g191,[-0.0010,-0.0010]); g192=move(g192,[-0.0010,-0.0010]); g193=move(g193,[-0.0010,-0.0010]); g194=move(g194,[-0.0010,-0.0010]); g195=move(g195,[-0.0010,-0.0010]); g196=move(g196,[-0.0010,-0.0010]); g197=move(g197,[-0.0010,-0.0010]); g198=move(g198,[-0.0010,-0.0010]); g199=move(g199,[-0.0010,-0.0010]); g2=move(g2,[0,0.0010]); g4=move(g4,[0,0.0010]); g5=move(g5,[0,0.0010]); g9=move(g9,[0,0.0010]); g10=move(g10,[0,0.0010]); g11=move(g11,[0,0.0010]); g18=move(g18,[0,0.0010]); g19=move(g19,[0,0.0010]); g20=move(g20,[0,0.0010]); g21=move(g21,[0,0.0010]); g22=move(g22,[0,0.0010]); g23=move(g23,[0,0.0010]); g36=move(g36,[0,0.0010]); g37=move(g37,[0,0.0010]); g38=move(g38,[0,0.0010]); g39=move(g39,[0,0.0010]); g40=move(g40,[0,0.0010]); g41=move(g41,[0,0.0010]); g42=move(g42,[0,0.0010]); g43=move(g43,[0,0.0010]); g44=move(g44,[0,0.0010]); g45=move(g45,[0,0.0010]); g46=move(g46,[0,0.0010]); g47=move(g47,[0,0.0010]); g72=move(g72,[0,0.0010]); g73=move(g73,[0,0.0010]); g74=move(g74,[0,0.0010]); g75=move(g75,[0,0.0010]); g76=move(g76,[0,0.0010]); g77=move(g77,[0,0.0010]); g78=move(g78,[0,0.0010]); g79=move(g79,[0,0.0010]); g80=move(g80,[0,0.0010]); g81=move(g81,[0,0.0010]); g82=move(g82,[0,0.0010]); g83=move(g83,[0,0.0010]); g84=move(g84,[0,0.0010]); g85=move(g85,[0,0.0010]); g86=move(g86,[0,0.0010]); g87=move(g87,[0,0.0010]); g88=move(g88,[0,0.0010]); g89=move(g89,[0,0.0010]); g90=move(g90,[0,0.0010]); g91=move(g91,[0,0.0010]); g92=move(g92,[0,0.0010]); g93=move(g93,[0,0.0010]); g94=move(g94,[0,0.0010]); g95=move(g95,[0,0.0010]); g128=move(g128,[0,0.0010]); g129=move(g129,[0,0.0010]); g130=move(g130,[0,0.0010]); g131=move(g131,[0,0.0010]); g132=move(g132,[0,0.0010]); g133=move(g133,[0,0.0010]); g134=move(g134,[0,0.0010]); g135=move(g135,[0,0.0010]); g136=move(g136,[0,0.0010]); g137=move(g137,[0,0.0010]); g138=move(g138,[0,0.0010]); g139=move(g139,[0,0.0010]); g140=move(g140,[0,0.0010]); g141=move(g141,[0,0.0010]); g142=move(g142,[0,0.0010]); g143=move(g143,[0,0.0010]); g144=move(g144,[0,0.0010]); g145=move(g145,[0,0.0010]); g146=move(g146,[0,0.0010]); g147=move(g147,[0,0.0010]); g148=move(g148,[0,0.0010]); g149=move(g149,[0,0.0010]); g150=move(g150,[0,0.0010]); g151=move(g151,[0,0.0010]); g152=move(g152,[0,0.0010]); g153=move(g153,[0,0.0010]); g154=move(g154,[0,0.0010]); g155=move(g155,[0,0.0010]); g156=move(g156,[0,0.0010]); g157=move(g157,[0,0.0010]); g158=move(g158,[0,0.0010]); g159=move(g159,[0,0.0010]); g180=move(g180,[0,0.0010]); g181=move(g181,[0,0.0010]); g182=move(g182,[0,0.0010]); g183=move(g183,[0,0.0010]); g184=move(g184,[0,0.0010]); g185=move(g185,[0,0.0010]); g186=move(g186,[0,0.0010]); g187=move(g187,[0,0.0010]); g188=move(g188,[0,0.0010]); g189=move(g189,[0,0.0010]); g190=move(g190,[0,0.0010]); g191=move(g191,[0,0.0010]); g192=move(g192,[0,0.0010]); g193=move(g193,[0,0.0010]); g194=move(g194,[0,0.0010]); g195=move(g195,[0,0.0010]); g196=move(g196,[0,0.0010]); g197=move(g197,[0,0.0010]); g198=move(g198,[0,0.0010]); g199=move(g199,[0,0.0010]); g2=move(g2,[0,5.0E-4]); g4=move(g4,[0,5.0E-4]); g5=move(g5,[0,5.0E-4]); g9=move(g9,[0,5.0E-4]); g10=move(g10,[0,5.0E-4]); g11=move(g11,[0,5.0E-4]); g18=move(g18,[0,5.0E-4]); g19=move(g19,[0,5.0E-4]); g20=move(g20,[0,5.0E-4]); g21=move(g21,[0,5.0E-4]); g22=move(g22,[0,5.0E-4]); g23=move(g23,[0,5.0E-4]); g36=move(g36,[0,5.0E-4]); g37=move(g37,[0,5.0E-4]); g38=move(g38,[0,5.0E-4]); g39=move(g39,[0,5.0E-4]); g40=move(g40,[0,5.0E-4]); g41=move(g41,[0,5.0E-4]); g42=move(g42,[0,5.0E-4]); g43=move(g43,[0,5.0E-4]); g44=move(g44,[0,5.0E-4]); g45=move(g45,[0,5.0E-4]); g46=move(g46,[0,5.0E-4]); g47=move(g47,[0,5.0E-4]); g72=move(g72,[0,5.0E-4]); g73=move(g73,[0,5.0E-4]); g74=move(g74,[0,5.0E-4]); g75=move(g75,[0,5.0E-4]); g76=move(g76,[0,5.0E-4]); g77=move(g77,[0,5.0E-4]); g78=move(g78,[0,5.0E-4]); g79=move(g79,[0,5.0E-4]); g80=move(g80,[0,5.0E-4]); g81=move(g81,[0,5.0E-4]); g82=move(g82,[0,5.0E-4]); g83=move(g83,[0,5.0E-4]); g84=move(g84,[0,5.0E-4]); g85=move(g85,[0,5.0E-4]); g86=move(g86,[0,5.0E-4]); g87=move(g87,[0,5.0E-4]); g88=move(g88,[0,5.0E-4]); g89=move(g89,[0,5.0E-4]); g90=move(g90,[0,5.0E-4]); g91=move(g91,[0,5.0E-4]); g92=move(g92,[0,5.0E-4]); g93=move(g93,[0,5.0E-4]); g94=move(g94,[0,5.0E-4]); g95=move(g95,[0,5.0E-4]); g128=move(g128,[0,5.0E-4]); g129=move(g129,[0,5.0E-4]); g130=move(g130,[0,5.0E-4]); g131=move(g131,[0,5.0E-4]); g132=move(g132,[0,5.0E-4]); g133=move(g133,[0,5.0E-4]); g134=move(g134,[0,5.0E-4]); g135=move(g135,[0,5.0E-4]); g136=move(g136,[0,5.0E-4]); g137=move(g137,[0,5.0E-4]); g138=move(g138,[0,5.0E-4]); g139=move(g139,[0,5.0E-4]); g140=move(g140,[0,5.0E-4]); g141=move(g141,[0,5.0E-4]); g142=move(g142,[0,5.0E-4]); g143=move(g143,[0,5.0E-4]); g144=move(g144,[0,5.0E-4]); g145=move(g145,[0,5.0E-4]); g146=move(g146,[0,5.0E-4]); g147=move(g147,[0,5.0E-4]); g148=move(g148,[0,5.0E-4]); g149=move(g149,[0,5.0E-4]); g150=move(g150,[0,5.0E-4]); g151=move(g151,[0,5.0E-4]); g152=move(g152,[0,5.0E-4]); g153=move(g153,[0,5.0E-4]); g154=move(g154,[0,5.0E-4]); g155=move(g155,[0,5.0E-4]); g156=move(g156,[0,5.0E-4]); g157=move(g157,[0,5.0E-4]); g158=move(g158,[0,5.0E-4]); g159=move(g159,[0,5.0E-4]); g180=move(g180,[0,5.0E-4]); g181=move(g181,[0,5.0E-4]); g182=move(g182,[0,5.0E-4]); g183=move(g183,[0,5.0E-4]); g184=move(g184,[0,5.0E-4]); g185=move(g185,[0,5.0E-4]); g186=move(g186,[0,5.0E-4]); g187=move(g187,[0,5.0E-4]); g188=move(g188,[0,5.0E-4]); g189=move(g189,[0,5.0E-4]); g190=move(g190,[0,5.0E-4]); g191=move(g191,[0,5.0E-4]); g192=move(g192,[0,5.0E-4]); g193=move(g193,[0,5.0E-4]); g194=move(g194,[0,5.0E-4]); g195=move(g195,[0,5.0E-4]); g196=move(g196,[0,5.0E-4]); g197=move(g197,[0,5.0E-4]); g198=move(g198,[0,5.0E-4]); g199=move(g199,[0,5.0E-4]); g2=move(g2,[-5.0E-4,0]); g4=move(g4,[-5.0E-4,0]); g5=move(g5,[-5.0E-4,0]); g9=move(g9,[-5.0E-4,0]); g10=move(g10,[-5.0E-4,0]); g11=move(g11,[-5.0E-4,0]); g18=move(g18,[-5.0E-4,0]); g19=move(g19,[-5.0E-4,0]); g20=move(g20,[-5.0E-4,0]); g21=move(g21,[-5.0E-4,0]); g22=move(g22,[-5.0E-4,0]); g23=move(g23,[-5.0E-4,0]); g36=move(g36,[-5.0E-4,0]); g37=move(g37,[-5.0E-4,0]); g38=move(g38,[-5.0E-4,0]); g39=move(g39,[-5.0E-4,0]); g40=move(g40,[-5.0E-4,0]); g41=move(g41,[-5.0E-4,0]); g42=move(g42,[-5.0E-4,0]); g43=move(g43,[-5.0E-4,0]); g44=move(g44,[-5.0E-4,0]); g45=move(g45,[-5.0E-4,0]); g46=move(g46,[-5.0E-4,0]); g47=move(g47,[-5.0E-4,0]); g72=move(g72,[-5.0E-4,0]); g73=move(g73,[-5.0E-4,0]); g74=move(g74,[-5.0E-4,0]); g75=move(g75,[-5.0E-4,0]); g76=move(g76,[-5.0E-4,0]); g77=move(g77,[-5.0E-4,0]); g78=move(g78,[-5.0E-4,0]); g79=move(g79,[-5.0E-4,0]); g80=move(g80,[-5.0E-4,0]); g81=move(g81,[-5.0E-4,0]); g82=move(g82,[-5.0E-4,0]); g83=move(g83,[-5.0E-4,0]); g84=move(g84,[-5.0E-4,0]); g85=move(g85,[-5.0E-4,0]); g86=move(g86,[-5.0E-4,0]); g87=move(g87,[-5.0E-4,0]); g88=move(g88,[-5.0E-4,0]); g89=move(g89,[-5.0E-4,0]); g90=move(g90,[-5.0E-4,0]); g91=move(g91,[-5.0E-4,0]); g92=move(g92,[-5.0E-4,0]); g93=move(g93,[-5.0E-4,0]); g94=move(g94,[-5.0E-4,0]); g95=move(g95,[-5.0E-4,0]); g128=move(g128,[-5.0E-4,0]); g129=move(g129,[-5.0E-4,0]); g130=move(g130,[-5.0E-4,0]); g131=move(g131,[-5.0E-4,0]); g132=move(g132,[-5.0E-4,0]); g133=move(g133,[-5.0E-4,0]); g134=move(g134,[-5.0E-4,0]); g135=move(g135,[-5.0E-4,0]); g136=move(g136,[-5.0E-4,0]); g137=move(g137,[-5.0E-4,0]); g138=move(g138,[-5.0E-4,0]); g139=move(g139,[-5.0E-4,0]); g140=move(g140,[-5.0E-4,0]); g141=move(g141,[-5.0E-4,0]); g142=move(g142,[-5.0E-4,0]); g143=move(g143,[-5.0E-4,0]); g144=move(g144,[-5.0E-4,0]); g145=move(g145,[-5.0E-4,0]); g146=move(g146,[-5.0E-4,0]); g147=move(g147,[-5.0E-4,0]); g148=move(g148,[-5.0E-4,0]); g149=move(g149,[-5.0E-4,0]); g150=move(g150,[-5.0E-4,0]); g151=move(g151,[-5.0E-4,0]); g152=move(g152,[-5.0E-4,0]); g153=move(g153,[-5.0E-4,0]); g154=move(g154,[-5.0E-4,0]); g155=move(g155,[-5.0E-4,0]); g156=move(g156,[-5.0E-4,0]); g157=move(g157,[-5.0E-4,0]); g158=move(g158,[-5.0E-4,0]); g159=move(g159,[-5.0E-4,0]); g180=move(g180,[-5.0E-4,0]); g181=move(g181,[-5.0E-4,0]); g182=move(g182,[-5.0E-4,0]); g183=move(g183,[-5.0E-4,0]); g184=move(g184,[-5.0E-4,0]); g185=move(g185,[-5.0E-4,0]); g186=move(g186,[-5.0E-4,0]); g187=move(g187,[-5.0E-4,0]); g188=move(g188,[-5.0E-4,0]); g189=move(g189,[-5.0E-4,0]); g190=move(g190,[-5.0E-4,0]); g191=move(g191,[-5.0E-4,0]); g192=move(g192,[-5.0E-4,0]); g193=move(g193,[-5.0E-4,0]); g194=move(g194,[-5.0E-4,0]); g195=move(g195,[-5.0E-4,0]); g196=move(g196,[-5.0E-4,0]); g197=move(g197,[-5.0E-4,0]); g198=move(g198,[-5.0E-4,0]); g199=move(g199,[-5.0E-4,0]); g200=rect2(0.01,0.01,'base','center','pos',[-0.02,0]); g201=rect2('.001','.001','base','center','pos',{'0','0'},'rot','0'); g202=rect2('.01','.01','base','center','pos',{'0','0'},'rot','0'); g203=rect2('.02','.02','base','center','pos',{'0','0'},'rot','0'); g204=rect2('.01','.01','base','center','pos',{'0','0'},'rot','0'); g205=rect2('.03','.03','base','center','pos',{'0','0'},'rot','0'); % Analyzed geometry clear s s.objs={g2,g4,g5,g9,g10,g11,g18,g19,g20,g21,g22,g23,g36,g37,g38,g39, ... g40,g41,g42,g43,g44,g45,g46,g47,g72,g73,g74,g75,g76,g77,g78,g79,g80, ... g81,g82,g83,g84,g85,g86,g87,g88,g89,g90,g91,g92,g93,g94,g95,g128,g129, ... g130,g131,g132,g133,g134,g135,g136,g137,g138,g139,g140,g141,g142,g143, ... g144,g145,g146,g147,g148,g149,g150,g151,g152,g153,g154,g155,g156,g157, ... g158,g159,g180,g181,g182,g183,g184,g185,g186,g187,g188,g189,g190,g191, ... g192,g193,g194,g195,g196,g197,g198,g199,g205}; s.name={'C1','C2','C3','C4','C5','C6','C7','C8','C9','C10','C11', ... 'C12','C13','C14','C15','C16','C17','C18','C19','C20','C21','C22','C23', ... 'C24','C25','C26','C27','C28','C29','C30','C31','C32','C33','C34','C35', ... 'C36','C37','C38','C39','C40','C41','C42','C43','C44','C45','C46','C47', ... 'C48','C49','C50','C51','C52','C53','C54','C55','C56','C57','C58','C59', ... 'C60','C61','C62','C63','C64','C65','C66','C67','C68','C69','C70','C71', ... 'C72','C73','C74','C75','C76','C77','C78','C79','C80','C81','C82','C83', ... 'C84','C85','C86','C87','C88','C89','C90','C91','C92','C93','C94','C95', ... 'C96','C97','C98','C99','C100','R1'}; s.tags={'g2','g4','g5','g9','g10','g11','g18','g19','g20','g21', ... 'g22','g23','g36','g37','g38','g39','g40','g41','g42','g43','g44','g45', ... 'g46','g47','g72','g73','g74','g75','g76','g77','g78','g79','g80','g81', ... 'g82','g83','g84','g85','g86','g87','g88','g89','g90','g91','g92','g93', ... 'g94','g95','g128','g129','g130','g131','g132','g133','g134','g135', ... 'g136','g137','g138','g139','g140','g141','g142','g143','g144','g145', ... 'g146','g147','g148','g149','g150','g151','g152','g153','g154','g155', ... 'g156','g157','g158','g159','g180','g181','g182','g183','g184','g185', ... 'g186','g187','g188','g189','g190','g191','g192','g193','g194','g195', ... 'g196','g197','g198','g199','g205'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); g205=scale(g205,1,.5,0,0); % Analyzed geometry clear s s.objs={g2,g4,g5,g9,g10,g11,g18,g19,g20,g21,g22,g23,g36,g37,g38,g39, ... g40,g41,g42,g43,g44,g45,g46,g47,g72,g73,g74,g75,g76,g77,g78,g79,g80, ... g81,g82,g83,g84,g85,g86,g87,g88,g89,g90,g91,g92,g93,g94,g95,g128,g129, ... g130,g131,g132,g133,g134,g135,g136,g137,g138,g139,g140,g141,g142,g143, ... g144,g145,g146,g147,g148,g149,g150,g151,g152,g153,g154,g155,g156,g157, ... g158,g159,g180,g181,g182,g183,g184,g185,g186,g187,g188,g189,g190,g191, ... g192,g193,g194,g195,g196,g197,g198,g199,g205}; s.name={'C1','C2','C3','C4','C5','C6','C7','C8','C9','C10','C11', ... 'C12','C13','C14','C15','C16','C17','C18','C19','C20','C21','C22','C23', ... 'C24','C25','C26','C27','C28','C29','C30','C31','C32','C33','C34','C35', ... 'C36','C37','C38','C39','C40','C41','C42','C43','C44','C45','C46','C47', ... 'C48','C49','C50','C51','C52','C53','C54','C55','C56','C57','C58','C59', ... 'C60','C61','C62','C63','C64','C65','C66','C67','C68','C69','C70','C71', ... 'C72','C73','C74','C75','C76','C77','C78','C79','C80','C81','C82','C83', ... 'C84','C85','C86','C87','C88','C89','C90','C91','C92','C93','C94','C95', ... 'C96','C97','C98','C99','C100','R1'}; s.tags={'g2','g4','g5','g9','g10','g11','g18','g19','g20','g21', ... 'g22','g23','g36','g37','g38','g39','g40','g41','g42','g43','g44','g45', ... 'g46','g47','g72','g73','g74','g75','g76','g77','g78','g79','g80','g81', ... 'g82','g83','g84','g85','g86','g87','g88','g89','g90','g91','g92','g93', ... 'g94','g95','g128','g129','g130','g131','g132','g133','g134','g135', ... 'g136','g137','g138','g139','g140','g141','g142','g143','g144','g145', ... 'g146','g147','g148','g149','g150','g151','g152','g153','g154','g155', ... 'g156','g157','g158','g159','g180','g181','g182','g183','g184','g185', ... 'g186','g187','g188','g189','g190','g191','g192','g193','g194','g195', ... 'g196','g197','g198','g199','g205'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); g205=scale(g205,1,2,0,0); % Geometry objects clear s s.objs={g2,g4,g5,g9,g10,g11,g18,g19,g20,g21,g22,g23,g36,g37,g38,g39, ... g40,g41,g42,g43,g44,g45,g46,g47,g72,g73,g74,g75,g76,g77,g78,g79,g80, ... g81,g82,g83,g84,g85,g86,g87,g88,g89,g90,g91,g92,g93,g94,g95,g128,g129, ... g130,g131,g132,g133,g134,g135,g136,g137,g138,g139,g140,g141,g142,g143, ... g144,g145,g146,g147,g148,g149,g150,g151,g152,g153,g154,g155,g156,g157, ... g158,g159,g180,g181,g182,g183,g184,g185,g186,g187,g188,g189,g190,g191, ... g192,g193,g194,g195,g196,g197,g198,g199,g205}; s.name={'C1','C2','C3','C4','C5','C6','C7','C8','C9','C10','C11', ... 'C12','C13','C14','C15','C16','C17','C18','C19','C20','C21','C22','C23', ... 'C24','C25','C26','C27','C28','C29','C30','C31','C32','C33','C34','C35', ... 'C36','C37','C38','C39','C40','C41','C42','C43','C44','C45','C46','C47', ... 'C48','C49','C50','C51','C52','C53','C54','C55','C56','C57','C58','C59', ... 'C60','C61','C62','C63','C64','C65','C66','C67','C68','C69','C70','C71', ... 'C72','C73','C74','C75','C76','C77','C78','C79','C80','C81','C82','C83', ... 'C84','C85','C86','C87','C88','C89','C90','C91','C92','C93','C94','C95', ... 'C96','C97','C98','C99','C100','R1'}; s.tags={'g2','g4','g5','g9','g10','g11','g18','g19','g20','g21', ... 'g22','g23','g36','g37','g38','g39','g40','g41','g42','g43','g44','g45', ... 'g46','g47','g72','g73','g74','g75','g76','g77','g78','g79','g80','g81', ... 'g82','g83','g84','g85','g86','g87','g88','g89','g90','g91','g92','g93', ... 'g94','g95','g128','g129','g130','g131','g132','g133','g134','g135', ... 'g136','g137','g138','g139','g140','g141','g142','g143','g144','g145', ... 'g146','g147','g148','g149','g150','g151','g152','g153','g154','g155', ... 'g156','g157','g158','g159','g180','g181','g182','g183','g184','g185', ... 'g186','g187','g188','g189','g190','g191','g192','g193','g194','g195', ... 'g196','g197','g198','g199','g205'}; fem.draw=struct('s',s); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry % Analyzed geometry clear s s.objs={g193,g192,g149,g22,g19,g74,g151,g132,g84,g184,g135,g152, ... g146,g83,g72,g183,g153,g80,g134,g81,g191,g36,g141,g181,g93,g20,g144, ... g47,g77,g23,g198,g73,g190,g199,g194,g145,g75,g131,g46,g78,g89,g82,g133, ... g150,g85,g41,g143,g86,g42,g128,g157,g44,g140,g196,g88,g9,g195,g40,g76, ... g130,g94,g10,g21,g129,g39,g197,g182,g185,g2,g43,g87,g180,g138,g188,g38, ... g142,g159,g91,g5,g186,g79,g205,g154,g95,g90,g147,g11,g4,g155,g18,g45, ... g156,g189,g92,g148,g137,g187,g139,g37,g158,g136}; s.name={'C94','C93','C70','C11','C8','C27','C72','C53','C37','C85', ... 'C56','C73','C67','C36','C25','C84','C74','C33','C55','C34','C92','C13', ... 'C62','C82','C46','C9','C65','C24','C30','C12','C99','C26','C91','C100', ... 'C95','C66','C28','C52','C23','C31','C42','C35','C54','C71','C38','C18', ... 'C64','C39','C19','C49','C78','C21','C61','C97','C41','C4','C96','C17', ... 'C29','C51','C47','C5','C10','C50','C16','C98','C83','C86','C1','C20', ... 'C40','C81','C59','C89','C15','C63','C80','C44','C3','C87','C32','R1', ... 'C75','C48','C43','C68','C6','C2','C76','C7','C22','C77','C90','C45', ... 'C69','C58','C88','C60','C14','C79','C57'}; s.tags={'g193','g192','g149','g22','g19','g74','g151','g132','g84', ... 'g184','g135','g152','g146','g83','g72','g183','g153','g80','g134','g81', ... 'g191','g36','g141','g181','g93','g20','g144','g47','g77','g23','g198', ... 'g73','g190','g199','g194','g145','g75','g131','g46','g78','g89','g82', ... 'g133','g150','g85','g41','g143','g86','g42','g128','g157','g44','g140', ... 'g196','g88','g9','g195','g40','g76','g130','g94','g10','g21','g129', ... 'g39','g197','g182','g185','g2','g43','g87','g180','g138','g188','g38', ... 'g142','g159','g91','g5','g186','g79','g205','g154','g95','g90','g147', ... 'g11','g4','g155','g18','g45','g156','g189','g92','g148','g137','g187', ... 'g139','g37','g158','g136'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(0.01,0.01,'base','center','pos',[0,0]); g3=rect2(0.0090,0.0090,'base','center','pos',[0,0]); % Analyzed geometry clear s s.objs={g193,g192,g149,g22,g19,g74,g151,g132,g84,g184,g152,g135, ... g146,g83,g72,g183,g153,g80,g134,g81,g191,g36,g141,g181,g93,g20,g144, ... g47,g77,g23,g198,g190,g73,g194,g199,g145,g131,g75,g46,g78,g82,g89,g150, ... g133,g85,g41,g143,g86,g42,g128,g157,g44,g140,g196,g88,g9,g195,g40,g76, ... g130,g94,g10,g21,g129,g39,g197,g182,g185,g2,g43,g87,g180,g138,g188,g38, ... g142,g159,g91,g5,g186,g79,g154,g95,g90,g147,g11,g4,g155,g18,g45,g156, ... g189,g92,g148,g137,g187,g139,g37,g158,g136,g3}; s.name={'C94','C93','C70','C11','C8','C27','C72','C53','C37','C85', ... 'C73','C56','C67','C36','C25','C84','C74','C33','C55','C34','C92','C13', ... 'C62','C82','C46','C9','C65','C24','C30','C12','C99','C91','C26','C95', ... 'C100','C66','C52','C28','C23','C31','C35','C42','C71','C54','C38','C18', ... 'C64','C39','C19','C49','C78','C21','C61','C97','C41','C4','C96','C17', ... 'C29','C51','C47','C5','C10','C50','C16','C98','C83','C86','C1','C20', ... 'C40','C81','C59','C89','C15','C63','C80','C44','C3','C87','C32','C75', ... 'C48','C43','C68','C6','C2','C76','C7','C22','C77','C90','C45','C69', ... 'C58','C88','C60','C14','C79','C57','R1'}; s.tags={'g193','g192','g149','g22','g19','g74','g151','g132','g84', ... 'g184','g152','g135','g146','g83','g72','g183','g153','g80','g134','g81', ... 'g191','g36','g141','g181','g93','g20','g144','g47','g77','g23','g198', ... 'g190','g73','g194','g199','g145','g131','g75','g46','g78','g82','g89', ... 'g150','g133','g85','g41','g143','g86','g42','g128','g157','g44','g140', ... 'g196','g88','g9','g195','g40','g76','g130','g94','g10','g21','g129', ... 'g39','g197','g182','g185','g2','g43','g87','g180','g138','g188','g38', ... 'g142','g159','g91','g5','g186','g79','g154','g95','g90','g147','g11', ... 'g4','g155','g18','g45','g156','g189','g92','g148','g137','g187','g139', ... 'g37','g158','g136','g3'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); g6=geomcomp({g3,g192,g149,g22,g19,g74,g151,g132,g84,g184,g152,g135,g146,g83,g72,g183,g153,g80,g134,g81,g191,g36,g141,g181,g93,g20,g144,g47,g77,g23,g198,g190,g73,g194,g199,g145,g131,g75,g46,g78,g82,g89,g150,g133,g85,g41,g143,g86,g42,g128,g157,g44,g140,g196,g88,g9,g195,g40,g76,g130,g94,g10,g21,g129,g39,g197,g182,g185,g2,g43,g87,g180,g138,g188,g38,g142,g159,g91,g5,g186,g79,g154,g95,g90,g147,g11,g4,g155,g18,g45,g156,g189,g92,g148,g137,g187,g139,g37,g158,g136,g193},'ns',{'g3','g192','g149','g22','g19','g74','g151','g132','g84','g184','g152','g135','g146','g83','g72','g183','g153','g80','g134','g81','g191','g36','g141','g181','g93','g20','g144','g47','g77','g23','g198','g190','g73','g194','g199','g145','g131','g75','g46','g78','g82','g89','g150','g133','g85','g41','g143','g86','g42','g128','g157','g44','g140','g196','g88','g9','g195','g40','g76','g130','g94','g10','g21','g129','g39','g197','g182','g185','g2','g43','g87','g180','g138','g188','g38','g142','g159','g91','g5','g186','g79','g154','g95','g90','g147','g11','g4','g155','g18','g45','g156','g189','g92','g148','g137','g187','g139','g37','g158','g136','g193'},'sf','g3-g192-g149-g22-g19-g74-g151-g132-g84-g184-g152-g135-g146-g83-g72-g183-g153-g80-g134-g81-g191-g36-g141-g181-g93-g20-g144-g47-g77-g23-g198-g190-g73-g194-g199-g145-g131-g75-g46-g78-g82-g89-g150-g133-g85-g41-g143-g86-g42-g128-g157-g44-g140-g196-g88-g9-g195-g40-g76-g130-g94-g10-g21-g129-g39-g197-g182-g185-g2-g43-g87-g180-g138-g188-g38-g142-g159-g91-g5-g186-g79-g154-g95-g90-g147-g11-g4-g155-g18-g45-g156-g189-g92-g148-g137-g187-g139-g37-g158-g136-g193','edge','none'); % Analyzed geometry clear s s.objs={g6}; s.name={'CO1'}; s.tags={'g6'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); g7=rect2(0.0020,0.0020,'base','corner','pos',[-0.0010,-0.0010]); [g8]=geomcopy({g7}); [g12,g13,g14,g15,g16]=geomcopy({g93,g94,g91,g90,g7}); [g17,g24,g25,g26,g27]=geomcopy({g12,g13,g14,g15,g16}); g17=move(g17,[0.0010,0]); g24=move(g24,[0.0010,0]); g25=move(g25,[0.0010,0]); g26=move(g26,[0.0010,0]); g27=move(g27,[0.0010,0]); g17=move(g17,[5.0E-4,-2.499999999999998E-4]); g24=move(g24,[5.0E-4,-2.499999999999998E-4]); g25=move(g25,[5.0E-4,-2.499999999999998E-4]); g26=move(g26,[5.0E-4,-2.499999999999998E-4]); g27=move(g27,[5.0E-4,-2.499999999999998E-4]); g17=move(g17,[5.0E-4,2.499999999999998E-4]); g24=move(g24,[5.0E-4,2.499999999999998E-4]); g25=move(g25,[5.0E-4,2.499999999999998E-4]); g26=move(g26,[5.0E-4,2.499999999999998E-4]); g27=move(g27,[5.0E-4,2.499999999999998E-4]); g17=move(g17,[-0.0020,0.0020]); g24=move(g24,[-0.0020,0.0020]); g25=move(g25,[-0.0020,0.0020]); g26=move(g26,[-0.0020,0.0020]); g27=move(g27,[-0.0020,0.0020]); g28=geomcomp({g7,g94,g91,g90,g93},'ns',{'g7','g94','g91','g90','g93'},'sf','g7-g94-g91-g90-g93','edge','none'); % Analyzed geometry clear s s.objs={g28}; s.name={'CO1'}; s.tags={'g28'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) for ii=1:n ii kx(ii)=(ii-1)/(n-1)*pi*1000/2 ky=0 % Constants fem.const = { %'d','5.3/1000', ... % 'r','0.8/1000', ... 'kx' ,num2str(kx(ii)), ...... 'ky',num2str(ky)} % 'k2','0', ... % 'a1x','d', ... % 'a1y','0', ... % 'a2x','0', ... % 'a2y','d'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.5'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 494; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2008/09/19 16:09:48 $'; fem.version = vrsn; % Initialize mesh fem.mesh=meshinit(fem, ... 'hmax',[500e-6], ... 'hmaxfact',0.8); % (Default values are not included) % Application mode 1 clear appl appl.mode.class = 'Acoustics'; appl.assignsuffix = '_aco'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.type = 'SH'; bnd.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.bnd = bnd; appl.var = {'freq','1.6e5'}; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1'}; src = cell(1,1); clear bnd bnd.expr = {{'p',{},{}},{{},{},'p'}}; bnd.map = {{'1','1','1'},{'1','1','1'}}; bnd.ind = {{'1'},{'2','4','5','6','7','8','9','10','11','12','13','14', ... '15','16','17','18','19','20'},{'3'}}; src{1} = {{},bnd,{}}; elem.src = src; geomdim = cell(1,1); clear bnd bnd.map = {{{},{},'2'},{{},'3',{}}}; bnd.ind = {{'1','3','5','6','7','8','9','10','11','12','13','14','15', ... '16','17','18','19','20'},{'2'},{'4'}}; geomdim{1} = {{},bnd,{}}; elem.geomdim = geomdim; elem.var = {'pconstr1','pconstr2'}; map = cell(1,3); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'19','20'}; submap.dg = '1'; submap.dv = {'1','2'}; map{2} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'1','19'}; submap.dg = '1'; submap.dv = {'2','20'}; map{3} = submap; elem.map = map; elemcpl{1} = elem; % Point constraint variables (used for periodic conditions) clear elem elem.elem = 'elpconstr'; elem.g = {'1'}; clear bnd bnd.constr = {{'0','pconstr2-(p)'},{'pconstr1-(p)','0'}}; bnd.cpoints = {{'2','2'},{'2','2'}}; bnd.ind = {{'2'},{'4'}}; elem.geomdim = {{{},bnd,{}}}; elemcpl{2} = elem; fem.elemcpl = elemcpl; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; fem.ode=ode; % Multiphysics fem=multiphysics(fem); % Extend mesh fem.xmesh=meshextend(fem); % Solve problem fem.sol=femeig(fem, ... 'solcomp',{'p'}, ... 'outcomp',{'p'}, ... 'blocksize','auto', ... 'neigs',100, ... 'shift',-9.424778*i, ... 'eigref','-i*2*pi*100'); % Save current fem structure for restart purposes fem0=fem; % Plot solution % postplot(fem, ... % 'tridata',{'p','cont','internal','unit','Pa'}, ... % 'trimap','jet(1024)', ... % 'solnum',1, ... % 'title','freq_aco(1)=0.037713-0.029703i Surface: Pressure [Pa]', ... % 'axis',[-0.001964162847596324,0.001964162847596324,-0.001373921991879247,0.0013739219918792467]); % (Default values are not included) % Application mode 1 clear appl appl.mode.class = 'Acoustics'; appl.assignsuffix = '_aco'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.type = 'SH'; bnd.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.bnd = bnd; appl.var = {'freq','1.6e5'}; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1'}; src = cell(1,1); clear bnd bnd.expr = {{'p',{},{}},{{},{},'p'}}; bnd.map = {{'1','1','1'},{'1','1','1'}}; bnd.ind = {{'1'},{'2','4','5','6','7','8','9','10','11','12','13','14', ... '15','16','17','18','19','20'},{'3'}}; src{1} = {{},bnd,{}}; elem.src = src; geomdim = cell(1,1); clear bnd bnd.map = {{{},{},'2'},{{},'3',{}}}; bnd.ind = {{'1','3','5','6','7','8','9','10','11','12','13','14','15', ... '16','17','18','19','20'},{'2'},{'4'}}; geomdim{1} = {{},bnd,{}}; elem.geomdim = geomdim; elem.var = {'pconstr1','pconstr2'}; map = cell(1,3); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'19','20'}; submap.dg = '1'; submap.dv = {'1','2'}; map{2} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'1','19'}; submap.dg = '1'; submap.dv = {'2','20'}; map{3} = submap; elem.map = map; elemcpl{1} = elem; % Point constraint variables (used for periodic conditions) clear elem elem.elem = 'elpconstr'; elem.g = {'1'}; clear bnd bnd.constr = {{'0','pconstr2-(p*exp(j*ky*y))'},{'pconstr1-(p*exp(j*kx*x))', ... '0'}}; bnd.cpoints = {{'2','2'},{'2','2'}}; bnd.ind = {{'2'},{'4'}}; elem.geomdim = {{{},bnd,{}}}; elemcpl{2} = elem; fem.elemcpl = elemcpl; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; fem.ode=ode; % Multiphysics fem=multiphysics(fem); freky(1:100,ii)=real(fem.sol.lambda(1:100)*i/2/pi) end %%% Save eigenfrequencies d=str2num(cell2mat(fem.const(2))); %%% Plot dispersion relation plot(kx,freky(1,:),'.')