ANSYS命令流学习笔记圆柱形shell单元的复合材料分析

ANSYS常用工具命令流简介：学习ANSYS有限元软件三年，积累了很多常用的命令流，适用于不同研究方向的建模，在建模过程中起到小工具的作用。

希望可以帮助到学习ANSYS 的人，之后会尽量持续更新。

为方便使用特地做了目录和简单的用途介绍如下：目录!1.查看面的法线方向：使用voffst命令面偏移创建体之前，查看面的法线方向。

!2.查看点和单元的最大编号。

!3.ansys窗口背景变白!4.ansys出结果图：所出的图位于ansys的计算文件夹中，为png文件。

!5.设置体的透明度：设置plotctrls→numbering仍看不多体单元编号时，通过设置体的透明度可看到位于体内部的编号。

!6.单位制对应表：建模之前先统一单元，局部模拟常用：mm、t、n、mpa。

!7.提取指定位置的节点的节点号：同样适用于面、体、单元等!8.提取某材料的体积：很方便的计算体的方法!9.后处理中显示塑性区!10.显示变形动画并保存:可以用来到处变形动画放在PPT中或者检查加载情况!11.设置变形放大系数：结果中的变形是夸张变形后的结果，利用此命令流可以将变形恢复，用于查看混凝土裂缝时的将变形复原。

!12.查看GUI操作的命令流：查看GUI操作对应命令流，将某些GUI操作保存为命令流可简化操作!13.查看开裂单元透明调整!14.检查形状不好的单元!15.关闭所有标识：关闭ansys右上角标识!16.单元轮廓不显示：应力云图中单元轮廓影响查看效果!17.显示壳单元厚度：壳单元在ansys中以面存在，不显示厚度。

!18.最大最小位移节点编号的提取!19.求解完后显示梁单元的截面!20.显示体、面、线、单元、点：GUI操作的简化!21. ANSYS结果窗格中DMX、SMX、SMN的含义!22.查看节点坐标系：在节点施加荷载之前必须先明确节点坐标系从而确定荷载的正负!23.查看切面应力：查看模型某个截面的应力情况!24.查看荷载列表!24.已选定节点个数统计!25.查看结构整体的质量、体积等!26.查看某个应力范围内的单元命令流!1.查看面的法线方向/PSYMB,adir,1/replot,all/PSYMB,adir,0/replot,allLDELE,215!2.查看点和单元的最大编号单元：esel,all*get,emax,elem,,num,max*get,emin,elem,,num,min点：*get,nmax,node,,num,max*get,nmax,kp,,num,max工具栏：list→ picked entities→ maximumlist→ picked entities→鼠标点选需要查看的点→OK !3.窗口背景变白/COLOR,DEFAULT/RGB,INDEX,100,100,100,0/RGB,INDEX,80,80,80,13/RGB,INDEX,60,60,60,14/RGB,INDEX,0,0,0,15/REPLOT!4.ansys出结果图!GET THE PNG IMAGE/SHOW,PNG,,0PNGR,COMP,1,-1PNGR,ORIENT,HORIZPNGR,COLOR,2PNGR,TMOD,1/GFILE,1500,!*/CMAP,_TEMPCMAP_,CMP,,SAVE/RGB,INDEX,100,100,100,0/RGB,INDEX,0,0,0,15/REPLOT/CMAP,_TEMPCMAP_,CMP/DELETE,_TEMPCMAP_,CMP/SHOW,CLOSE/DEVICE,VECTOR,0!*!5.体单元编号看不到时可设置体的透明度/trlcy,volu,1,all ! 修改所有体的透明度为100% /trlcy,defa!6.单位制对应表!查看单位/status,units!7.提取指定位置的节点的节点号nsel,s,loc,x,5nsel,r,loc,y,6nsel,r,loc,z,7*get,kcon,kp,,num,min也可用list→picked entities，这种方法更直观简便!8.提取某材料的体积vsel,s,mat, ,2vsum,all*get,v2,volu,0,volu!9.后处理中显示塑性区/post1plnsol,eppl,eqv,2!10.显示变形动画并保存GUI:Utility Menu>Plotctrls>Animate>Deformed Shape!11.设置变形放大系数/dscale,,0!12.查看GUI操作的命令流list→files→log file!13.查看开裂单元透明调整/TRLCY,elem,0.5,all!14.检查形状不好的单元check,esel,warn!15.关闭所有标识!16.单元轮廓不显示/GLINE,1,-1/GLINE,1,0 !显示!17.显示壳单元厚度主菜单——PlotCrls——Style——Size and Shape——第二项[ESHAPE]，将off 改为on（必须要在划分网格之后）!18.最大最小位移节点编号的提取allselnsort,u,sum,0,0,all !根据位移矢量大小按降序排列*get,max_u,sort,0,imaxallselnsort,u,sum,0,0,all !根据位移矢量大小按降序排列*get,max_u,sort,0,imin!19.求解完后显示梁单元的截面PlotCtrls-->Style-->Size and Shape-->勾选Shape!20.显示体、面、线、单元、点eplotvplotaplotnplot!21.ANSYS结果输出中DMX、SMX、SMN的含义DMX固定指：最大位移（Displacement Max）SMX：指定的、要查看的、那个Item解的最大值（Solution Max）SMN：指定的、要查看的、那个Item解的最小值（Solution Min）!21.ANSYS结果输出中DMX、SMX、SMN的含义DMX固定指：最大位移（Displacement Max）SMX：指定的、要查看的、那个Item解的最大值（Solution Max）SMN：指定的、要查看的、那个Item解的最小值（Solution Min）!22.查看节点坐标系plotctrls>symbols>nodal coordinate system来源：https:///view/b1a458fc5901020206409c40.html!23.查看切面应力wpstyl,defawprota,0,0,90/type,1,5/cplane,1来源：https:///view/0bab0662f12d2af90242e6c7.html!24.查看荷载列表flist!24.已选定节点个数统计nsel,s,loc,y,0*get,aaa,node,0,count !令aaa=节点数后续在parameters→scalar parameters中查看!25.查看结构整体的质量、体积等preprocessor→modeling→calc geom items→of geometry!26.查看某个应力范围内的单元/cont,1,9,-32, ,-20/replot。

