% COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5a (COMSOL 3.5.0.603, $Date: 2008/12/03 17:02:19 $) % Some geometry objects are stored in a separate file. % The name of this file is given by the variable 'flbinaryfile'. clear flclear fem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.5'; vrsn.ext = 'a'; vrsn.major = 0; vrsn.build = 603; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2008/12/03 17:02:19 $'; fem.version = vrsn; flbinaryfile='circ3_1inc.mphm'; % Geometry clear draw g10=flbinary('g10','draw',flbinaryfile); g11=flbinary('g11','draw',flbinaryfile); g12=flbinary('g12','draw',flbinaryfile); g13=flbinary('g13','draw',flbinaryfile); g14=flbinary('g14','draw',flbinaryfile); g15=flbinary('g15','draw',flbinaryfile); g4=flbinary('g4','draw',flbinaryfile); g5=flbinary('g5','draw',flbinaryfile); g6=flbinary('g6','draw',flbinaryfile); g7=flbinary('g7','draw',flbinaryfile); g8=flbinary('g8','draw',flbinaryfile); g9=flbinary('g9','draw',flbinaryfile); draw.p.objs = {g10,g11,g12,g13,g14,g15,g4,g5,g6,g7,g8,g9}; draw.p.name = {'PT1','PT10','PT11','PT12','PT2','PT3','PT4','PT5','PT6','PT7', ... 'PT8','PT9'}; draw.p.tags = {'g10','g11','g12','g13','g14','g15','g4','g5','g6','g7','g8', ... 'g9'}; g1=flbinary('g1','draw',flbinaryfile); g16=flbinary('g16','draw',flbinaryfile); g3=flbinary('g3','draw',flbinaryfile); draw.s.objs = {g1,g16,g3}; draw.s.name = {'C1','C2','CO1'}; draw.s.tags = {'g1','g16','g3'}; fem.draw = draw; fem.geom = geomcsg(fem); % Initialize mesh fem.mesh=meshinit(fem, ... 'hauto',5); % (Default values are not included) % Application mode 1 clear appl appl.mode.class = 'ConductiveMediaDC'; appl.assignsuffix = '_dc'; clear prop prop.elemdefault='Lag1'; appl.prop = prop; clear pnt pnt.name = {'','sink','source'}; pnt.Qj0 = {0,-1.0,1.0}; pnt.ind = [3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2]; appl.pnt = pnt; clear bnd bnd.type = {'nJ0','cont'}; bnd.ind = [1,1,1,1,2,2,2,2,1,1,2,2,2,2,1,1,1,2,2,2,2,1,1,1,1,1,1,1]; appl.bnd = bnd; clear equ equ.sigma = {1.0,10.0,1.0,10.0}; equ.ind = [1,2,3,4]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; fem.ode=ode; % Multiphysics fem=multiphysics(fem); % KJE added 2009Jul23 fem.xmesh = meshextend(fem); for isig = 1:5 sig2 = isig * 10; equ.sigma{2} = sig2; equ.sigma fem.sol = femstatic(fem); save(isig,1:539) = fem.sol.u; end