为了赶时间，开始研究comsol script了。手工算比较费时间，一个for循环解决很大问题，不过解的区分度不高，还得需要借助人工。网络上有两个文章较好，一个是“Primer on Scripting Comsol With Matlab”，一个是“Run Comsol Simulations from Matlab”。自己google搜一下吧，另外一些人的毕业论文里也有Comsol script的代码，自己找一找，虽然方向可能不一样，但是也会有所帮助。
20090423更新
comsol自己带了一个文件mlinterface.pdf，在comsol安装目录的doc/multiphysics文件夹下找，比较全面的介绍了script。
COMSOL Multiphysics允许用户自定义坐标系统，我们可以使用软件自带的函数进行坐标系统设定。例如：
1.3D球坐标系
选择Options > Expressions > Scalar Expressions，在对话框中输入：
r＝sqrt(x^2+y^2+z^2)
phi＝atan2(y,x)
theta＝acos(z/r)
2.2D极坐标
选择Options > Expressions > Scalar Expressions，在对话框中输入：
r＝sqrt(x^2+y^2)
phi＝atan2(y,x)
sqrt, atan2及acos 均为comsol內建函数
有人在comsol 交流群里发出来的，如果涉及版权问题，可以回复，我加上作者。
找到了，是在中仿科技的网站里的comsol faq ，http://www.cntech.com.cn/faq/comsol.html
在师姐的帮助下，终于计算出来了一种结构的一些解。原来这东西没有那么神秘，自己慢慢来。
Do it, it’s simple!!
COMSOL Multiphysics默认的长度单位是m，如果想在其他单位制下画图,可以选择Options>Preferences,打开Preferences对话框,在Modeling中的Default unit system中修改.
如果模型已经画后,想修改模型大小,可以选择draw>modify>scale,打开scale对话框,在Scale factor中输入要放缩的比例因子,在Scale base point中,输入放缩的基点.比如说,要把模型在坐标轴方向缩小10e3,可以在Scale base point中,输入坐标圆点位置,在Scale factor中每个坐标方向输入10e-3.
对有限元计算，无论是ansys，abaqus，msc还是comsol等，归结为一句话就是解微分方程。而解方程要有定解，就一定要引入条件， 这些附加条件称为定解条件。
定解条件的形式很多，只讨论最常见的两种——初始条件和边界条件。
初值和边值问题：
对一般的微分方程，求其定解，必须引入条件，这个条件大概分两类—初始条件和边界条件，如果方程要求未知量y(x)及其导数y′(x)在自变量的同一点x=x0 取给定的值，即y(x0 )=y0 ，y′(x0 )=y0′，则这种条件就称为初始条件，由方程和初始条件构成的问题就称为初值问题；
而在许多实际问题中，往往要求微分方程的解在在某个给定的区间a ≤ x ≤b的端点满足一定的条件，如y(a) = A , y(b) = B 则给出的在端点(边界点)的值的条件，称为边界条件，微分方程和边界条件构成数学模型就称为边值问题。
边值问题中的边界条件的形式多种多样，在端点处大体上可以写成这样的形式，Ay+By'=C，若B=0，A≠0,则称为第一类边界条件或狄里克莱(Dirichlet)条件；B≠0,A=0，称为第二类边界条件或诺依曼(Neumann)条件；A≠0,B≠0，则称为第三类边界条件或洛平(Robin)条件。
总体来说，
第一类边界条件：给出未知函数在边界上的数值；
第二类边界条件：给出未知函数在边界外法线的方向导数；
第三类边界条件：给出未知函数在边界上的函数值和外法向导数的线性组合。
对应于comsol，只有两种边界条件：
Dirichlet boundary（第一类边界条件）—在端点，待求变量的值被指定。
Neumann boundary（第二类边界条件）—待求变量边界外法线的方向导数被指定。