有限元分析软件ANSYS命令流中文说明4 4有限元分析软件ANSYS命令流中文说明4/42010-05-23 21：151设置分析类型ANTYPE，Antype,status,ldstep,action其中antype表示分析类型STATIC：静态分析MODAL：模态分析TRANS：瞬态分析SPECTR：谱分析2 KBC，KEY制定载荷为阶跃载荷还是递增载荷EKY=0递增方式KEY=1阶跃方式3 SOLVE开始一个求解运算4 LSSOLVE读入并求解多个载荷步5 TIME，time设置求解时间有时在分析中需要进入后处理，然后在保持进入后处理之前的状态的情况下接着算下去，可以使用以下的方法：PARSAV,ALL,PAR,TXT！PARSAV命令是储存ANSYS的参数，ALL代表所有参数，PAR是文件名，TXT是扩展名/SOLU ANTYPE,REST,CruStep-1,,CONTINUE！ANTYPE是定义分析类型的命令，REST代表重启动，CruStep代表本载荷步的编号PARRES,NEW,PAR,TXT！PARRES是恢复参数的命令，NEW表示参数是以刷新状态恢复，PAR和TXT 代表了储存了参数的文件名和扩展名如果有单元生死的问题，可以这样处理：ALLSEL,ALL*GET,E_SUM_MAX,ELEM,NUM,MAX！得到单元的最大编号，即单元的总数ESEL,S,LIVE！选中"生"的单元*GET,E_SUM_AL,ELEM,COUNT*DIM,E_POT_AL,E_SUM_MAX！单元选择的指示*DIM,E_NUM_AL,E_SUM_AL！单元编号的数组J=0！读出所选单元号*DO,I,1,E_SUM_MAX*VGET,E_POT_AL(I),ELEM,I,ESEL！对所有单元做循环，被选中的单元标志为"1"*IF,E_POT_AL(I),EQ,1,THEN J=J+1 E_NUM_AL(J)=I*ENDIF*ENDDO ALLSEL,ALL在重启动之后恢复单元生死状态*if,E_SUM_AL,ne,0,then*do,i,1,Num_Alive esel,a,E_NUM_AL(i)*enddo ealive,all allsel*endif/WINDOW,WN,XMIN,XMAX,YMIN,YMAX,NCOPY注意x的坐标是-1到1.67，y坐标是-1到1 Xmin=off on，FULL,LEFT,RIGH,TOP,BOT,LTOP,LBOT,RTOP,RBOT注意一个问题，除了1号窗口外，其他的不能用鼠标操作，只用先发/view 和/dist，然后用/replot。

