（必看）ANSYS转子动力学计算讨论(2)

判断方法很简单，只需判断错误文件file.err是否可写就可以了。

因为当ANSYS在运行时，file.err是不可写的，只有当它运行完毕，此文件才可写。

原作者webycn

VC中调用Ansys的方法如下，具体过程就不写了。

: 1 利用Ansys建立一个通用有限元模型，并进行计算求解。注意：根据模型中所有可变 的参数定义相应的参数，利用参数化建模。将生成的log 文件拷贝成另一文件，即宏命 令文件。当然可以直接用APDL编写此文件。 ：2 VC编写界面：

（1） 要实现前处理参数输入功能； (2） 根据用户输入的参数修改刚才得到的ansys宏命令文件中，参数对于的参数值 ； (3） 将生成的宏命令文件文件提交给Ansys程序进行批处理操作；

在VC＋＋ 中可利用CreateProcess函数创建一个进程，去执行其他程序，并且可以设置改进程的优先级。

Ansys提供了一种批处理方式的格式：\e. 其中InputName和OutName分别为输入输出文件名。 （4）程序需要判断Ansys的批处理操作何时结束。

在Window NT 操作系统中当一个进程运行完毕后，窗口的标题会有 \已完成\这几个字。如当Ansys批处理完成后，窗口标题会显示 \已完成”。可以利用这 一个特点来判断进程是否运行完毕。

： 3 利用结果文件（out文件，或者由你在宏命令中指定输出的文件）进行后处理

自评:其实这种方法不能是真正意义上的开发,只不过用vc创建一个进程,vc与ansys并没有真正意义上的结合.如果自己开发程序,你永远都无法脱离ansys这个软件环境.二次开发应该是DLL或ActiveX,或者其它的COM对象形式,可以离开开发软件的运行环境.不过,实现一定的批处理功能也不错. 希望多讨论

搞ansys二次开发的可以看看 这个帖子

【原创】自己做的转子临界转速计算的小程序,请高手指点

请指教为什么刚度的变化(即使在轴承的刚度选取范围内)会对临界转速的影响这么大,用matlab计算影响就没这么大.ps:如有错误请高手指出! 谢谢!! fini /clear r1=.020 r2=.030 r3=.225 r4=.030 r5=.025

r6=.020 r7=.0225 r8=.018 pi=acos(-1) /prep7

mp,ex,,2.1e11 mp,dens,,7850 mp,prxy,,0.3 !定义节点

n,1

n,5,.030 fill,1,5,3 n,13,.070 fill,5,13,7 n,17,.095 fill,13,17,3 n,21,.135 fill,17,21,3 n,31,.240 fill,21,31,9 n,35,.270 fill,31,35,3 n,43,.300 fill,35,43,7 n,47,.330 fill,43,47,3

n,48,.015,-.003 n,49,.255,-.003 !定义单元类型

et,1,4 !梁单元 et,2,14 !弹簧单元

et,3,21,,,4 !集中质量单元 !定义实常数

R,1,pi*R1**2,pi*R1**4/4,pi*R1**4/4,2*R1,2*R1, , RMORE, ,pi*R1**4/2, , , , ,

R,2,pi*R2**2,pi*R2**4/4,pi*R2**4/4,2*R2,2*R2, , RMORE, ,pi*R2**4/2, , , , ,

R,3,pi*R3**2,pi*R3**4/4,pi*R3**4/4,2*R3,2*R3, , RMORE, ,pi*R3**4/2, , , , ,

R,4,pi*R4**2,pi*R4**4/4,pi*R4**4/4,2*R4,2*R4, , RMORE, ,pi*R4**4/2, , , , ,

R,5,pi*R5**2,pi*R5**4/4,pi*R5**4/4,2*R5,2*R5, , RMORE, ,pi*R5**4/2, , , , ,

R,6,pi*R6**2,pi*R6**4/4,pi*R6**4/4,2*R6,2*R6, ,

RMORE, ,pi*R6**4/2, , , , ,

R,7,pi*R7**2,pi*R7**4/4,pi*R7**4/4,2*R7,2*R7, , RMORE, ,pi*R7**4/2, , , , ,

R,8,pi*R8**2,pi*R8**4/4,pi*R8**4/4,2*R8,2*R8, , RMORE, ,pi*R8**4/2, , , , , r,9,1.5e7,0.1 r,10,25.670 type,1 real,1 e,1,2 e,2,3 e,3,4 e,4,5 type,1 real,2 e,5,6 e,6,7 e,7,8 e,8,9 e,9,10 e,10,11 e,11,12 e,12,13 type,1 real,3 e,13,14 e,14,15 e,15,16 e,16,17 type,1 real,4 e,17,18 e,18,19 e,19,20 e,20,21 type,1 real,5 e,21,22 e,22,23 e,23,24 e,24,25 e,25,26 e,26,27 e,27,28

e,28,29 e,29,30 e,30,31 type,1 real,6 e,31,32 e,32,33 e,33,34 e,34,35 type,1 real,7 e,35,36 e,36,37 e,37,38 e,38,39 e,39,40 e,40,41 e,41,42 e,42,43 type,1 real,8 e,43,44 e,44,45 e,45,46 e,46,47 type,2 real,9 e,3,48 e,33,49 type,3 real,10 e,15 d,3,ux d,33,ux d,48,all d,49,all save finish /solu antype,2

MODOPT,lanb,40 EQSLV,FRONT MXPAND,40, , ,1 LUMPM,0

PSTRES,0

MODOPT,lanb,40,10,40000, ,OFF solve finish

1 转子临界转速概念

转子的固有频率除了与转子结构（和支承结构）参数有关外，它还随转子涡动转速和转子 自转转速的变化而变化。在转子不平衡力驱动下，转子一般作正同步涡动，当转子涡动转 速等于转子固有频率时，转子出现共振,相应转速就称为该转子的临界转速。

2 转子临界转速计算对程序的要求

计算转子临界转速必须能够考虑旋转结构涡动时产生的陀螺效应对转子临界转速的影响， 这是转子临界转速计算同其他非旋转结构固有频率计算的差异所在。一般有限元程序不具 备计算转子临界转速的功能。

3 ANSYS的临界转速计算功能 1) 计算转子临界转速可用单元 BEAM4； PIPE16。

COBIN14(用于模拟带阻尼的弹性支撑) 2) 单元特性及实常数 BEAM4和PIPE16:

Keyoption(7)=1

实常数Spin=转子自转角速度（ω） rad/s。 3) 特征值求解方法

选取DAMP方法求解特征值。 4) 计算结果处理

采用有限元方法计算转子临界转速时，转子会出现正进动和反进动。由于陀螺效应的作用 ，随着转子自转角速度的提高，反进动固有频率将降低，而正进动固有频率将提高。根据 临界转速的定义，应只对正进动固有频率（Ωc）进行分析。

在后处理中首先剔除负固有频率，然后分析各阶模态振型，确定同一阶振型的正进动和反 进动固有频率。

改变转子自转角速度（ω），计算出新的Ωc，最后画出Ωc～ω曲线，根据临界转速的定 义，当Ωc=ω时，Ωc即所求临界转速。需注意：由于Ωc的单位为Hz，而ω为rad/s，计算 时应转换单位。

4 算例

单转子结构如图所示，转子轴近似无质量，轮盘密度8*104Kg/m3,其余材料参数为： E=200Gpa μ=0.3

||