COMSOL4及链接Matlab使用心得！

kyolee · 发表于 2010-10-11 08:28:55

可能你的版本是4.0吧。升级到4.0a应该就有了

zhangkai0503 · 发表于 2010-10-25 19:39:27

好贴子，要收藏

anweihit · 发表于 2010-11-15 16:30:21

感谢楼主分享！不知道有没有与matlab一起配合的例子共享一下，那就更好了。

sjtueric · 发表于 2010-11-23 16:01:26

楼主,你好,我想请教一下关于comsol能否解决1个电路+多个电磁场模型?请联系我

Q385753429.谢谢了

ansys平安 · 发表于 2010-12-11 19:35:07

学习！学习！

petro2009 · 发表于 2010-12-16 11:28:12

感谢分享

sjtueric · 发表于 2010-12-16 21:16:33

不知道楼主做过comsol3.5连接simulink的模型没，我发现comsol中输出的结果正确，但是到了simlink中输出的结果就不一样了，不知道哪出错了，希望楼主指教，我的QQ号：1643026252

fn001cn · 发表于 2011-4-20 18:41:50

强人，加紧练习练习，看看怎么弄好他。

fn001cn · 发表于 2011-4-20 18:45:09

强人，加紧练习练习，看看怎么弄好他。

changyuhui · 发表于 2011-11-10 23:47:27

谢谢楼主的无私分享！！！

初学者，非常需要

noncall · 发表于 2011-11-11 18:22:06

学习了,谢谢无私分享

lshuai8002 · 发表于 2011-11-14 23:56:12

不错，共同学习

jj201314 · 发表于 2011-11-15 14:30:41

刚学习comsol，请问楼主COMSOL with Matlab 怎么打开？username和password输入什么？谢谢~~

zhouteng · 发表于 2011-11-18 20:52:52

感觉还是不如3.5script好用啊。大家有没讲comsoL with matlab的书籍或者视频啊？

complax · 发表于 2012-3-13 11:12:45

感谢分享体会

fessorji · 发表于 2012-3-13 21:48:25

学习学习不错

PursueDream · 发表于 2012-10-26 19:00:38

好贴学习了

jing_miao · 发表于 2012-11-1 16:32:19

楼主高手

zheng000012 · 发表于 2012-11-12 08:57:13

我下载了一段程序，在3.5中运行的，可是运行总是出错，程序如下，希望楼主能够帮着看看，谢谢

clear all, flclear fem, close all
disp('Fem2D_SymmetricBox.m');
N=2;
cWa=1483;
cSi=8490;
rhoWa=998;
rhoSi=2331;
alpha=2;
L=55.00e-3;
l=29.30e-3;
w=380e-6;
W=w*(1+cSi/cWa*alpha);
f0=cWa/2/w;
fs=f0;
points_mesh=60;
points=100;
y0=0.0;
paX=+1;
paY=+1;
if paX==-1
BCX=1;
else
BCX=0;
end
if paY==-1
BCY=1;
else
BCY=0;
end
const=struct('cWa',cWa,'cSi',cSi,'rhoWa',rhoWa,'rhoSi',rhoSi,...
'w',w,'l',l,'W',W,'L',L,'y0',y0,'BCX',BCX,'BCY',BCY);
fem.const=const;
g1=rect2(w/2,l/2,'base','corner','pos',[0 0]);
g2=rect2(W/2,L/2,'base','corner','pos',[0 0]);
s.objs={g1,g2};
fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
fem.equ.ind={1 2};
fem.bnd.ind={[4 6] [5 8] [1 3] [2 7]};
fem.sdim={'x','y'};
fem.dim={'p'};
fem.shape=2;
fem.form='general';
fem.equ.expr={...
'c',{'cWa','cSi'},...
'rho',{'rhoWa','rhoSi'},...
'Eaco',{'((abs(px)^2+abs(py)^2)/omega^2+abs(p)^2/c^2)/(4*rh0)'}};
fem.equ.shape={1 1};
fem.equ.ga={{'px/rho';'py/rho'}};
fem.equ.f={{'-omega^2*p/(rho*c^2)'}};
fem.equ.init={{'0'}};
fem.bnd.shape={1,1,1,1};
fem.bnd.r={...
{'0'},...
{'0-p'},...
{'0-BCX*p'},...
{'0-BCY*p'}};
fem.bnd.g={...
{'0'},...
{'0'},...
{'0'},...
{'0'}};
mSi=cSi/f0/points_mesh;
mWa=cWa/f0/points_mesn;
fem.mesh=meshinit(fen,'report','0ff','hmaxsub',[1 2;mWa mSi]);
fem.xmesh=meshextend(fem);
DOF=flngdof(fem);
omega_vec=fem.sol.lambda;
Ewa=zeros(1,length(omega_vec));
Esi=Ewa;
xvwa=(0+mWa/points):mWa/points:(w/2-mWa/points);
xvsi=(w/2+mSi/points):mSi/points:(W/2-mSi/points);
Edens_wa=zeros(length(omega_vec),length(xvwa));
Edens_si=zeros(length(omega_vec),lenghth(xvsi));
for j=1:length(omega_vec)
sfigure(3);
postplot(fem,'tridata','p','solnum',j)
xlabel('x'),ylabel('y')
axis equal
yv0wa=y0*ones(1,length(xvwa));
yv0si=y0*ones(1,length(xvsi));
Edens_wa(j,:)=postinterp(fem,'Eaco',[xvwa;yv0wa],...
'solnum',j,'D1',2,'edim',2);
Edens_si=postinerp(fem,'Eaco',[xvsi;yv0si],...
'solnum',j,'D1',2,'edim',2);
EmaxWa=max(Edens_wa(j,:));
EmaxSi=max(Edens_si(j,:));
Awa=w8l;
Asi=W*L-w*l;
Ewa(j)=1/(Awa/4)*postint(fem,'Eaco','D1',1,'edim',2,'solnum',j);
Esi(j)=1/(Asi/4)*postint(fem,'Eaco','D1',2,'edim',2,'solnum',j);
p_wa=postinterp(fem,'p',[xvwa;yv0wa],'solnum',j);
p_si=postinterp(fen,'p',[xvsi;yv0si],'solnum',j);
sfigure(4);
plot(xvwa/1e-3,p_wa/max(p_si),'-r',xvsi/1e-3,p_si/max(p_si),'-b');
xlabel('y [mm]'), ylabel('p/pmax [Pa]')
delta=cWa/(omega_vec(j)/2/pi)/2/w-1;
fprintf(1,'f0:%6.4f MH,DOF:%1,alpha:%2.3f\n',f0/1e6,DOF,alpha);
fprintf(1,['N0.:%i/%i,f:%6.4f MHz,delta:%6.8f,'...
'EmaxWa/EmaxSi:%2.3f,Ewa/Emaxsi:%2.3f\n'],...
j,length(omega_vec),omega_vec(j)/2/pi/1e6,delta,...
EmaxWa/EmaxSi,Ewa(j)/EmaxSi);
pause
end

复制代码

fengerccsyk · 发表于 2012-11-14 16:54:48

感觉两个联合计算效率好低

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[15.分享让我牛B] COMSOL4及链接Matlab使用心得！