请教sun老师（ansys模拟消声器）

古流新用

网格不好，可以自己设定lesize，细致划分，主要是几何模型要做好，分网就比较方便。在将模型时就要考虑到分网的方便。

  
确实，看了以上的命令流，感受了一下区别。有很大体会。
补充lesize命令流：
lesize,1,,,11,1
lesize,2,,,1,1
lesize,3,,,7,1
lesize,4,,,18,1
lesize,5,,,12,1
lesize,6,,,8,1
lesize,7,,,18,1
lesize,8,,,2,1
lesize,9,,,18,1
lesize,10,,,8,1
lesize,11,,,8,1
lesize,12,,,10,1
lesize,13,,,2,1
lesize,14,,,10,1
lesize,15,,,1,1
lesize,16,,,3,1
lesize,17,,,7,1
lesize,18,,,1,1
lesize,19,,,1,1
lesize,20,,,3,1
lesize,21,,,8,1
lesize,23,,,4,1
lesize,24,,,4,1
lesize,25,,,6,1
  
划完，网格如下附件

古流新用

下面定义fluid129单元
csys,1
nsel,s,loc,x,2
type,3
real,3
mat,2
esurf
alls
定义流体-结构接触面
这里怎么点选接触面的节点？
因为reflect的只是节点吧？所以开始我想选lines,再attached的方法是不行的
我现在想了一个办法，就是通过选type1，list其所有节点，查找xy坐标符合的，选出。这个方法，一、不一定对 ；二、比较费时。请教怎么点选？
还有就是不是命令流的方法了，选出type1，在选nodes，apply用鼠标一点一点选。
  
选那种分析？harmonic?这个问题问得我都……不要见笑呀……
解决完这个问题，我要加载入口声压，这个我看了一下书，是否就是在入口处的节点上设置压强呢？（有专门的声压这项吗？还是就是一般的压强）

lypzsl

     没细看你的结构，说的参考一下吧！
1。定义FSI可以直接定义在线上，不知你的模型可不可行！那样就不用选择节点了。
2。如果你想作声压图，当然是谐波分析（harmonic）。载荷就应该是Pressure，

lypzsl

还有一事：
  
你的模型应该可以建成二分之一吧？你现在的是一个整体的圆面，而且使用轴对称单元，我怀疑是否正确？有时间我试试！
看看sun老师给解释一下！

古流新用

谢谢！

古流新用

补充下面的命令流：（有些不会写，只会gui）
  
这步暂时没命令流，用gui，选element type1，选entitle,nodes apply,鼠标点选出腔体与流体的耦合面。
esel,s,type,,2
sf,all,fsi,1
alls
finish
  
/solu
antype,harmic
hropt,full
  
Main Menu>Solution>Define Loads>Apply>Fluid/ANSYS>Fluid-Struct>On Nodes    选出入口处的7个nodes（这步完成后，看到了绿色箭头指向腔内）
  
harfrq,0,3000
nsubst,100
kbc,1
solve
finish
  
