COMSOL with MATLAB 应用实例（1）【天乐树】

mxio

建立几何体，本例简单介绍如何运用CwM画几何实体
目的：在一个方形阵列中每个单元画三个半径随机的圆。
效果：

实现方法：

1. % This is a script for Comsol(4.x) with Matlab.
   
2. % It product random radius circles inside a square box array.
   
3. %
   
4. % Copyright to mxio at Shanghai University.
   
5. %
   
6. %                                                        mxio
   
7. %                                                10,July,2013
   
8. 
   
9. 
   
10. 
    
11. import com.comsol.model.*
    
12. import com.comsol.model.util.*
    
13. 
    
14. model = ModelUtil.create('Model');
    
15. model.modelNode.create('mod1');
    
16. model.geom.create('geom1', 2);
    
17. 
    
18. % Model Parameters
    
19. rd_seed = rng('shuffle');               % random number seed
    
20. sq_length = 1;                          % square length
    
21. sq_posx = 0;                            % square position x
    
22. sq_posy = 0;                            % square position y
    
23. sq_disx = 1.5;                          % square displacement x
    
24. sq_disy = 1.5;                          % square displacement y
    
25. sq_numx = 5;                            % square number x
    
26. sq_numy = 5;                            % square number y
    
27. sq_num = sq_numx * sq_numy;             % square number total
    
28. 
    
29. cl_num = 3;                             % circle number inside the square
    
30. cl_radius_normal = sq_length * 0.45;    % circle radius : half for the maximum
    
31. co_min = 0.2;                           % random number must be > co_min
    
32. %------------------------------------------------------------------------------------------------------------
    
33. idx = 1;                                % index for square
    
34. idx2 = 1;                               % index for circle
    
35. for i = 1: sq_numx
    
36.     for j = 1:sq_numy
    
37.         sq_name = ['sq',num2str(idx)];
    
38.         model.geom('geom1').feature.create(sq_name, 'Square');
    
39.         model.geom('geom1').feature(sq_name).set('base', 'center');
    
40.         model.geom('geom1').feature(sq_name).set('size', num2str(sq_length));
    
41.         model.geom('geom1').feature(sq_name).setIndex('pos', num2str(sq_posx + (i-1)*sq_disx), 0);
    
42.         model.geom('geom1').feature(sq_name).setIndex('pos', num2str(sq_posx + (j-1)*sq_disy), 1);
    
43.             for k = 1:cl_num
    
44.                 cl_name = ['c', num2str(idx2)];
    
45.                 co_rand = rand();
    
46.                 while co_rand < co_min
    
47.                     co_rand = rand();
    
48.                 end
    
49.                 cl_radius = cl_radius_normal * co_rand;
    
50.                 model.geom('geom1').feature.create(cl_name, 'Circle');
    
51.                 model.geom('geom1').feature(cl_name).set('r', num2str(cl_radius));
    
52.                 model.geom('geom1').feature(cl_name).setIndex('pos', num2str(sq_posx + (i-1)*sq_disx), 0);
    
53.                 model.geom('geom1').feature(cl_name).setIndex('pos', num2str(sq_posx + (j-1)*sq_disy), 1);
    
54.                 idx2 = idx2 +1;
    
55.             end
    
56.         idx = idx + 1;
    
57.     end
    
58. end
    
59. model.geom('geom1').runAll;
    
60. %------------------------------------------------------------------------------------------------------------
    
61. 
    
62. mphsave(model,'rand_circle')

复制代码

mixo
2013.7.10

edodo

好东西  不过好像有些只有两个圆，如果根据随机数判断下，剔除这些就更完美了。

mxio

都是三环，只是半径相差不大而已。不知道应用背景的条件下，没必要写太复杂，只是拿来当教学案例罢了。。。。

longwen36

mxio 发表于 2013-7-10 20:49
都是三环，只是半径相差不大而已。不知道应用背景的条件下，没必要写太复杂，只是拿来当教学案例罢了。。。 ...

太赞了！但是不知道效率如何，没有测试。我做复合材料模拟时需要很多小球（成千上万）随机分散在基体中，按此方法不知道会不会死掉。

hhwapp

matlab不太会用

hhwapp

matlab不太会用，

zhaokw

好东西，学习学习

s1302563090

谢谢分享

等待的巧

谢谢分享。

你的微笑

谢谢分享。

hhm258

感谢分享，不错的

chendingshi

不错，谢谢分享

luozhu3

好帖子a!