ansys网格模型导入SYSNOISE

ANSYS Workbench 模型导入与网格划分(吴淑芳)

u v

a1 a7

a2 a8
x x

a3 a9
y y

a4 x2 a5 xy a6 y2 a10 x2 a11xy a12 y 2

网格划分基础
1.网格类型及单元阶次
三维网格
网格划分基础
1.网格类型及单元阶次
网格划分基础
2.网格划分原则
• （1）网格数量 • （2）网格疏密 • （3）单元形状及评价 • （4）单元阶次 • （5）网格质量 • （6）网格分界面和分界点 • （7）位移协调性 • 网格• 划（分8）没有细定节式处，理只能根据经验划分网格，宽广的有限元知识和丰富的经验是保证
选用高阶单元可进步计算精度，由于高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界，所以当结构外形不规则、应力分布或变形很复杂时可以选用高阶单元。
但高阶单元的节点数较多，在网格数目相同的情况下由高阶单元组成的模型规模要大得多，因此在使用时应权衡考虑计算精度和时间。
网格划分基础
2.网格划分原则
（2）网格疏密划分疏密不同的网格主要用于应力分析(包括静应力和动应力)，而
计算固有特性时则趋于采用较均匀的钢格形式。这是由于固有频率和振型主要取决于结构质量分布和刚度分布，不存在类似应力集中的现象。同样，在结构温度场计算中也趋于采用均匀网格。
网格划分基础
2.网格划分原则
机械与动力工程学院CAD/CAM工程技术研究中心
ANSYS Workbench
模型导入与网格划分
主讲：吴淑芳
2014年11月
主要内容
一、模型导入
网格划分基础
二、网格划分
全局网格控制局部网格控制

ansys计算的位移结果导入sysnoise详细步骤

Ansys中计算的响应位移结构导入sysnoise（例子中用到的ansys为9.0版；sysnoise为5.5版）1、网格模型的导出将ANSYS中划好的网格模型，用cdwrite命令（命令格式如下）进行导出。

cdwrite,all,file-name,cbd如图4-1-4所示：图4-1-1其中：file-name为输出文件的文件名；cbd为输出的文件格式，无须改动。

运行命令后，ANSYS会在C:\Documents and Settings\用户名下生成两个如上文件名的文件，分别为cbd格式和iges格式，在此要用到的是cbd格式的文件。

2、响应结果的导出在对板件进行完谐响应分析后，将网格模型每个节点的响应结果输出，同样是在ANSYS软件界面上方Ansys Commond Prompt命令栏中输入如下命令流：/post1nfreq=(freqe-freqb)/df+1*get,nnod,node,,count*dim,r,array,nnod,7,nfreq*do,i,1,nfreq,1set,,i*get,nd,node,,num,min*do,j,1,nnod,1r(j,1,i)=ndr(j,2,i)=ux(nd)r(j,4,i)=uy(nd)r(j,6,i)=uz(nd)nd=ndnext(nd)*enddoset,,i,,1*get,nd,node,,num,min*do,j,1,nnod,1r(j,3,i)=ux(nd)r(j,5,i)=uy(nd)r(j,7,i)=uz(nd)nd=ndnext(nd)*enddo*enddo*cfopen,result,fre*do,i,1,nfreq,1*vwrite,'SYSNOISE',' DISPLAC','EMENTS','FILE' (a,a,a7,a4)*vwrite,'Rev 5.5 ','IBM P2ES','SL 20-AP','R-04' (a,a,a,a4)*vwrite,'Displace','_Structu','re'(a,a,a2)*vwrite,' 20-APR-','2004 10:','07:13'(a,a,a5)*vwrite,'FREQUENC','Y'(a,a1)freq=freqb+(i-1)*df*vwrite,freq(e20.8)*vwrite,'NODAL DI','SPLACEME','NT V ALUE','S' (a,a,a,a1)*do,j,1,nnod,1nseq=chrval(j)a=r(j,1,i)nd=chrval(a)uxr=r(j,2,i)uxi=r(j,3,i)uyr=r(j,4,i)uyi=r(j,5,i)uzr=r(j,6,i)uzi=r(j,7,i)*vwrite,' ',nseq,' ',nd,uxr,uxi,uyr(a2,a8,a2,a8,e20.8,e20.8,e20.8)*vwrite,'','','',uyi,uzr,uzi(a8,a8,a4,e20.8,e20.8,e20.8)*vwrite,'','','',0,0,0(a8,a8,a4,e20.8,e20.8,e20.8)*vwrite,'','','',0,0,0(a8,a8,a4,e20.8,e20.8,e20.8)*enddo*enddo*cfclos需要注意的是：nfreq=(freqe-freqb)/df+1一句中，freqe为要提取的截止频率，freqb为起始频率，两频率都不能超出之前在ansys中分析的频率。

基于ANSYS和SYSNOISE的弹性结构噪声分析

株洲４２０１０８）
ｔ
（湖南工业大学机械学院，湖南
要：对弹性结构的噪声分析过程中，结构表面的声学量很难精确求解，结构的形貌往往需要优化．利
用声学分析软件ｓＮｏｓ可以计算结构表面及任意点的辐射声场，提出了将有限元分析软件ＡＹＳＹｓＩＥＮＳ得到的
学一ｙ
院
璐
ＡｂｔａｔＩｅｃｕｓｆｈｏｓｎｌｓｓｎｔｅｅａｔｃｓｒｃｕｅｔｅａｏｓｉｏｕｆｈｔｕｔｒｌｓｒｃ：ｎｔｏｒｅｏｅｎｉｅａａｙｉｌｓｉｔｕｔｒ．ｈｃｕｔｖｌｍｅｏｅｓｃｕａｈｔｏｈｃｔｒ
中图分类号：Ｔ５２Ｈ１２Ｂ３，Ｔ２文献标识码：Ａ文章编号：１７— ３４２１）ｌ０４一４６２７０（０１Ｏ－０ｌ０
城锄
市Ｃ
ＮｏｓｉｅＡｎａｙｉｆＥｌｓｉｔｕｔｅｂａｅｎＡＮＳＹＳａｄＹＳＯＩｌｓｓｏａｔｃＳｒｃｕｒｓｄＯｎＳＮＳＥ
ＬＩＳｕＸＩｈｎ，ＯＮＧｏｇ－ｎＹｎｇａｇ
（ＣｏｌｇｆＭｅｈｎｉｓＨｕｎＵｎｖｒｉｆｅｈｏｏｙＺｚｏ，ｎｎ４１０８Ｃｈｉａｌｅｅｏｃａｃ，ｎａｉｅｓｔｏＴｃｎｌｇ，ｈｕｈｕＨｕａ２０，ｙｎ）
随着人们生活水平的提高，噪声污染越来越严重，特别是中、低频噪声倍受关注．弹性结构

