本文介绍了轮毂的ANSYS APDL命令流建模及模态分析实例相关内容。
ANSYS命令流及注释
五个辐条的轮毂
!
!初始化ANSYS环境
!
FINISH
/CLEAR !清空内存
/FILNAM,WHEEL5 !文件名
/TITILE,WHEEL5 PARAMETER MODELING !工作名
!
!定义几何尺寸参数
!
R1=180
R2=157
R3=75
R4=75
R5=30
R6=28
R7=20
R8=90
R9=60
S_HOLE=5
TH1=48
TH2=23
TH3=11
TH4=180
TH5=40
TH6=45
TH7=105
TH8=25
TH9=15
TH10=25
TH11=13
/VIEW,1,1,1,1 !改变视图/ANG,1
/PNUM,LINE,1
/PNUM,AREA,1
/PNUM,VOLU,1
/NUMBER,1
!
!关键点
!
/PREP7
k,1,r5,r7,0
k,2,r4-ky(1),ky(1),0
k,3,r4,0,0
k,4,r1,0,0
k,5,kx(4),th5-th9,0
k,6,r1-th8,ky(5),0
k,7,kx(6),th4/2,0
k,8,kx(7)+th11,ky(7)+th10,0 k,9,kx(8),th4-th3,0
k,10,kx(4),ky(9),0
k,11,kx(4),th4,0
k,12,r2,ky(11),0
k,13,kx(12),ky(8),0
k,14,kx(7)-th3,ky(7),0
k,15,kx(14),th5,0
k,16,r3+r6,ky(15),0
k,17,kx(3),r7+th1,0
k,18,kx(1),ky(17),0
k,19,kx(16),ky(17),0
k,20,kx(2),0,0
k,21,0,0,0
k,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,4
lstr,4,5
lstr,5,6
lstr,6,7
lstr,7,8
lstr,8,9
lstr,9,10
lstr,10,11
lstr,11,12
lstr,12,13
lstr,13,14
lstr,14,15
lstr,15,16
larc,16,17,19,r6
lstr,17,18
lstr,18,1
al,all
cm,an-all,area !形成组件
!
!创建实体模型
!
allsel,all
vrotat,an-all,,,,,,21,22,360,S_hole, !旋转拉伸形成体
cm,v-an-all,volu
!
!减去孔洞
!
vsel,none
wpro,,-90, !绕Y轴转动工作平面
cswpla,11,1,1,1
csys,11
wpoff,r8*sin(180/s_hole),r8*cos(180/s_hole)
RPR4,3,-th5,th5/2,r9,, !创建三角形
adele,96
LFILLT,182,181,10, , !在直线182，181间形成半径10的圆角LFILLT,182,183,10, ,
LFILLT,183,181,10, ,
LARC,98,100,21,144,
ldele,182
asel,none
al,181,184,187,185,183,186 !连接各线形成面
cm,sanjiao_hole,area
vext,sanjiao_hole,,,0,0,th5,,,, !以th5为厚度形成体
cm,v_hole,volu
vgen,s_hole,all,,,,360/s_hole,,,0 !旋转拉伸形成s_hole个体cm,v-hole,volu
vsel,all
vsbv,v-an-all,v-hole !布尔运算减去体，形成孔洞
cm,v-an-all,volu
ALLSEL,ALL
!
!定义单元属性
!
et,1,solid45
mp,ex,1,71000 !铝合金材料特性mp,nuxy,1,0.33
mp,dens,1,2720
!
!划分单元创建网格模型
!
SMRT,5 !自由网格划分MSHAPE,1,3D
MSHKEY,0
FLST,5,5,6,ORDE,2
FITEM,5,11
FITEM,5,-15
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
finish
!
!保存
!
save
APLOT
/SOLU
FLST,2,5,5,ORDE,5 !约束固定
FITEM,2,19
FITEM,2,38
FITEM,2,57
FITEM,2,76
FITEM,2,95
/GO
DA,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,37
FITEM,2,54
FITEM,2,354
/GO
F,P51X,FX, %JIAZAI%
/STA TUS,SOLU !求解
SOLVE
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
/SOLU
ANTYPE,2 !模态求解
MSA VE,0
MODOPT,LANB,10
EQSLV,SPAR
MXPAND,10, , ,1
LUMPM,0
PSTRES,0
MODOPT,LANB,10,0,0, ,OFF
/STA TUS,SOLU
SOLVE
FINISH
Save
模型图
网格划分
位移图
应变图
应力图
应力模态(其中之一)
全面讲解ANSYS高速旋转轮盘模态分析讲解 高速旋转轮盘模态分析 在进行高速旋转机械的转子系统动力设计时，需要对转动部件进行模态分析，求解出其固有频率和相应的模态振型。通过合理的设计使其工作转速尽量远离转子系统的固有频率。而对于高速部件，工作时由于受到离心力的影响，其固有频率跟静止时相比会有一定的变化。为此 ，在进行模态分析时需要考虑离心力的影响。我通过该例子学习到了如何用ANSYS进行有预应力的结构的模态分析。 一．例子描述 本例子是对某高速旋转轮盘进行考虑离心载荷引起的预应力的模态分析，求解出该轮盘的 前10阶固有频率及其对应的模态振型。轮盘截面形状如图1所示，该轮盘安装在某转轴上以120 00转/分的速度高速旋转。相关参数为：弹性模量EX＝2.1E5Mpa，泊松比PRXY＝0.3，密度DE NS＝7.8E-9T/mm3。 图1、轮盘截面图 1-5关键点坐标： 1(-10, 150, 0) 2(-10, 140, 0) 3(-3, 140, 0) 4(-4, 55, 0)
5(-15, 40, 0) L=15 RS=5 二．A nsys求解的具体步骤 1.启动ansys，定义工作名、工作标题 ①定义工作名：Example of dynamic ②工作标题：dynamic analysis of a disc 2、选择单元类型 本例将选用六面体结构实体单元来分析，但在建模过程中需要使用四边形平面单元，所有需要定义两种单元类型：PLANE42和SOLID45，设置完成后，如图2，在Element Types (单元类型定义)对话框的列表框中将会列出刚定义的两种单元类型：PLANE42、 SOLID45， 图2、定义单元类型 3、设置材料属性 由于要进行的是考虑离心力引起的预应力作用下的轮盘的模态分析，材料的弹性模量EX 和密度DENS必须定义。 ①定义材料的弹性模量EX 弹性模量 EX＝2.1E5 泊松比 PRXY＝0.3 ②定义材料的密度DENS DENS =7.8E-9 4、建立实体模型 对于本实例的有限元模型，首先需要建立轮盘的截面几何模型，然后对其进行网格划分，最后通过截面的有限元网格扫描出整个轮盘的有限元模型。具体的操作过程如下。
/PREP7 !进入前处理 ANTYPE,STATIC !设置分析类型为静力结构分析 PSTRES,ON !用于后面的模态分析中考虑预应力(该开关不影响静力分析) ET,1,LINK10 !选取单元类型1(单向杆单元) KEYOPT,1,3,0 !设置仅承受拉应力，KEYOPT(3)=0 R,1,306796E-8,543248E-8 !设置实常数，包括绳索截面积(306796E-8)，初始应变(543248E-8) MP,EX,1,30E6 !定义材料的弹性模量(1号材料) MP,DENS,1,73E-5 !定义材料的密度(1号材料) N,1 ! 定义第1号节点 N,14,100 ! 定义第14号节点 FILL ! 均分填满第2号至第13号节点 E,1,2 !由节点1及节点2生成单元 EGEN,13,1,1 !依序复制生成13个单元 D,ALL,ALL ! 对所有节点施加固定约束 FINISH ! 前处理结束 /SOLU ! 进入求解模块，求解预应力引起的应力状态 SOLVE ! 求解 FINISH ! 退出求解模块 /POST1 ! 进入一般的后处理 ETABLE,STRS,LS,1 !针对LINK10单元,建立单元列表STRS,通过LS及特征号1来获得单元的轴向应力 *GET,STRSS,ELEM,13,ETAB,STRS !针对单元列表STRS, 提取13号单元的应力 FINISH ! 后处理结束 /POST26 ! 进入时间历程后处理，处理支反力 RFORCE,2,1,F,X !将1号节点上的x方向支反力提取，并存储到2号变量中 STORE ! 存储 *GET,FORCE,V ARI,2,EXTREM,VMAX !将2号变量的最大值赋给参数FORCE /SOLU ! 再次进入求解模块，模态分析 ANTYPE,MODAL ! 模态分析 MODOPT,SUBSP,3 ! 选择子空间迭代法，求3阶模态 MXPAND,3 ! 设定3阶模态扩展 PSTRES,ON ! 用于在模态分析中考虑预应力(还需在前面的静力分析中也同时打开) DDELE,2,UX,13 ! 删除从2号节点到13号节点上的UX约束 DDELE,2,UY,13 !删除从2号节点到13号节点上的UY约束 SOLVE !求解 *GET,FREQ1,MODE,1,FREQ ! 提取第1阶模态共振频率，并赋值给参数FREQ1 *GET,FREQ2,MODE,2,FREQ ! 提取第2阶模态共振频率，并赋值给参数FREQ2 *GET,FREQ3,MODE,3,FREQ ! 提取第3阶模态共振频率，并赋值给参数FREQ3 *STATUS !列出所有参数的实际内容
压电变换器的自振频率分析及详细过程 1.模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型)，即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS提供的模态提取方法有:子空间法(subspace)、分块法(block lancets),缩减法(reduced/householder)、动态提取法(power dynamics)、非对称法(unsymmetric),阻尼法(damped), QR阻尼法(QR damped)等，大多数分析都可使用子空间法、分块法、缩减法。 ANSYS的模态分析是线形分析，任何非线性特性，例如塑性、接触单元等，即使被定义了也将被忽略。 2.模态分析操作过程 一个典型的模态分析过程主要包括建模、模态求解、扩展模态以及观察结果四个步骤。 (1).建模 模态分析的建模过程与其他分析类型的建模过程是类似的，主要包括定义单元类型、单元实常数、材料性质、建立几何模型以及划分有限元网格等基本步骤。 (2).施加载荷和求解 包括指定分析类型、指定分析选项、施加约束、设置载荷选项，并进行固有频率的求解等。 指定分析类型，Main Menu- Solution-Analysis Type-New Analysis,选择Modal。 指定分析选项，Main Menu-Solution-Analysis Type-Analysis Options，选择MODOPT(模态提取方法〕，设置模态提取数量MXPAND. 定义主自由度，仅缩减法使用。 施加约束，Main Menu-Solution-Define Loads-Apply-Structural-Displacement。 求解，Main Menu-Solution-Solve-Current LS。 (3).扩展模态 如果要在POSTI中观察结果，必须先扩展模态，即将振型写入结果文件。过程包括重新进入求解器、激话扩展处理及其选项、指定载荷步选项、扩展处理等。 激活扩展处理及其选项，Main Menu-Solution-Load Step Opts-Expansionpass-Single Expand-Expand modes。 指定载荷步选项。 扩展处理，Main Menu-solution-Solve-Current LS。 注意:扩展模态可以如前述办法单独进行，也可以在施加载荷和求解阶段同时进行。本例即采用了后面的方法 (4).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等
Ansys命令流 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,,
K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面 ...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识 /SOLU !进入求解模块标识 !施加约束和载荷 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放 /DEVICE !设置图例的显示,如:风格,字体等 /REPLOT !重新显示当前图例 /RESET !恢复缺省的图形设置 /VIEW !设置观察方向 /ZOOM !对图形显示窗口的某一区域进行缩放
/FILNAME,Elastic-Plasitc,1 /TITLE, Elastic-Plasitc Analysis !前处理。 /PREP7 !**定义梁单元189。 ET,1,BEAM189 !定义单元。 !**梁截面1。 SECTYPE, 1, BEAM, HREC, , 0 !定义梁截面。SECOFFSET, CENT SECDATA,50,100,6,6,6,6,0,0,0,0 !定义梁截面完成。 !**定义材料。 MPTEMP,,,,,,,, !定义弹塑性材料模型。MPTEMP,1,0 MPDATA,EX,1,,2.05e5 MPDATA,PRXY,1,,0.3 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,150,18600,,,, !定义弹塑性材料模型。!**建立几何模型。 K,1, , , , K,2 ,900, K,3 ,,50 LSTR, 1, 2 !**网格划分。 FLST,5,1,4,ORDE,1 !定义网格密度。FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义网格密度完成。CM,_Y,LINE !网格划分。 LSEL, , , , 1 CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , 3, ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 LMESH, 1 !网格划分完成。 !施加载荷及求解。 FINISH /SOL
!**施加约束。 FLST,2,1,3,ORDE,1 !施加约束。FITEM,2,1 /GO DK,P51X, , , ,0,UX,UY,UZ,ROTX, , , FLST,2,1,3,ORDE,1 FITEM,2,2 /GO DK,P51X, , , ,0,UY,UZ,ROTX, , , , !施加约束完成。 !**加载。 FLST,2,50,2,ORDE,2 FITEM,2,1 FITEM,2,-50 SFBEAM,P51X,1,PRES,100, , , , , , LSWRITE,1, !定义载荷步1完成。FLST,2,50,2,ORDE,2 !定义载荷步2。FITEM,2,1 FITEM,2,-50 SFEDELE,P51X,1,PRES LSWRITE,2, !定义载荷步2完成。!设定求解步并求解。 LSSOLVE,1,2,1,
高速旋转轮盘模态分析 在进行高速旋转机械的转子系统动力设计时，需要对转动部件进行模态分析，求解出其固有频率和相应的模态振型。通过合理的设计使其工作转速尽量远离转子系统的固有频率。而对于高速部件，工作时由于受到离心力的影响，其固有频率跟静止时相比会有一定的变化。为此，在进行模态分析时需要考虑离心力的影响。通过该实验掌握如何用ANSYS进行有预应力的结构的模态分析。 一．问题描述 本实验是对某高速旋转轮盘进行考虑离心载荷引起的预应力的模态分析，求解出该轮盘的前5阶固有频率及其对应的模态振型。轮盘截面形状如图所示，该轮盘安装在某转轴上以12000转/分的速度高速旋转。相关参数为：弹性模量EX＝2.1E5Mpa，泊松比PRXY＝0.3， 密度DENS＝7.8E-9Tn/mm 3。 1-5关键点坐标： 1(-10, 150, 0) 2(-10, 140, 0) 3(-3, 140, 0) 4(-4, 55, 0) 5(-15, 40, 0) L=10+（学号×0.1） RS=5 二．分析具体步骤 1.定义工作名、工作标题、过滤参数 ①定义工作名：Utility menu > File > Jobname ②工作标题：Utility menu > File > Change Title（个人学号） 2.选择单元类型 本实验将选用六面体结构实体单元来分析，但在建模过程中需要使用四边形平面单元，所有需要定义两种单元类型：PLANE42和SOLID45，具体操作如下： Main Menu >Preprocessor > Element Type > Add/Edit/Delete
①“ Structural Solid”→“ Quad 4node 42” →Apply（添加PLANE42为1号单元） ②“ Structural Solid”→“ Quad 8node 45” →ok（添加六面体单元SOLID45为2号单元） 在Element Types (单元类型定义)对话框的列表框中将会列出刚定义的两种单元类型：PLANE42、SOLID45，关闭Element Types (单元类型定义)对话框，完成单元类型的定义。 3.设置材料属性 由于要进行的是考虑离心力引起的预应力作用下的轮盘的模态分析，材料的弹性模量EX 和密度DENS必须定义。 ①定义材料的弹性模量EX Main Menu >Preprocessor > Material Props > Material Models> Structural > Linear > Elastic >Isotropic 弹性模量EX＝2.1E5 泊松比PRXY＝0.3 ②定义材料的密度DENS Main Menu >Preprocessor > Material Props > Material Models>density DENS =7.8E-9 4.实体建模 对于本实例的有限元模型，首先需要建立轮盘的截面几何模型，然后对其进行网格划分，最后通过截面的有限元网格扫描出整个轮盘的有限元模型。具体的操作过程如下。 ①创建关键点操作：Main Menu > Preprocessor > Modeling > Create > Keypoints > In Active CS 列出各点坐标值Utility menu >List > Keypoints >Coordinate only
ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT
1 旋转摩擦 (1) 2. 电磁三d命令流实例（论坛看到） (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19) 5. 如何施加恒定的角速度？Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子（二维旋转） (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0
模态分析步骤 第1步： 载入模型Plot>Volumes 第2步： 指定分析标题并设置分析范畴 1设置标题等Utility Menu>File>Change Title Utility Menu>File> Change Jobname Utility Menu>File>Change Directory 2选取菜单途径MainMenu>Preference ,单击Structure,单击OK第3步： 定义单元类型 MainMenu>Preprocessor>ElementType>Add/Edit/Delete,出现Element Types 对话框,单击Add出现Library of Element Types对话框,选择Structural Solid,再右滚动栏选择Brick 20node 95,然后单击OK，单击Element Types对话框中的Close 按钮就完成这项设置了。 第4步： 指定材料性能 选取菜单途径MainMenu>Preprocessor>MaterialProps>MaterialModels。出现DefineMaterialModelBehavior对话框,在右侧Structural>Linear>Elastic>Isotropic,指定材料的弹性模量和泊松系数，Structural>Density指定材料的密度，完成后退出即可。 第5步： 划分网格
选取菜单途径Main Menu>Preprocessor>Meshing>MeshTool,出现MeshTool 对话框，一般采用只能划分网格，点击SmartSize,下面可选择网格的相对大小（太小的计算比较复杂，不一定能产生好的效果，一般做两三组进行比较），保留其他选项，单击Mesh出现Mesh Volumes对话框，其他保持不变单击Pick All，完成网格划分。 第6步： 进入求解器并指定分析类型和选项 选取菜单途径Main Menu>Solution>Analysis Type>New Analysis，将出现New Analysis对话框,选择Modal单击OK。 选取Main Menu>Solution> Analysis Type>Analysis Options，将出现Modal Analysis对话框，选中Subspace模态提取法，在Number ofmodes to extract处输入相应的值（一般为5或10，如果想要看更多的可以选择相应的数字），单击OK，出现Subspace Model Analysis对话框，选择频率的起始值，其他保持不变，单击OK。 第7步： 施加边界条件.选取 MainMenu>Solution>Defineloads>Apply>Structural>Displacement,出现 ApplyU,ROTonKPS对话框，选择在点、线或面上施加位移约束,单击OK会打开约束种类对话框，选择（AllDOF,UX,UY,UZ）相应的约束,单击apply或OK即可。 第8步： 指定要扩展的模态数。选取菜单途径 MainMenu>Solution>LoadStepOpts>ExpansionPass>ExpandModes,出现Expand Modes对话框,在number of modes to expand处输入第6步相应的数字，单击OK 即可。（当选取MainMenu>Solution>AnalysisType>AnalysisOptions，将出现ModalAnalysis对话框，选中Subspace模态提取法，在Number of modes to extract处输入相应的值（一般为5或10，如果想要看更多的可以选择相应的数字），同时选择number of modes to expand输入相应值时，这步可以省略）
/FILNAME, Beam,1 ！定义工作文件名。 /TITLE, Beam Analysis ！定义工作标题。/PREP7 !定义单元。 ET,1,BEAM188 !定义材料属性。 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.1e5 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7.9e-6 ! 定义杆件截面■200。 SECTYPE, 1, BEAM, RECT, , 0 SECOFFSET, CENT SECDATA,10,10,0,0,0,0,0,0,0,0 !建立几何模型。 K,1, ,, , K,2,350,, , !生成立柱。 LSTR, 1, 2 !以上完成几何模型。 !以下进行网格划分。 FLST,5,1,4,ORDE,1 FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义单元大小。!分配、划分平板结构。 LMESH, 1 !分析类型施加载荷并求解。 FINISH /SOLUTION ANTYPE,2 ！定义分析类型及求解设置。MSAVE,0 ！模态提取方法。
MODOPT,LANB,10 EQSLV,SPAR MXPAND,10, , ,0 ！模态扩展设置。 LUMPM,0 PSTRES,0 MODOPT,LANB,10,0,0, ,OFF MXPAND,10,0,0,1,0.001, !施加约束。 FLST,2,2,3,ORDE,2 FITEM,2,1 FITEM,2,-2 /GO DK,P51X, , , ,0,ALL, , , , , , !求解。 FINISH /SOL /STATUS,SOLU SOLVE !定义谱分析。 FINISH /SOLUTION ANTYPE,8 SPOPT,PSD,10,1 PSDUNIT,1,DISP,386.4, PSDFRQ,1, ,13.8,40,50.6,73,120 !定义谱—频率表。PSDFRQ,1, ,134,178,233, , PSDV AL,1,1,4,0.6,3,5 PSDV AL,1,6,2,6, , FLST,2,2,1,ORDE,2 !施加谱。 FITEM,2,1 FITEM,2,-2 D,P51X,UX,1.0 PFACT,1,BASE, !计算PSD激励参与系数。PSDRES,DISP,REL !设置输出选项。PSDRES,VELO,OFF PSDRES,ACEL,OFF
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
三 生成关键点和线部分 1.生成关键点 K,关键点编号,X坐标,Y坐标,Z坐标 例:K,1,0,0,0 2.在激活坐标系生成直线 LSTR,关键点P1,关键点P2 例LSTR,1,2 3.在两个关键点之间连线 L,关键点P1,关键点P2 例L,1,2 注:此命令会随当前的激活坐标系不同而生成直线或弧线 4.由三个关键点生成弧线 LARC,关键点P1,关键点P2,关键点PC,半径RAD 例LARC,1,3,2,0.05 注:关键点PC是用来控制弧线的凹向 5.通过圆心半径生成圆弧
CIRCLE,关键点圆心,半径RAD,,,,圆弧段数NSEG 例:CIRCLE,1,0.05,,,,4 6.通过关键点生成样条线 BSPLIN,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6 例:BSPLIN,1,2,3,4,5,6 7.生成倒角线 LFILLT,线NL1,线NL2,倒角半径RAD 例LFILLT,1,2,0.005 8.通过关键点生成面 A,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6,P7,P8... 例:A,1,2,3,4 9.通过线生成面 AL,线L1,线L2,线L3,线L4,线L5,线L6,线L7,线L8,线L9,线L10 例:AL,5,6,7,8 10.通过线的滑移生成面
ASKIN,线NL1,线NL2,线NL3,线NL4,线NL5,线NL6,线NL7,线NL8,线NL9 例:ASKIN,1,4,5,6,7,8 注:线1为滑移的导向线 四 目标:掌握常用的实体-面的生成 生成矩形面 1.通过矩形角上定位点生成面 BLC4,定位点X方向坐标XCORNER,定位点Y方向坐标YCORNER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 例:BLC4,0,0,5,3,0 2.通过矩形中心定位点生成面 BLC5,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 注:与上条命令的不同就在于矩形的定位点不一样 例:BLC5,2.5,1.5,5,3,0 3.通过在工作平面定义矩形X.Y坐标生成面 RECTNG,矩形左边界X坐标X1,矩形右边界X坐标X2,矩形下边界Y
大桥全长2996.8m，其中主桥采用跨度为101.5+188.5+580+217.5+159.5+116m的钢桁梁斜拉桥；非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁；东引桥采用16孔梁长32.6m预应力混凝土简支箱梁；跨大堤桥采用48.9+86+48.8m预应力混凝土连续箱梁；西引桥采用15孔梁长32.6m预应力混凝土简支梁及2孔梁长24.6m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁。 主桥采用103+188.5+580+217.5+159.5+117.5m两塔钢桁斜拉桥方案，全长1366m。主梁为三片主桁钢桁梁，桁间距2x14m，节间长14.5m，桁高15m。主塔为钢筋混凝土结构，塔顶高程+204.00m，塔底高程-6.00m，斜拉索为空间三索面，立面上每塔两侧共18对索，全桥216根斜拉索。所有桥墩上均设竖向和横向约束，4#塔与主梁之间设纵向水平约束，3#塔与梁间使用带限位功能的粘滞阻尼器。主梁为”N”字型桁式，横向采用三片桁结构，主桁的横向中心距各为14m，桁高15m，节间距14.5m[2]。 结构构造 主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面，上弦杆内高1000mm，内宽1200mm，板厚20～48mm。下弦杆内高1400mm，宽1200mm，板厚20～56mm。下弦杆顶板向桁内侧加宽700mm与整体桥面板焊接。腹杆主要采用H型截面。H型杆件宽1200mm，高720和760mm，板厚20～48mm。根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。主桁采用焊接杆件，整体节点。在节点外以高强度螺栓拼接的结构形式，上下弦杆四面等强对接拼装。H型腹杆采用插入式连接。箱型腹杆采用四面与主桁节点对拼的连接形式。主桁拼接采用M30高强螺栓。
第一章模态分析 §模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性（固有频率和振型），即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性，如塑性和接触（间隙）单元，即使定义了也将被忽略。ANSYS提供了七种模态提取方法，它们分别是子空间法、分块Lanczos法、PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼。后面将详细介绍模态提取方法。 §模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样，无论进行何种类型的分析，均可从用户图形界面（GUI）上选择等效于命令的菜单选项来建模和求解问题。 后面的“模态分析实例（命令流或批处理方式）”将给出进行该实例模态分析时要输入的命令（手工或以批处理方式运行ANSYS时）。而“模态分析实例（GUI方式）” 则给出了以从ANSYS GUI中选择菜单选项方式进行同一实例分析的步骤。（要想了解如何使用命令和GUI选项建模，请参阅<
?ANSYS动力学分析指南 作者: 安世亚太 第一章模态分析 §1.1模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性（固有频率和振型），即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性，如塑性和接触（间隙）单元，即使定义了也将被忽略。ANSYS提供了七种模态提取方法，它们分别是子空间法、分块Lanczos法、PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼。后面将详细介绍模态提取方法。 §1.2模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样，无论进行何种类型的分析，均可从用户图形界面（GUI）上选择等效于命令的菜单选项来建模和求解问题。 后面的“模态分析实例（命令流或批处理方式）”将给出进行该实例模态分析时要输入的命令（手工或以批处理方式运行ANSYS时）。而“模态分析实例（GUI方式）” 则给出了以从ANSYS GUI中选择菜单选项方式进行同一实例分析的步骤。（要想了解如何使用命令和GUI选项建模，请参阅<
/FILNAME,SolidBeam ,1 !定义工作文件名。/TITLE,SolidBeam Analysis !定义工作标题。/PREP7 ET,1,SOLID95 !定义材料属性。 MPTEMP,,,,,,,, !定义材料属性。MPTEMP,1,0 MPDATA,EX,1,,2.06e5 MPDATA,PRXY,1,,0.3 !建立几何模型 K,1,,,, K,2,450,,, K,3,450,-55,, K,4,,-100,, FLST,2,4,3 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 A,P51X VOFFST,1,45, , !网格划分。 FLST,5,4,4,ORDE,2 FITEM,5,9 FITEM,5,-12 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,3, , , , ,1 !* FLST,5,4,4,ORDE,4
FITEM,5,2 FITEM,5,4 FITEM,5,6 FITEM,5,8 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,6, , , , ,1 !* FLST,5,4,4,ORDE,4 FITEM,5,1 FITEM,5,3 FITEM,5,5 FITEM,5,7 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,30, , , , ,1 !* CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !加载。
ANSYS命令流总结(全) 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 -指定阶段状或者用跳板装载里面一个负荷步骤。