1、适用PLANE182和SHELL181建立单元时，旋转成体以后要手动将两个单元删除。

而是用MESH200时，系统会在运算时自动删除或关闭该单元状态。

2、EXPOPT体扫掠相关信息的定义，可用于旋转产生体，定义产生体的单元属性、单元尺寸、是否删除源面上的网格、是否自动选择源面和目标面等。

3、施加周向和径向约束需要在柱坐标系下进行。

4、柱坐标系的平面为XOY平面，在切换到柱坐标系时，要注意工作平面是哪个平面，如果图形不满足XOY平面，应建立局部坐标系，调整平面。

5.局部坐标系中蓝色代表Z轴，白色代表X轴，黄色代表Y轴。

6、转速的施加一般是按照弧度进行施加。

7、考虑预应力模态分析的求解，需要先进性预应力效应打开时静力分析，求解一次，然后点击FINISH，之后再进入求解，选择模态分析进行求解。

一、二1.1、采用三维单元分析平面问题时需要约束其Z向的全部位移。

三1、梁柱铰接可以通过两种形式设置（350）：（1）梁柱连接处共用一个节点，使用梁单元弯矩释放功能，将转动自由度和平动自由度释放，可近似实现铰接。

（2）在建立模型是，梁柱不共用节点，通过约束方程，耦合节点位移实现铰接。

使用CP 命令。

2、BEAM188单元画内力图时，设置KEYOPT,1,3,3设置三次形函数可消除弯矩图出现锯齿状。

3、施加重力加速度时，方向与实际方向相反。

4、单元表中坐标系以单元坐标系为准。

四1、求解塑性极限荷载时，荷载大小未知：（1）可以通过理论求解，将大概的理论之求解出来，通过加载距离理论值相近的值进行求解，选取求解不收敛失败前最后一个荷载值作为极限荷载。

（2）通过不断试错，加一个很大的荷载，加到知道计算不收敛，将最大的荷载乘以不收敛的时刻点作为屈服极限荷载。

2、非线性求解中，时间点（TIME）可以写10也可以写1，最终数值要乘以时间点，最好填写1，方便计算。

求解完成后，打开Results Summary中TIME/FREQ列乘以施加力，可得到每一荷载子步施加的力的大小。

ANSYS结构分析单元功能与特性/可以组成一一些命令，一般是一种总体命令（session),三十也有特殊，比如是处理/POST1! 是注释说明符号，，与其他软件的说明是一样的，ansys不作为命令读取，* 此符号一般是APDL的标识符，也就是ansys的参数化语言，如*do ,,,*enddo等等NSEL的意思是node select，即选择节点。

s就是select，选择。

DIM是定义数组的意思。

array 数组。

MP命令用来定义材料参数。

K是建立关键点命令。

K,关键点编号,x坐标,y坐标，z坐标。

K, NPT, X, Y, Z是定义关键点，K是命令，NPT是关键点编号，XYZ是坐标。

NUMMRG, keypoint 用这个命令，要保证关键点的位置完全一样，只是关键点号不一样的才行。

这个命令对于重复的线面都可以用。

这个很简单，压缩关键。

Ngen 复制节点e,节点号码：这个命令式通过节点来形成单元NUMCMP,ALL：压缩所有编号，这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50：通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧，此法会在一段区间内，以一定的步长逐步搜索根，相比常用的牛顿迭代法所要耗费的计算量大得多，但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。

LNSRCH激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。