扩张式消声器消声特性理论研究和实验分析

ｉｓｆｏｕｎｄｔｈａｔｔｈｅｓｏｕｎｄｗａｖｅｏｆｔｈｅｍｕｆｆｌｅｒｉｓｎｏｔｓｔｒｉｃｔｌｙｓｐｒｅａｄａｃｃｏｒｄｉｎｇｔｏｔｈｅｕｎｉｄｉｍｅｎｓｉｏｎａｌｐｌａｎｅ
ｗａｖｅｆｏｒｍ．Ｔｈｅｎｏｉｓｅａｔｔｅｎｕａｔｉｏｎｐｅｒｆｏｒｍａｎｃｅｏｆｔｈｅｍｕｆｆｌｅｒｉｓｉｍｐｒｏｖｅｄｂｙｍｅｎｄｉｎｇｔｈｅｍｕｆｆｌｅｒ．Ｋｅｙｗｏｒｄｓ：ａｃｏｕｓｔｉｃｓ；Ｓｙｓｎｏｉｓｅ；ｆｉｎｉｔｅｅｌｅｍｅｎｔｓｉｍｕｌａｔｉｏｎ；ｅｘｐａｎｓｉｏｎｍｕｆｆｌｅｒ；ｉｎｓｅｒｔｉｏｎｌｏｓｓ
３）对原有结构改进的有限元分析中可知，改进后消声器的消声特性明显优于原有结构。从而表明双扩张式消声器消声特性优于单扩张式的；并得出消声器的扩张比是影响消声器消声特性的主要参数，而消声器扩张腔长度、直径大小、入口腔和扩张腔间圆角大小、人口和出口轴线的错开多少都不同程度地影响消声效果，为进一部改进设计奠定了基础。
第１期
蹒１单扩张式消声器漓声特性的瑷论分析曲线豳
出现后，消声腔消声作用已大大降低，易形成径向平面波，一维平西波翡分拼结果就不准确，使消声器消声分析在高频失效。
２基于有限元的仿真分析
由于以上理论分析方面存在一定缺陷，而有限元三维数值仿真解法只要处理好边值问题，在求解域蠹可以较耩灌地算出复杂清声器声学模型在任俺
扩张式消声器消声特性理论研究和实验分析
１０５
文章编号：１００６—１３５５（２００８）０１－０１０５－０３
扩张式消声器消声特性理论研究和实验分析
张晓龙，李功宇
（昆明理工大学机电工程学院，昆明６５００９３）
摘要：对某工业消声器进行了理论研究和实验分析，通过分析试验得到这类消声器的消声特性，指出在实际

ANSYSIMPORT

ANSYSIMPORT
通过UG可以用多种文件格式导入ANSYS，如：prt、model、iges、x_t等，不同格式导入的模型时，可能在ANSYS系统内修复的工作不尽相同，不能说那种格式一定就好，所以多试试不同的接口，但IGES通常问题较多，推荐不采用。

复杂模型导入时，如果一次性导入全部实体，困难较大，而且划分网格也很困难，所以推荐拆分后导入，如果导入实体有困难可以考虑先导入面模型，然后再在ANSYS里面想办法生成体。

但这里面也有问题，
分几次导入的实体，要想在ANSYS里面进行布尔操作，经常会遇到问题，不如CAD系统里面那么容易，所以模型怎么拆分导入，导面还是导体，这个就需要自己去权衡。

总之，复杂模型导入ANSYS时，花在前处理的时间较多，可能需对模型进行修复和重构，不过复杂模型多做几个的话，你的前处理能力将肯定会有大大的提高，祝你顺利！
直接将UG倒入也不很不错。

不过导入后，记住存一下盘，这样会加快速度。

另外导入后，采用/facet,normal，可以将面和体正常显示出来
我做复杂模型时多用solidworks做，然后用sat格式导入ansys，目前还未发现有数据丢失的问题。

利用ANSYS谐响应分析结果导入LMS Virtual lab中进行声学分析步骤

1.前期用ANSYS对模型进行动力学分析，然后保存结果文件.rst格式的，然后导入到Vrituallab12中进行声学分析，可能步骤有些长，大家尽量慢慢看，如果有不明白的，或者我的步骤有错误的，大家可以指正，还有我的VL版本是12的，12的版本和以前的微有不同，在后边大家会发现的。

我的Q1728993717.2.进入声学模块：开始—Acoustics—Acoustics Harmonic BEM ；3.导入Ansys分析结果文件.rst格式：文件—Import—默认即可，看好单位，与模型统一；4.更改文件名称，便于后续操作：在特征树中点开Nodes and Elements—右键点其子选项（就是带有齿轮标志那个）—属性—特征属性—更改名称—StructuresMesh.5.提取声学面网格：开始—Structures—Cavity Meshing—插入—Pre/Acoustics Meshers—Pre/Acoustics Meshers—Skin Meshers，出现一下图框，在Grid to Skin 区域选择结构网格即：StructuresMesh，其余都默认不用改，之后点击应用，Close。