solve完成后，想出等值声压图，但作不出来，如下：
plnsol,pres(异或使用gui，都出某种freedom ,not avail的提示语）
我推测的原因，我看了节点求解结果的主要分类，有dof solution,stress,strain-total,nodal force data,energy,等。没有对应的压强呀。
请教，解法！谢谢！

lypzsl

声压在dof solution－pres，就是点dof solution后，右侧菜单的最下面二行。
  你的freedom not avail 的提示语，问题可能还出在流体的建模上。

古流新用

solve 时出waring.
coefficient ratio exceeds 1.0e8,check the results.
  
后处理时，用plnsol或按以上的gui方式（dof solution--press,确实是倒数第二项）。得到warning.
The degree of freedom solution is not avaiable .The PLNSOL command is ignored.
  
依以上warning，我自己觉得：
提示DOF的计算没有，证明算选的计算没计算声压，是不是计算类型错了？
  
还有，我看云图没出，我就想用list列表的方式试试，结果还是没有。

古流新用

是不是算得不收敛呀？
我施加的dof press，却得不出dof press的等值图，所以应该是算的问题，
如果是不收敛，怎么改？改变边界条件，改模型，还有什么方法？

古流新用

请教，如上！
还有，SF,all,FSI,1  
其中，最后的1,我看英文帮助，是interface number，是指面的数量（比如，1个面），还是面的编码（编码为一的面）？

sunhaolan

搞得好麻烦，索性全部给你做完。
载荷消声器入口处声压100dB，也就是2Pa。结果显示为扩展为3/4整体模型80Hz分析得到的分贝图（dB）。
/prep7
et,1,plane42  
et,2,fluid29  
et,3,fluid129  
r,2,2e-5
r,3,2,0,0  
et,4,fluid29,,1,0  
  
mp,ex,1,2.10e11  
mp,dens,1,7600  
mp,nuxy,1,0  
mp,dens,2,1.22  
mp,sonc,2,1460  
  
k,1,0,0.55
k,2,0.9,0.55
k,3,1,0.6
k, 4,,0.6
k,5,2  
k,6,1.5
k,7,1.5,0.2
k,8,1,0.2
k,9,0,2
k,10,0.9,0.6
k,12,0.9,0.15
k,13,1,0.15
k,15
k,16,0.9
k,17,1.5,0.15
  
l,1,15   
l,1,4   
l,4,9   
l,15,16  
l,16,6   
l,6,5   
l,4,10   
l,10,3   
l,1,2   
l,2,12   
l,3,8   
l,8,7   
l,12,13  
l,13,17  
l,7,17   
l,17,6   
larc,5,9,3,2
l,2,10   
l,8,13   
l,12,16
  
al,1,4,20,10,9   
al,20,5,16,14,13
al,2,9,18,7  
al,18,10,13,19,11,8  
al,14,15,12,19  
  
k,100,2,2
l,3,100
k,101,2,0.2
l,7,101
lsel,s,,,21,22
CM,sub,LINE
alls
lsbl,17,sub,,delete,keep
lsbl,21,25,,delete,keep
lsbl,22,25,,delete,keep
ldele,26,27,1,1
  
al,3,7,8,17,24   
al,17,25,21,12,11   
al,15,16,6,23,21
  
lesize,1,,,9,1  
lesize,2,,,1,1  
lesize,3,,,11,1  
lesize,4,,,10,1  
lesize,5,,,8,1  
lesize,6,,,5,1  
lesize,7,,,10,1  
lesize,8,,,2,1  
lesize,9,,,10,1  
lesize,10,,,6,1  
lesize,11,,,6,1  
lesize,12,,,6,1  
lesize,13,,,2,1  
lesize,14,,,6,1  
lesize,15,,,1,1  
lesize,16,,,3,1  
lesize,17,,,11,1  
lesize,18,,,1,1  
lesize,19,,,1,1  
lesize,20,,,3,1  
lesize,21,,,5,1  
lesize,23,,,4,1  
lesize,24,,,6,1  
lesize,25,,,6,1  
  
asel,s,,,1,2
asel,a,,,6,8
aatt,2,2,2   
LCCAT,7,8
MSHKEY,1
AMESH,6
  
LCCAT,12,21
MSHKEY,1
AMESH,7
  
LCCAT,15,16
MSHKEY,1
AMESH,8
  
LCCAT,10,20
MSHKEY,1
AMESH,1
  
LCCAT,13,14
MSHKEY,1
AMESH,2
alls
ldele,26,29  
ldele,22
  
asel,s,,,3,5
aatt,1,,1   
mshkey,1
amesh,3
amesh,5
lccat,10,18
lccat,11,19
mshkey,1
amesh,4
alls
  
esel,s,type,,1  
nsle,s  
esln,s,0  
nsle,s  
esel,inve  
nsle,s  
emodif,all,type,4
alls
  
wprot,45
CSYS,4   
ARSYM,Y,all, , , ,0,1   
wprot,-45
  
nsym,y,1000,all  
esym,,1000,all  
nummerg,all  
  
csys,1  
nsel,s,loc,x,2  
type,3  
real,3  
mat,2  
esurf  
alls  
csys,0
  
nsel,s,loc,y,-1.5,-.9   
nsel,r,loc,x,.15
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-1.5,-1   
nsel,r,loc,x,.2
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-.9,.9   
nsel,r,loc,x,.55
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-1,1   
nsel,r,loc,x,.6
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,.9,1.5   
nsel,r,loc,x,.15
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,1,1.5  
nsel,r,loc,x,.2
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-1.5  
nsel,r,loc,x,.15,.2
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-1  
nsel,r,loc,x,.2,.6
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-.9  
nsel,r,loc,x,.15,.55
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,.9  
nsel,r,loc,x,.15,.55
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,1  
nsel,r,loc,x,.2,.6
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,1.5  
nsel,r,loc,x,.15,.2
esel,s,type,,2   
sf,all,fsi,1
alls
  
/solu  
antype,harmic  
hropt,full  
nsel,s,loc,y,-1.5   
nsel,r,loc,x,,.15   
d,all,pres,2
alls
HARF,80,80
SOLVE
  
/post1
SET,FIRST
esel,s,type,,2   
esel,a,type,,4   
PLESOL, NMISC,4, 0,1.0
/EXPAND,27,AXIS,,,10
  
其中注意模型轴对称方向错误，已经修正。

sunhaolan

古流新用 wrote:
solve 时出waring.
coefficient ratio exceeds 1.0e8,check the results.
  
后处理时，用plnsol或按以上的gui方式（dof solution--press,确实是倒数第二项）。得到warning.
The degree of freedom solution is not avaiable .The PLNSOL command is ignored.
  
依以上warning，我自己觉得：
提示DOF的计算没有，证明算选的计算没计算声压，是不是计算类型错了？
  
还有，我看云图没出，我就想用list列表的方式试试，结果还是没有。

没有结果是因为FSI耦合载荷没有加上，我的经验是：注意FSI载荷加在几何模型上是不行的，起不了作用。
一定要选中耦合流体单元和节点，再加FSI载荷。
其实选择的问题不复杂，也不难，从坐标选取就可以了。
coefficient ratio exceeds 1.0e8,check the results.
上面这个错误是正常的，提示你不同材料之间的物性相差太大，流体与结构之间的物性相差很大是正常的。不是错误，只是一个提示。

sunhaolan

古流新用 wrote:
请教，如上！
还有，SF,all,FSI,1  
其中，最后的1,我看英文帮助，是interface number，是指面的数量（比如，1个面），还是面的编码（编码为一的面）？

1是指界面的数量。不是面的编码，也不是数量。
我们所用的就是一层界面，所以永远都是1。

古流新用

sun老师，首先表示深深的谢意，占用你的时间，实在是……
  你对我问题的回应，让我明白了很多，上面做的命令流，我从中也受益匪浅，再次感谢！
  
上面那条r,2,2e-5  是做什么用的？
还有wprot,45  
CSYS,4  
ARSYM,Y,all, , , ,0,1  
wprot,-45  
变换wp，是为什么？
此外，我下面的是一个文献上的成图，请教该图是不是没设耦合呀？
怎么外形这么齐整？后处理的命令我不通，是不是处理方法不同呀？

sunhaolan

古流新用 wrote:
上面那条r,2,2e-5  是做什么用的？

指定声流体单元实常数（空气中的基准声压），便于我们在后处理中看dB值的结果。

还有
wprot,45   !旋转工作平面45度，也就是y轴逆时针旋转45度。
CSYS,4      !指定当前Active坐标为工作平面坐标。
ARSYM,Y,all, , , ,0,1    !对y轴做映射，实现模型的逆时针旋转45度。
wprot,-45    !工作平面还原为初始状态。
变换wp，是为什么？   因为我们原来建的模型是对于x轴对称的，而ansys中的轴对称是对于y轴对称的，所以我们要把原来建的模型给他转过来，才能实现周对称模型的求解。
此外，我下面的是一个文献上的成图，请教该图是不是没设耦合呀？
怎么外形这么齐整？后处理的命令我不通，是不是处理方法不同呀？
这个图实在太模糊了，看不出是什么结果。

古流新用

sunhaolan wrote:
  
nsel,s,loc,y,-1.5,-.9   
nsel,r,loc,x,.15
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-1.5,-1   
nsel,r,loc,x,.2
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-.9,.9   
nsel,r,loc,x,.55
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-1,1   
nsel,r,loc,x,.6
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,.9,1.5   
nsel,r,loc,x,.15
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,1,1.5  
nsel,r,loc,x,.2
esel,s,type,,2   
sf,all,fsi,1
alls
!不好意思，红色的命令重复了，删除掉就行了。
nsel,s,loc,y,-1.5,-.1   
nsel,r,loc,x,.15
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-.9,.9   
nsel,r,loc,x,.55
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-1,1   
nsel,r,loc,x,.6
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,.9,1.5   
nsel,r,loc,x,.15
esel,s,type,,2   
sf,all,fsi,1
alls

nsel,s,loc,y,-1.5  
nsel,r,loc,x,.15,.2
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-1  
nsel,r,loc,x,.2,.6
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,-.9  
nsel,r,loc,x,.15,.55
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,.9  
nsel,r,loc,x,.15,.55
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,1  
nsel,r,loc,x,.2,.6
esel,s,type,,2   
sf,all,fsi,1
alls
  
nsel,s,loc,y,1.5  
nsel,r,loc,x,.15,.2
esel,s,type,,2   
sf,all,fsi,1
alls
  
  

sun老师，为什么是16个施加fsi的命令，那个图形只有12条边，上面是有重叠的吧？
没错，是12条。

d,all,pres,2
  

这句pre的方向是什么？
我照作了一次，为什么方向如下
方向没有问题，它的显示就是这样的。
另外你图里的模型好像不对阿。应该是这样的。

古流新用

sun老师，你说“对不起“，我可受不起，真的！你指导我这么多，我说谢谢，都是无法表达的。
  
sun老师，我想做成以x轴为对称轴显示云图，我怎么把显示好的图，改过来？
  
还有，我还是不放心压力，那个确实是指向腔内吗？我有点嘀咕，那个显示真的是向右指的，只是显示问题吗？最后云图的显示为什么给人感觉，耦合造成的变形这么大？

sunhaolan

古流新用 wrote:
sun老师，你说“对不起“，我可受不起，真的！你指导我这么多，我说谢谢，都是无法表达的。
  
sun老师，我想做成以x轴为对称轴显示云图，我怎么把显示好的图，改过来？
  
还有，我还是不放心压力，那个确实是指向腔内吗？我有点嘀咕，那个显示真的是向右指的，只是显示问题吗？最后云图的显示为什么给人感觉，耦合造成的变形这么大？

不用客气。
想把显示好的图改过来？恐怕不好改，你可以把坐标给隐藏掉，然后存成图片后，用acdsee啊什么的看图软件给他旋转90度，掉过来。
声压是标量，不是矢量，没有方向，无所谓指向哪里，所以箭头只是一个标志性的显示，不用管。
云图给的显示是为了让你更直观的看到，你可以看看结构上及耦合流体边界上的节点的变形，数值是很小的。

skyhawkwu

不对称结构的消声器可以用二维简化嘛？

mz6868

如果要是非对称结构的复杂空间，并伴有阻尼和非均匀介质，想求在一定条件下的声压级怎么办好，用sysnoice好还是ansys好？谢谢