SPLINE:P1，P2，P3，P4，P5，P6，XV1，YV1，ZV1，XV6，YV6，ZV6（生成分段样条曲线）*DIM，Par，Type，IMAX，JMAX，KMAX，Var1，Var2，Var3（定义载荷数组的名称）【注】Par: 数组名Type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableIMAX，JMAX，KMAX各维的最大下标号Var1，Var2，Var3 各维变量名，缺省为row,column,plane(当type为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, V ALUELab为参数名称value为参数值例如：/config，MXEL，10000的意思是最大单元数为10000杆单元:LINK1、8、10、11、180梁单元：BEAM3、4、23、24，44，54，188，189管单元:PIPE16，17，18，20，59，602D实体元:PLANE2，25，42，82，83，145，146，182，1833D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209弹簧单元:COMBIN7，14，37，39，40质量单元:MASS21接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178矩阵单元:MATRIX27，50表面效应元:SURF153，154粘弹实体元:VISCO88，89，106，107，108， 超弹实体元:HYPER56，58，74，84，86，158耦合场单元:SOLID5，PLANE13，FLUID29，30，38，SOLID62，FLUID79，FLUID80，81， SOLID98，FLUID129，INFIN110，111，FLUID116，130 界面单元:INTER192，193，194，195 显式动力分析单元:LINK160，BEAM161，PLANE162，SHELL163，SOLID164，COMBI16杆单元(Large deflection)，F-大应变(Large strain)或有限应变(Finite strain)，B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应下降(Adaptive descent )等。

要计算固有频率和模态，就必须选择模态分析；可以进行下列类型的分析：结构静力分析、结构动力分析、结构屈曲分析、结构非线性分析、热力学分析、电磁场分析、声场分析、压电分析、流体动力分析。

载荷步与子步定义单位制：/UNITS定义单元类型：ET,1,BEAM4；ET,2,SHELL91定义材料属性：MP,EX,1,4.45E10（材料参考号为1的材料X方向的杨氏模量为4.45E10；MP,DENS,2,7.8E3（材料参考号为2的材料密度为7.8E3）。

若加惯性载荷（如重力），必须定义能求出质量的参数，如密度DENS若施加热载荷，必须定义温度膨胀系数ALPX进入求解器：/SOLU结构分析中，可以将随时间变化的有关变量定义为一位数组，时间作为基本变量，表格的定义：*DIM；Parameters→Array Parameters→Define/Edit位移：UX、UY、UZ、ROTX、ROTY、ROTZ集中力（FX、FY、FZ）和力矩（MX、MY、MZ）表面压力：PRES；温度载荷：TEMP；能量密度：FLUE保存：SA VE；开始求解：SOLVE；推出求解器：FINISH应力：SX、SY、SZ；应变：EPELX、EPEL Y、EPELZ表面载荷：在结构分析中，指施加的压力；体载荷：在结构分析中，有温度和流场两种；惯性载荷有加速度、角加速度、角速度等，惯性载荷只有在模型具有质量时才有效。

→Linearized Strs：用来显示薄膜单元的线性化应力。

Stress→von Mises SEQV：第四强度理论应力子模型：可以让用户把模型的一部分截取后作为一个子模型，重新细分网格，进一步分析。

General Postproc→Submodeling一般单元中的节点以字母I、J、K等表示。

结构分析可进行：静力分析、模态分析、谐波分析、瞬态动力分析、谱分析、屈曲分析、显式动力分析、断裂力学分析、复合材料分析、疲劳分析和p－Method方法。

ANSYS结构解析单元功能与特征/POST1/可以构成一一些命令，一般是一种整体命令（ session),三十也有特别，比方是办理 ! 是说明说明符号，，与其余软件的说明是相同的， ansys 不作为命令读取，*此符号一般是 APDL 的表记符，也就是 ansys 的参数化语言，如 *do ,,,*enddo 等等NSEL 的意思是node select，即选择节点。

s 就是 select，选择。

DIM是定义数组的意思。

array 数组。

MP 命令用来定义资料参数。

K 是建立要点点命令。

K, 要点点编号 ,x 坐标 ,y 坐标， z 坐标。