6.在次回到声学模块：开始—Acoustics—Acoustics Harmonic BEM ；7.命名声学网格：点开特征树中的Nodes and Elements—右键Skin Meshpar1.—属性—特征属性—改名称—AcousticsMesh；到这步之后为了方便起见，可以将结构网格StructuresMesh隐藏：右键StructuresMesh—Hide/Show；8.设定分析类型：工具—Edit the Model Type Definitions—点击“是”出现对话框如下：按照图所示设置即可；9.设置网格类型：工具—Set Mesh parts Type：之后，在左边选中StructuresMesh，然后点右边的Set as Structures；同理，选中AcousticsMesh点击右边Set as Acoustics；然后确定即可；10.声学网格前处理：插入—Acoustic Mesh Prepocessing set 出现如下：在Mesh Parts 中选声学网格AcousticsMesh—确定即可；11.定义材料：插入—Materials—New Materials—New Fluid Materials按下图选着填写即可：其实就更改个Materials ID为Air 其余就都是默认即可，不用更改什么，然后点击确定。

SolidWorks与ANSYS之间的数据交换方法研究

4结语利用宏程序编写非圆曲线程序，不但节省了使用软件造型、自动编程、程序传输等方面的时间，而且使程序的可读性、简洁性和合理性大大增加，减少了加工的准备工作时间，从而提高了生产效率。

宏程序在一些三维曲面的编程中也得到了广泛的应用，对于不同的数控系统，宏程序的编写指令和格式有所差异，但编写的总体方法和思路基本相同。

随着数控系统的发展，二次开发功能将不断得到加强和完善，宏程序功能也会得到进一步的加强。

参考文献：［1］陈银清．宏程序编程在数控加工中的应用研究［J ］.机床与液压，2009（5）：42-45.［2］王荣兴．加工中心培训教程［K ］.北京：机械工业出版社，2006.作者简介：倪贵华（1980-），江苏南通人，南京理工大学，机械电子工程专业，实验师，数控铣床高级技师，加工中心高级技师，从事数控技术领域编程加工与数控技术维修，电子信箱：04sk323@.责任编辑：武伟民收稿日期：2011－02－01煤矿机械Coal Mine MachineryVol.32No.09Sep.2011第32卷第09期2011年09月!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!引言SolidWorks 软件是世界上第1个基于Windows 开发的三维CAD 系统，其功能强大，用途广泛，是一款智能型高级CAD/CAE/CAM 组合软件。

ANSYS 软件是集结构、流体、电场、磁场、声场、热分析于一体的大型通用有限元分析软件，由世界上最大的有限元分析软件公司之一的美国ANSYS 公司开发，其完善的有限元分析功能已为全球业界所广泛接受。

因此，如何把2个软件结合起来，将在SolidWorks 中建立的复杂模型完整地导入到ANSYS 中去，进而完成所需的有限元分析则成为问题的关键所在。

1SolidWorks 与ANSYS 之间的数据交换ANSYS 提供的与多数CAD 软件的专用接口，实现数据的共享和交换，如Pro/E,NASTRAN,Alo -gor,I-DEAS,AutoCAD 等，而SolidWorks 又能输出IGES 、Parasolid 等格式的文件，这使两软件之间的数据交换成为可能。

ANSYS的基本步骤讲解

ANSYS的基本步骤讲解1.问题定义：在使用ANSYS之前，您需要明确要解决的问题。

定义问题包括确定您要分析的物理现象，所使用的材料，边界条件和所需的结果等。

2.创建几何模型：根据问题定义，您需要创建一个几何模型来表示分析的对象。

ANSYS 提供了各种建模工具，可以用于创建二维和三维的几何形状。

您可以使用ANSYS自带的CAD工具或导入其他CAD软件创建的模型。

3.划分网格：对几何模型进行网格划分是进行仿真分析的关键步骤。

ANSYS提供了各种网格划分工具，可以根据需要选择不同的划分技术和网格密度。

良好的网格划分可以提高仿真的准确性和效率。

4.材料属性定义：根据问题定义，您需要为模型中的不同部分定义适当的材料属性。

ANSYS提供了一个材料库，可以选择多种不同的材料，并根据需要定义其属性参数，如弹性模量，热导率等。

5.添加约束条件：在仿真中，机械结构通常受到约束条件的限制。

您需要添加适当的约束条件来代表物理世界中的限制。

ANSYS提供了各种约束选项，包括固定支撑、自由支撑、弹簧等。

通过添加这些约束条件，您可以更准确地模拟实际场景。

6.应用载荷：在仿真中，您需要明确地定义作用在模型上的载荷。

这可能是一个力，一个压力，一个温度或者其他物理效应。

您需要在模型的相关位置上添加适当的载荷。

ANSYS提供了各种载荷选项，可以精确描述应用的载荷。

7.设置仿真参数：在进行仿真之前，您需要设置一些模拟参数。

这些参数可以控制仿真过程和结果的精度。

例如，您可以设置时间步长、迭代次数、收敛标准等。

8.运行仿真：配置完所有参数后，可以开始运行仿真。

ANSYS将根据所设置的仿真参数对模型进行计算。

这可能需要一段时间，具体取决于模型的大小和复杂程度。

9.结果分析：仿真运行结束后，可以进行结果分析。

ANSYS提供了丰富的结果可视化和后处理工具，可以帮助您更好地理解和解释结果。

您可以查看模型的位移、应力、应变、温度分布等结果。

10.结果验证：最后，您需要对仿真结果进行验证。

[原创]图文并茂详解如何在ANSYS 12 Workbench中导入ICEM网格进行分析--ripperjack

ansys11 workbench如何使用icem生成网格大家也许比较熟悉了，但是12里面如何导入icem 的网格一直没有一个详细的教程，本人折腾了很久终于成功搞定，特此分享一下。

首先新建一个static strutural 分析，在geometry上依次右键transfer data from new-finite element modeler(FEM)。

