日志
通过一个例子来谈论关于模型耦合的问题
天气: 晴朗
心情: 高兴
通过一个例子来谈论关于模型耦合的问题在一个隔声量仿真的过程中，对SYSNOISE的模型偶合有了更加深刻的感悟，下面我把我试的过程中的心得说一下：
假设是在一个长1m，直径50mm的驻波管中间位置放置一个隔声板（建模时用的是壳单元），板的直径和管的直径相同。在管的最左端有一个声源，在管的最右端是全吸声的劈尖。管其他部位都是全反射的刚性边界。
就这个模型，目的是想计算板的隔声量。
下边依次列出我所采用的模型及出现的问题，以及出现问题的原因。
一。FEM FLUID +FEM STRUCTRUE +FEM FLUID
起初我认为中间的FEM STRUCTRUE 可以同时和两边的FEM FLUID耦合，但是失败了，SYSNOISE根本就无法运算。后来思考了一下，在耦合中，由于FEM FLUID和FEM STRUCTRUE同时计算，而且共用一个矩阵，在FEM STRUCTRUE 与其中一个 FEM FLUID耦合后，再与另外一个FEM FLUID耦合，势必它的结果会被覆盖，不知道这样的考虑是否正确，但是的确是运行不了。
二。DBEM FLUID COUPLED+FEM STRUCTRUE +DBEM FLUID COUPLED
接下来将FEM FLUID换成DBEM FLUID，同样是不能计算，原因我想应该和上边的一样，这里不在重复说了。三。DBEM FLUID 。
在尝试了一个FEM STRUCTRUE 的壳单元模型不能同时和两个流体单元耦合后，我放弃了这个念头。开始思考分步操作，先提取没有板时板所在位置的声压数据，再将此数据双倍施加在板和后半段管耦合模型的板上。
三。DBEM、FEM STRUCTRUE+DBEM FLUID COUPLED
先建立整个管的DBEM模型（之所以采用DBEM模型而不用FEM模型，是因为在用FEM模型时，加面声源计算出来的结果很离谱，听说是SYSNOISE的BUG吧），提取中间管所在位置的声压数据，存成RESULT文件。然后将板的FEM STRUCTRUE与管的后半段的空气模型DBEM FLUID COUPLED进行耦合，将得到的RESULT文件作为压力边界施加在板模型上，在DBEM FLUID COUPLED和FEM STRUCTRUE耦合时，发现下面的问题。
其他所有的耦合都可以通过LINK完成，但是唯独DBEM FLUID 和FEM STRUCTURE的耦合要通过在MODAL\OPTION里选择COUPLE来耦合，因为在DBEM FLUID 和FEM STRUCTURE的耦合中，DBEM FLUID 和FEM STRUCTURE共用一个模型（或者是模型中的一部分），也就是说你只需要导进一个模型就可以了，SYSNOISE会认为在同一个模型中，既有结构模型又有流体模型，并且两者在所有面上是耦合关系。若要模型在部分面上耦合，则需要定义所要耦合的面，通过GEOMETRY\INTERFACE定义耦合面。但是不管我怎样设置，依然不能计算，现在还没弄清楚原因所在。
四。DBEM、FEM STRUCTRUE+FEM FLUID
于是将DBEM FLUID COUPLED换成了FEM FLUID，这样得到的结果还比较理想。但是有个问题，FEM计算出来的结果通过COLOR MAP看后半段空气柱表面的声压云图时，按理论上来说，空气柱表面没有任何边界条件设置，默认的应该是刚性的全反射边界，表面振速法向振速应该为零，但是在云图观察空气柱表面时，空气柱还存在驻波现象，但是提取内部的节点声压绘成曲线时，所有的点的变化非常小，可能这也是FEM的一个BUG吧。
五。FEM STRUCTRUE+IBEM BAFFLED
利用这样的模型的耦合计算出来的结果和四中计算出来的结果非常吻合。
说了好多，可能看的人也会觉得累，想说的太多了，有点乱，呵呵！
