ANSYS第一类完全椭圆积分宏命令，替代查表法

在受压柱大挠度理论分析中，要用到第一类完全椭圆积分计算，可以查表得到其数值，但实在太麻烦了。于是发挥ANSYS之APDL特长，特编制了宏命令可供应用，从而方便计算。

计算方法采用级数法实现，其计算公式如下：

宏命令及其应用如下：

!---------------------------------

!第一类完全椭圆积分宏命令FAN.MAC

!采用级数近似解

!使用方法：FAN,K,N

!K---积分式中的常数(ARG1)

!N---级数项数（缺省为600）(ARG2)

!FANV---返回计算结果

!在宏中K用ARG1表示

!循环变量用ARG3

!-------------------------------------

*IF,ARG2,LT,600,THEN

ARG2=600

*ENDIF

FZM=1/2

FANV=1 ARG1*ARG1/4

*DO,ARG3,2,ARG2

FZM=FZM*(2*ARG3-1)/(2*ARG3)

FANV=FANV (FZM)**2*ARG1**(2*ARG3)

*ENDDO

FANV=FANV*ACOS(-1)/2

FZM=

!==========================

!宏命令使用举例

*DIM,BG,,18

*DO,I,1,18

CITA=(I-1)*10/180*ACOS(-1)

AI=SIN(CITA/2)

FAN,AI,800!调用宏命令计算

BG(I)=FANV!保存数值

*ENDDO!

*STAT,BG

未去查表，只与某参考资料的计算表格对比如上表。

从计算结果可知，宏命令的计算结果完全满足要求，且极其方便（程序计算必然较差别方便了）。该FAN.MAC编制有两个技巧，一是项数从2开始，借第一项形成阶乘；第二是边计算边利用，不做阶乘形成大数值，例如分子或分母的阶乘，否则会超过APDL对大数值的限制。

该FAN.MAC宏命令放在ANSYS工作目录，可像ANSYS内部命令一样任意使用，也可放在宏库文件中使用。

因此，某些计算完全可以用APDL去完成，毕竟APDL与高级语言程序FORTRAN类似，也是一种编程语言，如能灵活运用APDL就可提高效率。当然这是在使用ANSYS软件的基础上说的，此类计算也可采用其他编程语言或工具完成，但需要熟悉其他语言或工具，与其去学习其他语言或工具，倒不如直接采用APDL实现，还可提高APDL编程水平和技巧。

---------------------------------------------------------------------------------------------

版权声明：

原创文章，来源铁大王新敏，本文已经授权，欢迎分享，如需转载请联系作者。