如下图：双击新生成的FEM，选择需要导入的.uns 网格文件后进入FEM，如下：在FEM中我们可以在component里面看到ICEM中定义好的各个面，如下图中的IN、OUT、WALL三个面接着，右键geomery synthesis-insert-initial geometry ，如下图但是结果却令人失望，似乎wall这个面有些问题，如下图所示继续，在intial geometry 右键选covert to parasolid我们可以看到，新生成的parasolid geometry只有两个面，周围的面丢失了。

为什么会这样本人也不清楚，具体原因有待高手解答，既然导入ICEM的component不行那么就试试看skin detection tools。

首先，先把刚才生成的parasolid 和initial geometry删除，并且同时删除wall、in、out这三个ICEM中定义好的component，接着对skin detection tools 右键－create skin components，如下图所示：FEM自动检测到的面如下：接下来，重复上面的步骤，insert initial geometry，并且生成parasolid geometry，这时候的结果如下图，我们发现，parasolid里面全是面，没有实体。

接着，对parasolid geometry右键-add a sew tool选中这个立方体的6个面apply这时候我们已经把这六个面缝合成一个实体了保存FEM返回Workbench，点updata project，workbench会自动把FEM中的几何与网格导入static strutural导入后如下图：static strutural中的setup查看导入的几何与网格，OK！。

ansys与sysnoise转换问题

ansys与sysnoise转换问题ansys与sysnoise转换问题看到大家都比较关心这个问题，为此提供一种采用ansys命令流直接生成sysnoise的fre文件（振动模态表面位移），以供参考。

