COMSOL Multiphysics 脚本功能详述.ppt

1、COMSOL Multiphysics 脚本功能详述,中仿科技,议程安排,COMSOL vs. MATLAB 脚本相互调用 二次开发功能 LiveLink for MATLAB JAVA 示例 讨论,Livelink for Pro/ENGINEER ,Livelink for Autodesk Inventor,CAD Import Module,Livelink for MATLAB,优化 Optimization Module,材料库 Material Library,Livelink for SolidWorks ,COMSOL Multiphysics,低频电磁场 AC/DC Mod

2、ule,高频电磁场 RF Module,微机电 MEMS Module,等离子体 PLASMA MODULE,结构力学 Structural Mechanics Module,声学 Acoustics Module,传热 Heat Transfer Module,计算流体力学 CFD MODULE,化学反应工程 CHEMICAL REACTION ENGINEERING MODULE,电池与燃料电池 BATTERIES model.save(test,M);,从模型对象中提取信息到命令行 查阅COMSOL Multiphysics API参考指南，了解针对模型对象的完整函数列表 MPHMATR

3、IX集成刚度矩阵,从模型对象中提取信息,Demo: 在命令行中建模,在命令行中建立一个3D台阶流动模型,首先建立一个模型对象 添加一个3D几何 在几何中添加一个圆柱,在MATLAB命令行中设定模型,设定物理场 创建网格,在MATLAB命令行中建模,求解 创建绘图组 可以在MATLAB中绘图 或在服务器中绘图（需要带选项启动：-graphics）,在MATLAB命令行中建模,修改脚本进行高级求解运算,控制流命令 if条件分支语句 switch多重条件分支 for循环 while条件循环语句 中断命令 break跳出循环 continue继续下一步循环 return从函数返回 可能用到的函数 sp

4、rintf、eval等格式转换函数,Demo：生成随机图形,利用MATLAB脚本生成随机图形 运行Demo脚本 model = randcir; mphgeom(model,geom1),Demo: 在COMSOL GUI中调用MATLAB脚本,在材料属性、边界条件等中使用MATLAB脚本函数 不需要启动COMSOL with MATLAB 被调用的脚本必须能够被MATLAB直接使用（注意路径） 推荐加入到目录列表中,Demo：MATLAB中后处理,计算网格节点上的结果 mpheval pd=mpheval(model,T,Edim,1); 各种后处理积分 mphint val=mphint(

5、model,T),计算点的结果 mphinterp T = mphinterp(model,T,coord,0.3 0.4;0.5 0.5),Demo:获得总刚矩阵,mphmatrix model = mphload(point_source.mph) MA=mphmatrix(model,sol1,out,K,L); spy(MA.K),Demo：二次开发界面,推荐步骤 从GUI建立简单模型 另存为M文件 在MATLAB中开发界面 将M文件拆分，插入以MATLAB界面命令脚本,兼容运行模式,启动COMSOL with MATLAB (3.5a compatibility) 以FEMGUI类型

6、设定FEM结构体 运行3.5a M文件，调用 MPHV4 将 FEMGUI 结构体转换成 COMSOL 4.0 模型对象,MPHPLOT 不支持在 MATLAB 绘制所有的绘图组，可使用服务器方式 comsolserver matlab -graphics 在COMSOL GUI 和 MATLAB 之间交换数据只能采用保存文件后的导入/导出方式,局限性,JAVA API,提供应用程序接口（API） 采用JAVA编译器 多种求解运行方式 更加灵活的开发和应用 并不仅仅局限于MATLAB,JAVA API使用方法,编译成Java的库 comsolcompile -jdkroot 运行求解 COMS

7、OL桌面 菜单项文件打开，选择Model Class File (*.class) COMSOL桌面中的批处理任务 在工作序列中右键添加外部类 命令行方式 comsolbatch -inputfile -outputfile -batchlog ,格式化字符串的应用,sprintf-生成包含指令的字符串 eval-执行上述字符串 Demo str=sprintf(model.mesh(mesh1).feature(size).set(hauto, %d); , 10); eval(str);,参考资料,LiveLink for MATLAB User Guide of LiveLink for MATLAB JAVA COMS