ANSYS房屋框架结构模型命令流APDL
ANSYS 房屋 框架结构模型 命令流 APDL
!********************提示********************************************!第一步:建立框架结构模型!第二步:框架结构施加重力荷载!第三步:框架结构施加活荷载!第四步:框架结构施加风荷载!第五步:框架结构荷载组合运算!第六步:框架结构模态分析!第七步:框架结构屈曲特征值分析!第八步:框架结构地震时程弹性分析!第九步:结束!采用地震波时间文件为:TIME.TXT!采用地震波文件分别为:AC_X.TXT/AC_Y.TXTFINISH/CLEAR/FILENAME,FRAME_SHEAR_WALL/TITLE, FRAME_SHEAR_WALL!************************************************************************** !************************第一步******************************************** !********************框架结构建模型***************************************/PREP7!采用单位为N/m/S 国际单位制!定义单元类型ET,1,BEAM4ET,2,SHELL63!定义实常数R,1,0.25,0.0052,0.0052,0.5,0.5R,2,0.2025,0.0034,0.0034,0.45,0.45R,3,0.06,0.0002,0.00045,0.3,0.2R,4,0.25,0.25,0.25,0.25R,5,0.1,0.1,0.1,0.1!定义材料参数MP,EX,1,3.0E10MP,NUXY,1,0.2MP,DENS,1,2500!关键点*DO,II,1,1K, 1+(II-1)*20, 0, 0, (II-1)*6K, 2+(II-1)*20, 6, 0, (II-1)*6K, 3+(II-1)*20, 12, 0, (II-1)*6K, 4+(II-1)*20, 17, 0, (II-1)*6K, 5+(II-1)*20, 0, 5, (II-1)*6K, 6+(II-1)*20, 6, 5, (II-1)*6K, 7+(II-1)*20, 12, 5, (II-1)*6K, 8+(II-1)*20, 17, 5, (II-1)*6K, 9+(II-1)*20, 0, 10, (II-1)*6K, 10+(II-1)*20, 6, 10, (II-1)*6K, 11+(II-1)*20, 12, 10, (II-1)*6K, 12+(II-1)*20, 17, 10, (II-1)*6K, 14+(II-1)*20, 6, 1.875, (II-1)*6K, 16+(II-1)*20, 6, 3.125, (II-1)*6K, 18+(II-1)*20, 6, 1.875, 3.9K, 20+(II-1)*20, 6, 3.125, 3.9*ENDDO*DO,II,1,10K, 1+II*20, 0, 0, (II-1)*3+6K, 2+II*20, 6, 0, (II-1)*3+6K, 3+II*20, 12, 0, (II-1)*3+6K, 4+II*20, 17, 0, (II-1)*3+6K, 5+II*20, 0, 5, (II-1)*3+6K, 6+II*20, 6, 5, (II-1)*3+6K, 7+II*20, 12, 5, (II-1)*3+6K, 8+II*20, 17, 5, (II-1)*3+6K, 9+II*20, 0, 10, (II-1)*3+6K, 10+II*20, 6, 10, (II-1)*3+6K, 11+II*20, 12, 10, (II-1)*3+6K, 12+II*20, 17, 10, (II-1)*3+6K, 13+II*20, 0, 1.875, (II-1)*3+6K, 14+II*20, 6, 1.875, (II-1)*3+6K, 15+II*20, 0, 3.125, (II-1)*3+6K, 16+II*20, 6, 3.125, (II-1)*3+6K, 17+II*20, 0, 1.875, (II-1)*3+6+1.5 K, 18+II*20, 6, 1.875, (II-1)*3+6+2 K, 19+II*20, 0, 3.125, (II-1)*3+6+1.5 K, 20+II*20, 6, 3.125, (II-1)*3+6+2 *ENDDO*DO,II,11,11K, 1+II*20, 0, 0, (II-1)*3+6K, 2+II*20, 6, 0, (II-1)*3+6K, 3+II*20, 12, 0, (II-1)*3+6K, 4+II*20, 17, 0, (II-1)*3+6K, 5+II*20, 0, 5, (II-1)*3+6K, 6+II*20, 6, 5, (II-1)*3+6K, 7+II*20, 12, 5, (II-1)*3+6K, 8+II*20, 17, 5, (II-1)*3+6K, 9+II*20, 0, 10, (II-1)*3+6K, 10+II*20, 6, 10, (II-1)*3+6K, 11+II*20, 12, 10, (II-1)*3+6K, 12+II*20, 17, 10, (II-1)*3+6K, 13+II*20, 0, 1.875, (II-1)*3+6K, 14+II*20, 6, 1.875, (II-1)*3+6K, 15+II*20, 0, 3.125, (II-1)*3+6K, 16+II*20, 6, 3.125, (II-1)*3+6K, 17+II*20, 0, 1.875, (II-1)*3+6+1.8 K, 18+II*20, 6, 1.875, (II-1)*3+6+2.5 K, 19+II*20, 0, 3.125, (II-1)*3+6+1.8 K, 20+II*20, 6, 3.125, (II-1)*3+6+2.5 *ENDDO*DO,II,12,12K, 1+II*20, 0, 0, 39.5K, 2+II*20, 6, 0, 39.5K, 3+II*20, 12, 0, 39.5K, 4+II*20, 17, 0, 39.5K, 5+II*20, 0, 5, 39.5K, 6+II*20, 6, 5, 39.5K, 7+II*20, 12, 5, 39.5K, 8+II*20, 17, 5, 39.5K, 9+II*20, 0, 10, 39.5K, 10+II*20, 6, 10, 39.5K, 11+II*20, 12, 10, 39.5K, 12+II*20, 17, 10, 39.5*ENDDO*DO,II,2,13L, (II-1)*20+1, (II-1)*20+2L, (II-1)*20+2, (II-1)*20+3L, (II-1)*20+3, (II-1)*20+4L, (II-1)*20+5, (II-1)*20+6L, (II-1)*20+6, (II-1)*20+7L, (II-1)*20+7, (II-1)*20+8L, (II-1)*20+9, (II-1)*20+10L, (II-1)*20+10, (II-1)*20+11L, (II-1)*20+11, (II-1)*20+12L, (II-1)*20+5, (II-1)*20+9L, (II-1)*20+6, (II-1)*20+10L, (II-1)*20+3, (II-1)*20+7L, (II-1)*20+7, (II-1)*20+11L, (II-1)*20+4, (II-1)*20+8L, (II-1)*20+8, (II-1)*20+12*ENDDO*DO,II,1,12L, (II-1)*20+1, II*20+1L, (II-1)*20+2, II*20+2L, (II-1)*20+3, II*20+3L, (II-1)*20+4, II*20+4L, (II-1)*20+5, II*20+5L, (II-1)*20+6, II*20+6L, (II-1)*20+7, II*20+7L, (II-1)*20+8, II*20+8L, (II-1)*20+9, II*20+9L, (II-1)*20+10, II*20+10L, (II-1)*20+11, II*20+11L, (II-1)*20+12, II*20+12*ENDDO*DO,II,1,1L, (II-1)*20+14, (II-1)*20+18L, (II-1)*20+16, (II-1)*20+20*ENDDO*DO,II,2,12L, (II-1)*20+13, (II-1)*20+17L, (II-1)*20+15, (II-1)*20+19L, (II-1)*20+14, (II-1)*20+18L, (II-1)*20+16, (II-1)*20+20*ENDDOA, 1, 21, 25, 5*DO,II,2,12A, (II-1)*20+1, II*20+1, II*20+5, (II-1)*20+5, (II-1)*20+15, (II-1)*20+19, (II-1)*20+17, (II-1)*20+13 *ENDDO*DO,II,1,12A, (II-1)*20+2, II*20+2, II*20+6, (II-1)*20+6, (II-1)*20+16, (II-1)*20+20, (II-1)*20+18, (II-1)*20+14 *ENDDO*DO,II,1,12A, II*20+1, II*20+2, II*20+6, II*20+5A, II*20+2, II*20+3, II*20+7, II*20+6A, II*20+3, II*20+4, II*20+8, II*20+7A, II*20+5, II*20+6, II*20+10,II*20+9A, II*20+6, II*20+7, II*20+11,II*20+10A, II*20+7, II*20+8, II*20+12,II*20+11*ENDDO!立柱网格划分LSEL,S,,,181,192 !底层立柱LATT,1,1,1LESIZE,ALL,,,10LMESH,ALLLSEL,S,,,192,324 !2~12层立柱LATT,1,2,1LESIZE,ALL,,,5LMESH,ALL!梁网格划分LSEL,S,,,1,180LATT,1,3,1,LESIZE,ALL,,,8LMESH,ALL!剪力墙网格划分LSEL,S,,,372 !外剪力墙两端、内剪力墙顶端LSEL,A,,,413LSEL,A,,,461LESIZE,ALL,,,8LSEL,s,,,418 !内剪力墙底端LSEL,A,,,420LESIZE,ALL,,,3LSEL,s,,,325 !内剪力墙底端内侧LSEL,A,,,326LESIZE,ALL,,,6LSEL,S,,,375 !剪力墙空洞顶部LSEL,A,,,379LSEL,A,,,383LSEL,A,,,387LSEL,A,,,391LSEL,A,,,395LSEL,A,,,399LSEL,A,,,403LSEL,A,,,407LSEL,A,,,411LSEL,A,,,415LSEL,A,,,419LSEL,A,,,423LSEL,A,,,427LSEL,A,,,431LSEL,A,,,435LSEL,A,,,439LSEL,A,,,443LSEL,A,,,447LSEL,A,,,451LSEL,A,,,455LSEL,A,,,459LSEL,A,,,463LESIZE,ALL,,,2!剪力墙网格划分ASEL,S,,,1,24AATT,1,4,2AMESH,ALL!楼板网格划分ASEL,S,,,25,96AATT,1,5,2AMESH,ALLNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点/eshape,1.0/VIEW,1,1,1,1/ANG,1,270,XM,0/REPlotFINISH!************************************************************************** !************************第二步******************************************** !********************框架结构施加重力荷载********************************* !框架结构施加重力荷载ANTYPE,STATICNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点ACEL,0,0,9.8SOLVEFINISH/POST1SET,FIRSTPLNSOL,U,Z,0,1FINISH!************************************************************************** !************************第三步******************************************** !*************框架结构施加楼面活荷载D=3KN/M^2***************************** !框架结构施加楼面活荷载D=3KN/M^2FINISH/SOLUANTYPE,STATICNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点ACEL,0,0,9.8SOLVEFINISH/POST1SET,FIRSTPLNSOL,U,Z,0,1FINISH!************************************************************************** !************************第四步******************************************** !***********框架结构施加风荷载(基本风压=0.25KN/M^2,体形系数为1.0********** !框架结构施加风荷载(先Y方向,后X方向)分两次分别施加ANTYPE,STATIC*DIM,LOAD_1,ARRAY,12LOAD_1(1)=3.78,2.16,2.39,2.57,2.72,2.84,2.95,3.17,3.20,3.29,3.39,3.51*DIM,LOAD_2,ARRAY,12LOAD_2(1)=7.56,4.32,4.78,5.14,5.44,5.68,5.90,6.34,6.40,6.58,6.76,7.02*DIM,LOAD_3,ARRAY,12LOAD_3(1)=6.93,3.96,4.38,4.71,4.99,5.21,5.41,5.81,5.87,6.03,6.20,6.44*DIM,LOAD_4,ARRAY,12LOAD_4(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93*DIM,LOAD_A,ARRAY,12LOAD_A(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93*DIM,LOAD_B,ARRAY,12LOAD_B(1)=6.30,3.60,3.98,4.28,4.54,4.74,4.92,5.28,5.34,5.48,5.64,5.86*DIM,LOAD_C,ARRAY,12LOAD_C(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93*DO,II,1,12FK,1+20*II,FY,LOAD_1(II)*ENDDO*DO,II,1,12FK,2+20*II,FY,LOAD_2(II)*ENDDO*DO,II,1,12FK,3+20*II,FY,LOAD_3(II)*ENDDO*DO,II,1,12FK,4+20*II,FY,LOAD_4(II)*ENDDOSOLVEFKDELE,ALL,ALL !第二次施加风荷载,删除第一次的。
ansys workbench中apdl的用法
ANSYS Workbench中的APDL(ANSYS Parametric Design Language)是一种参数化设计语言,用于在ANSYS软件中自动化建模和求解过程。
以下是APDL的一些用法:
1. 创建模型:使用APDL可以创建各种类型的模型,包括结构、流体动力学、电磁等。
在创建模型时,可以通过定义参数、约束条件和载荷等来自动化建模过程。
2. 优化设计:APDL可以用于优化设计,通过调整参数、约束条件和载荷等,获得最佳的设计方案。
3. 自动化求解:使用APDL可以自动化求解过程,包括网格划分、求解设置、结果后处理等。
4. 批处理操作:通过APDL,可以对一组模型进行批处理操作,例如批量分析、批量结果后处理等。
5. 自定义功能:使用APDL可以自定义功能,例如创建自定义的命令流、宏等,扩展ANSYS软件的功能。
在使用APDL时,需要注意以下几点:
1. 学习APDL需要一定的编程基础和数学知识。
2. 在使用APDL之前,需要了解ANSYS软件的基本操作和功能。
3. 在编写APDL脚本时,需要注意语法错误和逻辑错误,并进
行充分的测试和验证。
4. 在使用APDL进行复杂模型的分析时,需要注意计算资源和内存的分配,以确保计算过程的稳定性和效率。
ANSYS中的APDL命令
结合自身经验,谈ANSYS中的APDL命令(二)发表时间:2009-5-10 作者: 倪欣来源: e-works关键字: ANSYS APDL 命令流在ANSYS中,命令流是由一条条ANSYS的命令组成的一个命令组合,这些命令按照一定顺序排布,能够完成一定的ANSYS功能,本文是作者结合自身经验所总结的一些命令。
1.1 /prep7(进入前处理)定义几何图形:关键点、线、面、体(1).csys,kcnkcn , 0 迪卡尔坐标系1 柱坐标2 球4 工作平面5 柱坐标系(以Y轴为轴心)n 已定义的局部坐标系(2).numstr, label, value 设置以下项目编号的开始nodeelemkplineareavolu注意:vclear, aclear, lclear, kclear 将自动设置节点、单元开始号为最高号,这时如需要自定义起始号,重发numstr(3).K, npt, x,y,z, 定义关键点Npt:关键点号,如果赋0,则分配给最小号(4).Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imoveItime:拷贝份数Np1,Np2,Ninc:所选关键点Dx,Dy,Dz:偏移坐标Kinc:每份之间节点号增量noelem: “0” 如果附有节点及单元,则一起拷贝。
“1”不拷贝节点和单元imove:“0” 生成拷贝“1”移动原关键点至新位置,并保持号码,此时(itime,kinc,noelem)被忽略注意:MAT,REAL,TYPE 将一起拷贝,不是当前的MAT,REAL,TYPE(5).A, P1, P2, ……… P18 由关键点生成面(6).AL, L1,L2, ……,L10 由线生成面面的法向由L1按右手法则决定,如果L1为负号,则反向。
(线需在某一平面内坐标值固定的面内)(7).vsba, nv, na, sep0,keep1,keep2 用面分体(8).vdele, nv1, nv2, ninc, kswp 删除体kswp: 0 只删除体1 删除体及面、关键点(非公用)(9).vgen, itime, nv1, nv2, ninc, dx, dy, dz, kinc, noelem, imove 移动或拷贝体itime: 份数nv1, nv2, ninc:拷贝对象编号dx, dy, dz :位移增量kinc: 对应关键点号增量noelem,:0:同时拷贝节点及单元1:不拷贝节点及单元imove:0:拷贝体1:移动体(10).cm, cname, entity 定义组元,将几何元素分组形成组元cname: 由字母数字组成的组元名entity: 组元的类型(volu, area, line, kp, elem, node)(11).cmgrp, aname, cname1, ……,cname8 将组元分组形成组元集合aname: 组元集名称cname1……cname8: 已定义的组元或组元集名称1.2 定义几个所关心的节点,以备后处理时调用节点号。
ANSYS计算结果无难事,APDL经典命令让你的模型舞起来
ANSYS计算结果无难事,APDL经典命令让你的模型舞起来1、让你的ANSYS模型'舞'起来ANSYS计算结果的动画可采用ANTIME、ANMODE、ANCNTR、ANHARM等自动生成动画(AVI格式),使结果展示更加生动直观,相信使用ANSYS的都会制作。
然而,几何模型或有限元模型则无动画显示功能,有时为展示模型本身,会从多个角度截取图片。
那么,模型能否也可制作动画呢?答案是肯定的。
利用ANSYS的图形存储命令/SEG可以实现此功能,让你的模型动起来。
具体过程详见命令流中及其注释,动画上传总是失败,自己生成不要观看吧。
Finish$/clear$/prep7!简单的创建几何模型以减少篇幅blc4,0,0,4,2,5cyl4,2,4,1,,2,,4!关闭图例信息/plopts,info,off!以下开始制作模型动画!删除当前储存的图形/seg,dele/seg,multi,jhdh,1 !独立存储且不覆盖,文件名为jhdh/auto,1 !自动计算与图形区合适显示方式!正视/view,1,0,0,1$vplot!侧视/view,1,1$vplot!俯视/view,1,,1$vplot!D视图/view,1,1,1,1$vplot!循环36次,每次改变10度视角*do,i,1,36$/ang,1,10,ys,1$/replot$*enddo!关闭图形存储操作,保存为jhdh.avi文件/seg,off$/anfile,save,jhdh,avi其实比较简单,一旦进入模型动画制作过程,所有的xPLOT(x=KLAVNE)绘制的图形都将进入动画序列,按显示过程形成一部连续的动画。
2、用一个命令解决ANSYS数据列表分页早年初学ANSYS时,经常用到xLIST(如NLIST、ELIST、KLIST、LLIST、ALIST、VLIST等命令)和PRxSOL(如PRNSOL、PRESOL、PRRSOL、PRETAB、PRPATH)等列表命令,并希望将这些内容保存到TXT文件中,然后再导入EXCEL中处理。
ANSYS中的APDL命令总结
在ANSYS中,命令流是由一条条ANSYS的命令组成的一个命令组合,这些命令按照一定顺序排布,能够完成一定的ANSYS功能,这些功能一般来说通过菜单操作也能够实现(而那些命令流能够实现,菜单操作实现不了的单个命令比较少见)。
以下命令是结合我自身经验,和前辈们的一些经验而总结出来的,希望对大家有帮助。
(1).Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元(2).Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点type: S: 选择一组新节点(缺省)R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0”使用正负号“1”仅用绝对值(3).Esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元type: S: 选择一组单元(缺省)R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item:Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号(4). mp, lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens)ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号(缺省为当前材料号)c : 材料特性值,或材料之特性,温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数(5). 定义DP材料:首先要定义EX和泊松比:MP,EX,MAT,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MAT进入单元表并编辑添加单元表:TBDATA,1,CTBDATA,2,ψTBDATA,3,……如定义:EX=1E8,NUXY=0.3,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,0.3TB,DP,1TBDATA,1,27TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg(6). 根据需要耦合某些节点自由度cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotz ,allnode1-node17: 待耦合的节点号。
ANSYS APDL命令流详解-9网格划分技术
3.1 定义单元属性 —实常数
2. 变厚度壳实常数定义
命令:RTHICK,Par,ILOC,JLOC,KLOC,LLOC
Par---节点厚度的数组参数(以节点号引用),如mythick(19)
表示在节点19的壳体厚度。
ILOC---单元I节点的厚度在实常数组中的位置,缺省为1。
JLOC---单元J节点的厚度在实常数组中的位置,缺省为2。
D=0.1 PI=acos(-1) a0=pi*d*d/4 I0=pi*D**4/64 IX=pi*D**4/32 R,3,a0,i0,i0,d,d,0 !定义第3组实常数的AREA、IZZ、IYY、TKZ、TKY、THETA Rmore,0,ix,0,0,0,2.0!定义第3组实常数的其它实常数值
et,1,beam4
!定义BEAM4单元的参考号为1
et,3,beam189
!定义BEAM189单元的参考号为3
keyopt,1,2,1
!BEAM4单元考虑应力刚度时关闭一致切线刚度矩阵
keyopt,3,1,1
!考虑BEAM189的第7个自由度,即翘曲自由度
!当然这些参数也可在ET命令中一并定义,如上述四条命令
*DO,i,1,MXNODE !以节点号循环对厚度数组赋值
THICK(i) = 0.5 + 0.2*NX(i) + 0.02*NY(i)**2
与下列两条命令等效:
et,1,beam4,,1
et,3,beam189,1
3.1 定义单元属性 —单元类型
3. 自由度集 命令:DOF,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6,Lab7,Lab8,Lab9,Lab10
4.改变单元类型 命 令:ETCHG,Cnv
ANSYS+APDL+命令集合
ANSYS APDL 命令集合(2)2007-11-23 16:50:04/post1中的几个命令:set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入的数据lstep :荷载步数sbstep:子步数,缺省为最后一步time:时间点(如果弧长法则不用)nset: data set numberdscale, wn, dmult 显示变形比例wn: 窗口号(或all),缺省为1dmult, 0或auto : 自动将最大变形图画为构件长的5%pldisp, kund 显示变形的结构kund: 0 仅显示变形后的结构1 显示变形前和变形后的结构2 显示变形结构和未变形结构的边缘PRETAB,LAB1, LAB2, ……LAB9 沿线单元长度方向绘单元表数据LABn : 空:所有ETABLE命令指定的列名列名:任何ETABLE命令指定的列名PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据LABI:节点I的单元表列名LABJ:节点J的单元表列名FACT: 显示比例,缺省为1kund: 0 不显示未变形的结构1 变形和未变形重叠2 变形轮廓和未变形边缘etable, lab,item,comp将单元的某项结果制作成表格,以供pretable命令输出,lab: 字段名称,自己指定item: 结果的顶目名称,在每个单元的说明中有(在单元说明表中冒号左边的comp, 结果项目名称的分量,在单元说明表中冒号右边的比如将plane42单元的x应力分量制成表etable,sx,x,xLACAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2定义区域坐标系统,该命令执行后,ANSYS坐标系统自动更改为新建立的坐标系统,故可以定义许多区域坐标系统,以辅助有限元模型的建立。
KCN:该区域坐标系统的确定代号,大于10的任何一个号码都可以。
如何学习ANSYS命令流及APDL简解
ANSYS没有像其他软件一样的工具图标,但是可以把一些经常 用的功能通过自定义按钮的方式添加到Toolbar里面。
主要用到了*ABBR(Abbreviation的缩写)命令,其格式为 *ABBR,按钮上的标签名称,按钮调用的命令或宏的名称
命令*ABBR的使用格式如下: *ABBR,Abbr,string abbr:显示在工具条上的缩写名,最多可包含8个字符. string:Abbr所代表的ANSYS命令,系统UIDL函数名或宏名的缩写字符串.如果*ABBR命令 的值域string是一个宏名,那么该宏必须放在指定的宏搜索路径下.
9
Dare Design 2014
上篇
简介及准备
如何学习ALeabharlann SYS命令和APDL >ANSYS命令按照功能可分为三个大类:前 处理命令、后处理命令、和结果查看命令, 每个大类有自己对应的处理器,/PREP7处 理器,后处理/POST1、/POST26等。 ANSYS有超过1000条命令,很难有人把 这些完完全全记住,我建议先学习APDL语 法及规则,记住常用的关键词,配合这些关 键词套用需要的命令。然后了解常用的 ANSYS命令。对于ANSYS常用命令的学习 ,网上资料很多,更详细的用法可以在 ANSYS主菜单-help-help topic中查找。
10
Dare Design 2014
上篇
简介及准备
如何学习ANSYS命令和APDL >
图1 APDL在ANSYS14.5 帮助中的位置
11
Dare Design 2014
工欲善其事必先利其
器
(编辑器选择)
多数命令流可以直接拷贝到ANSYS命令提示符栏中运行,但部 分命令不支持这种方式,而且执行速度比较慢。更常用的方法是 将命令流在文本文档中整理好,然后通过/input命令或Utility Menu>File>Read Input From运行。 最常见的文本编辑器就是WINDOWS自带的记事本了,具有 系统自带、功能简单,使用方便的特点,能够满足基本的需求。
ANSYS APDL命令流详解-6几何建模的其它常用命令
2.4.3 选择与组件
图素选择是建模和结果处理的重要手段,而在选择中 使 用组件将更加方便。本章前述也介绍了一些选择命令,选 择命 令如表所示。
图 素
命令
NSEL
NSLE
节 点
NSLK
NSLL
NSLA
NSLV
关 键
KSEL
点 KSLN
KSLL
ASEL 面 ASLL
ASLV
ALLSEL
功能说明
选择节点的基本命令 选择与所选单元相关的节点 选择与所选关键点相关节点 选择与所选线相关的节点 选择与所选面相关的节点 选择与所选体相关的节点 选择关键点的基本命令 选择与所选节点相关关键点 选择与所选线相关的关键点 选择面的基本命令 选择与所选线相关的面 选择与所选体相关的面 选择所有图素
图 素
命令
ESEL
单 ESLN 元 ESLL
ESLA
ESLV
LSEL 线 LSLA
LSLK
体 VSEL
VSLA
组 件
CMSEL
DOFSEL
功能说明 选择单元的基本命令 选择与所选节点相关单元 选择与所选线相关的单元 选择与所选面相关的单元 选择与所选体相关的单元 选择线的基本命令 选择与所选面相关的线 选择与所选关键点相关的 线选择体的基本命令 选择与所选面相关的体 选择组件
如为SOLU则保存模型数据和后处理数据。 该命令将当前的数据信息保存到文件中,而且可产生备份(
扩展名为DBB),以防系统崩溃或操作失误时可恢复上次保存
的数据库。事实上,当采用命令流时不必担心这个问题,也不
必经常保存数据库,除非求解或后处理花费时间较长;或者在
编写命令流时需要通过GUI方式得到验证,可不断保存和恢复 数据库,以加快编写速度和质量。
ANSYS APDL命令流详解-13通用与时间历程后处理技术
5.1 通用后处理—图形显示结果
命令
功能
PLDISP 显示结构变形图
PLNSOL 显示节点结果图
PLESOL 显示单元结果图
PLVECT
PLCRACK
以矢量方式显示 结 果图
显示裂缝或压碎图
备注
结果内容很多,根据需要选用Item及 Comp 结果内容很多,根据需要选用Item及 Comp 如U、S等
Lab2,Lab3---用户定义的分矢量标识符,如Item为预定标识符,必须为空。 LabP---合成矢量标识符,缺省为Item。 Mode---显示方式控制。
如为空,则采用/DEVICE中的KEY参数指定的方式; 如MODE=RAST则采用光栅模式; 如MODE=VECT则采用矢量模式显示。 Loc---显示单元场结果的矢量位置。
5.1 通用后处理—结果输出控制选项
用于图形显示和列表显示,如导出结果的方式和显示比例等设置。
命令
功能
备注
AVPRIN
定义矢量和主轴的计算方法
用于计算主应力或主应变等时
AVRES
定义结果数据平均处理
仅适用于PowerGraphics模式
/EFACET
设置单元每边的分段数目
见2.4.1中
/DSCALE
第5章 通用与时间历程后处理技术
ANSYS两个后处理器: ★通用后处理器POST1:
查看整个模型在各个时间点上的结果. ★时间历程后处理器POST26:
查看整个模型上的某一点结果随时间的变化曲线。 后处理可在求解完后直接进入,也可在重新进入 ANSYS后读入文件进入后处理。
5.1 通用后处理—读入结果文件
ANSYS APDL命令汇总
ANSYS APDL命令汇总AA,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点(P1~P9)来定义面(Area),最少使用三个点才能围成面,同时产生转围绕些面的线。
点要依次序输入,输入的顺序会决定面的法线方向。
如果超过四个点,则这些点必须在同一个平面上。
Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPsABBR*ABBR,Abbr,String--定义一个缩略语.Abbr:用来表示字符串"String"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.ABBRESABBRES,Lab,Fname,Ext-从一个编码文件中读出缩略语.Lab:指定读操作的标题,NEW:用这些读出的缩略语重新取代当前的缩略语(默认)CHANGE:将读出的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语.Ext:如果"Fname"是空的,则缺省的扩展命是"ABBR".ABBSAVABBSAV,Lab,Fname,Ext-将当前的缩略语写入一个文本文件里Lab:指定写操作的标题,若为ALL,表示将所有的缩略语都写入文件(默认ADDadd, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic:变量号name: 变量的名称ADELEAdele,na1,na2,ninc,kswp !kswp=0时只删除掉面积本身,=1时低单元点一并删除。
ADRAGAdrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 !面积的建立,沿某组线段路径,拉伸而成AFILLTAfillt,na1,na2,rad !建立圆角面积,在两相交平面间产生曲面,rad为半径。
ANSYS 分析的命令流(APDL语言)
ANSYS 分析的命令流(APDL语言)前面分析的过程都是基于ANSYS用户图形界面(GUI命令流)形式,GUI 形式非常的直观明了,分析过程中通常会采用此种的分析方式。
但是这种分析方式会有一个缺点,就是当我们操作失误时候,没有后退的功能,因此我们不得不重新操作,这样就会给分析带来很多的麻烦。
但如果以命令流的形式进行计算的话就能很容易的减少这些麻烦。
下面就是本次点机前处理和分析计算过程的命令流/UNITS,SI !采用国际单位制/TITLE,2D DIANJI Static Analysis !定义分析名称KEYW,MAGNOD,1*AFUN,DEG !指定角度单位为度/PREP7 !进入前处理器ET,1,PLANE53 !设定单元类型为plane53EMUNIT,MKS !电磁单位LOCAL,12,0,0,0,0,90 !定义局部坐标系12--17LOCAL,13,0,0,0,0,30LOCAL,14,0,0,0,0,330LOCAL,15,0,0,0,0,270LOCAL,16,0,0,0,0,210LOCAL,17,0,0,0,0,150MP,MURX,1,1.0 ! 定义第一种材料的相对磁导率MP,MURX,2,1.0 ! 定义第二种材料的相对磁导率TB,BH,3TBPT,,35.03,0.1 ! 第三种材料的B-H磁化特性TBPT,,46.97,0.2TBPT,,57.32,0.3TBPT,,66.08,0.4TBPT,,74.04,0.5TBPT,,82.01,0.6TBPT,,91.56,0.7TBPT,,103.5,0.8TBPT,,117.83,0.9TBPT,,136.15,1TBPT,,160.83,1.1TBPT,,203.03,1.2TBPT,,286.62,1.3TBPT,,461.78,1.4TBPT,,955.41,1.5TBPT,,2547.8,1.6!TBPLOT,BH,3 !绘制B-H曲线TB,BH,4TBPT,,130,0.1 ! 定义第四种材料的B-H磁化特性TBPT,,170,0.2TBPT,,197,0.3TBPT,,218,0.4TBPT,,250,0.5TBPT,,290,0.6TBPT,,338,0.7TBPT,,400,0.8TBPT,,472,0.9TBPT,,570,1TBPT,,682,1.1TBPT,,810,1.2TBPT,,970,1.3TBPT,,1600,1.4TBPT,,2520,1.5TBPT,,3520,1.6!TBPLOT,BH,4 !绘制B-H曲线MP,MURX,5,1.17 ! 定义第五种材料的相对磁导率MP,MGXX,5,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,6,1.17 ! 定义第六种材料的相对磁导率MP,MGXX,6,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,7,1.17 ! 定义第七种材料的相对磁导率MP,MGXX,7,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,8,1.17 ! 定义第八种材料的相对磁导率MP,MGXX,8,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,9,1.17 ! 定义第九种材料的相对磁导率MP,MGXX,9,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,10,1.17 ! 定义第十种材料的相对磁导率MP,MGXX,10,4.4E+005 ! 定义矫顽力矢量为X方向分量RECTNG,-0.00275,0.00275,0.0485,0.062 ! 创建一个矩形RECTNG,-0.00025,0.00025,0.048,0.0485 ! 再创建一个矩形AADD,ALL ! 将两个矩形融合在一起,形成槽结构CSYS,1 !转动工作平面到柱坐标下AGEN,36,3,38,1,,10,0 ! 复制槽/PNUM,AREA,ALL ! 给所有槽编号PCIRC,0.048,0.074 ! 创建定子圆环AOVLAP,ALL !对定子圆环和槽进行交迭布尔操作使得在连接线上共节点PCIRC,0.037,0.047,78,102 ! 创建永磁体极CSYS,1 !转到柱坐标下AGEN,6,3,8,1,,60,0 ! 复制6个磁极PCIRC,0.023,0.047 ! 创建转子圆环AOVLAP,3,10,45,51,57,63,69 ! 对转子圆环与磁极进行交迭操作PCIRC,0.047,0.0475 !创建两层气隙PCIRC,0.0475,0.048AOVLAP,ALLAGLUE,ALL ! 将所有交接面黏在一起NUMCMP,AREA !压缩面编号ASEL,S,AREA,,1,35 !选中槽ASEL,A,AREA,,42AATT,1 !材料定义为1号材料ASEL,S,AREA,,45ASEL,A,AREA,,43 !选中气隙AATT,2 !材料定义为2号材料ASEL,S,AREA,,46 !选中定子铁芯AATT,3 !材料定义为3号材料ASEL,S,AREA,,44 !选中转子磁轭AATT,4 !材料定义为4号材料ASEL,S,AREA,,40 !选中第1-6个永磁体磁极AATT,5,,,12 !材料定义为5-10号材料ASEL,S,AREA,,39AATT,6,,,13ASEL,S,AREA,,38AATT,7,,,14ASEL,S,AREA,,41AATT,8,,,15ASEL,S,AREA,,37AATT,9,,,16ASEL,S,AREA,,36AATT,10,,,17/PNUM,MAT,1 !打开材料编号APLOT !重新显示ALLSEL,ALL !选择所有实体SMRTSIZE,1 !指定智能网格划分等级L AMESH,ALL !智能网格划分!ESEL,S,MAT,,1!CM,ARM,ELEM!FMAGBC,'ARM'!SA VE,DIANJI_2D_MESH.db !保存!SA VEFINISH !完成网格划分/SOLU !进入求解器ANTYPE,STATIC !选择静态磁场分析NROPT,AUTO !采用牛顿---拉夫森迭代方法!ASEL,S,AREA,,1,35 !选中槽!ASEL,A,AREA,,42!BFE,ALL,JS,1,,,15000000 !施加电流密度LOCAL,11,1 ! 定义柱坐标系FLST,2,8,4,ORDE,4FITEM,2,289FITEM,2,-292FITEM,2,561FITEM,2,-564DL,P51X, ,AZ,0, !加边界条件ALLSEL,ALLMAGSOLV!SA VE,EMAGE_2D_MESH.db !保存!SA VEFINISH/POST1PLF2D,27,0,10,1 !显示磁力线图FMAGSUM,'ARM' !对电磁力求和plvect,B,,,,VECT,ELEM,ON,0 !显示磁感应强度矢量plvect,H,,,,VECT,ELEM,ON,0 !显示磁场强度矢量/GRAPHICS,POWERA VRES,2PLNSOL,B,SUM !显示磁通密度等值云图flishESEL,S,MAT,,2CM,ARM,ELEM!生成一个组件FMAGSUM,'ARM'!以下命令流是保存图片时使用,使其背景为白,结合PlotCtrlsHard ﹥Copy命令可设置格式/GRAPHICS,POWER/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 !背景反白显示/SHOW,JPEG,,0JPEG,QUAL,75,JPEG , ORIENT, HORIZJPEG , Color, 2JPEG , TMOD, 1/GFILE,600,EPLOT/SHOW,CLOSE/DEVICE,VECTOR,0。
ANSYS APDL命令流详解-9网格划分技术
面映射网格划分/体映射网格划分 扫掠生成体网格/单元有效性检查/网格修改
3.4 网格划分实例
基本模型的网格划分/复杂面模型的网格划分 复杂体模型的网格划分
创建几何模型后,必须生成有限元模型才能 分析计算,生成有限元模型的方法就是对几何 模型进行网格划分,网格划分主要过程包括三 个步骤: ⑴ 定义单元属性
3.1 定义单元属性 —材料属性
每一组材料属性有一个材料参考号,用于识别各个材料特
性组。一个模型中可有多种材料特性组。 1. 定义线性材料属性 命令:MP,Lab,MAT,C0,C1,C2,C3,C4 Lab---材料性能标识,其值可取:
EX:弹性模量(也可为EY、EZ)。 ALPX:线膨胀系数(也可为ALPY、ALPZ)。 PRXY:主泊松比(也可为PRYZ、PRXZ)。 NUXY:次泊松比(也可为NUYZ、NUXZ)。 GXY:剪切模量(也可为GYZ、GXZ)。 DAMP:用于阻尼的K矩阵乘子,即阻尼比。 DMPR:均质材料阻尼系数。 MU:摩擦系数 。 DENS:质量密度。
与下列两条命令等效:
et,1,beam4,,1
et,3,beam189,1
3.1 定义单元属性 —单元类型
3. 自由度集 命令:DOF,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6,Lab7,Lab8,Lab9,Lab10
4.改变单元类型 命 令:ETCHG,Cnv
5. 单元类型的删除与列表 删除命令:ETDELE,ITYP1,ITYP2,INC 列表命令:ETLIST,ITYP1,ITYP2,INC
单元属性包括单元类型、实常数、材料特性、单元坐标系和截 面号等。
⑵ 定义网格控制选项
★对几何图素边界划分网格的大小和数目进行设置; ★没有固定的网格密度可供参考; ★可通过评估结果来评价网格的密度是否合理。
APDL命令流:将ansys分析结果输出为tecplot格式
APDL命令流:将ansys分析结果输出为tecplot格式/post1file,Jobname,rst !指明从哪⼀个结果⽂件中读取数据*get,NodeNum,NODE,0,COUNT !得到模型的所有节点数*get,nd,NODE,0,NUM,MIN !得到模型的最⼩节点编号*dim,nodes,array,NodeNum !定义⼀个存储节点的数组*dim,xyz,array,NodeNum,3 !定义⼀个存储节点坐标的数组*dim,NodeTemp,array,NodeNum !定义⼀个存储节点的数组set,4,1 !读⼊第四载荷步、第⼀个⼦步的结果*do,i,1,NodeNum,1nodes(i)=nd !将节点编号存储在nodes数组中xyz(i,1)=NX(nd) !在xyz(i,1)中存储节点的X坐标xyz(i,2)=NY(nd) !在xyz(i,1)中存储节点的Y坐标xyz(i,3)=NZ(nd) !在xyz(i,1)中存储节点的Z坐标*Get,NodeTemp(i),NODE,nodes(i),TEMPnd=NDNEXT(nd) !取得下⼀个节点编号*enddo*get,ElementNum,ELEM,0,COUNT !得到模型的所有单元数*get,el,ELEM,0,NUM,MIN !得到模型的最⼩单元编号*dim,elems,array,ElementNum !定义⼀个存储单元的数组*dim,ndlst,array,ElementNum,8 !定义⼀个存储节点坐标的数组*do,i,1,ElementNum,1elems(i)=el !将单元编号存储在elems数组中*do,j,1,8,1ndlst(i,j)=NELEM(el,j) !将单元对应的节点编号存储在ndlst数组中*enddoel=ELNEXT(el) !取得下⼀个单元编号*enddo*cfopen,mesh,dat !输出的⽂件名为:mesh.dat!写TecPlot⽂件的⽂件头*vwrite(’TITLE=Ansys Temperature Analysis’)!写变量名*vwrite(’VARIABLES=”X”,”Y”,”Z”,”Temp”‘)!写节点数和单元数*vwrite,NodeNum,ElementNum(’ZONE n=’f6.1,’ e=’f6.1,’ f=fepoint’,' et=brick’)!写节点的坐标和温度值*vwrite,xyz(1,1),xyz(1,2),xyz(1,3),NodeTemp(1)(4f12.6)!写⼀个空⾏*vwrite(’ ‘)!写各单元对应的节点号*vwrite,ndlst(1,1),ndlst(1,2),ndlst(1,3),ndlst(1,4),ndlst(1,5),ndlst(1,6),ndlst(1,7),ndlst(1,8)(8f6.0)!关闭⽂件*CFCLOSEfinish。
ANSYS软件APDL命令流建模的体会
ANSYS软件APDL命令流建模的体会首先申明,本人学习AN SYS基本上是靠自己一点一点琢磨出来的,由于本人喜欢用APD L命令流,故总结出来的几点经验也就比较适合用AP DL命令的朋友。
1、多看help,ANSYS的help为我们提供了很强大的功能,我最喜欢的是其中对各个命令有关参数的说明和解释部分,不管是建模、加载、后处理等,都可以通过a pdl命令来实现。
只要你知道命令,如“a att”,在hel p搜索栏输入“a at t”,回车,弹出aa tt的有关页码,一般其中有一个只有“aa tt”的一项,确认,即可看到你要查询的aatt命令的有关参数意义,本人常用的命令有:e t---定义单元类型mp---定义材料属性k----建关键点,l----建线条a---由关键点建立面al---由线建立面v----由关键点建立体v l---由线建立体va--由面建立体lsel---在很多很多线中选择你需要的目标线,数量可以无限多……a sel---在很多很多面中选择你需要的目标面,数量也可以无限多……v sel---在很多很多体中选择你需要的目标体,数量也可以无限多……l att----给选中的线按材料编号赋属性(前提是首先已定义好材料)a att---给选中的面按材料编号赋属性vatt-----给选中的体按材料编号赋属性ac el---按坐标轴赋体积力,lmesh,ames h,vme sh---对线、面、体进行剖分d---在节点上加约束边界dl---在线上加载约束边界da----在面上加载约束边界2、以上只是列出了常见的几个命令,但是a nsys提供的命令是很多的,我们不可能都记得,计算记得,也不知道其有关参数是如何定义的,那不要紧,我们可以与界面操作结合起来学习。
ansysapdl命令流输入方法
ansysapdl命令流输入方法ANSYS APDL命令流输入方法概述ANSYS APDL(ANSYS Parametric Design Language)是ANSYS软件中的一种命令流输入方法,它允许用户通过输入一系列的命令来定义模型、设置分析条件、运行分析和获取结果。
本文将介绍ANSYS APDL的基本输入方法和常用命令,帮助读者了解如何使用APDL进行工程分析。
ANSYS APDL使用命令流输入方法,即通过输入一系列的命令来完成模型定义、分析设置和结果获取。
用户可以将这些命令保存在一个文本文件中,然后通过ANSYS界面中的命令输入窗口或者批处理脚本来执行这些命令。
下面是一个简单的命令流输入示例:! 定义节点n,1,0,0n,2,1,0n,3,1,1n,4,0,1! 定义单元et,1,PLANE42r,1,0.1type,1mat,1real,1secnum,1secdata,1emodif,1,1e7! 定义边界条件d,1,1,ux,0d,2,2,ux,0d,3,2,uy,0d,4,1,uy,0! 施加荷载f,2,fx,1000! 求解solve! 输出结果finish在这个示例中,首先通过n命令定义了4个节点,然后通过et、r、type、mat、real等命令定义了一个平面应力单元,并给定了单元的材料属性和几何参数。
接着使用d命令定义了边界条件,限制了节点1和节点2的x方向位移为0,节点2和节点3的y方向位移为0。
然后使用f命令施加了一个1000N的水平力。
最后使用solve 命令求解模型,并使用finish命令结束分析。
常用命令ANSYS APDL提供了丰富的命令用于定义模型、设置分析条件和获取结果。
下面是一些常用的命令:1. n:定义节点,可以指定节点的坐标;2. et:定义单元类型,可以选择不同类型的单元,如平面应力单元、平面应变单元等;3. r:定义单元属性,如单元的厚度、截面积等;4. type:定义单元类型属性,如材料特性、几何特性等;5. mat:定义材料属性,如杨氏模量、泊松比等;6. real:定义实体特性,如单元类型、材料类型等;7. secnum:定义截面类型,如圆形、矩形等;8. secdata:定义截面参数,如半径、宽度、高度等;9. emodif:修改单元材料属性;10. d:定义边界条件,如位移约束、固支约束等;11. f:定义荷载,如力、压力等;12. solve:求解模型;13. finish:结束分析。
ANSYS APDL命令流详解-11荷载及其施加
固定约束、支座 沉 降等
力、力矩等
施加在模型面上的分布力
压力、线荷载
施加体积荷载或场荷载
温度
施加物理惯性引起的荷载
从一种分析得到的结果,作 为另一种分析的荷载
重力加速度、角 速 度角加速度等
热分析的温度等
4.1 荷载及其施加---荷载
★荷载即可施加在几何模型(关键点、硬点、线、面、 体)上,也可施加在有限元模型(节点、单元)上, 或者二者混合使用。
2. 关键点自由度约束及相关命令
命令:DK,KPOI,Lab,VALUE,VALUE2,KEXPND,Lab2,Lab3,Lab4,Lab5,Lab6
KPOI---关键点编号,也可取ALL或元件名。
KEXPND---扩展控制参数。如为0则仅施加约束到关键点上的节
点;如为1则扩展到关键点之间(两关键点所连线)
★施加在几何模型上的荷载独立于有限元网格,不必 为修改网格而重新加载;
★施加在有限元模型上且要修改网格,则必须先删除 荷载再修改网格,然后重新施加荷载。
★不管施加到何种模型上,在求解时荷载全部转换 (自动或人工)到有限元模型上。
4.1 荷载及其施加---施加自由度约束
在结构分析中自由度共有7个,自由度的方向均依从节点坐 标系。约束可施加在节点、关键点、线和面上。
NODE---拟施加约束的节点号,其值可取ALL、组件名。 Lab---自由度标识符,如UX、ROTZ等。如为ALL,则为所有适宜的自由度。 VALUE---自由度约束位移值或表式边界条件的表格名称。 VALUE2---约束位移值的第二个数,如为复数输入时,VALUE为实部,而VALUE2为虚部。 NEND,NINC---节点编号范围和编号增量,缺省时NEND=NODE,NINC=1。 Lab2,Lab3,Lab4,Lab5,Lab6---其它自由度标识符,VALUE对这些自由度也有效。 各自由度的方向用节点坐标系确定,转角约束位移用弧度输入
ANSYS APDL命令流详解-16结构弹性稳定分析
7 . 1 特征值屈曲分析的步骤--获得静力解
注意几个问题:
⑴ 必须激活预应力效应。
命令PSTRES设为ON便可考虑预应力效应。
⑵ 由屈曲分析所得到的特征值是屈曲荷载系数,而屈 曲荷载等于该系数乘以所施加的荷载。若施加单位荷 载,则该屈曲荷载系数就是屈曲荷载;若施加了多种 不同类型的荷载,则将所有荷载按该系数缩放即为屈 曲荷载。
正常求解:屈曲荷载=屈曲荷载系数×(恒载+活载) 实际要求:屈曲荷载=1.0×(恒载+K×活 载 ) 其实现方法是通过调整所施加的活载大小(例如放大K 倍),然后进行屈曲分析,如果所求得的屈曲荷载系数不等 于 1.0,则继续修改K值重新分析,直到屈曲荷载系数为1.0为 止。 K的初值通常可采用第一次的屈曲荷载系数,然后调整3
~ 4次即可达到要求。
⑸ 非零约束。如同静力分析一样,可以施加非零约束。同样以 屈曲荷载系数对非零约束进行缩放得到屈曲荷载。
⑹ 静力求解完成后,退出求解层。
7 . 1 特征值屈曲分析的步骤--获得特征值屈曲解
该过程需要静力分析中得到的.EMAT和.ESAV文件,且数据库 中包含有模型数据,以备需要时恢复。如下步骤: ⑴ 进入求解层 命
第7章 结构弹性稳定分析
1. 特征值屈曲分析的步骤 2. 构件的特征值屈曲分析 3. 结构的特征值屈曲分析
第7章 结构弹性稳定分析
结构失稳或结构屈曲: 当结构所受载荷达到某一值时,若增加一微小的
增量,则结构的平衡位形将发生很大的改变,这种现 象叫做结构失稳或结构屈曲。 结构稳定问题一般分为两类: ★ 第一类失稳:又称平衡分岔失稳、分枝点失稳、特 征值屈曲分析。结构失稳时相应的荷载可称为屈曲荷 载、临界荷载、压屈荷载或平衡分枝荷载。 ★ 第二类失稳:结构失稳时,平衡状态不发生质变, 也称极值点失稳。结构失稳时相应的荷载称为极限荷 载或压溃荷载。 ●跳跃失稳:当荷载达到某值时,结构平衡状态发生 一明显的跳跃,突然过渡到非邻近的另一具有较大位 移的平衡状态。可归入第二类失稳。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
!********************提示********************************************!第一步:建立框架结构模型!第二步:框架结构施加重力荷载!第三步:框架结构施加活荷载!第四步:框架结构施加风荷载!第五步:框架结构荷载组合运算!第六步:框架结构模态分析!第七步:框架结构屈曲特征值分析!第八步:框架结构地震时程弹性分析!第九步:结束!采用地震波时间文件为:TIME.TXT!采用地震波文件分别为:AC_X.TXT/AC_Y.TXTFINISH/CLEAR/FILENAME,FRAME_SHEAR_WALL/TITLE, FRAME_SHEAR_WALL!************************************************************************** !************************第一步******************************************** !********************框架结构建模型***************************************/PREP7!采用单位为N/m/S 国际单位制!定义单元类型ET,1,BEAM4ET,2,SHELL63!定义实常数R,1,0.25,0.0052,0.0052,0.5,0.5R,2,0.2025,0.0034,0.0034,0.45,0.45R,3,0.06,0.0002,0.00045,0.3,0.2R,4,0.25,0.25,0.25,0.25R,5,0.1,0.1,0.1,0.1!定义材料参数MP,EX,1,3.0E10MP,NUXY,1,0.2MP,DENS,1,2500!关键点*DO,II,1,1K, 1+(II-1)*20, 0, 0, (II-1)*6K, 2+(II-1)*20, 6, 0, (II-1)*6K, 3+(II-1)*20, 12, 0, (II-1)*6K, 4+(II-1)*20, 17, 0, (II-1)*6K, 5+(II-1)*20, 0, 5, (II-1)*6K, 6+(II-1)*20, 6, 5, (II-1)*6K, 7+(II-1)*20, 12, 5, (II-1)*6K, 8+(II-1)*20, 17, 5, (II-1)*6K, 9+(II-1)*20, 0, 10, (II-1)*6K, 10+(II-1)*20, 6, 10, (II-1)*6K, 11+(II-1)*20, 12, 10, (II-1)*6K, 12+(II-1)*20, 17, 10, (II-1)*6K, 14+(II-1)*20, 6, 1.875, (II-1)*6K, 16+(II-1)*20, 6, 3.125, (II-1)*6K, 18+(II-1)*20, 6, 1.875, 3.9K, 20+(II-1)*20, 6, 3.125, 3.9*ENDDO*DO,II,1,10K, 1+II*20, 0, 0, (II-1)*3+6K, 2+II*20, 6, 0, (II-1)*3+6K, 3+II*20, 12, 0, (II-1)*3+6K, 4+II*20, 17, 0, (II-1)*3+6K, 5+II*20, 0, 5, (II-1)*3+6K, 6+II*20, 6, 5, (II-1)*3+6K, 7+II*20, 12, 5, (II-1)*3+6K, 8+II*20, 17, 5, (II-1)*3+6K, 9+II*20, 0, 10, (II-1)*3+6K, 10+II*20, 6, 10, (II-1)*3+6K, 11+II*20, 12, 10, (II-1)*3+6K, 12+II*20, 17, 10, (II-1)*3+6K, 13+II*20, 0, 1.875, (II-1)*3+6K, 14+II*20, 6, 1.875, (II-1)*3+6K, 15+II*20, 0, 3.125, (II-1)*3+6K, 16+II*20, 6, 3.125, (II-1)*3+6K, 17+II*20, 0, 1.875, (II-1)*3+6+1.5 K, 18+II*20, 6, 1.875, (II-1)*3+6+2 K, 19+II*20, 0, 3.125, (II-1)*3+6+1.5 K, 20+II*20, 6, 3.125, (II-1)*3+6+2 *ENDDO*DO,II,11,11K, 1+II*20, 0, 0, (II-1)*3+6K, 2+II*20, 6, 0, (II-1)*3+6K, 3+II*20, 12, 0, (II-1)*3+6K, 4+II*20, 17, 0, (II-1)*3+6K, 5+II*20, 0, 5, (II-1)*3+6K, 6+II*20, 6, 5, (II-1)*3+6K, 7+II*20, 12, 5, (II-1)*3+6K, 8+II*20, 17, 5, (II-1)*3+6K, 9+II*20, 0, 10, (II-1)*3+6K, 10+II*20, 6, 10, (II-1)*3+6K, 11+II*20, 12, 10, (II-1)*3+6K, 12+II*20, 17, 10, (II-1)*3+6K, 13+II*20, 0, 1.875, (II-1)*3+6K, 14+II*20, 6, 1.875, (II-1)*3+6K, 15+II*20, 0, 3.125, (II-1)*3+6K, 16+II*20, 6, 3.125, (II-1)*3+6K, 17+II*20, 0, 1.875, (II-1)*3+6+1.8 K, 18+II*20, 6, 1.875, (II-1)*3+6+2.5 K, 19+II*20, 0, 3.125, (II-1)*3+6+1.8 K, 20+II*20, 6, 3.125, (II-1)*3+6+2.5 *ENDDO*DO,II,12,12K, 1+II*20, 0, 0, 39.5K, 2+II*20, 6, 0, 39.5K, 3+II*20, 12, 0, 39.5K, 4+II*20, 17, 0, 39.5K, 5+II*20, 0, 5, 39.5K, 6+II*20, 6, 5, 39.5K, 7+II*20, 12, 5, 39.5K, 8+II*20, 17, 5, 39.5K, 9+II*20, 0, 10, 39.5K, 10+II*20, 6, 10, 39.5K, 11+II*20, 12, 10, 39.5K, 12+II*20, 17, 10, 39.5*ENDDO*DO,II,2,13L, (II-1)*20+1, (II-1)*20+2L, (II-1)*20+2, (II-1)*20+3L, (II-1)*20+3, (II-1)*20+4L, (II-1)*20+5, (II-1)*20+6L, (II-1)*20+6, (II-1)*20+7L, (II-1)*20+7, (II-1)*20+8L, (II-1)*20+9, (II-1)*20+10L, (II-1)*20+10, (II-1)*20+11L, (II-1)*20+11, (II-1)*20+12L, (II-1)*20+5, (II-1)*20+9L, (II-1)*20+6, (II-1)*20+10L, (II-1)*20+3, (II-1)*20+7L, (II-1)*20+7, (II-1)*20+11L, (II-1)*20+4, (II-1)*20+8L, (II-1)*20+8, (II-1)*20+12*ENDDO*DO,II,1,12L, (II-1)*20+1, II*20+1L, (II-1)*20+2, II*20+2L, (II-1)*20+3, II*20+3L, (II-1)*20+4, II*20+4L, (II-1)*20+5, II*20+5L, (II-1)*20+6, II*20+6L, (II-1)*20+7, II*20+7L, (II-1)*20+8, II*20+8L, (II-1)*20+9, II*20+9L, (II-1)*20+10, II*20+10L, (II-1)*20+11, II*20+11L, (II-1)*20+12, II*20+12*ENDDO*DO,II,1,1L, (II-1)*20+14, (II-1)*20+18L, (II-1)*20+16, (II-1)*20+20*ENDDO*DO,II,2,12L, (II-1)*20+13, (II-1)*20+17L, (II-1)*20+15, (II-1)*20+19L, (II-1)*20+14, (II-1)*20+18L, (II-1)*20+16, (II-1)*20+20*ENDDOA, 1, 21, 25, 5*DO,II,2,12A, (II-1)*20+1, II*20+1, II*20+5, (II-1)*20+5, (II-1)*20+15, (II-1)*20+19, (II-1)*20+17, (II-1)*20+13 *ENDDO*DO,II,1,12A, (II-1)*20+2, II*20+2, II*20+6, (II-1)*20+6, (II-1)*20+16, (II-1)*20+20, (II-1)*20+18, (II-1)*20+14 *ENDDO*DO,II,1,12A, II*20+1, II*20+2, II*20+6, II*20+5A, II*20+2, II*20+3, II*20+7, II*20+6A, II*20+3, II*20+4, II*20+8, II*20+7A, II*20+5, II*20+6, II*20+10,II*20+9A, II*20+6, II*20+7, II*20+11,II*20+10A, II*20+7, II*20+8, II*20+12,II*20+11*ENDDO!立柱网格划分LSEL,S,,,181,192 !底层立柱LATT,1,1,1LESIZE,ALL,,,10LMESH,ALLLSEL,S,,,192,324 !2~12层立柱LATT,1,2,1LESIZE,ALL,,,5LMESH,ALL!梁网格划分LSEL,S,,,1,180LATT,1,3,1,LESIZE,ALL,,,8LMESH,ALL!剪力墙网格划分LSEL,S,,,372 !外剪力墙两端、内剪力墙顶端LSEL,A,,,413LSEL,A,,,461LESIZE,ALL,,,8LSEL,s,,,418 !内剪力墙底端LSEL,A,,,420LESIZE,ALL,,,3LSEL,s,,,325 !内剪力墙底端内侧LSEL,A,,,326LESIZE,ALL,,,6LSEL,S,,,375 !剪力墙空洞顶部LSEL,A,,,379LSEL,A,,,383LSEL,A,,,387LSEL,A,,,391LSEL,A,,,395LSEL,A,,,399LSEL,A,,,403LSEL,A,,,407LSEL,A,,,411LSEL,A,,,415LSEL,A,,,419LSEL,A,,,423LSEL,A,,,427LSEL,A,,,431LSEL,A,,,435LSEL,A,,,439LSEL,A,,,443LSEL,A,,,447LSEL,A,,,451LSEL,A,,,455LSEL,A,,,459LSEL,A,,,463LESIZE,ALL,,,2!剪力墙网格划分ASEL,S,,,1,24AATT,1,4,2AMESH,ALL!楼板网格划分ASEL,S,,,25,96AATT,1,5,2AMESH,ALLNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点/eshape,1.0/VIEW,1,1,1,1/ANG,1,270,XM,0/REPlotFINISH!************************************************************************** !************************第二步******************************************** !********************框架结构施加重力荷载********************************* !框架结构施加重力荷载ANTYPE,STATICNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点ACEL,0,0,9.8SOLVEFINISH/POST1SET,FIRSTPLNSOL,U,Z,0,1FINISH!************************************************************************** !************************第三步******************************************** !*************框架结构施加楼面活荷载D=3KN/M^2***************************** !框架结构施加楼面活荷载D=3KN/M^2FINISH/SOLUANTYPE,STATICNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点ACEL,0,0,9.8SOLVEFINISH/POST1SET,FIRSTPLNSOL,U,Z,0,1FINISH!************************************************************************** !************************第四步******************************************** !***********框架结构施加风荷载(基本风压=0.25KN/M^2,体形系数为1.0********** !框架结构施加风荷载(先Y方向,后X方向)分两次分别施加ANTYPE,STATIC*DIM,LOAD_1,ARRAY,12LOAD_1(1)=3.78,2.16,2.39,2.57,2.72,2.84,2.95,3.17,3.20,3.29,3.39,3.51*DIM,LOAD_2,ARRAY,12LOAD_2(1)=7.56,4.32,4.78,5.14,5.44,5.68,5.90,6.34,6.40,6.58,6.76,7.02*DIM,LOAD_3,ARRAY,12LOAD_3(1)=6.93,3.96,4.38,4.71,4.99,5.21,5.41,5.81,5.87,6.03,6.20,6.44*DIM,LOAD_4,ARRAY,12LOAD_4(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93*DIM,LOAD_A,ARRAY,12LOAD_A(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93*DIM,LOAD_B,ARRAY,12LOAD_B(1)=6.30,3.60,3.98,4.28,4.54,4.74,4.92,5.28,5.34,5.48,5.64,5.86*DIM,LOAD_C,ARRAY,12LOAD_C(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93*DO,II,1,12FK,1+20*II,FY,LOAD_1(II)*ENDDO*DO,II,1,12FK,2+20*II,FY,LOAD_2(II)*ENDDO*DO,II,1,12FK,3+20*II,FY,LOAD_3(II)*ENDDO*DO,II,1,12FK,4+20*II,FY,LOAD_4(II)*ENDDOSOLVEFKDELE,ALL,ALL !第二次施加风荷载,删除第一次的。