!从ANSYS中导出模态频率及振型数据/Post1*get,nodenum,node,,num,max !读取节点数*set,modenum,20 !指定模态的阶数*Cfopen,modefile_ansys,fre !打开文件*Vwrite('SYSNOISE MODES STRUCTURE FILE') !头文件*Vwrite('Rev 5.5 Windows NT 30NOV2000')*Vwrite('SYSNOISE Default Model')*Vwrite('23-JUN-2008 10:22:59')*Do,i,1,modenum,1*Vwrite('REAL MODE FREQUENCY EIGENVALUE')*Get,modefq,mode,i,freqbb='%i%'modeanglefq = modefq*modefq*2*2*3.141593*3.141593 *Vwrite, bb, modefq, modeanglefq(A10,2e20.8)*Vwrite('DISPLACEMENT VALUES')set,1,i*Do,ii,1,nodenum,1*Get,modeshpx,node,ii,u,x*Get,modeshpy,node,ii,u,y*Get,modeshpz,node,ii,u,z*Get,modeshprx,node,ii,rot,x*Get,modeshpry,node,ii,rot,y*Get,modeshprz,node,ii,rot,zbb1 ='%ii%'bb2 =' '*vwrite,bb1,bb1,modeshpx,modeshpy,modeshpz(2A10,3e20.8)*vwrite,bb2,modeshprx,modeshpry,modeshprz(A20,3e20.8)*Enddo*Enddo*cfclosFINISH或者参考haohaosir的原创帖/doc/976957018.html,/forum/thread-10954-1-1.html/POST26NSOL,2,10,U,Y,UY_2 ! 定义第二个变量为UY_2,值为10号节点Y 方向的位移,节点可任选,但要保证其值非0 XVAR,1 ! 定义时间变量为坐标横轴PLVAR,2*GET,num_var,VARI,0, NSETS ! 将变量长度值赋给变量num_var k=num_var*DIM,SYSNOISE_TITLE,CHAR,5,4SYSNOISE_TITLE(1,1)='SYSNOISE 'SYSNOISE_TITLE(1,2)=' ACCELER 'SYSNOISE_TITLE(1,3)='ATIONS 'SYSNOISE_TITLE(1,4)=' FILE'SYSNOISE_TITLE(2,1)='Rev 5.5 'SYSNOISE_TITLE(2,2)=' IBM P2E 'SYSNOISE_TITLE(2,3)='SSL 11'SYSNOISE_TITLE(2,4)='-AUG-02 'SYSNOISE_TITLE(3,1)='ACCELERA'SYSNOISE_TITLE(3,2)='TION_St'SYSNOISE_TITLE(3,3)='ructure 'SYSNOISE_TITLE(4,1)='11-AUG-2002 'SYSNOISE_TITLE(4,2)=' 10:07 'SYSNOISE_TITLE(4,3)=':13 'SYSNOISE_TITLE(5,1)='TIME'*CFOPEN,ACCE,fre ! 定义数据结果文件，ACCE.fre*VWRITE,SYSNOISE_TITLE(1,1) ,SYSNOISE_TITLE(1,2) ,SYSNO ISE_TITLE(1,3) ,SYSNOISE_TITLE(1,4)%C%C%C%C*get,nmax,node,,num,max, !得到节点编号的最大数*DIM,T_U,ARRAY,k,4! 定义保存时间历程结果的矩阵T_U，有num_var行，4列VGET,T_U(1,1,1),1 ! 矩阵的第一列保存第一个时间历程变量（默认为TIME）*dim,num_t,,1,5 !定义一个数组*DIM,SYSNOISE_SUBT1,CHAR,1,4SYSNOISE_SUBT1(1,1)='NODAL AC'SYSNOISE_SUBT1(1,2)='CELERATI 'SYSNOISE_SUBT1(1,3)='ON VALU'SYSNOISE_SUBT1(1,4)='ES'*do,i,1,k,1num_t(1,1)=T_U(i,1,1) !按次序取TIME值*VWRITE,num_t(1,1)%20.8E*VWRITE,SYSNOISE_SUBT1(1,1) ,SYSNOISE_SUBT1(1,2) ,SYS NOISE_SUBT1(1,3) ,SYSNOISE_SUBT1(1,4) %C%C%C%C *do,j,1,nmax,1NSOL,2,j,U,X,UX ! 定义第二个变量为UX,值为j号节点X方向的位移DERIV,3,2,1,,,,,1 ! 将位移对时间求导，得到速度，存为变量3DERIV,4,3,1,,,,,1 ! 将速度对时间求导，得到加速度，存为变量4 VGET,T_U(1,2,1),4 ! 矩阵的第2列保存第四个时间历程变量，j号节点X方向的加速度时间历程结果num_t(1,2)=T_U(i,2,1) !将与i对应的j节点的X加速度值赋给num_t的第2列NSOL,5,j,U,Y,UY ! 定义第5个变量为UY,值为j号节点Y方向的位移DERIV,6,5,1,,,,,1 ! 将位移对时间求导，得到速度，存为变量6DERIV,7,6,1,,,,,1 ! 将速度对时间求导，得到加速度，存为变量7 VGET,T_U(1,3,1),7 ! 矩阵的第3列保存第7个时间历程变量，j号节点Y方向的加速度时间历程结果num_t(1,3)=T_U(i,3,1) !将与i对应的j节点的Y加速度值赋给num_t的第3列NSOL,8,j,U,Z,UZ ! 定义第8个变量为UZ,值为j号节点Z方向的位移DERIV,9,8,1,,,,,1 ! 将位移对时间求导，得到速度，存为变量9DERIV,10,9,1,,,,,1 ! 将速度对时间求导，得到加速度，存为变量10 VGET,T_U(1,4,1),10 ! 矩阵的第4列保存第10个时间历程变量，j 号节点Z方向的加速度时间历程结果num_t(1,4)=T_U(i,4,1) !将与i对应的j节点的Z加速度值赋给num_t的第4列num_t(1,5)=j*VWRITE,num_t(1,5),num_t(1,5),num_t(1,2),0,num_t(1,3) ! 将时间历程结果数据写到数据结果文件中,0为虚部及旋转加速度值%10I%10I%20.8E%20.8E%20.8E*VWRITE,0,num_t(1,4),0(20X,E20.8,E20.8,E20.8)*VWRITE,0,0,0(20X,E20.8,E20.8,E20.8)*VWRITE,0,0,0(20X,E20.8,E20.8,E20.8)*enddo*IF,i,EQ,k,THEN*EXIT*ENDIF*VWRITE,SYSNOISE_TITLE(1,1) ,SYSNOISE_TITLE(1,2) ,SYSNO ISE_TITLE(1,3) ,SYSNOISE_TITLE(1,4)%C%C%C%C*enddo*CFCLOSFINISH原创]第2种方法将Ansys计算结果(比如加速度值)写为fre文件的命令流/POST26NSOL,2,231,U,Z,UY_2 ! 定义第二个变量为UY_2,值为10号节点Y方向的位移,节点可任选,但要保证其值非0XVAR,1 ! 定义时间变量为坐标横轴PLVAR,2*GET,num_var,VARI,0, NSETS ! 将变量长度值赋给变量num_var k=num_var*DIM,SYSNOISE_TITLE,CHAR,5,4SYSNOISE_TITLE(1,1)='SYSNOISE 'SYSNOISE_TITLE(1,2)=' ACCELER 'SYSNOISE_TITLE(1,3)='ATIONS 'SYSNOISE_TITLE(1,4)=' FILE'SYSNOISE_TITLE(2,1)='Rev 5.5 'SYSNOISE_TITLE(2,2)=' IBM P2E 'SYSNOISE_TITLE(2,3)='SSL 11'SYSNOISE_TITLE(2,4)='-AUG-02 'SYSNOISE_TITLE(3,1)='ACCELERA'SYSNOISE_TITLE(3,2)='TION_St'SYSNOISE_TITLE(3,3)='ructure 'SYSNOISE_TITLE(4,1)='11-AUG-2002 'SYSNOISE_TITLE(4,2)=' 10:07 'SYSNOISE_TITLE(4,3)=':13 'SYSNOISE_TITLE(5,1)='TIME'*CFOPEN,ACCE1,fre ! 定义数据结果文件，ACCE.fre*VWRITE,SYSNOISE_TITLE(1,1) ,SYSNOISE_TITLE(1,2) ,SYSNO ISE_TITLE(1,3) ,SYSNOISE_TITLE(1,4) %C%C%C%C*get,nmax,node,,num,max, !得到节点编号的最大数kk=nmax*3+1*DIM,T_U,ARRAY,k,kk! 定义保存时间历程结果的矩阵T_U，有num_var行，kk列VGET,T_U(1,1,1),1 ! 矩阵的第一列保存第一个时间历程变量（默认为TIME）*dim,num_t,,1,5 !定义一个数组*DIM,SYSNOISE_SUBT1,CHAR,1,4SYSNOISE_SUBT1(1,1)='NODAL AC'SYSNOISE_SUBT1(1,2)='CELERATI 'SYSNOISE_SUBT1(1,3)='ON VALU'SYSNOISE_SUBT1(1,4)='ES's=1*do,j,1,nmax,1NSOL,2,j,U,X,UX ! 定义第二个变量为UX,值为j号节点X方向的位移DERIV,3,2,1,,,,,1 ! 将位移对时间求导，得到速度，存为变量3DERIV,4,3,1,,,,,1 ! 将速度对时间求导，得到加速度，存为变量4 s=s+1VGET,T_U(1,s,1),4 ! 矩阵的第s列保存第四个时间历程变量，j号节点X方向的加速度时间历程结果NSOL,5,j,U,Y,UY ! 定义第5个变量为UY,值为j号节点Y方向的位移DERIV,6,5,1,,,,,1 ! 将位移对时间求导，得到速度，存为变量6DERIV,7,6,1,,,,,1 ! 将速度对时间求导，得到加速度，存为变量7 s=s+1VGET,T_U(1,s,1),7 ! 矩阵的第s列保存第7个时间历程变量，j号节点Y方向的加速度时间历程结果NSOL,8,j,U,Z,UZ ! 定义第8个变量为UZ,值为j号节点Z方向的位移DERIV,9,8,1,,,,,1 ! 将位移对时间求导，得到速度，存为变量9DERIV,10,9,1,,,,,1 ! 将速度对时间求导，得到加速度，存为变量10 s=s+1VGET,T_U(1,s,1),10 ! 矩阵的第s列保存第10个时间历程变量，j 号节点Z方向的加速度时间历程结果*enddo*do,i,1,k,1num_t(1,1)=T_U(i,1,1) !按次序取TIME值*VWRITE,num_t(1,1)%20.8E*VWRITE,SYSNOISE_SUBT1(1,1) ,SYSNOISE_SUBT1(1,2) ,SYS NOISE_SUBT1(1,3) ,SYSNOISE_SUBT1(1,4) %C%C%C%C s=1*do,j,1,nmaxs=s+1num_t(1,2)=T_U(i,s,1) !将与i对应的j节点的X加速度值赋给num_t的第2列s=s+1num_t(1,3)=T_U(i,s,1) !将与i对应的j节点的Y加速度值赋给num_t的第3列s=s+1num_t(1,4)=T_U(i,s,1) !将与i对应的j节点的Z加速度值赋给num_t的第4列num_t(1,5)=j*VWRITE,num_t(1,5),num_t(1,5),num_t(1,2),0,num_t(1,3) ! 将时间历程结果数据写到数据结果文件中,0为虚部及旋转加速度值%10I%10I%20.8E%20.8E%20.8E*VWRITE,0,num_t(1,4),0(20X,E20.8,E20.8,E20.8)*VWRITE,0,0,0(20X,E20.8,E20.8,E20.8)*VWRITE,0,0,0(20X,E20.8,E20.8,E20.8)*enddo*IF,i,EQ,k,THEN*EXIT*ENDIF*VWRITE,SYSNOISE_TITLE(1,1) ,SYSNOISE_TITLE(1,2) ,SYSNO ISE_TITLE(1,3) ,SYSNOISE_TITLE(1,4)%C%C%C%C*enddo*CFCLOSFINISH该方法与上次传的相比：计算量大大减少、因此速度要快！因为第一种方法每写一步数据，变量需要计算一次。

