COMSOL与MATLAB.pdf

COMSOL与MATLAB.pdf

《COMSOL与MATLAB.pdf》由会员分享，可在线阅读，更多相关《COMSOL与MATLAB.pdf（54页珍藏版）》请在冰豆网上搜索。

VERSION4.4LiveLinkforMATLABTM简介?

Partnumber:

CM020010LiveLinkforMATLAB?

20092013COMSOL?

COMSOL4.4|i?

ii|?

|1?

2|?

COMSOLwithMATLAB?

COMSOLwithMATLAB?

comsolservermatlab?

COMSOLwithMATLAB|3?

“?

”?

4|?

2a?

1?

2b?

“?

”?

“?

”?

|5?

blk=geom.feature.create（blk,Block）;blk.setIndex（size,2,0）;blk.setIndex（size,3,1）;?

geom.feature.create（blk,Block）;geom.feature（blk）.setIndex（size,2,0）;geom.feature（blk）.setIndex（size,3,1）;?

1?

“?

COMSOLwithMATLAB”?

2?

model=ModelUtil.create（Model）;?

6|?

“?

”?

1?

2?

3?

ModelUtil.setServerBusyHandler（ServerBusyHandler

（2）;?

|71?

model.param.set（L,9cm,Lengthofthebusbar）;model.param.set（rad_1,6mm,Radiusofthefillet）;model.param.set（tbb,5mm,Thickness）;model.param.set（wbb,5cm,Width）;model.param.set（mh,6mm,Maximumelementsize）;model.param.set（htc,5W/m2/K,Heattransfercoefficient）;model.param.set（Vtot,20mV,Appliedelectricpotential）;?

1?

geom1=model.geom.create（geom1,3）;?

2?

wp1=geom1.feature.create（wp1,WorkPlane）;wp1.set（quickplane,xz）;3?

r1=wp1.geom.feature.create（r1,Rectangle）;r1.set（size,L+2*tbb0.1）;?

4?

r2=wp1.geom.feature.create（r2,Rectangle）;r2.set（size,L+tbb0.1-tbb）;r2.set（pos,0tbb）;5?

dif=wp1.geom.feature.create（dif,Difference）;dif.selection（input）.set（r1）;dif.selection（input2）.set（r2）;8|?

geom1.run;?

6?

fil1=wp1.geom.feature.create（fil1,Fillet）;fil1.selection（point）.set（dif

（1）,3）;fil1.set（radius,tbb）;7?

fil2=wp1.geom.feature.create（fil2,Fillet）;fil2.selection（point）.set（fil1

（1）,6）;fil2.set（radius,2*tbb）;?

|98?

ext1=geom1.feature.create（ext1,Extrude）;ext1.selection（input）.set（wp1）;ext1.set（distance,wbb）;?

mphgeom（model）?

9?

wp2=geom1.feature.create（wp2,WorkPlane）;wp2.set（planetype,faceparallel）;wp2.selection（face）.set（ext1

（1）,8）;10?

c1=wp2.geom.feature.create（c1,Circle）;c1.set（r,rad_1）;11?

ext2=geom1.feature.create（ext2,Extrude）;ext2.selection（input）.set（wp2）;ext2.set（distance,-2*tbb）;10|?

12?

wp3=geom1.feature.create（wp3,WorkPlane）;wp3.set（planetype,faceparallel）;wp3.selection（face）.set（ext1

（1）,4）;13?

c2=wp3.geom.feature.create（c2,Circle）;c2.set（r,rad_1）;c2.set（pos,-L/2+1.5e-2-wbb/4）;14?

copy=wp3.geom.feature.create（copy,Copy）;copy.selection（input）.set（c2）;copy.set（disply,wbb/2）;15?

ext3=geom1.feature.create（ext3,Extrude）;ext3.selection（input）.set（wp3.c2wp3.copy）;ext3.set（distance,-2*tbb）;16?

geom1.run;?

|11?

1?

sel1=model.selection.create（sel1）;sel1.set（234567）;sel1.name（Tibolts）;2?

