function out = model % % ichichaa.m % % Model exported on Nov 16 2015, 19:13 by COMSOL 4.4.0.150. import com.comsol.model.* import com.comsol.model.util.* model = ModelUtil.create('Model'); model.modelPath('C:\Users\Annisanurf\Documents\SIMULASI\Latihan'); model.modelNode.create('comp1'); model.geom.create('geom1', 3); model.mesh.create('mesh1', 'geom1'); model.physics.create('mf', 'InductionCurrents', 'geom1'); model.physics.create('cir', 'Circuit', 'geom1'); model.study.create('std1'); model.study('std1').feature.create('time', 'Transient'); model.study('std1').feature('time').activate('mf', true); model.study('std1').feature('time').activate('cir', true); model.geom('geom1').lengthUnit('mm'); model.geom('geom1').scaleUnitValue(true); model.geom('geom1').feature.create('imp1', 'Import'); model.geom('geom1').feature('imp1').set('filename', 'C:\Users\Annisanurf\Documents\SIMULASI\Latihan\icichaa.mphbin'); model.geom('geom1').feature('imp1').importData; model.geom('geom1').run('imp1'); model.geom('geom1').feature.create('sph1', 'Sphere'); model.geom('geom1').feature('sph1').set('r', '20'); model.view('view1').hideObjects.create('hide1'); model.view('view1').hideObjects('hide1').init; model.view('view1').hideObjects('hide1').add({'imp1(6)'}); model.geom('geom1').runPre('fin'); model.geom('geom1').run; model.material.create('mat1'); model.material('mat1').name('Air'); model.material('mat1').set('family', 'air'); model.material('mat1').propertyGroup('def').set('relpermeability', '1'); model.material('mat1').propertyGroup('def').set('relpermittivity', '1'); model.material('mat1').propertyGroup('def').set('dynamicviscosity', 'eta(T[1/K])[Pa*s]'); model.material('mat1').propertyGroup('def').set('ratioofspecificheat', '1.4'); model.material('mat1').propertyGroup('def').set('electricconductivity', '0[S/m]'); model.material('mat1').propertyGroup('def').set('heatcapacity', 'Cp(T[1/K])[J/(kg*K)]'); model.material('mat1').propertyGroup('def').set('density', 'rho(pA[1/Pa],T[1/K])[kg/m^3]'); model.material('mat1').propertyGroup('def').set('thermalconductivity', 'k(T[1/K])[W/(m*K)]'); model.material('mat1').propertyGroup('def').set('soundspeed', 'cs(T[1/K])[m/s]'); model.material('mat1').propertyGroup('def').func.create('eta', 'Piecewise'); model.material('mat1').propertyGroup('def').func('eta').set('funcname', 'eta'); model.material('mat1').propertyGroup('def').func('eta').set('arg', 'T'); model.material('mat1').propertyGroup('def').func('eta').set('extrap', 'constant'); model.material('mat1').propertyGroup('def').func('eta').set('pieces', {'200.0' '1600.0' '-8.38278E-7+8.35717342E-8*T^1-7.69429583E-11*T^2+4.6437266E-14*T^3-1.06585607E-17*T^4'}); model.material('mat1').propertyGroup('def').func.create('Cp', 'Piecewise'); model.material('mat1').propertyGroup('def').func('Cp').set('funcname', 'Cp'); model.material('mat1').propertyGroup('def').func('Cp').set('arg', 'T'); model.material('mat1').propertyGroup('def').func('Cp').set('extrap', 'constant'); model.material('mat1').propertyGroup('def').func('Cp').set('pieces', {'200.0' '1600.0' '1047.63657-0.372589265*T^1+9.45304214E-4*T^2-6.02409443E-7*T^3+1.2858961E-10*T^4'}); model.material('mat1').propertyGroup('def').func.create('rho', 'Analytic'); model.material('mat1').propertyGroup('def').func('rho').set('funcname', 'rho'); model.material('mat1').propertyGroup('def').func('rho').set('args', {'pA' 'T'}); model.material('mat1').propertyGroup('def').func('rho').set('expr', 'pA*0.02897/8.314/T'); model.material('mat1').propertyGroup('def').func('rho').set('dermethod', 'manual'); model.material('mat1').propertyGroup('def').func('rho').set('argders', {'pA' 'd(pA*0.02897/8.314/T,pA)'; 'T' 'd(pA*0.02897/8.314/T,T)'}); model.material('mat1').propertyGroup('def').func.create('k', 'Piecewise'); model.material('mat1').propertyGroup('def').func('k').set('funcname', 'k'); model.material('mat1').propertyGroup('def').func('k').set('arg', 'T'); model.material('mat1').propertyGroup('def').func('k').set('extrap', 'constant'); model.material('mat1').propertyGroup('def').func('k').set('pieces', {'200.0' '1600.0' '-0.00227583562+1.15480022E-4*T^1-7.90252856E-8*T^2+4.11702505E-11*T^3-7.43864331E-15*T^4'}); model.material('mat1').propertyGroup('def').func.create('cs', 'Analytic'); model.material('mat1').propertyGroup('def').func('cs').set('funcname', 'cs'); model.material('mat1').propertyGroup('def').func('cs').set('args', {'T'}); model.material('mat1').propertyGroup('def').func('cs').set('expr', 'sqrt(1.4*287*T)'); model.material('mat1').propertyGroup('def').func('cs').set('dermethod', 'manual'); model.material('mat1').propertyGroup('def').func('cs').set('argders', {'T' 'd(sqrt(1.4*287*T),T)'}); model.material('mat1').propertyGroup('def').addInput('temperature'); model.material('mat1').propertyGroup('def').addInput('pressure'); model.material('mat1').propertyGroup.create('RefractiveIndex', 'Refractive index'); model.material('mat1').propertyGroup('RefractiveIndex').set('n', '1'); model.material('mat1').set('family', 'air'); model.material('mat1').selection.set([1]); model.material.create('mat2'); model.material('mat2').name('Air (2)'); model.material('mat2').set('family', 'air'); model.material('mat2').propertyGroup('def').set('relpermeability', '1'); model.material('mat2').propertyGroup('def').set('relpermittivity', '1'); model.material('mat2').propertyGroup('def').set('dynamicviscosity', 'eta(T[1/K])[Pa*s]'); model.material('mat2').propertyGroup('def').set('ratioofspecificheat', '1.4'); model.material('mat2').propertyGroup('def').set('electricconductivity', '0[S/m]'); model.material('mat2').propertyGroup('def').set('heatcapacity', 'Cp(T[1/K])[J/(kg*K)]'); model.material('mat2').propertyGroup('def').set('density', 'rho(pA[1/Pa],T[1/K])[kg/m^3]'); model.material('mat2').propertyGroup('def').set('thermalconductivity', 'k(T[1/K])[W/(m*K)]'); model.material('mat2').propertyGroup('def').set('soundspeed', 'cs(T[1/K])[m/s]'); model.material('mat2').propertyGroup('def').func.create('eta', 'Piecewise'); model.material('mat2').propertyGroup('def').func('eta').set('funcname', 'eta'); model.material('mat2').propertyGroup('def').func('eta').set('arg', 'T'); model.material('mat2').propertyGroup('def').func('eta').set('extrap', 'constant'); model.material('mat2').propertyGroup('def').func('eta').set('pieces', {'200.0' '1600.0' '-8.38278E-7+8.35717342E-8*T^1-7.69429583E-11*T^2+4.6437266E-14*T^3-1.06585607E-17*T^4'}); model.material('mat2').propertyGroup('def').func.create('Cp', 'Piecewise'); model.material('mat2').propertyGroup('def').func('Cp').set('funcname', 'Cp'); model.material('mat2').propertyGroup('def').func('Cp').set('arg', 'T'); model.material('mat2').propertyGroup('def').func('Cp').set('extrap', 'constant'); model.material('mat2').propertyGroup('def').func('Cp').set('pieces', {'200.0' '1600.0' '1047.63657-0.372589265*T^1+9.45304214E-4*T^2-6.02409443E-7*T^3+1.2858961E-10*T^4'}); model.material('mat2').propertyGroup('def').func.create('rho', 'Analytic'); model.material('mat2').propertyGroup('def').func('rho').set('funcname', 'rho'); model.material('mat2').propertyGroup('def').func('rho').set('args', {'pA' 'T'}); model.material('mat2').propertyGroup('def').func('rho').set('expr', 'pA*0.02897/8.314/T'); model.material('mat2').propertyGroup('def').func('rho').set('dermethod', 'manual'); model.material('mat2').propertyGroup('def').func('rho').set('argders', {'pA' 'd(pA*0.02897/8.314/T,pA)'; 'T' 'd(pA*0.02897/8.314/T,T)'}); model.material('mat2').propertyGroup('def').func.create('k', 'Piecewise'); model.material('mat2').propertyGroup('def').func('k').set('funcname', 'k'); model.material('mat2').propertyGroup('def').func('k').set('arg', 'T'); model.material('mat2').propertyGroup('def').func('k').set('extrap', 'constant'); model.material('mat2').propertyGroup('def').func('k').set('pieces', {'200.0' '1600.0' '-0.00227583562+1.15480022E-4*T^1-7.90252856E-8*T^2+4.11702505E-11*T^3-7.43864331E-15*T^4'}); model.material('mat2').propertyGroup('def').func.create('cs', 'Analytic'); model.material('mat2').propertyGroup('def').func('cs').set('funcname', 'cs'); model.material('mat2').propertyGroup('def').func('cs').set('args', {'T'}); model.material('mat2').propertyGroup('def').func('cs').set('expr', 'sqrt(1.4*287*T)'); model.material('mat2').propertyGroup('def').func('cs').set('dermethod', 'manual'); model.material('mat2').propertyGroup('def').func('cs').set('argders', {'T' 'd(sqrt(1.4*287*T),T)'}); model.material('mat2').propertyGroup('def').addInput('temperature'); model.material('mat2').propertyGroup('def').addInput('pressure'); model.material('mat2').propertyGroup.create('RefractiveIndex', 'Refractive index'); model.material('mat2').propertyGroup('RefractiveIndex').set('n', '1'); model.material('mat2').set('family', 'air'); model.material.create('mat3'); model.material('mat3').name('Copper'); model.material('mat3').set('family', 'copper'); model.material('mat3').propertyGroup('def').set('relpermeability', '1'); model.material('mat3').propertyGroup('def').set('electricconductivity', '5.998e7[S/m]'); model.material('mat3').propertyGroup('def').set('thermalexpansioncoefficient', '17e-6[1/K]'); model.material('mat3').propertyGroup('def').set('heatcapacity', '385[J/(kg*K)]'); model.material('mat3').propertyGroup('def').set('relpermittivity', '1'); model.material('mat3').propertyGroup('def').set('density', '8700[kg/m^3]'); model.material('mat3').propertyGroup('def').set('thermalconductivity', '400[W/(m*K)]'); model.material('mat3').propertyGroup.create('Enu', 'Young''s modulus and Poisson''s ratio'); model.material('mat3').propertyGroup('Enu').set('poissonsratio', '0.35'); model.material('mat3').propertyGroup('Enu').set('youngsmodulus', '110e9[Pa]'); model.material('mat3').propertyGroup.create('linzRes', 'Linearized resistivity'); model.material('mat3').propertyGroup('linzRes').set('alpha', '0.0039[1/K]'); model.material('mat3').propertyGroup('linzRes').set('rho0', '1.72e-8[ohm*m]'); model.material('mat3').propertyGroup('linzRes').set('Tref', '298[K]'); model.material('mat3').set('family', 'copper'); model.material('mat2').selection.set([]); model.material('mat2').selection.all; model.material('mat2').selection.set([2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19]); model.material('mat3').selection.all; model.material('mat3').selection.set([2 3 4 5 7 8 10 11 12 13 14 15 16 17 18 19]); model.material('mat2').selection.set([9]); model.material('mat1').selection.set([1 6]); model.view('view1').hideEntities.create('hide1'); model.view('view1').hideEntities('hide1').geom(3); model.view('view1').hideEntities('hide1').add([1]); model.view('view1').hideEntities.clear; model.material('mat1').propertyGroup('def').set('electricconductivity', {'10[S/m]'}); model.material('mat2').propertyGroup('def').set('electricconductivity', {'10[S/m]'}); model.physics('mf').prop('ShapeProperty').set('order_magneticvectorpotential', 1, '1'); model.physics('mf').feature.create('mtcd1', 'MultiTurnCoilDomain', 3); model.view('view1').hideEntities.create('hide1'); model.view('view1').hideEntities('hide1').geom(3); model.view('view1').hideEntities('hide1').add([1]); model.physics('mf').feature('mtcd1').selection.set([2 3 4 5]); model.physics('mf').feature('mtcd1').set('CoilType', 1, 'Numeric'); model.physics('mf').feature('mtcd1').set('N', 1, '30'); model.physics('mf').feature('mtcd1').set('CoilExcitation', 1, 'CircuitCurrent'); model.physics('mf').feature('mtcd1').feature.create('ccc1', 'CoilCurrentCalculation', 3); model.physics('mf').feature('mtcd1').feature('ccc1').selection.set([2 3 4 5]); model.physics('mf').feature('mtcd1').feature('ccc1').feature.create('ci1', 'CoilInsulation', 2); model.physics('mf').feature('mtcd1').feature('ccc1').feature('ci1').selection.set([6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25]); model.physics('mf').feature('mtcd1').feature('ccc1').feature.create('ct1', 'CoilTerminal', 2); model.physics('mf').feature('mtcd1').feature('ccc1').feature('ct1').selection.set([6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 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 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100]); model.physics('mf').feature('mtcd1').feature('ccc1').feature('ci1').selection.set([6 7 9 10 11 12 13 15 16 17 19 20 22 23 24 25]); model.physics('mf').feature('mtcd1').feature('ccc1').feature('ct1').selection.set([14]); model.physics('mf').feature.create('gfa1', 'GaugeFixingA', 3); model.physics('mf').feature('gfa1').selection.set([2 3 4 5]); model.physics('mf').feature('gfa1').set('SetPointConstraint', 1, '1'); model.physics('mf').feature.create('mtcd2', 'MultiTurnCoilDomain', 3); model.physics('mf').feature('mtcd2').selection.set([7 8 12 13]); model.physics('mf').feature('mtcd2').set('CoilType', 1, 'Numeric'); model.physics('mf').feature('mtcd2').set('N', 1, '30'); model.physics('mf').feature('mtcd2').set('CoilExcitation', 1, 'CircuitCurrent'); model.physics('mf').feature('mtcd2').feature.create('ccc1', 'CoilCurrentCalculation', 3); model.physics('mf').feature('mtcd2').feature('ccc1').selection.set([7 8 12 13]); model.physics('mf').feature('mtcd2').feature('ccc1').feature.create('ci1', 'CoilInsulation', 2); model.physics('mf').feature('mtcd2').feature('ccc1').feature('ci1').selection.set([31 32 33 34 35 36 48 49 53 54 55 57 58 59 60 61]); model.physics('mf').feature('mtcd2').feature('ccc1').feature.create('ct1', 'CoilTerminal', 2); model.physics('mf').feature('mtcd2').feature('ccc1').feature('ct1').selection.set([52]); model.physics('mf').feature.create('gfa2', 'GaugeFixingA', 3); model.physics('mf').feature('gfa2').selection.set([7 8 12 13]); model.physics('mf').feature('gfa2').set('SetPointConstraint', 1, '1'); model.physics('mf').feature.create('mtcd3', 'MultiTurnCoilDomain', 3); model.physics('mf').feature('mtcd3').selection.set([16 17 18 19]); model.physics('mf').feature('mtcd3').set('CoilType', 1, 'Numeric'); model.physics('mf').feature('mtcd3').set('N', 1, '30'); model.physics('mf').feature('mtcd3').set('CoilExcitation', 1, 'CircuitCurrent'); model.physics('mf').feature('mtcd3').feature.create('ccc1', 'CoilCurrentCalculation', 3); model.physics('mf').feature('mtcd3').feature('ccc1').selection.set([16 17 18 19]); model.physics('mf').feature('mtcd3').feature('ccc1').feature.create('ci1', 'CoilInsulation', 2); model.physics('mf').feature('mtcd3').feature('ccc1').feature('ci1').selection.set([80 81 83 84 85 86 87 89 90 91 93 94 96 97 98 99]); model.physics('mf').feature('mtcd3').feature('ccc1').feature.create('ct1', 'CoilTerminal', 2); model.physics('mf').feature('mtcd3').feature('ccc1').feature('ct1').selection.set([88]); model.physics('mf').feature.create('gfa3', 'GaugeFixingA', 3); model.physics('mf').feature('gfa3').selection.set([16 17 18 19]); model.physics('mf').feature('gfa3').set('SetPointConstraint', 1, '1'); model.physics('mf').feature.create('mtcd4', 'MultiTurnCoilDomain', 3); model.physics('mf').feature('mtcd4').selection.set([10 11 14 15]); model.physics('mf').feature('mtcd4').set('CoilType', 1, 'Numeric'); model.physics('mf').feature('mtcd4').set('N', 1, '30'); model.physics('mf').feature('mtcd4').set('CoilExcitation', 1, 'CircuitCurrent'); model.physics('mf').feature('mtcd4').feature.create('ccc1', 'CoilCurrentCalculation', 3); model.physics('mf').feature('mtcd4').feature('ccc1').selection.set([10 11 14 15]); model.physics('mf').feature('mtcd4').feature('ccc1').feature.create('ci1', 'CoilInsulation', 2); model.physics('mf').feature('mtcd4').feature('ccc1').feature('ci1').selection.set([42 43 44 45 46 47 50 51 69 70 71 73 74 75 76 77]); model.physics('mf').feature('mtcd4').feature('ccc1').feature.create('ct1', 'CoilTerminal', 2); model.physics('mf').feature('mtcd4').feature('ccc1').feature('ct1').selection.set([68]); model.physics('mf').feature.create('gfa4', 'GaugeFixingA', 3); model.physics('mf').feature('gfa4').selection.set([10 11 14 15]); model.physics('mf').feature('gfa4').set('SetPointConstraint', 1, '1'); model.physics('cir').feature.create('V1', 'VoltageSource', -1); model.physics('cir').feature('V1').set('Connections', 1, 1, '5'); model.physics('cir').feature('V1').set('sourceType', 1, 'SineSource'); model.physics('cir').feature('V1').set('value', 1, '12.5[V]'); model.physics('cir').feature('V1').set('freq', 1, '100[kHz]'); model.physics('cir').feature.create('IvsU1', 'ModelDeviceIV', -1); model.physics('cir').feature('IvsU1').set('V_src', 1, 'root.comp1.mf.VCoil_1'); model.physics('cir').feature.create('R1', 'Resistor', -1); model.physics('cir').feature('R1').set('Connections', 1, 1, '0'); model.physics('cir').feature('R1').set('Connections', 2, 1, '1'); model.physics('cir').feature('R1').set('R', 1, ['50[' 'ohm' ']']); model.physics('cir').feature.create('IvsU2', 'ModelDeviceIV', -1); model.physics('cir').feature('IvsU2').set('Connections', 1, 1, '2'); model.physics('cir').feature('IvsU2').set('V_src', 1, 'root.comp1.mf.VCoil_2'); model.physics('cir').feature.create('R2', 'Resistor', -1); model.physics('cir').feature('R2').set('Connections', 1, 1, '0'); model.physics('cir').feature('R2').set('Connections', 2, 1, '2'); model.physics('cir').feature('R2').set('R', 1, ['50[' 'ohm' ']']); model.physics('cir').feature.create('IvsU3', 'ModelDeviceIV', -1); model.physics('cir').feature('IvsU3').set('Connections', 1, 1, '3'); model.physics('cir').feature('IvsU3').set('V_src', 1, 'root.comp1.mf.VCoil_3'); model.physics('cir').feature.create('R3', 'Resistor', -1); model.physics('cir').feature('R3').set('Connections', 1, 1, '0'); model.physics('cir').feature('R3').set('Connections', 2, 1, '3'); model.physics('cir').feature('R3').set('R', 1, ['50[' 'ohm' ']']); model.physics('cir').feature.create('IvsU4', 'ModelDeviceIV', -1); model.physics('cir').feature('IvsU4').set('Connections', 1, 1, '4'); model.physics('cir').feature('IvsU4').set('V_src', 1, 'root.comp1.mf.VCoil_4'); model.physics('cir').feature.create('R4', 'Resistor', -1); model.physics('cir').feature('R4').set('Connections', 1, 1, '0'); model.physics('cir').feature('R4').set('Connections', 2, 1, '4'); model.physics('cir').feature('R4').set('R', 1, ['50[' 'ohm' ']']); model.physics('cir').feature('V1').set('Connections', 2, 1, '1'); model.physics('cir').feature('V1').set('Connections', 2, 1, '0'); model.physics('cir').feature('IvsU1').set('Connections', 2, 1, '5'); model.mesh('mesh1').run; model.view('view1').hideEntities.clear; model.study('std1').feature.create('ccc', 'CoilCurrentCalculation'); model.study('std1').feature.create('ccc2', 'CoilCurrentCalculation'); model.study('std1').feature('ccc2').set('CoilName', '2'); model.study('std1').feature.create('ccc3', 'CoilCurrentCalculation'); model.study('std1').feature('ccc3').set('CoilName', '3'); model.study('std1').feature.create('ccc4', 'CoilCurrentCalculation'); model.study('std1').feature('ccc4').set('CoilName', '4'); model.study('std1').feature.move('time', 4); model.study('std1').feature('time').set('tlist', 'range(0,5.0e-8,3.0e-5)'); model.study('std1').feature('time').set('rtolactive', 'on'); model.study('std1').feature('time').set('rtol', '0.001'); model.sol.create('sol1'); model.sol('sol1').study('std1'); model.sol('sol1').feature.create('st1', 'StudyStep'); model.sol('sol1').feature('st1').set('study', 'std1'); model.sol('sol1').feature('st1').set('studystep', 'ccc'); model.sol('sol1').feature.create('v1', 'Variables'); model.sol('sol1').feature('v1').set('control', 'ccc'); model.sol('sol1').feature.create('e1', 'Eigenvalue'); model.sol('sol1').feature('e1').set('shift', '0'); model.sol('sol1').feature('e1').set('neigs', 1); model.sol('sol1').feature('e1').set('control', 'ccc'); model.sol('sol1').feature.create('su1', 'StoreSolution'); model.sol('sol1').feature.create('st2', 'StudyStep'); model.sol('sol1').feature('st2').set('study', 'std1'); model.sol('sol1').feature('st2').set('studystep', 'ccc2'); model.sol('sol1').feature.create('v2', 'Variables'); model.sol('sol1').feature('v2').set('initmethod', 'sol'); model.sol('sol1').feature('v2').set('initsol', 'sol1'); model.sol('sol1').feature('v2').set('notsolmethod', 'sol'); model.sol('sol1').feature('v2').set('notsol', 'sol1'); model.sol('sol1').feature('v2').set('control', 'ccc2'); model.sol('sol1').feature.create('e2', 'Eigenvalue'); model.sol('sol1').feature('e2').set('shift', '0'); model.sol('sol1').feature('e2').set('neigs', 1); model.sol('sol1').feature('e2').set('control', 'ccc2'); model.sol('sol1').feature.create('su2', 'StoreSolution'); model.sol('sol1').feature.create('st3', 'StudyStep'); model.sol('sol1').feature('st3').set('study', 'std1'); model.sol('sol1').feature('st3').set('studystep', 'ccc3'); model.sol('sol1').feature.create('v3', 'Variables'); model.sol('sol1').feature('v3').set('initmethod', 'sol'); model.sol('sol1').feature('v3').set('initsol', 'sol1'); model.sol('sol1').feature('v3').set('notsolmethod', 'sol'); model.sol('sol1').feature('v3').set('notsol', 'sol1'); model.sol('sol1').feature('v3').set('control', 'ccc3'); model.sol('sol1').feature.create('e3', 'Eigenvalue'); model.sol('sol1').feature('e3').set('shift', '0'); model.sol('sol1').feature('e3').set('neigs', 1); model.sol('sol1').feature('e3').set('control', 'ccc3'); model.sol('sol1').feature.create('su3', 'StoreSolution'); model.sol('sol1').feature.create('st4', 'StudyStep'); model.sol('sol1').feature('st4').set('study', 'std1'); model.sol('sol1').feature('st4').set('studystep', 'ccc4'); model.sol('sol1').feature.create('v4', 'Variables'); model.sol('sol1').feature('v4').set('initmethod', 'sol'); model.sol('sol1').feature('v4').set('initsol', 'sol1'); model.sol('sol1').feature('v4').set('notsolmethod', 'sol'); model.sol('sol1').feature('v4').set('notsol', 'sol1'); model.sol('sol1').feature('v4').set('control', 'ccc4'); model.sol('sol1').feature.create('e4', 'Eigenvalue'); model.sol('sol1').feature('e4').set('shift', '0'); model.sol('sol1').feature('e4').set('neigs', 1); model.sol('sol1').feature('e4').set('control', 'ccc4'); model.sol('sol1').feature.create('su4', 'StoreSolution'); model.sol('sol1').feature.create('st5', 'StudyStep'); model.sol('sol1').feature('st5').set('study', 'std1'); model.sol('sol1').feature('st5').set('studystep', 'time'); model.sol('sol1').feature.create('v5', 'Variables'); model.sol('sol1').feature('v5').set('initmethod', 'sol'); model.sol('sol1').feature('v5').set('initsol', 'sol1'); model.sol('sol1').feature('v5').set('notsolmethod', 'sol'); model.sol('sol1').feature('v5').set('notsol', 'sol1'); model.sol('sol1').feature('v5').set('control', 'time'); model.shape('shape1').feature('shfun1'); model.sol('sol1').feature.create('t1', 'Time'); model.sol('sol1').feature('t1').set('tlist', 'range(0,5.0e-8,3.0e-5)'); model.sol('sol1').feature('t1').set('plot', 'off'); model.sol('sol1').feature('t1').set('plotfreq', 'tout'); model.sol('sol1').feature('t1').set('probesel', 'all'); model.sol('sol1').feature('t1').set('probes', {}); model.sol('sol1').feature('t1').set('probefreq', 'tsteps'); model.sol('sol1').feature('t1').set('control', 'time'); model.sol('sol1').feature('t1').feature.create('se1', 'Segregated'); model.sol('sol1').feature('t1').feature('se1').feature.remove('ssDef'); model.sol('sol1').feature('t1').feature('se1').feature.create('ss1', 'SegregatedStep'); model.sol('sol1').feature('t1').feature('se1').feature('ss1').set('segvar', {'comp1_currents'}); model.sol('sol1').feature('t1').feature('se1').feature('ss1').set('linsolver', 'dDef'); model.sol('sol1').feature('t1').feature('se1').feature.create('ss2', 'SegregatedStep'); model.sol('sol1').feature('t1').feature('se1').feature('ss2').set('segvar', {'comp1_A' 'comp1_mf_psi'}); model.sol('sol1').feature('t1').feature.create('i1', 'Iterative'); model.sol('sol1').feature('t1').feature('i1').set('linsolver', 'gmres'); model.sol('sol1').feature('t1').feature('i1').set('prefuntype', 'right'); model.sol('sol1').feature('t1').feature('i1').feature.create('mg1', 'Multigrid'); model.sol('sol1').feature('t1').feature('i1').feature('mg1').feature('pr').feature.create('va1', 'Vanka'); model.sol('sol1').feature('t1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'comp1_mf_psi'}); model.sol('sol1').feature('t1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka'); model.sol('sol1').feature('t1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'comp1_mf_psi'}); model.sol('sol1').feature('t1').feature('se1').feature('ss2').set('linsolver', 'i1'); model.sol('sol1').feature('t1').feature.remove('fcDef'); model.sol('sol1').feature('v2').set('notsolnum', 'auto'); model.sol('sol1').feature('v2').set('notsolvertype', 'solnum'); model.sol('sol1').feature('v2').set('solnum', 'auto'); model.sol('sol1').feature('v2').set('solvertype', 'solnum'); model.sol('sol1').feature('v3').set('notsolnum', 'auto'); model.sol('sol1').feature('v3').set('notsolvertype', 'solnum'); model.sol('sol1').feature('v3').set('solnum', 'auto'); model.sol('sol1').feature('v3').set('solvertype', 'solnum'); model.sol('sol1').feature('v4').set('notsolnum', 'auto'); model.sol('sol1').feature('v4').set('notsolvertype', 'solnum'); model.sol('sol1').feature('v4').set('solnum', 'auto'); model.sol('sol1').feature('v4').set('solvertype', 'solnum'); model.sol('sol1').feature('v5').set('notsolnum', 'auto'); model.sol('sol1').feature('v5').set('notsolvertype', 'solnum'); model.sol('sol1').feature('v5').set('solnum', 'auto'); model.sol('sol1').feature('v5').set('solvertype', 'solnum'); model.sol('sol1').attach('std1'); model.sol('sol1').feature('t1').feature.create('fc1', 'FullyCoupled'); model.sol('sol1').feature('t1').feature('fc1').set('linsolver', 'dDef'); model.sol('sol1').feature('t1').feature('fc1').set('jtech', 'onevery'); model.sol('sol1').feature('t1').feature('fc1').set('maxiter', '25'); model.result.create('pg1', 'PlotGroup3D'); model.result('pg1').name('Magnetic Flux Density Norm (mf)'); model.result('pg1').set('oldanalysistype', 'noneavailable'); model.result('pg1').set('data', 'dset1'); model.result('pg1').feature.create('mslc1', 'Multislice'); model.result('pg1').feature('mslc1').set('oldanalysistype', 'noneavailable'); model.result('pg1').feature('mslc1').set('data', 'parent'); model.sol('sol1').runAll; model.result('pg1').run; model.result.numerical.create('gev1', 'EvalGlobal'); model.result.numerical('gev1').set('expr', 'mf.VCoil_1'); model.result.numerical('gev1').set('descr', 'Coil voltage'); model.result.table.create('tbl1', 'Table'); model.result.table('tbl1').comments('Global Evaluation 1 (mf.VCoil_1)'); model.result.numerical('gev1').set('table', 'tbl1'); model.result.numerical('gev1').setResult; model.result.numerical('gev1').set('dataseries', 'maximum'); model.result.table.create('tbl2', 'Table'); model.result.table('tbl2').comments('Global Evaluation 1 (mf.VCoil_1)'); model.result.numerical('gev1').set('table', 'tbl2'); model.result.numerical('gev1').setResult; model.result.table('tbl1').clearTableData; model.result.table('tbl2').clearTableData; model.result.table.remove('tbl1'); model.result.table.remove('tbl2'); model.result.numerical.remove('gev1'); model.result.export.create('data1', 'Data'); model.result.export('data1').set('solnum', {'114'}); model.result.export('data1').set('expr', {'mf.Bx'}); model.result.export('data1').set('descr', {'Magnetic flux density, x component'}); model.result.export('data1').set('unit', {'T'}); model.result.export('data1').set('expr', {'mf.Bx' 'mf.By'}); model.result.export('data1').set('descr', {'Magnetic flux density, x component' 'Magnetic flux density, y component'}); model.result.export('data1').set('expr', {'mf.Bx' 'mf.By' 'mf.Bz'}); model.result.export('data1').set('descr', {'Magnetic flux density, x component' 'Magnetic flux density, y component' 'Magnetic flux density, z component'}); model.result.export('data1').set('filename', 'C:\Users\Annisanurf\Documents\SIMULASI\Latihan\crop.txt'); model.result.export('data1').run; model.result.export('data1').set('location', 'regulargrid'); model.result.export('data1').set('regulargridx3', '32'); model.result.export('data1').set('regulargridy3', '32'); model.result.export('data1').set('regulargridz3', '32'); model.result.export('data1').set('filename', 'C:\Users\Annisanurf\Documents\SIMULASI\Latihan\medan_port_1.txt'); model.result.export('data1').run; model.result.numerical.create('gev1', 'EvalGlobal'); model.result.numerical('gev1').set('descr', ' Coil voltage'); model.result.numerical('gev1').set('expr', 'mf.VCoil_2'); model.result.numerical('gev1').set('dataseries', 'maximum'); model.result.table.create('tbl1', 'Table'); model.result.table('tbl1').comments('Global Evaluation 1 (mf.VCoil_2)'); model.result.numerical('gev1').set('table', 'tbl1'); model.result.numerical('gev1').setResult; model.result.numerical('gev1').set('expr', 'mf.VCoil_3'); model.result.numerical('gev1').set('table', 'tbl1'); model.result.numerical('gev1').appendResult; model.result.numerical('gev1').set('expr', 'mf.VCoil_4'); model.result.numerical('gev1').set('table', 'tbl1'); model.result.numerical('gev1').appendResult; model.result.export.create('tbl1', 'Table'); model.result.export('tbl1').set('filename', 'C:\Users\Annisanurf\Documents\SIMULASI\Latihan\konduktivitas_udara_port_1.txt'); model.result.export('tbl1').run; out = model;