在ansys中导入自定义本构模型的算例

在ansys中导入自定义本构模型的算例邓肯-张模型的关键点是材料的弹性模量随大小主应力差及小主应力(围压)的变化而变化，用APDL实现之的基本思路是：给每个单元定义一个材料号，分级施加荷载，在每个荷载步结束时提取出各单元的大小主应力，据此计算出下个荷载步的弹性模量Et,修改各单元之MP，用于下一步计算。

以下是一个简单算例，copy出去可直接运行。

!!!常规三轴试验模拟!!!by taomingxing,NWPU!!!2003.7.16FINISH/CLEAR/TITLE,Numerical Simulation of three axes testing of soils/PREP7*dim,SUy,array,50 !Settlement records*dim,MaxPs,array,120 !Max history p1-p3*dim,MaxDs,array,120 !Max history Ds!*dim,EEt,array,50 !Et of elememt!!!Duncan-Chang Model!!!Symbols：c-粘滞力，Fai-内摩擦角，Sf-破坏强度(p1-p3)f，Ds-应力水平，Pa-大气压,P3-围压*CREATE,Duncan-Chang !Creat Macro file*afun,deg !Unit of angle*set,Pa,1e5*set,P1,-ArrS3(i) !注意：岩土工程中应力为拉负压正*set,P3,-ArrS1(i)*if,P3,LT,0.1*Pa,thenP3=0.1*Pa!围压最小取值*endifSf=2*(c*cos(Fai)+P3*sin(Fai))/(1-sin(Fai)) !Mohr-Coulomb破坏强度(p1-p3)f Ds=(P1-P3)/Sf !应力水平，*if,Ds,GT,0.95,thenDs=0.95!应力水平最大取值*endif!判断加卸荷，如果(P1-P3)小于历史最大值视为卸荷－再加荷过程*if,MaxPs(i),LT,P1-P3,thenEi=k*Pa*(P3/Pa)**nEt=Ei*(1-Rf*Ds)**2!加荷情况的切线模量MaxPs(i)=P1-P3 !保存历史最大应力*elseif,MaxPs(i),GE,P1-P3Et=Kur*Pa*(P3/Pa)**n!卸荷模量*endifmp,ex,i,Et !修改单元i的Et mp,nuxy,i,Mu*END!!!单元类型et,1,42!平面四节点单元KEYOPT,1,3,2!平面应变!!!以下定义材料初始模量mp,ex,1,3.728e7 !砂土的弹性模量mp,nuxy,1,0.33mp,dens,1,1800!!!建立几何模型blc4,0,0,0.08,0.15!8cm X 15cm /PNUM,AREA,1/REPLOT!!!网格划分aesize,all,0.01mat,1amesh,all!!!边界条件nsel,s,loc,y,0d,all,Uy !底边界竖向约束nsel,s,loc,x,0d,all,Ux !左侧边界水平向约束nsel,all/replotfini/SOLUtime,0.01 !施加围压sfl,all,pres,2e5 !200kPasolve!!!分级施加荷载，实现非线性计算!!!荷载增量10kPa,共50级*DO,ti,1,50!!!取出计算结果，修改弹性模量/POST1*get,SUy(ti),node,29,u,y !Settlement record of time ti ETABLE,EtabS1,S,1 !取各单元第一主应力ETABLE,EtabS3,S,3!取各单元第三主应力*dim,ArrS1,array,120*dim,ArrS3,array,120*do,Num,1,120 !Num为单元编号*get,ArrS1(Num),elem,Num,etab,EtabS1!将单元结果存入数组*get,ArrS3(Num),elem,Num,etab,EtabS3*enddo/PREP7!^^修改砂土单元的Et,单元号1-120c=0 $Fai=35 $Rf=0.7 $k=400 $n=0.6 $Mu=0.33 $Kur=326.7*do,i,1,120 !各单元循环计算*use,Duncan-Chang,c,Fai,Rf,k,n,Mu,Kur !调用Duncan-Chang宏文件*enddo!EEt(ti)=ET !保存第120单元之ET/SOLUtime,tisfl,3,pres,2e5+1e4*ti !施加荷载，增量1e4solve !对ti级荷载情况求解*ENDDO。

