时隔一年……又要开始零基础学习comsol了……
今天学习的case:在不同位置加热物体的实现方法(谁知道这个case是哪年的呀
除了域的激活和去激活之外,这个例子还展示了其他
重要的建模技术,即
•如何使用框选择移动边界条件
•如何使用之前的解决方案作为初始条件在选定的领域。
一个简便的建模方式是首先在comsol中设置好相关传热entire,然后保存。mph文件,最后在matlab中进行加载。
在matlab程序中将会使用以下内容:
mphopen
mphglobal
mphplot
(主要是后处理部分的学习)
model = mphopen('domain_activation_llmatlab');
domInd = [2,3,5,4];
ht = model.physics('ht');
for i = 1:8
k = mod(i,4);
if k == 0
k = 4;
end
ht.selection.set([1 domInd(k)]);
ht.feature('init2').selection.set(domInd(k));
model.study('std1').run;
%画点
if i==1
cpt1 = model.result.dataset.create('cpt1', 'CutPoint3D');
cpt1.set('pointx', '0 L/2 L');
cpt1.set('pointy', '0 L/2 L');
cpt1.set('pointz', 'L/10');
pg1 = model.result.create('pg1', 'PlotGroup1D');
pg1.set('data', 'cpt1');
ptgr1 = pg1.feature.create('ptgr1', 'PointGraph');
ptgr1.set('legend', 'on');
pg2 = model.result.create('pg2', 'PlotGroup3D');
surf1 = pg2.feature.create('surf1', 'Surface');
surf1.set('rangecoloractive', 'on');
surf1.set('rangecolormax', '336');
surf1.set('rangecolormin', '293.15');
ht.feature('init1').set('T', 1, 'T');
v1 = model.sol('sol1').feature('v1');
v1.set('initsol', 'sol1');
end