mphviewselection（model,sel1）;12|?

5.998e7S/m7.407e5S/m?

385J/（kg*K）710J/（kg*K）?

11?

8700kg/m34940kg/m3?

400W/（m*K）7.5W/（m*K）1?

mat1=model.material.create（mat1）;?

COPPERTITANIUM?

|132?

mat1.materialModel（def）.set（electricconductivity,5.998e7S/m）;mat1.materialModel（def）.set（heatcapacity,385J/（kg*K）;mat1.materialModel（def）.set（relpermittivity,1）;mat1.materialModel（def）.set（density,8700kg/m3）;mat1.materialModel（def）.set（thermalconductivity,400W/（m*K）;3?

mat1.name（Copper）;4?

mat2=model.material.create（mat2）;mat2.materialModel（def）.set（electricconductivity,7.407e5S/m）;mat2.materialModel（def）.set（heatcapacity,710J/（kg*K）;mat2.materialModel（def）.set（relpermittivity,1）;mat2.materialModel（def）.set（density,4940kg/m3）;mat2.materialModel（def）.set（thermalconductivity,7.5W/（m*K）;mat2.name（Titanium）;5?

mat2.selection.named（sel1）;?

1?

ht=model.physics.create（ht,HeatTransfer,geom1）;14|?

2?

hf1=ht.feature.create（hf1,HeatFluxBoundary,2）;?

3?

hf1.set（HeatFluxType,InwardHeatFlux）;hf1.selection.set（1:

79:

1416:

42）;hf1.set（h,htc）;4?

mphgeom（model,geom1,facemode,off,facelabels,on）?

5?

ec=model.physics.create（ec,ConductiveMedia,geom1）;6?

pot1=ec.feature.create（pot1,ElectricPotential,2）;pot1.selection.set（43）;pot1.set（V0,Vtot）;7?

gnd1=ec.feature.create（gnd1,Ground,2）;gnd1.selection.set（815）;?

|15?

1?

emh=model.multiphysics.create（emh,ElectromagneticHeatSource,geom1,3）;emh.selection.all;2?

emh.set（EMHeat_physics,ec）;emh.set（Heat_physics,ht）;?

1?

mesh=model.mesh.create（mesh,geom1）;2?

size=mesh.feature（size）;size.set（hmax,mh）;size.set（hmin,mh-mh/3）;size.set（hcurve,0.2）;3?

ftet=mesh.feature.create（ftet,FreeTet）;16|?

4?

mesh.run;?

5?

mphmesh（model）?

1?

std=model.study.create（std）;stat=std.feature.create（stat,Stationary）;?

|172?

ModelUtil.showProgress（true）;?

3?

std.run;?

18|?

1?

pg=model.result.create（pg,PlotGroup3D）;2?

surf=pg.feature.create（surf,Surface）;surf.set（expr,T）;?

3?

surf.set（rangecoloractive,on）;surf.set（rangecolormin,322.6）;surf.set（rangecolormax,323.3）;4?

mphplot（model,pg,rangenum,1）?

|19?

“?

”?

1?

data=model.result.export.create（data,Data）;data.setIndex（expr,T,0）;data.set（filename,filepathTemperature.txt）;data.run;20|?

%Model:

%Version:

COMSOL4.4.0%Date:

Sep302013,11:

25%Dimension:

3%Nodes:

1700%Expressions:

1%Description:

Temperature%Lengthunit:

m%xyzT（K）0.09999999999999999-0.0038380818225469850.09639433064874708322.949889247876060.09749606929798824-0.0069997502975423650.1322.949261374022460.09999999999999996-0.0055555555555555560.1322.94824689779085?

1?

mphsave（model,/busbar）;?

mphsave（model,/busbar.m）;?

|21?

MATLAB?

“?

”?

1?

2?

“?

COMSOLwithMATLAB”?

3?

22|?

4?

-?

COMSOL4.4startedlisteningonport2036-?

5?

6?

model=ModelUtil.model（Model2）;?

ModelUtil.tags7?

mphmesh（model）8?

model.result.remove（pg2）;?

|