sysnoise介绍

LMS公司的sysnoise软件是首家将边界元（BEM）应用到声学领域的商业软件。

一般来说用sysnoise基本上都是：通过结构的FEM得到结构表面的边界条件（振速、模态向量等），然后作为声场边界元的载荷输入条件，再用BEM 计算辐射声场。

也就是所谓的间接边界元（Indirect BEM）方法。

边界元的方法实际上也存在对高频信号处理能力的问题，基本上也只能做到500-1kHz左右。

sysnoise目前已经整合到LMS公司的整体CAE仿真解决方案b产品中了,如今其最新版本已出到了b Revision10。

对于中高频的声学问题，目前来说比较好的解决方法只有一种，也就是统计能量法（SEA, Statistical Energy Analysis）。

原理是从能量守恒的角度去进行分析，将整个模型分为若干个子系统，考虑到各个子系统内部的能量损耗因素，以及子系统之间耦合的能量损耗因素，形成一组能量守恒方程，来求解各个子系统的能量分布情况，而声能量的表征也就是要求解的声学响应。

但是统计能量法也有其缺点：无法对结构或声场的细节进行分析，因为它是基于能量及统计学的方法进行的，如果子空间太小，模态密度不够高，就符合不了统计方法的假设，结果也就不可信了。

LMS Virtual. Lab 软件的声学模块包括边界元法、有限元法无限元法技术，能够有效地模拟结构声学耦合现象，处理结构内部外部的辐射噪声问题，并无缝链接多种CAE软件，如ANSYS，MSC，NASTRAN与ABQAQUS等。

对结构声学的数值分析方法可以分为两大类：离散方法和能量方法。

能量方法主要是指统计能量分析（Statistical Energy Analysis(SEA)）和能量有限元法（Energy Finite Element Analysis(EFEA)）。

能量方法使用于中高频激励作用下模态密集结构振动与声的计算分析。

对中低频区的振动声辐射问题，通常采用单元离散方法。

基于SYSNOISE软件的薄板振动声辐射数值仿真研究

基于SYSNOISE软件的薄板振动声辐射数值仿真研究作者：刘先锋薛伟飞陈进来源：LMS1.引言车身壁板结构厚度小、质量轻，特别容易产生振动并且辐射噪声，因此计算汽车车身薄板受到外部激励时的振动和声辐射特性是十分必要的，然而在大多数工程实际问题中，结构振动引起的声辐射常常是无法用解析解的形式予以解决，SYSNOISE是国际著名振动和声学测试分析软件公司LMS（Leuven Measurement System International）研发的大型声学计算分析软件，能快速地进行声学计算分析。

本文利用振动声学软件SYSNOISE，对不同约束和不同激励点位置条件下薄板声辐射进行研究，比较了加筋板与板的声辐射功率级，得到一些抑制薄板结构振动和声辐射的方法，从而为汽车车身设计及减振降噪提供了可靠的预报和措施。

2.薄板结构振动和声辐射的理论2.1 声辐射功率的计算机械噪声大部分是由结构振动而辐射的，结构声辐射功率表示了声辐射系统向外辐射噪声的能力，它不仅与振动弹性物体固有的物理特性有关，还与激励力大小、频率以及辐射声环境有关。

根据空气介质的连续条件，认为邻近振动表面一层的振动速度就是振动表面的速度。

振动表面任一点的振动速度为V(x,ω)，如果振动表面为平面，设振动表面任意一点的振动声压为p(y,ω)，则由瑞利积分得：式中r为表面上任意两点x和y之间的距离。

于是可以得到振动表面声强为：由表面声强就可以得到结构辐射功率为：3.仿真研究以受垂向激励的矩形金属薄板振动声辐射为例进行仿真计算。

板长2287mm，宽600mm，厚1.2mm，材料密度ρs=7850kg/m3，弹性模量E=2.1×1011N/m2，泊松比μ=0.3，空气密度ρ=1.21kg/m3，空气中声速c=343m/s。

加强筋与薄板材料相同。

考虑到在实际工程中，薄板受到正弦激励是不常见的，因此选择随机激励，激励的方向垂直于薄板表面，激励带宽为50～1000Hz。

ANSYS Workbench 模型导入与网格划分-文档资料

ANSYS Workbench
模型导入与网格划分
主讲：吴淑芳
2014年11月
1
主要内容
一、模型导入
网格划分基础
全局网格控制二、网格划分局部网格控制
虚拟拓扑
网格检查
网格划分实例
2020/7/2
22
主要内容
一、模型导入
网格划分基础
全局网格控制二、网格划分局部网格控制
虚拟拓扑
网格检查
网格划分实例
（8）细节处理
34
对称性模型
35
对称性模型
36
对称性模型
37
对称性模型
38
网格划分基础
3.ANSYS Workbench网格划分界面
• ANSYS 网格划分不能单独启动，只能在 Workbench 中调用分析系统或【Mesh】组件启动，如图所示。
39
网格划分基础
3.ANSYS Workbench网格划分界面
2.网格划分原则
（4）单元阶次
图中是一悬臂梁分别用线性和二次三角形单元离散时，其顶端位移随网格数目的收敛情况。可以看出，但网格数目较少时，两种单元的计算精度相差很大，这时采用低阶单元是不合适的。当网格数目较多时，两种单元的精度相差并不很大，这时采用高阶单元并不经济。例如在离散细节时，由于细节尺寸限制，要求细节四周的网格划分很密，这时采用线性单元更合适。
31
网格划分基础
2.网格划分原则
（6）网格分界面和分界点
• 结构中的一些特殊界面和特殊点应分为网格边界或节点以便定义材料特性、物理特性、载荷和位移约束条件。
• 即应使网格形式满足边界条件特点，而不应让边界条件来适应网格。
• 常见的特殊界面和特殊点有材料分界面、几何尺寸突变面、分布载荷分界线(点)、集中载荷作用点和位移约束作用点等。具有上述几种界面的结构及其网格划分形式。