K, NPT, X, Y , Z 是定义要点点， K 是命令， NPT 是要点点编号， XYZ 是坐标。

NUMMRG , keypoint 用这个命令，要保证要点点的地点完整相同，不过要点点号不一样样的才行。

这个命令关于重复的线面都可以用。

这个很简单，压缩要点。

Ngen 复制节点e,节点号码：这个命令式经过节点来形成单元NUMCMP,ALL ：压缩所有编号，这样你所有的线都会挨次次重新编号 ~你若是需要固定的线固定的标号NSUBST,100,500,50 ：经过指定子步数来设置载荷步的子步LNSRCH 线性搜寻是求解非线性代数方程组的一种技巧，此法会在一段区间内，以必定的步长逐渐搜寻根，对比常用的牛顿迭代法所要耗费的计算量大得多，但它可以防备在一些状况下牛顿迭代法出现的跳跃现象。

LNSRCH激活线性搜寻PRED 激活自由度求解展望NEQIT 指定一个荷载步中的最大子步数AUTOTS自动求解控制打开自动时间步长.KBC -指定阶段状也许用跳板装载里面一个负荷步骤。

SPLINE :P1， P2， P3，P4， P5， P6， XV1 ， YV1 ， ZV1 ， XV6 ，YV6 ， ZV6 （生成分段样条曲线）*DIM ， Par，Type ，IMAX ，JMAX ， KMAX ， Var1，Var2， Var3（定义载荷数组的名称）【注】 Par: 数组名Type： array 数组，仿佛fortran, 下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8 个字符）tableIMAX ， JMAX ， KMAX各维的最大下标号Var1， Var2，Var3 各维变量名，缺省为row,column,plane( 当 type 为 table 时 )/config 是设置 ansys 配置参数的命令格式为 /CONFIG, Lab, V ALUELab 为参数名称value 为参数值比方： /config ， MXEL ，10000 的意思是最大单元数为10000杆单元 : LINK1、 8、 10、 11、 180梁单元： BEAM3、 4、 23、 24，44， 54， 188， 189管单元 : PIPE16， 17， 18， 20， 59， 602D实体元 : PLANE2， 25， 42， 82， 83， 145，146， 182， 1833D实体元 : SOLID45， 46， 64，65， 72， 73，92， 95， 147，148， 185， 186，187， 191壳单元 : SHELL28， 41， 43， 51， 61， 63， 91， 93， 99， 143， 150， 181，208， 209弹簧单元 : COMBIN7， 14， 37，39， 40质量单元 : MASS21接触单元 : CONTAC12， 52， TARGE169， 170， CONTA171， 172， 173， 174， 175， 178矩阵单元 : MATRIX27， 50表面效应元 : SURF153， 154粘弹实体元 : VISCO88， 89， 106， 107， 108，超弹实体元 : HYPER56， 58， 74， 84， 86， 158耦合场单元 : SOLID5， PLANE13， FLUID29， 30，38， SOLID62， FLUID79， FLUID80，81，SOLID98， FLUID129， INFIN110 ， 111， FLUID116，130界面单元 : INTER192， 193， 194， 195显式动力解析单元 : LINK160， BEAM161， PLANE162， SHELL163， SOLID164， COMBI16杆单元单元名称简称节点数节点自由度特征备注LINK12D杆2Ux,Uy EPCSDGB常用杆元LINK83D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉EDGB模拟缆索的废弛及或仅受压杆缝隙LINK113D线性调理EGB模拟液压缸和大转器动LINK1803D有限应变杆EPCDFGB另可考虑粘弹塑性E- 弹性 (Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度deflection)， F- 大应变 (Large strain)或有限应变(Finite strain)，B-单元存亡(Birth and dead),G-化 (Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应降落(Adaptive descent)等。

ANSYS分析实例与工程应用命令流学习笔记1大纲静力分析：2杆、3梁、5薄膜和板壳、4实体单元梁单元：简化计算，结构总体受力情况实体单元：较复杂的结构，局部细节的受力情况稳定性分析：6振动、模态分析：7简单振动和梁的振动、8膜板和实体振动2杆系结构的静力分析2.1铰接杆在外力作用下的变形二维杆单元LINK1*AFUN,DEG：三角函数默认为弧度，改为角度后处理：结构变形图、显示节点位移和杆件应力2.2人字形屋架的静力分析后处理：杆单元的轴力、轴向应力、轴向应变2.3超静定拉压杆的反力计算后处理：节点反力2.4平行杆件与刚性梁连接的热应力问题定义3点的UY为耦合自由度，即三者的UY位移相等温度（增量）后处理：寻找特定位置的节点和单元，并从单元表中提取它们的内力2.5端部有间隙的杆的热膨胀二维带厚度的平面应力单元PLANE42、二维接触单元CONTACT26温度（始、末）后处理：定义水平应力和铅直应力单元表，并提取3号单元的应力结果*Status,ParmFINISH定义数组变量，将计算结果通过数组变量输出到文件3梁的弯曲静力分析3.1单跨等截面超静定梁的平面弯曲二维弹性梁单元BEAM3后处理：定义以两端弯矩和剪力的单元表，并列出单元表数据并用单元表数据绘制剪力图和弯矩图更细的节点划分方案，更精细3.2四跨连续梁的内力计算体素建模：keypoint, line, area, volume便于细分单元3.3七层框架结构计算3.4工字形截面外伸梁的平面弯曲3.5矩形截面梁的纵横弯曲分析考虑应力强化效应后处理：迭代过程3.6空间刚架静力分析三维梁单元BEAM43.7悬臂梁的双向弯曲三维8节点耦合场实体单元SOLID5三维20节点固体单元SOLID92三维10节点耦合场实体单元SOLID98三维结构实体自适应单元SOLID147定义宏程序，对应四种工况，各种结果差别不大3.8圆形截面悬臂杆的弯扭组合变形三维直管单元PIPE16（只定义外直径，不定义内直径）3.9悬臂等强度梁的弯曲四边形壳单元SHELL63（这里用退化的三角形单元，并使用节点耦合自由度保证模型的对称变形）三维非对称锥形梁单元BEAM44（定义横截面主轴，单元宽度线性变化）计算结果都很好，但壳体单元更能模拟出等强度梁的实际几何形状，更直观，截面定义更简单。

本文介绍了轮毂的ANSYS APDL命令流建模及模态分析实例相关内容。

ANSYS命令流及注释五个辐条的轮毂!!初始化ANSYS环境!FINISH/CLEAR !清空内存/FILNAM,WHEEL5 !文件名/TITILE,WHEEL5 PARAMETER MODELING !工作名!!定义几何尺寸参数!R1=180R2=157R3=75R4=75R5=30R6=28R7=20R8=90R9=60S_HOLE=5TH1=48TH2=23TH3=11TH4=180TH5=40TH6=45TH7=105TH8=25TH9=15TH10=25TH11=13/VIEW,1,1,1,1 !改变视图/ANG,1/PNUM,LINE,1/PNUM,AREA,1/PNUM,VOLU,1/NUMBER,1!!关键点!/PREP7k,1,r5,r7,0k,2,r4-ky(1),ky(1),0k,3,r4,0,0k,4,r1,0,0k,5,kx(4),th5-th9,0k,6,r1-th8,ky(5),0k,7,kx(6),th4/2,0k,8,kx(7)+th11,ky(7)+th10,0 k,9,kx(8),th4-th3,0k,10,kx(4),ky(9),0k,11,kx(4),th4,0k,12,r2,ky(11),0k,13,kx(12),ky(8),0k,14,kx(7)-th3,ky(7),0k,15,kx(14),th5,0k,16,r3+r6,ky(15),0k,17,kx(3),r7+th1,0k,18,kx(1),ky(17),0k,19,kx(16),ky(17),0k,20,kx(2),0,0k,21,0,0,0k,22,0,th1+r7,0*ask,s_hole,'the number of hole',5 !宏!!创建轮毂面!lstr,1,2 !连接1，2关键点，形成直线larc,2,3,20,r7 !以20点为圆心r7为半径，2，3点为端点作弧线lstr,3,4lstr,4,5lstr,5,6lstr,6,7lstr,7,8lstr,8,9lstr,9,10lstr,10,11lstr,11,12lstr,12,13lstr,13,14lstr,14,15lstr,15,16larc,16,17,19,r6lstr,17,18lstr,18,1al,allcm,an-all,area !形成组件!!创建实体模型!allsel,allvrotat,an-all,,,,,,21,22,360,S_hole, !旋转拉伸形成体cm,v-an-all,volu!!减去孔洞!vsel,nonewpro,,-90, !绕Y轴转动工作平面cswpla,11,1,1,1csys,11wpoff,r8*sin(180/s_hole),r8*cos(180/s_hole)RPR4,3,-th5,th5/2,r9,, !创建三角形adele,96LFILLT,182,181,10, , !在直线182，181间形成半径10的圆角LFILLT,182,183,10, ,LFILLT,183,181,10, ,LARC,98,100,21,144,ldele,182asel,noneal,181,184,187,185,183,186 !连接各线形成面cm,sanjiao_hole,areavext,sanjiao_hole,,,0,0,th5,,,, !以th5为厚度形成体cm,v_hole,voluvgen,s_hole,all,,,,360/s_hole,,,0 !旋转拉伸形成s_hole个体cm,v-hole,voluvsel,allvsbv,v-an-all,v-hole !布尔运算减去体，形成孔洞cm,v-an-all,voluALLSEL,ALL!!定义单元属性!et,1,solid45mp,ex,1,71000 !铝合金材料特性mp,nuxy,1,0.33mp,dens,1,2720!!划分单元创建网格模型!SMRT,5 !自由网格划分MSHAPE,1,3DMSHKEY,0FLST,5,5,6,ORDE,2FITEM,5,11FITEM,5,-15CM,_Y,VOLUVSEL, , , ,P51XCM,_Y1,VOLUCHKMSH,'VOLU'CMSEL,S,_YVMESH,_Y1finish!!保存!saveAPLOT/SOLUFLST,2,5,5,ORDE,5 !约束固定FITEM,2,19FITEM,2,38FITEM,2,57FITEM,2,76FITEM,2,95/GODA,P51X,ALL,*DEL,_FNCNAME !函数加载*DEL,_FNCMTID*DEL,_FNC_C1*DEL,_FNCCSYS*SET,_FNCNAME,'jiazai'*DIM,_FNC_C1,,1*SET,_FNC_C1(1),5*SET,_FNCCSYS,11! /INPUT,111.func,,,1*DIM,%_FNCNAME%,TABLE,6,7,1,,,,%_FNCCSYS% !! Begin of equation: 1000*{X}/cos(180/s_hole)*SET,%_FNCNAME%(0,0,1), 0.0, -999*SET,%_FNCNAME%(2,0,1), 0.0*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%*SET,%_FNCNAME%(4,0,1), 0.0*SET,%_FNCNAME%(5,0,1), 0.0*SET,%_FNCNAME%(6,0,1), 0.0*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 1000, 0, 0, 2 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 2*SET,%_FNCNAME%(0,3,1), 0, -1, 0, 180, 0, 0, 17*SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, -1, 4, 17 *SET,%_FNCNAME%(0,5,1), 0.0, -1, 10, 1, -3, 0, 0 *SET,%_FNCNAME%(0,6,1), 0.0, -3, 0, 1, -2, 4, -1 *SET,%_FNCNAME%(0,7,1), 0.0, 99, 0, 1, -3, 0, 0 ! End of equation: 1000*{X}/cos(180/s_hole) FLST,2,3,1,ORDE,3 !确定加载点位置FITEM,2,37FITEM,2,54FITEM,2,354/GOF,P51X,FX, %JIAZAI%/STA TUS,SOLU !求解SOLVE/VIEW,1,1,1,1/ANG,1/REP,FAST/SOLUANTYPE,2 !模态求解MSA VE,0MODOPT,LANB,10EQSLV,SPARMXPAND,10, , ,1LUMPM,0PSTRES,0MODOPT,LANB,10,0,0, ,OFF/STA TUS,SOLUSOLVEFINISHSave模型图网格划分位移图应变图应力图应力模态(其中之一)。