在计算正方形晶格排列的二维声子晶体的能带结构的时候,为什么使用脚本求解出来的固有频率和在同样参数下运行一次COMSOL显示的固有频率不一样。脚本里的计算是f=i*lambda/(2*pi),在MT方向上kx=ky=0.25时,结果如下:
1.0e+005 *
0.2660 0.4348 0.4891 0.8547 0.9897 1.0394 1.1663 1.3375
而在同样的kx=ky=0.25,运行一次COMSOL,在“后处理--绘图参数--特征频率”下显示的频率却为:
说明一下,我用脚本画出来的能带图是与别人用平面波计算出来的完全相符。不光是在kx=ky=0.25,在其他的好些点,COMSOL运行一次显示的结果都与脚本计算的不一样,不知道为什么,求高手相助!
如果有同类问题的,也可以讨论一下!
————————————————————————————————————————————————————————
【解决办法】说起来这只是一个小的是设置问题造成的,不过也困扰了我一两天。出现这个问题需要选中"求解--求解器参数--高级--使用约束矩阵厄米特变换和对称检测"。
https://blog.sciencenet.cn/blog-506324-505767.html
上一篇:
“红月亮”图解下一篇:
COMSOL使用问题2:散射声场【已解决】