ANSYS12 ICEM导入Workbench步骤

ANSYS12 ICEM导入Workbench步骤
ANSYS11，和ANSYS13都可以直接调用ICEM进行网格划分，单独12不可以，期间的导入只能靠手工操作。

1，在ICEM中生成网格后，通过转换成无结构或者结构网格，生成.uns文件并保存。

2，启动Workbench，在点选
右击model——import mesh，选择第一步中生成的.uns文件。

3，右击edit，进入编辑，删除part_1，右击
，insert---initial geometry
initial geometry 右键选covert to parasolid后，生成
4，此时生成的几何中只是面体，没有体积。

需要对生成的Parasolid进行缝合，右键parasolid geometry -选择add a sew tool
单击黄色区域出现
，
选择
所有的body，然后Apply。

出现
，在Sew1右击，选择。

(注：若在这一部中出现警告，或者不成功，可能是容差的问题，删除刚生成的Sew1,重新来过第4步，在
之前，将容差变大，
，例如将.01变成.1,然后再
）
5，保存，关闭
6，拖动到workbench，但不要与
进行连接。

8、拖动A2到B4
（一定是中的与
中的相连接，
9，点击，就完成了输入
注意问题
1，第3步中，一定要删除，不然在initial geometry时出现错误
此时生成的模型与原来的模型有很大区别
2，一定要在前5步完成之后再与其他模块连接，否则影响最后生成网格的质量
3，一定是model与model连接，不然不会生成网格。

!从ANSYS中导出模态频率及振型数据

!从ANSYS中导出模态频率及振型数据!从ANSYS中导出模态频率及振型数据/Post1*get,nodenum,node,,num,max !读取节点数*set,modenum,20 !指定模态的阶数*Cfopen,modefile_ansys,fre !打开文件*Vwrite('SYSNOISE MODES STRUCTURE FILE') !头文件*Vwrite('Rev 5.5 Windows NT 30NOV2000')*Vwrite('SYSNOISE Default Model')*Vwrite('23-JUN-2008 10:22:59')*Do,i,1,modenum,1*Vwrite('REAL MODE FREQUENCY EIGENVALUE')*Get,modefq,mode,i,freqbb='%i%'modeanglefq = modefq*modefq*2*2*3.141593*3.141593 *Vwrite, bb, modefq, modeanglefq(A10,2e20.8)*Vwrite('DISPLACEMENT VALUES')set,1,i*Do,ii,1,nodenum,1*Get,modeshpx,node,ii,u,x*Get,modeshpy,node,ii,u,y*Get,modeshpz,node,ii,u,z*Get,modeshprx,node,ii,rot,x*Get,modeshpry,node,ii,rot,y*Get,modeshprz,node,ii,rot,zbb1 ='%ii%'bb2 =' '*vwrite,bb1,bb1,modeshpx,modeshpy,modeshpz(2A10,3e20.8)*vwrite,bb2,modeshprx,modeshpry,modeshprz(A20,3e20.8)*Enddo*Enddo*cfclosFINISH*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,3X,F10.3,5X,F8.3,5X,F8.3)其中1X表示一个空格数，A8是输出字符型LABEL，占8位（列），F10.3表示浮点型数值，占据10列，小数3位，F8.3表示浮点型数值，占据8列，小数3位。

利用ANSYS谐响应分析结果导入LMS Virtual lab中进行声学分析步骤

1.前期用ANSYS对模型进行动力学分析，然后保存结果文件.rst格式的，然后导入到Vrituallab12中进行声学分析，可能步骤有些长，大家尽量慢慢看，如果有不明白的，或者我的步骤有错误的，大家可以指正，还有我的VL版本是12的，12的版本和以前的微有不同，在后边大家会发现的。

我的Q1728993717.2.进入声学模块：开始—Acoustics—Acoustics Harmonic BEM ；3.导入Ansys分析结果文件.rst格式：文件—Import—默认即可，看好单位，与模型统一；4.更改文件名称，便于后续操作：在特征树中点开Nodes and Elements—右键点其子选项（就是带有齿轮标志那个）—属性—特征属性—更改名称—StructuresMesh.5.提取声学面网格：开始—Structures—Cavity Meshing—插入—Pre/Acoustics Meshers—Pre/Acoustics Meshers—Skin Meshers，出现一下图框，在Grid to Skin 区域选择结构网格即：StructuresMesh，其余都默认不用改，之后点击应用，Close。

6.在次回到声学模块：开始—Acoustics—Acoustics Harmonic BEM ；7.命名声学网格：点开特征树中的Nodes and Elements—右键Skin Meshpar1.—属性—特征属性—改名称—AcousticsMesh；到这步之后为了方便起见，可以将结构网格StructuresMesh隐藏：右键StructuresMesh—Hide/Show；8.设定分析类型：工具—Edit the Model Type Definitions—点击“是”出现对话框如下：按照图所示设置即可；9.设置网格类型：工具—Set Mesh parts Type：之后，在左边选中StructuresMesh，然后点右边的Set as Structures；同理，选中AcousticsMesh点击右边Set as Acoustics；然后确定即可；10.声学网格前处理：插入—Acoustic Mesh Prepocessing set 出现如下：在Mesh Parts 中选声学网格AcousticsMesh—确定即可；11.定义材料：插入—Materials—New Materials—New Fluid Materials按下图选着填写即可：其实就更改个Materials ID为Air 其余就都是默认即可，不用更改什么，然后点击确定。