ANSYS工程结构数值分析命令流1
ANSYS工程结构数值分析第一章及第二章部分命令流
ANSYS工程结构数值分析第一章及第二章部分命令流!e1.1--平面桁架分析!----------前处理------------------ /filname,truss,1!定义工作标题/title,The Analysis of Plane truss /replot/prep7k,1k,2,6k,3,12k,4,18k,5,24k,6,6,8k,7,12,8k,8,18,8l,1,2l,2,3l,3,4l,4,5l,1,6l,6,7l,7,8l,5,8l,2,6l,3,7l,4,8l,6,3l,3,8/pnum,kp,1/pnum,line,1lplotet,1,link1r,1,0.01r,2,0.005r,3,0.0125mp,ex,1,2.1e11mp,prxy,1,0.3lsel,s,tan1,ylatt,1,1,1lsel,s,loc,y,1,7lsel,r,loc,x,6,18latt,1,2,1lsel,s,loc,y,1,7lsel,u,loc,x,6,18latt,1,3,1allsel,all/pnum,real,1/replotlesize,all,,,1lmesh,allfinish!----------加载和求解---------------------------- /solu antype,0lplotdk,1,ux,,,,uydk,5,uyfk,6,fy,-200000fk,8,fy,-200000fk,2,fy,-400000fk,3,fy,-400000fk,4,fy,-400000/pnum,real,0/pnum,line,0/pnum,kp,0/pbc,all,,1/psymb,ldiv,-1lplotsolvefinish!----------后处理------------------------------- /post1pldispetable,mforce,smisc,1etable,mstress,ls,1plls,mforce,mforce,1plls,mstress,mstress,1prrsolfinish!/exit,nosav!e1.2--平面桁架分析!----------前处理-------------------------------------- /filname,truss,1!定义工作标题/title,The Analysis of Plane truss/prep7et,1,link1r,1,0.01r,2,0.005r,3,0.0125mp,ex,1,2.1e11mp,prxy,1,0.3n,1n,2,6n,3,12n,4,18n,5,24n,6,6,8 n,7,12,8 n,8,18,8 type,1 mat,1 real,1e,1,2e,2,3e,3,4e,4,5e,6,7e,7,8 real,2e,2,6e,3,7e,4,8e,6,3e,3,8 real,3e,1,6e,8,5 finish/solu antype,0 d,1,ux,,,,,uyf,6,fy,-200000f,8,fy,-200000f,2,fy,-400000f,3,fy,-400000f,4,fy,-400000solvefinish/post1pldispetable,mforce,smisc,1 etable,mstress,ls,1 plls,mforce,mforce,1 plls,mstress,mstress,1 prrsolfinish!/exit,nosavFINISH/CLEAR/PREP7!1.定义弹簧参数D=4C=8N=10DZ=C*DT=DZ/2.5*IF,T,LT,D,THENT=D*ENDIF!2.创建一圈螺旋线CSYS,1K,1,DZ/2,0,-T/2K,2,DZ/2,180CSYS,0LSYMM,Z,1LSYMM,Y,2,,,,,1NUMMRG,ALLCM,L1,LINE!3.在螺旋线端部创建簧丝截面KWPA VE,1WPROTA,,90CYL4,,,D/2!4.沿L1路径拖拉圆面创建体、复制体等VDRAG,1,,,,,,L1 VGEN,N,ALL,,,,,TNUMMRG,KPWPCSYSFINISH/CLEAR/PREP7!1.定义弹簧参数D=4C=8N=10DZ=C*DT=DZ/2.5*IF,T,LT,D,THENT=D*ENDIF!2.创建一圈螺旋线CSYS,1K,1,DZ/2,0,-T/2K,2,DZ/2,180CSYS,0LSYMM,Z,1LSYMM,Y,2,,,,,1NUMMRG,ALLCM,L1,LINE!3.在螺旋线端部创建簧丝截面KWPA VE,1WPROTA,,90CYL4,,,D/2!4.沿L1路径拖拉圆面创建体、复制体等VDRAG,1,,,,,,L1 VGEN,N,ALL,,,,,TNUMMRG,KPWPCSYSFINISH/CLEAR/PREP7!1.定义弹簧参数D=4C=8N=10DZ=C*DT=DZ/2.5*IF,T,LT,D,THENT=D*ENDIF!2.创建一圈螺旋线CSYS,1K,1,DZ/2,0,-T/2K,2,DZ/2,180CSYS,0LSYMM,Z,1LSYMM,Y,2,,,,,1NUMMRG,ALLCM,L1,LINE!3.在螺旋线端部创建簧丝截面KWPA VE,1WPROTA,,90CYL4,,,D/2!4.沿L1路径拖拉圆面创建体、复制体等VDRAG,1,,,,,,L1 VGEN,N,ALL,,,,,TNUMMRG,KPWPCSYSFINISH/CLEAR/PREP7!1.定义弹簧参数D=4C=8N=10DZ=C*DT=DZ/2.5*IF,T,LT,D,THENT=D*ENDIF!2.创建一圈螺旋线CSYS,1K,1,DZ/2,0,-T/2K,2,DZ/2,180CSYS,0LSYMM,Z,1LSYMM,Y,2,,,,,1NUMMRG,ALLCM,L1,LINE!3.在螺旋线端部创建簧丝截面KWPA VE,1WPROTA,,90CYL4,,,D/2!4.沿L1路径拖拉圆面创建体、复制体等VDRAG,1,,,,,,L1 VGEN,N,ALL,,,,,TNUMMRG,KPWPCSYSFINISH/CLEAR/PREP7!1.定义弹簧参数D=4C=8N=10DZ=C*DT=DZ/2.5*IF,T,LT,D,THENT=D*ENDIF!2.创建一圈螺旋线CSYS,1K,1,DZ/2,0,-T/2K,2,DZ/2,180CSYS,0LSYMM,Z,1LSYMM,Y,2,,,,,1NUMMRG,ALLCM,L1,LINE!3.在螺旋线端部创建簧丝截面KWPA VE,1WPROTA,,90CYL4,,,D/2!4.沿L1路径拖拉圆面创建体、复制体等VDRAG,1,,,,,,L1 VGEN,N,ALL,,,,,TNUMMRG,KPWPCSYSFINISH/CLEAR/PREP7!1.定义弹簧参数D=4C=8N=10DZ=C*DT=DZ/2.5*IF,T,LT,D,THENT=D*ENDIF!2.创建一圈螺旋线CSYS,1K,1,DZ/2,0,-T/2K,2,DZ/2,180CSYS,0LSYMM,Z,1LSYMM,Y,2,,,,,1NUMMRG,ALLCM,L1,LINE!3.在螺旋线端部创建簧丝截面KWPA VE,1WPROTA,,90CYL4,,,D/2!4.沿L1路径拖拉圆面创建体、复制体等VDRAG,1,,,,,,L1 VGEN,N,ALL,,,,,TNUMMRG,KPWPCSYSFINISH/CLEAR/PREP7!1.定义弹簧参数D=4C=8N=10DZ=C*DT=DZ/2.5*IF,T,LT,D,THENT=D*ENDIF!2.创建一圈螺旋线CSYS,1K,1,DZ/2,0,-T/2K,2,DZ/2,180CSYS,0LSYMM,Z,1LSYMM,Y,2,,,,,1NUMMRG,ALLCM,L1,LINE!3.在螺旋线端部创建簧丝截面KWPA VE,1WPROTA,,90CYL4,,,D/2!4.沿L1路径拖拉圆面创建体、复制体等VDRAG,1,,,,,,L1 VGEN,N,ALL,,,,,TNUMMRG,KPWPCSYSFINISH/CLEAR/PREP7!1.定义弹簧参数D=4C=8N=10DZ=C*DT=DZ/2.5*IF,T,LT,D,THENT=D*ENDIF!2.创建一圈螺旋线CSYS,1K,1,DZ/2,0,-T/2K,2,DZ/2,180CSYS,0LSYMM,Z,1LSYMM,Y,2,,,,,1NUMMRG,ALLCM,L1,LINE!3.在螺旋线端部创建簧丝截面KWPA VE,1WPROTA,,90CYL4,,,D/2!4.沿L1路径拖拉圆面创建体、复制体等VDRAG,1,,,,,,L1 VGEN,N,ALL,,,,,TNUMMRG,KPWPCSYS。
(完整版)ANSYS最常用命令流+中文注释(超级大全)
(完整版)ANSYS最常用命令流+中文注释(超级大全)ANSYS最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes,用于2个solid相减操作,最终目的是要nv1-nv2=?通过后面的参数设置,可以得到很多种情况:sepo项是2个体的边界情况,当缺省的时候,是表示2个体相减后,其边界是公用的,当为sepo的时候,表示相减后,2个体有各自的独立边界。
keep1与keep2是询问相减后,保留哪个体?当第一个为keep时,保留nv1,都缺省的时候,操作结果最终只有一个体,比如:vsbv,1,2,sepo,,keep,表示执行1-2的操作,结果是保留体2,体1被删除,还有一个1-2的结果体,现在一共是2个体(即1-2与2),且都各自有自己的边界。
如vsbv,1,2,,keep,,则为1-2后,剩下体1和体1-2,且2个体在边界处公用。
同理,将v换成a 及l是对面和线进行减操作!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: 材料编号(缺省为当前材料号)co: 材料特性值,或材料之特性,温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数定义DP材料:首先要定义EX和泊松比:MP,EX,MA T,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MA T进入单元表并编辑添加单元表: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,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type,是选择的方式,有选择(s),补选(a),不选(u),全选(all)、反选(inv)等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择,loc 就是根据坐标选取,它的comp就可以是实体的某方向坐标!其余还有材料类型、实常数等MIN, VMAX, VINC,这个就不必说了吧!,例:vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14,17,19,21,23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令,输入后直接回车,就可以显示刚刚选择了的体、面或节点,很实用的哦!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”仅用绝对值下面是单元生死第一个载荷步中命令输入示例:!第一个载荷步TIME,... !设定时间值(静力分析选项)NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿-拉夫森选项ESTIF,... !设定非缺省缩减因子(可选)ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)D,ALL,ALL,0 !约束所有不活动的结点自由度(可选)NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。
Ansys命令流大全
Ansys命令流大全ANSYS是一款广泛应用于工程领域的仿真软件,它能够对复杂工程问题进行建模、分析和优化。
本文将提供一个包含常用ANSYS命令的大全,帮助读者快速了解和掌握ANSYS软件的使用。
一、前言ANSYS是一款功能强大的工程仿真软件,它提供了丰富的建模和分析工具,适用于多个领域的工程问题。
掌握ANSYS的命令流能够有效提高工程师的工作效率,快速完成复杂问题的仿真和分析。
二、ANSYS常用命令1. 创建几何模型由于ANSYS提供了多种创建几何模型的工具,我们可以使用命令流来进行几何模型的创建和编辑。
以下是一些常用的几何模型命令:(1)BLOCK:创建矩形或立方体体素模型。
(2)CYLIND:创建圆柱体模型。
(3)SWEEP:创建沿路径扫掠的模型。
2. 定义材料属性在进行仿真分析之前,需要定义材料的物理属性。
以下是一些常用的材料属性命令:(1)MP: 定义材料的参数,如密度、弹性模量、泊松比等。
(2)EX: 定义材料的弹性模量。
(3)DENS: 定义材料的密度。
3. 设定网格划分网格划分对于仿真分析的准确性和计算效率非常重要。
以下是一些常用的网格划分命令:(1)SIZE:设定初始网格尺寸。
(2)MESH:进行自动的网格划分。
(3)ESIZE:设定特定区域的网格尺寸。
4. 定义边界条件在进行仿真分析之前,需要定义边界条件以模拟实际工程环境。
以下是一些常用的边界条件命令:(1)D:定义位移边界条件。
(2)S:定义约束条件。
(3)F:定义外部力或施加力。
5. 设置分析类型ANSYS提供了多种分析类型,如结构分析、热分析、流体分析等。
以下是一些常用的分析类型命令:(1)SOLVE:执行数值分析求解。
(2)ANTYPE:设定分析类型。
(3)FILE:设置解算文件名和保存路径。
6. 查看和后处理结果分析完成后,我们需要查看和后处理结果。
以下是一些常用的结果查看和后处理命令:(1)PLOT:绘制结果曲线或图像。
ansys结构分析入门-命令流
1 ANSYS概述1.1 ANSYS简介ANSYS是一种广泛的商业套装工程分析软件。
所谓工程分析软件,主要是在机械结构系统受到外力负载所出现的反应,例如应力、位移、温度等,根据该反应可知道机械结构系统受到外力负载后的状态,进而判断是否符合设计要求。
一般机械结构系统的几何结构相当复杂,受的负载也相当多,理论分析往往无法进行。
想要解答,必须先简化结构,采用数值模拟方法分析。
由于计算机行业的发展,相应的软件也应运而生,ANSYS软件在工程上应用相当广泛,在机械、电机、土木、电子及航空等领域的使用,都能达到某种程度的可信度,颇获各界好评。
使用该软件,能够降低设计成本,缩短设计时间。
到80年代初期,国际上较大型的面向工程的有限元通用软件主要有:ANSYS, NASTRAN, ASKA, ADINA, SAP等。
以ANSYS为代表的工程数值模拟软件,是一个多用途的有限元法分析软件,它从1971年的2.0版本与今天的5.7版本已有很大的不同,起初它仅提供结构线性分析和热分析,现在可用来求结构、流体、电力、电磁场及碰撞等问题的解答。
它包含了前置处理、解题程序以及后置处理,将有限元分析、计算机图形学和优化技术相结合,已成为现代工程学问题必不可少的有力工具。
1.2ANSYS软件主要功能ANSYS软件是融结构、热、流体、电磁、声学于一体的大型通用有限元软件,可广泛的用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、生物医学、水利、日用家电等一般工业及科学研究。
该软件提供了不断改进的功能清单,具体包括:结构高度非线性分析、电磁分析、计算流体力学分析、设计优化、接触分析、自适应网格划分及利用ANSYS参数设计语言扩展宏命令功能。
1.3ANSYS软件主要特点主要技术特点:∙唯一能实现多场及多场耦合分析的软件∙唯一实现前后处理、求解及多场分析统一数据库的一体化大型FEA软件∙唯一具有多物理场优化功能的FEA软件∙唯一具有中文界面的大型通用有限元软件∙强大的非线性分析功能∙多种求解器分别适用于不同的问题及不同的硬件配置∙支持异种、异构平台的网络浮动,在异种、异构平台上用户界面统一、数据文件全部兼容∙强大的并行计算功能支持分布式并行及共享内存式并行∙多种自动网格划分技术∙良好的用户开发环境支持的图形传递标准:∙SAT∙Parasolid∙STEP与CAD软件的接口∙Unigraphics∙Pro/ENGINEER∙I-Deas∙Catia∙CADDS∙SolidEdge∙SolidWorks1.4 运行环境(ANSYS5.7)Computer: Pentium-class systemMemory (RAM): 64 MB以上Hard Disk: 500MB以上自由空间.Operating System: Microsoft Windows 2000, Windows NT 4.0 (SP 5 or higher) or Windows 98 Graphics: A Windows 2000/NT 4.0 or 98 supported Graphics Card, capable of 1024x768 in High Color (16-bit).A 17 inch monitor (or larger) compatible with the above mentioned card is recommended.2ANSYS 的基本使用2.1 ANSYS环境简介ANSYS有两种模式:一种是交互模式(Interactive Mode),另一个是非交互模式(Batch Mode)。
ansys命令流使用方法
ansys命令流使用方法
在ANSYS中,命令流是一种用于执行特定操作的自动化工具。
以下是ANSYS命令流使用的一般步骤:
1. 打开ANSYS软件并加载您要使用的工程文件。
2. 在ANSYS Graphical User Interface (GUI) 中,将鼠标指针放
在工具栏上。
在“Run”下拉菜单中选择“Command Line”。
3. 在命令行窗口中,输入和编辑您想执行的命令。
您可以使用ANSYS的命令语言以及相关命令进行模型操作、网格生成、
求解等。
4. 您可以通过多种方式输入命令:直接在命令行中输入、从脚本文件中读取、从ANSYS GUI中的日志文件中复制粘贴等。
5. 您可以使用命令流中的参数和变量来进行自动化操作。
使用“!VARIABLE”语句定义变量,并通过“!VARIABLE = value”语
句赋值。
6. 使用ANSYS的各种功能命令对模型进行操作。
例如,在预
处理阶段,您可以使用命令生成几何体、定义材料属性、设定网格、添加边界条件等。
7. 在求解阶段,使用命令启动求解器,设置求解器选项,运行求解器,并监视求解器的输出。
8. 在结果后处理阶段,使用命令读取并处理结果数据,生成图形、报告等。
9. 执行命令流,您可以一次性执行整个命令流,或者逐个执行命令。
10. 您还可以将命令流保存为脚本文件,以便将来再次使用。
以上是ANSYS命令流的一般用法,具体的命令和语法取决于您的特定需求和ANSYS的版本。
建议您参考ANSYS的官方文档和教程,以获得更详细和准确的使用说明。
ANSYS最全命令流解释大全
ANSYS最全命令流解释大全一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目(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次项的系数二、定义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三、单元生死载荷步!第一个载荷步TIME,... !设定时间值(静力分析选项)NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿-拉夫森选项ESTIF,... !设定非缺省缩减因子(可选)ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)D,ALL,ALL,0 !约束所有不活动的结点自由度(可选)NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
ANSYS工程结构数值分析命令流大纲
第二章几何建模技术与技巧
图形控制命令:
显示风格设置命令:
生成动画命令:
组件命令:
第三章网格划分技术及技巧
梁截面:(BEAM188/BEAM189)
内部网格划分控制:
面映射网格划分:
网格修改:
第四章加载与求解技术
施加集中荷载及其相关命令:
施加惯性荷载:
施加耦合场荷载:
非线性选项:
静态分析求解控制:
求解代价估计:
第五章通用与时间历程后处理技术
列表显示结果及其相关命令:
节点特殊结果计算命令:
路径及操作:
变量定义及其相关命令:
变量运算命令:。
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 -指定阶段状或者用跳板装载里面一个负荷步骤。
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, VALUELab为参数名称 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杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK183D有限应EPCDFG另可考虑粘0 变杆 B 弹塑性E-弹性(Elasticity),P-塑性(Plasticity),C-蠕变(Creep),S-膨胀(Swelling),D-大变形或大挠度(Large deflection),F-大应变(Large strain)或有限应变(Finite strain),B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening),A-自适应下降(Adaptive descent)等。
ANSYS单元特性之命令流算例
EX1.1 (LINK1)(1)进入后处理模块,显示节点位移和杆件内力MID_NO DE = NODE (A/2,-B,0 )! 寻找距离位置(A/2,-B,0)最近的点,存入MID_NODE*GET,DISP,NODE,MID_NO DE,U,Y!提取节点MI D_NOD E上的位移UY,若果已知要求的节点,直接提取即可。
LEFT_E L = ENEARN (MID_NO DE)! 需找距离节点MID_N ODE最近的单元,存入LEFT_ELETABLE,STRS,L S,1! 用轴向应力S AXL的编号“LS,1”定义单元表S TRS*GET,STRSS,ELEM,LEFT_E L,ETAB,STRS! 从单元表STRS中提取LEFT_EL单元的应力结果,存入变量STRSS。
注意:提取的轴向应力结果具体到指定的单元。
(2)申明数组,提取计算结果,并比较计算误差*DIM,LABEL,CHAR,2!定义2个元素的字符型数组LABE L*DIM,VALUE,,2,3!定义2*3的数值型数组V ALU ELABEL(1) = 'STRS_M Pa','DEF_mm' ! 给字符型数组的第1个元素赋值*VFILL,VALUE(1,1),DATA,1,-0.05498! 给其他数值型数组中的元素赋值*VFILL,VALUE(1,2),DATA,STRSS,DISP*VFILL,VALUE(1,3),DATA,ABS(STRSS/1 ) ,ABS( DISP /0.05498)/OUT,EX1_1,out !将输出内容重定向到文件EX1_1.out/COM ! 以注释形式输出内容/COM,------------------- EX1.1 RESULT S COMPAR ISON---------------------/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,----------------------------------------------------------------/OUT ! 结束数据重定向,关闭输出文件FINISH*LIST,EX1_1,out ! 列表显示文件EX1_1.out的内容EX1.2 (LINK1)/PNUM, NODE,1!打开节点编号显示/NUMBER, 2!只显示编号,不使用色彩列表显示节点位移和单元的计算结果PRDISP! 列表显示节点位移值计算结果ETABLE, MFORX,SMISC,1!以杆单元的轴力为内容,建立单元表M F ORXETABLE, SAXL, LS, 1 !以杆的轴向应力为内容,建立单元表S A XLETABLE, EPELAX L, LEPEL, 1! 以杆单元的轴向应变为内容,建立单元表E PELAX LPRETAB! 显示单元表中的计算结果/NUMBER, 0!显示编号,并使用彩色PLETAB, MFORX!用色度图显示单元表MF ORX中杆件轴力图EX1.3 (LINK1)NSEL,S,LOC,Y,1.0 !选择所有位于Y=1.0位置上的节点FSUM!累计叠加选择集中所有节点上的反力*GET,REAC_1,FSUM,,ITEM,FY !将累加结果中的FY(Y方向的力)保存到变量R E AC_1中EX1.4 (LINK1)R,1,65e-6!定义第1类实常数,杆件截面面积为65mm^2,在转化为国际单位制时操作TREF,70 ! 设定参考温度为70度BFUNIF,TEMP,80 ! 温度从原来的70度均匀上升到80度(TREF+10)EX1.5 (PLANE42 AND CONT AC26)ETABLE,STRSX,S,X!定义X方向的应力为单元表S TRS X*GET,STRSSX,ELEM,3,ETAB,STRSX!从单元表STRS X中提取3号单元的X向应力,存入STRS S X。
(完整版)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 -指定阶段状或者用跳板装载里面一个负荷步骤。
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 )等。
Ansys命令流大全(整理)
Ansys命令流大全(整理)1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点(P1~P9)来定义面(Area),最少使用三个点才能围成面,同时产生转围绕些面的线。
点要依次序输入,输入的顺序会决定面的法线方向。
如果超过四个点,则这些点必须在同一个平面上。
Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR,Abbr,String--定义一个缩略语.Abbr:用来表示字符串"String"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.3、ABBRES,Lab,Fname,Ext-从一个编码文件中读出缩略语.Lab:指定读操作的标题,NEW:用这些读出的缩略语重新取代当前的缩略语(默认)CHANGE:将读出的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语.Ext:如果"Fname"是空的,则缺省的扩展命是"ABBR".4、ABBSA V,Lab,Fname,Ext-将当前的缩略语写入一个文本文件里Lab:指定写操作的标题,若为ALL,表示将所有的缩略语都写入文件(默认)5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic:变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp !kswp=0时只删除掉面积本身,=1时低单元点一并删除。
7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 !面积的建立,沿某组线段路径,拉伸而成。
8、Afillt,na1,na2,rad !建立圆角面积,在两相交平面间产生曲面,rad为半径。
ansys命令流中文说明
ansys命令流中文说明展开全文KB、KE: 待划分线的定向关键点起始、终止号SECNUM: 截面类型号u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分SECID:由SECTYPE命令分配的截面编号MESHKEY:0:不显示网格划分1:显示网格划分u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元SCALE: 0:简单显示线、面单元1:使用实常数显示单元形状u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元xnode: 仅为产生surf151 或surf152单元时使用tlab: 仅用来生成接触元或目标元top 产生单元且法线方向与所覆盖的单元相同,仅对梁或壳有效,对实体单元无效Bottom产生单元且法线方向与所覆盖的单元相反,仅对梁或壳有效,对实体单元无效Reverse 将已产生单元反向Shape: 空与所覆盖单元形状相同Tri 产生三角形表面的目标元注意:选中的单元是由所选节点决定的,而不是选单元,如同将压力加在节点上而不是单元上u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的itemlabel: 要合并的项目node: 节点, Elem,单元,kp: 关键点(也合并线,面及点)mat: 材料,type: 单元类型,Real: 实常数cp:耦合项,CE:约束项,CE: 约束方程,All:所有项toler: 公差Gtoler:实体公差Action: sele 仅选择不合并空合并switch: 较低号还是较高号被保留(low, high)注意:可以先选择一部分项目,再执行合并。
如果多次发生合并命令,一定要先合并节点,再合并关键点。
合并节点后,实体荷载不能转化到单元,此时可合并关键点解决问题。
u 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 选择线及相关关键点、节点和单元u 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”仅用绝对值u NSLL,type, nkey 选择与所选线相联系的节点u nsla, type, nkey: 选择与选中面相关的节点type:s 选一套新节点r 从已选节点中再选a 附加一部分节点到已选节点u 从已选节点中去除一部分nkey: 0 仅选面内的节点1 选所有和面相联系的节点(如面内线,关键点处的节点)u esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元Type: S: 选择一组单元(缺省)R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组(?)Stat: 显示当前选择状态Item: Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号u ALLSEL, LABT, ENTITY 选中所有项目LABT: ALL: 选所有项目及其低级项目BELOW: 选指定项目的直接下属及更低级项目ENTITY: ALL: 所有项目(缺省)VOLU:体高级AREA:面LINE :线KP:关键点ELEM:单元NODE:节点低级u Tshap,shape 定义接触目标面为2D、3D的简单图形Shape: line:直线Arc:顺时针弧Tria:3点三角形Quad:4点四边形………….2.6 根据需要耦合某些节点自由度u cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotznode1-node17: 待耦合的节点号。
ANSYS最全命令流解释大全(DOC)
一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目(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次项的系数二、定义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三、单元生死载荷步!第一个载荷步TIME,... !设定时间值(静力分析选项)NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿-拉夫森选项ESTIF,... !设定非缺省缩减因子(可选)ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)D,ALL,ALL,0 !约束所有不活动的结点自由度(可选)NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
ANSYS一段命令流解释
ANSYS命令流解释/PREP7 ! 进入前处理器ET,1,BEAM4 ! 定义单元beam4,是个三维梁单元KEYOPT,1,2,0 ! 定义单元的关键选项,如果后面是0,代表默认的,可以先不用理解,KEYOPT,1,6,0KEYOPT,1,7,1 ! 这里定义第七个关键选项,定义为编号1,也就是计算陀螺阻尼矩阵方程,要求! IYY 等于IZZ,也就是两个转动惯量要相等,这两个量要在实常数中定义! 也就是下面的命令RKEYOPT,1,9,0KEYOPT,1,10,0*SET,p,acos(-1) ! 定义三个参数,分别是派(4.1315)、第一个半径R1,第二个半径r2*SET,R1,5*SET,R2,60 ! 半径的单位! 定义单元的实常数,有两个,因为有两个半径,分别就是下面的R,1! 和r,2! r命令的定义中需要根据使用单元beam4来一一对应,不同的单元R命令! 定义的意思是不一样的,具体每个意思,下图看,一定要一一对应! 这里就先是截面积,Z向转动惯量,Y像转动惯量,是一样的,从上面的! 关键选项定义中可以看的出来,KEYOPT,1,7,1 R,1,p*R1**2,p*R1**4/4,p*R1**4/4,2*R1,2*R1, ,RMORE, ,p*R1**4/2, , ,2175, , ! 这个也是定义实常数,因为命令只能定义6个数,从第七个! 就要使用这个命令,编号可以从下图中对应下R,2,p*R2**2,p*R2**4/4,p*R2**4/4,2*R2,2*R2, , ! 定义第二个实常数R,2 RMORE, ,p*R2**4/2, , ,2175, ,MPTEMP,,,,,,,, ! 这里来定义材料属性,目前楼主提供的是GUI操作以后的! ,这个需要相应简化,因为GUi操作中一般包括了温度的考虑,所以使用了几个温度的命令! 目前看来只是有一个温度,所以不用考虑温度,简化后的,在后面的注释中,注释开始是!(感叹号)MPTEMP,1,0MPDATA,EX,1,,2e5MPDATA,PRXY,1,,.3MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,1e-10! 简化后的材料属性! mp,ex,1,2e5 ! 分别定义弹性模量,泊松比,密度,其中的1代表第一个材料属性! mp,PRXY,1,0.3! mp,DENS,1,1e-10MPTEMP,,,,,,,, ! 这里就是定义第二个材料属性了!同样包括弹性模量,泊松比,密度MPTEMP,1,0MPDATA,EX,2,,2E5MPDATA,PRXY,2,,.3MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,8E-8K, ,,,, ! 定义两个关键点,第一个是0,第二个坐标是:(100,0,0) K, ,100,,,TYPE, 1 ! 选择第一个单元,也就是beam4MAT, 1 ! 选择第一个材料属性REAL, 1 ! 选择第一个实常数ESYS, 0 ! 选择默认的单元坐标系统LSTR, 1, 2 ! 用刚才生成的两个关键点,建立一条直线LESIZE,ALL, , ,200, ,1, , ,1, ! 把建立的这条直线划分为200个LMESH, 1 ! 根据划分的线的段数,网格这条线,从而生成有限元模型D,1,UX ! 对于第一个节点(也就是原点位置的那个节点)进行约束,三个平移方向都约束,但是转动方向不约束,! 因为单元beam4是有六个自由度的D,1,UYD,1,UZD,102,UY ! 对于第102个节点,这样来算他的位置,总共长度是100,划分为200分,那么102是那个位置,楼主应该知道了吧! 选择了Y方向和Z方向的平移约束D,102,UZFLST,2,1,2,ORDE,1 ! 后面的两个命令结合起来就是选择了第200个节点,但是没有具体操作,! 只是有了拾取选择这个点的操作,这两个命令一般是GUI直接生成的,写命令流的一般不使用这两个操作! 所以这两个命令其实什么任务都没有做目前FITEM,2,200对上面的命令一一般是需要简化的,简化后的就是:/PREP7ET,1,BEAM4KEYOPT,1,7,1*SET,p,acos(-1)*SET,R1,5*SET,R2,60R,1,p*R1**2,p*R1**4/4,p*R1**4/4,2*R1,2*R1, ,RMORE, ,p*R1**4/2, , ,2175, ,R,2,p*R2**2,p*R2**4/4,p*R2**4/4,2*R2,2*R2, ,RMORE, ,p*R2**4/2, , ,2175, ,mp,ex,1,2e5mp,PRXY,1,0.3mp,DENS,1,1e-10mp,ex,2,2e5mp,PRXY,2,0.3mp,DENS,2,8e-8K,1,,,,K,2,100,,,MAT, 1REAL, 1LSTR, 1, 2LESIZE,ALL, , ,200LMESH, 1D,1,UXD,1,UYD,1,UZD,102,UYD,102,UZFLST,2,1,2,ORDE,1FITEM,2,200简化后的命令楼主结合本人说的,可以理解下这段话因为只是在建模,所以没有执行什么计算的,思路其实我已经说的很明确了:首先定义单元,指定一个关键选项包括了陀螺阻尼效果,然后再定义两个实常数,然后定义了两个材料属性,然后建立一个一个线,使用第一个实常数,和第一个材料属性对这个线进行了网格!其次再对这个线进行了约束的施加,原点位置三个平移方向全约束,对102个节点位置,进行了两个方向的约束。
ANSYS命令流学习笔记1-link单元的静力分析
前言刚学习命令流,作为学习的记录。
接触ANSYS已经4年了,开始用了一年多的经典界面。
作为一名工程师,自从用了WorkBench,基本上放弃了经典界面。
对于工程问题上,有前辈表示WorkBench也足够了,但是有时候大家也会写一些APDL命令,毕竟WB的友好界面很难实现底层功能。
所以开始学习如何运用命令流,从写简单的命令开始。
通过学习的过程也强化自己的一些认识,笔记也能让自己的思考更严谨。
至于以后能不能写复杂了,那就看命运吧。
主要参考来自王新敏老师的《ANSYS工程结构的数值分析》,Saeed Moaveni的《有限元分析-ANSYS理论与应用》,以及各网站帖子。
link单元的静力分析1. 命令流内容/clear/filname,linkf,1/title, beam static analysis/com,structural !定义一些文件属性/prep7 !进入前处理et,1,link180 !三维桁架单元r,1,10 !实常数,面积mp,ex,1,21e10mp,nuxy,1,0.3n,1,0,0n,2,100,0n,3,50,100/pnum,node,1nplote,1,2e,2,3e,1,3/pnum,elem,1 !显示编码eplotfinish !结束前处理/solu !进入求解d,1,all,0 !约束d,2,all,0 !约束f,3,fy,-100 !加力/pbc,all,1 !显示BCeplot !显示solvefinish !求解结束/post1 !进入后处理etable,axforce,smisc,1 !创建轴向力单元表etable,axstress,ls,1 !创建轴向应力单元表/pbc,all,1pldisp,1 !变形,图形形式pletab,axstressprnsol,u,comp !节点位移,表格形式prrsol !反作用力,表格形式输出反作用力pretab !单元表,输出之前定义的单元表格,轴向力及轴向应力finish !后处理结束EOF !end of file总结工程上来说1.2节点形成的杆1是没用的。
ansys命令流实例
ansys命令流实例1.结构静力分析--梁分析/PREP7 !进入前处理器ET,1,BEAM3 !定义单元类型R,1,0.25,0.0052,0.5 !定义实常数MP,EX,1,210E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,1N,3,5N,4,7N,5,7N,6,9N,7,11N,8,11N,9,13N,10,14 !定义节点E,1,2E,2,3E,3,4E,5,6E,6,7E,8,9E,9,10 !生成单元CP,1,UX,4,5CP,2,UY,4,5CP,3,UX,7,8CP,4,UY,7,8 !耦合节点FINISH/SOLU !进入求解器D,2,UXD,2,UYD,3,UYD,6,UYD,9,UY !施加位移约束F,10,FY,-4 !施加集中约束SFBEAM,1,1,PRES,4,4SFBEAM,2,1,PRES,4,4 !施加均布力SOLVE !求解/POST1 !进入后处理器PLDISP !绘制结构变形图PRDISP !列出各节点的位移ETABLE,IMOMENT,SMISC,6ETABLE,JMOMENT,SMISC,12ETABLE,ISHEAR,SMISC,2ETABLE,JSHEAR,SMISC,8 !将节点弯矩、剪力制表PRETAB !列表显示单元的弯矩、剪力/TITLE,SHEAR FORCE DISTRIBUTION !设置剪力分布图的标题PLLS,ISHEAR,JSHEAR !绘制剪力分布图/TITLE,BENDING MOMENT IDSTRIBUTION !设置弯矩分布图的标题PLLS,IMOMENT,JMOMENT !绘制弯矩分布图流固耦合实例RAD=0.8 !底面半径H=1G=9.8OMEGAR=2ROU=1000 !定义参数变量/PREP7 !进入前处理器ET,1,FLUID79 !选择单元类型KEYOPT,1,3,1 !设置单元关键字MP,EX,1,2E9 !设置杨氏模量MP,DENS,1,ROU !设置材料密度K,1K,2,RADK,3,RAD,HK,4,,H !生成关键点A,1,2,3,4 !连接关键点生成面积LESIZE,ALL,,,10 !设置网格划分精度AMESH,ALL !将面积划分网络/SOLU !进入求解器DL,2,,UXDL,1,,UYNSEL,S,LOC,XDSYM,SYMM,XD,ALL,UXNSEL,ALL !施加位移约束ACEL,,GOMEGA,,OMEGAR !施加惯性力SOLVE !求解/POST1SET,LAST !进入通用后处理器PLNSOL,U,X,0,1 !绘制应力云图UCENT=UY(22)UEDGE=UY(12)UELEV=UEDGE-UCENT !提取节点位移结构静力分析--壳结构内力分析LENGTH=100YOUNG=200000THICKNESS=2FORCE=1000DENSITY=9E-6 !将材质、载荷、板的几何尺寸等参数化/PREP7 !进入前处理器MP,EX,1,YOUNGMP,NUXY,1,0.3MP,DENS,1,DENSITY !定义材质ET,1,SHELL63 !定义单元类型R,1,THICKNESS,THICKNESS,THICKNESS,THICKNESS !定义实常数!构建结构的几何模型K,1,0,0K,2,LENGTH,0K,3,LENGTH,LENGTHK,4,0,LENGTH !定义关键点A,1,2,3,4!通过关键点生成面LSEL,ALLLESIZE,ALL,,,16AMESH,ALL !设定网格划分参数,划分网格FINISH/SOLU !进入求解器NSEL,S,LOC,X,0,0D,ALL,ALL,0 !选择X=0的节点将其固定NSEL,S,LOC,X,LENGTH,LENGTHD,ALL,ALL,0 !选择X=LENGTH的节点将其固定NSEL,S,LOC,X,0.5*LENGTH,0.5*LENGTHNSEL,R,LOC,Y,0.5*LENGTH,0.5*LENGTHF,ALL,FZ,FORCEALLSEL !捕捉板的中心点并在中心点处施加集中力荷载SOLVE !求解FINISH/POST1 !进入后处理器/DSC,,10PLNSOL,U,Z,0,1 !绘图显示板的竖向变形NSEL,ALL !提取板的最大竖向变形NSORT,U,Z,1,1 !将节点的位移绝对值以升序排序*GET,MAXDEFLECTION,SORT,0,MAX !提取位移最大值并赋给变量NSEL,S,LOC,X,0NSEL,A,LOC,X,LENGTH,LENGTH !选择固定边节点NSORT,S,EQV,1,1 !将节点等效应力的绝对值以升序排序*GET,MAXSTRESS,SORT,0,MAX !提取等效应力最大值并赋给变量结构静力分析--平面桁架分析/PREP7 !进入前处理器ET,1,LINK1 !选择单元R,1,0.1 !定义实常数MP,EX,1,30E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,4N,3,8N,4,12N,5,0,3N,6,4,3N,7,8,3 !生成节点E,1,2E,2,3E,3,4E,4,7E,3,7E,2,7E,2,6E,2,5E,1,5E,5,6E,6,7 !生成单元FINISH/SOLU !进入求解器D,1,UXD,1,UYD,5,UX !实加位移约束F,2,FY,-15F,3,FY,-15F,4,FY,-15 !施加集中力SOLVE !求解/POST1 !进入能用后处理器PRESOL,FORC !列表显示反力。
ANSYS命令流详解(超全)
ANSYS命令流详解(超全)一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目(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次项的系数二、定义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三、单元生死载荷步!第一个载荷步TIME,... !设定时间值(静力分析选项)NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿-拉夫森选项ESTIF,... !设定非缺省缩减因子(可选)ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)D,ALL,ALL,0 !约束所有不活动的结点自由度(可选)NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
ansys命令流语法
ansys命令流语法ANSYS命令流语法是使用ANSYS软件进行仿真分析的关键部分。
它是一种将命令以特定顺序组合在一起的方式,以实现特定的分析目标。
本文将介绍ANSYS命令流语法的基本语法规则和常用命令,以及如何使用它们进行仿真分析。
一、ANSYS命令流语法的基本语法规则1. 命令的基本格式:命令[选项] [参数1, 参数2, ...]2. 命令的执行顺序:ANSYS命令流是按照命令的顺序逐条执行的。
如果需要改变执行顺序,可以使用条件语句、循环语句等控制结构。
3. 注释:可以在命令流中添加注释,以"!"开头。
注释部分不会被执行,可以用于解释命令的用途或添加说明。
4. 变量和参数:可以使用变量和参数来存储和传递数据。
变量以"$"开头,参数以"%"开头。
二、常用命令1. Preprocessor命令:用于定义和准备分析模型的预处理操作。
- /PREP7:进入预处理器界面。
- ET,MP,REAL等:定义单元类型、材料属性、实数等。
- K,L,A等:创建节点、单元、区域等。
2. Solution命令:用于设置和运行分析求解器。
- /SOLU:进入求解器界面。
- SOLVE,ANTYPE等:设置分析类型、求解选项等。
- D,S等:定义边界条件、加载条件等。
3. Postprocessor命令:用于后处理和分析结果的可视化。
- /POST1:进入后处理器界面。
- PLOT,PDEF等:绘制图形、定义图形属性等。
- PRINT,*VWRITE等:输出结果数据。
三、使用ANSYS命令流语法进行仿真分析使用ANSYS命令流语法进行仿真分析的一般步骤如下:1. 导入几何模型:使用CAD软件创建几何模型,并将其导入ANSYS 中。
2. 定义材料属性:根据实际材料的物理特性,使用MP命令定义材料属性。
3. 网格划分:使用网格划分命令划分几何模型,生成有限元网格。
ANSYS工程结构数值分析
作者简介
这是《ANSYS工程结构数值分析》的读书笔记,暂无该书作者的介绍。
谢谢观看
பைடு நூலகம்
目录分析
在当今的工程领域,数值分析已经成为了一个不可或缺的工具。其中, ANSYS是一款广泛使用的工程分析软件,能够帮助工程师和设计师对复杂工程结 构进行有效的数值分析。为了更好地理解和使用ANSYS,一本名为《ANSYS工程结 构数值分析》的书籍应运而生。本书将对这本书的目录进行分析,探讨其内容构 成及特点。
工程结构数值分析是一门涵盖广泛领域的学科,它涵盖了结构分析、力学、 计算机科学等多个方面。在众多工程结构数值分析的书籍中,《ANSYS工程结构 数值分析》以其独特的视角和深入浅出的讲解方式受到了广大读者的喜爱。这本 书不仅提供了对工程结构数值分析的基本理解,还通过丰富的案例和实用的技巧 展示了如何利用ANSYS软件进行实际操作。
从整体来看,这本书的目录结构非常清晰,内容由浅入深,循序渐进。全书 共分为十章,涵盖了工程结构数值分析的基本原理、ANSYS软件的基础知识、各 种工程结构的数值分析方法以及实例分析等内容。这样的编排使得读者能够很容 易地找到自己需要的知识点,并进行针对性的学习。
每个章节的内容都十分详尽。作者在编写过程中,对每个知识点都进行了详 细的讲解,使得读者能够全面了解和掌握工程结构数值分析的基本概念和方法。 同时,书中还提供了大量的实例和练习题,这有助于读者更好地理解和应用所学 知识,提高分析问题和解决问题的能力。
这本书还注重理论与实践相结合。在介绍了基本的数值分析原理和ANSYS软 件知识后,作者给出了各种工程结构的数值分析实例。这些实例涉及到了桥梁、 建筑、机械等领域,具有很强的实用性。通过这些实例的学习,读者可以更好地 理解数值分析方法在工程实践中的应用,提高自己的实际操作能力。
ANSYS工程结构数值分析---命令流实例
!------------------------------------------------------!EX8.26 钢筋混凝土简支梁数值分析!分离式模型,关闭压碎,keyopt(1)=0,keyopt(7)=1 !力加载,位移收敛准则,误差1.5%,1/4模型分析!--------------------------------------------finish/clear/config,nres,2000/prep7!1.定义单元与材料性质--------------------et,1,solid65,,,,,,,1 !K1=0,k7=1et,2,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,-1tb,miso,1,,11tbpt,,0.0002,fc*0.19tbpt,,0.0004,fc*0.36tbpt,,0.0006,fc*0.51tbpt,,0.0008,fc*0.64tbpt,,0.001,fc*0.75tbpt,,0.0012,fc*0.84tbpt,,0.0014,fc*0.91tbpt,,0.0016,fc*0.96tbpt,,0.0018,fc*0.99tbpt,,0.002,fctbpt,,0.0033,fc*0.85mp,ex,2,2.0e5mp,prxy,2,0.3tb,biso,2tbdata,,300,0pi=acos(-1)r,1,0.25*pi*22*22r,2,0.25*pi*22*22/2r,3,0.25*pi*10*10r,4,0.25*pi*10*10/2!2.创建几何模型blc4,,,150/2,300,2000/2*do,i,1,9wpoff,,,100vsbw,all*enddo wpcsys,-1 wpoff,,,50 vsbw,all wpcsys,-1 wprota,,-90 wpoff,,,30 vsbw,all wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1!3.划分钢筋网格elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,linelatt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,tan1,z lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2lesize,all,elemsizlsel,allcmsel,s,zjcmsel,a,zjbcmsel,a,jljcmsel,a,gjcmsel,a,gjbcm,gj,linelmesh,alllsel,all!4.划分混凝土网格vatt,1,,1mshkey,1esize,elemsizvmesh,allallsel,all!5.施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symm!-----------------p0=180000q0=p0/150/100asel,s,loc,z,0,50asel,r,loc,y,300sfa,all,1,pres,q0allsel,all!6.求解控制设置/soluantype,0nsubst,80outres,all,allautos,onneqit,50cnvtol,u,,0.015solve!7.进入POST1查看结果/post1set,lastesel,s,type,,2etable,saxl,ls,1plls,saxl,saxlesel,s,type,,1/device,vector,onplcrackeall!8.进入时程后处理/post26nsol,2,205,u,yprod,3,1,,,,,,p0/1000prod,4,2,,,,,,-1xvar,4plvar,3!EX2.6D 60kg/m钢轨实体建模FINISH/CLEAR/PREP7!1.定义参数H1=30.5H2=48.5H3=48.5B1=45.75B2=29.25B3=16.5B4=36.5B5=10XD1=1/3XD2=1/9XD3=1/3XD4=1/20R0=400R1=2R2=4R3=40R4=20R5=25R7=13R8=80R9=300R10=2FAI=43H=H1+2*H2+H3!2.创建轨底主要组成部分的线K,1K,2,B1+B2K,3,B1+B2,H1-B1*XD1-B2*XD2K,4,B1,H1-B1*XD1K,5,0,H1*DO,I,1,4L,I,I+1*ENDDO!3.计算确定轨头下三个圆弧的位置CTA=A TAN(XD3)BTA=(90-FAI)*ACOS(-1)/180-CTAA=1/TAN(CTA)C1=B3/2+R0-R6*SIN(CTA)-(R5-R6)*COS(BTA)C2=R5*SIN(BTA)+R6*COS(CTA)-R6*SIN(BTA)-H2C=(C1-C2/TAN(CTA))/(R0-R5)C3=A*C-SQRT(A*A-C*C+1)C3=C3/(A*A+1)REFA=ASIN(C3)C1=R6*SIN(CTA)+(R5-R6)*COS(BTA)+(R0-R5)*COS(REFA) !4.创建轨腰和轨头结构线ROUX=B3/2+R0-C1ROUY=ROUX*TAN(CTA)ANGI=ASIN((H1+H2)/R0)K,6,B3/2+R0-R0*COS(ANGI)K,7,B3/2+R0-R0*COS(REFA),H1+H2+R0*SIN(REFA) LARC,6,7,2,R0DX1=R6*(COS(BTA)-SIN(CTA))DY1=R6*(COS(CTA)-SIN(BTA))K,8,ROUX-DX1,H1+2*H2+ROUY-DY1K,9,ROUX,H1+2*H2+ROUYLARC,7,8,2,R5LARC,8,9,2,R6K,10,B4,H1+2*H2+B4*XD3K,11,B4-(H3-B4*XD3)*XD4,HL,9,10L,10,11ANGI=ASIN(B5/R9)K,12,0,HK,13,B5,H-(R9-R9*COS(ANGI)) LARC,12,13,1,R9DX1=B4-B5+R8*SIN(ANGI)ANGI1=ASIN(DX1/R8)Y1=R8*COS(ANGI)-R8*COS(ANGI1) K,14,B4,KY(13)-Y1LARC,13,14,1,R8!5.倒角或弧线连接LFILLT,1,2,R1LFILLT,2,3,R2LFILLT,3,4,R3LPTN,4,5LFILLT,17,18,R4LFILLT,8,9,R10LPTN,9,11LFILLT,20,22,R7LDELE,15,16,1,1LDELE,19,21,2,1LSYMM,X,ALLNUMMRG,ALLNUMCMP,ALLAL,ALLASUMVOFFST,1,300!EX2.9B 斜向圆台建模finish/clear/prep7r0=3r1=2h=4cyl4,,,r0wpoff,r1,,hcyl4,,,r1*do,i,1,4askin,i,i+4*enddova,all!斜向延伸创建圆台vgen,1,1,,,,3*r0,,,,1wpcsyscyl4,,,r0vext,7,,,r1,,h,r1/r0,r1/r0!EX2.9B 已知函数方程时的建模!玫瑰线1==================================== !极坐标方程为ρ=asin(4Φ/3)finish/clear/prep7a=20csys,1*do,i,0,360*3fei=i*acos(-1)/180r0=a*sin(fei*4/3)k,,r0,i*enddo*do,i,1,360*3l,i,i+1*enddo!玫瑰线2=================================== !极坐标方程为ρ=asin(4Φ)finish/clear/prep7a=20csys,1*do,i,0,360fei=i*acos(-1)/180r0=a*sin(fei*4)k,,r0,i*enddo*do,i,1,360l,i,i+1*enddo!长辐圆内旋轮线============================ !参数方程x=(a-b)cost+λcos(a-b)/bt!参数方程y=(a-b)sint-λsin(a-b)/btfinish/clear/prep7a=50b=5lmda=15*afun,deg*do,i,0,360x1=(a-b)*cos(i)+lmda*cos((a-b)/b*i)y1=(a-b)*sin(i)-lmda*sin((a-b)/b*i)k,,x1,y1*enddo*do,i,1,360l,i,i+1*enddo!椭圆抛物面================================ !方程为z=x*x/a/a+y*y/b/b!使用蒙皮方法finish/clear/prep7a=20b=10n=20*do,i,1,niz=i/10lsel,nonewpoff,,,izcyl4,,,1adele,alllsscale,all,,,a*sqrt(iz),b*sqrt(iz),,,,1*enddoallsel*do,i,1,4*(n-1),4askin,i,i+4askin,i+1,i+5askin,i+2,i+6askin,i+3,i+7*enddo!EX2.9C 旋转图素建模!类花瓶建模--------------------finish/clear/prep7K,1K,2,100K,3,300,400K,4,150,550k,5,130,880K,6,300,1150k,7,400,1100k,8,500,1150bsplin,allk,1001,,1000arotat,all,,,,,,1,1001,360,5!类轮建模----------------------finish/clear/prep7!先创建类轮结构的一个断面,然后旋转之blc4,,,4,6cyl4,2,6,1.5asba,1,2k,100,-3k,101,-3,10vrotat,all,,,,,,100,101!EX3.1 自定义箱形截面finish/clear/prep7k,1k,2,2k,3,2,2.2k,4,3,2.3k,5,3,2.5k,6,0,2.5a,1,2,3,4,5,6blc4,,0.2,1.7,2asba,1,2wpoff,1.7wprota,,,90asbw,allwpoff,,,0.3asbw,allwpoff,,0.2wprota,,90asbw,allwpoff,,,-2asbw,allarsym,x,allwpcsysaglue,allet,1,plane82lesize,all,,,1amesh,allsecwrite,myboxfinish/clear/prep7et,1,beam189mp,ex,1,3.0e10mp,prxy,1,0.167sectype,1,beam,mesh secread,mybox,,,meshk,1k,2,10k,100,5,5l,1,2lesize,all,,,20latt,1,,1,,100,,1lmesh,all/eshape,1Eplot!EX3.2 自定义多种材料截面finish/clear/prep7Ro=1.5Ri=1.0csys,1cyl4,,,ricyl4,,,roaptn,allwprota,,90asbwa,allwprota,,,90asbw,allwpcsyset,1,plane82mymat1=4mymat2=7mp,ex,mymat1,1.0mp,ex,mymat2,2.0asel,s,loc,x,0,riaatt,mymat1,,1asel,s,loc,x,ri,roaatt,mymat2,,1allselesize,0.25mshape,0,2dmshkey,1amesh,allsecwrite,mycsolid,sectfinish/clear/prep7et,1,beam189mym1=4mym2=7mp,ex,mym1,3.0e10mp,prxy,mym1,0.167 mp,ex,mym2,2.1e11mp,prxy,mym2,0.3 sectype,1,beam,mesh secread,mycsolid,sect,,mesh k,1k,2,,,10l,1,2lesize,all,,,20latt,,,1,,,,1lmesh,all/eshape,1/pnum,mat,1eplot!EX3.3A 单个方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8K,1K,2,,,1000L,1,2K,100,,500,500LATT,1,,1,,100,,1LGEN,4,1,,,500LESIZE,ALL,,,10LMESH,ALL/ESHAPE,1EPLOT!EX3.3B 同时使用KB和KE方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0L,1,2K,100,,DLK,200,DXC,-DLK,300,2*DXC,DLK,301,2*DXC+DLK,400K,500,8*DXCLGEN,5,1,,,DXCLATT,1,,1,,100,,1LSEL,S,,,2LATT,1,,1,,200,,1LSEL,S,,,3LATT,1,,1,,300,301,1LSEL,S,,,4LATT,1,,1,,400,,1LSEL,S,,,5LATT,1,,1,,500,,1LSEL,ALLLESIZE,ALL,,,50LMESH,ALL/ESHAPE,1EPLOT!EX3.3C 双方位关键点翻转示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,I SECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0K,3,DXCK,4,DXC,,L0/2K,5,DXC,,L0L,1,2L,3,4L,4,5K,100,,DLk,101,,-DLK,200,DXC,DLk,201,2*DXCK,202,DXC,-DLLATT,1,,1,,100,101,1lesize,all,,,100LSEL,S,,,2LATT,1,,1,,200,201,1lesize,all,,,50LSEL,S,,,3LATT,1,,1,,201,202,1lesize,all,,,50LSEL,ALLLMESH,ALL/ESHAPE,1EPLOT!EX3.4A 两种单元形状和两种网格划分比较finish/clear/prep7et,1,plane82k,1k,2,8k,3,7,6k,4,1,6a,1,2,3,4esize,1mshape,0mshkey,0!mshape,0!mshkey,1!mshape,1!mshkey,0!mshape,1!mshkey,1amesh,all!EX3.4B 中间节点位置控制网格划分比较finish/clear/prep7cyl4,,,4,,8,60lesize,all,,,2mshape,0mshkey,1mshmid,0!mshmid,1!mshmid,2amesh,all!EX3.5 线上单元尺寸设置示例!下边密上边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10lsel,s,loc,x,0lesize,all,,,9,1/8lsel,s,loc,x,10lesize,all,,,9,8lsel,allmshape,0mshkey,1amesh,all!中间密外边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10,-1/5lsel,s,tan1,xlesize,all,,,9,-1/8lsel,allmshape,0mshkey,1amesh,all!EX3.6 合并线和连接线以进行映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,3mshape,0mshkey,1lccat,1,2lcomb,4,5amesh,all!EX3.7 合并线或连接线的网格划分设置finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,10lesize,6,,,8lesize,1,,,4lesize,2,,,3lesize,4,,,5lesize,5,,,2mshape,0mshkey,1lcomb,1,2!lccat,1,2!如采用lcatt,1,2则采用7个划分数lcomb,4,5amesh,all!EX3.7B 简化映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,6mshape,0amap,1,2,5,3,4!EX3.8A 过渡四边形映射网格finish/clear/prep7k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,8lesize,3,,,3lesize,4,,,7lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.8B 过渡四边形映射网格finish/clear/prep7et,1,plane42k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,11lesize,3,,,3lesize,4,,,2lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.9A 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,3lesize,5,,,15allselmshape,0,3dvmesh,all!EX3.9B 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,4lesize,7,,,12mshape,0,3dmshkey,1vmesh,all!EX3.9C 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lsel,s,,,1,8,7lesize,all,,,6lsel,s,,,2,7,5lesize,all,,,3lsel,s,,,4,5lesize,all,,,7lsel,s,,,3,6,3lesize,all,,,2lsel,s,,,9,12,1lesize,all,,,5allselmshape,0,3dmshkey,1vmesh,all!EX3.9D 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,2,,,2lesize,4,,,2lesize,5,,,2lesize,1,,,3lesize,3,,,3lesize,7,,,6lesize,6,,,7lesize,11,,,7lesize,9,,,3lesize,10,,,3lesize,12,,,3allselmshape,0,3dmshkey,1vmesh,all!EX3.10 带两孔的长方体的扫掠网格划分finish/clear/prep7a=10r=2et,1,mesh200,6et,2,solid45blc4,,,2*a,a,acyl4,a/2,a/2,r,,,,awprota,,90cyl4,1.5*a,a/2,r,,,,-avsbv,1,2vsbv,4,3wprota,,,90wpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpcsyswpoff,,a/2,a/2vsbw,allwprota,,90vsbw,allwpcsysesize,1amap,105,15,16,26,63amap,107,16,13,60,26amap,108,13,14,28,60amap,103,15,14,28,63lesize,94,,,4lesize,79,,,5vsel,s,loc,x,0,avsweep,allasel,s,loc,y,aasel,r,loc,x,a,2*alccat,2,45lccat,57,71lccat,68,78lccat,65,72mshape,0,2dmshkey,1amesh,allvsel,s,loc,x,a,2*avsweep,allallsel/view,1,1,2,3Eplot!EX3.11 拉伸类命令生成体单元网格finish/clear/prep7et,1,82et,2,95blc4,,,4,4blc4,6,,4,4blc4,12,,4,4blc4,18,,4,4esize,1amesh,allesize,,8vrotat,1,,,,,,1,4,90vext,2,,,,,10,0.5,0.5voffst,3,10vdrag,4,,,,,,35!EX3.12 圆的网格划分finish/prep7et,1,plane82r0=10cyl4,,,r0cyl4,3*r0,,,,r0 wprota,,90asbw,allwprota,,,90asbw,allwpoff,,,3*r0asbw,allwpcsys,-1asel,s,loc,x,-r0,r0lsla,slesize,all,,,8mshape,0,2dmshkey,1amesh,allallselasel,s,loc,x,2*r0,4*r0 lesize,all,,,8lsel,r,length,,r0lesize,all,,,8,0.1,1 amesh,allallsel!EX3.13 圆环的网格划分finish/clear/prep7et,1,plane82r0=10cyl4,,,r0/3,,r0,90cyl4,2*r0,,r0/10,,r0,90 asel,s,loc,x,-r0,r0lsla,slesize,all,,,8lsel,r,length,,r0*2/3 lesize,all,,,3,,1 mshape,0,2dmshkey,1amesh,allasel,s,loc,x,2*r0,4*r0lesize,5,,,12lesize,7,,,6lsel,s,length,,r0*9/10lesize,all,,,7amesh,all!EX3.14 圆柱面和圆柱体的网格划分!圆柱面finish/clear/prep7r0=10h0=50et,1,shell63cyl4,,,r0adele,1cm,l1cm,linek,50k,51,,,h0l,50,51adrag,l1cm,,,,,,5lsel,s,loc,z,0lesize,all,,,6lsel,s,length,,h0lesize,all,,,8mshape,0,2dmshkey,1amesh,all!圆柱体finish/clear/prep7r0=10h0=50et,1,solid95cyl4,,,r0,,,,h0wprota,,90vsbw,allwprota,,,90vsbw,allmshape,0,3dmshkey,1lsel,s,loc,z,0lsel,s,length,,h0lesize,all,,,8vmesh,all!EX3.15 圆锥的网格划分!扫掠网格划分finish/clear/prep7cone,10,,,15,,90et,1,200,7et,2,95lsel,s,,,5,6lesize,all,,,12,0.5lsel,alllesize,3,,,6mshape,0,2dmshkey,1amesh,3vsweep,1,3,4vsymm,x,allvsymm,y,allvglue,all!六面体映射网格划分finish/clear/prep7cone,10,,,15,,90et,2,95lesize,all,,,8mshape,0,3dmshkey,1vmesh,allvsymm,x,allvsymm,y,allvglue,all!EX3.16 正多边形面的通用网格划分finish/clear/prep7ns=11rpr4,ns,,,10kp0=100+nsk,kp0*if,mod(ns,2),eq,0,then*do,i,1,ns/2l,kp0,2*i-1*enddo*else*do,i,1,nsl,kp0,i*enddo*endiflsel,s,,,ns+1,2*nscm,l1cm,lineallselasbl,1,l1cm*if,mod(ns,2),eq,0,then cmsel,s,l1cmlesize,all,,,8lsel,invelesize,all,,,4*elselesize,all,,,8*endifallselmshape,0mshkey,1amesh,all!EX3.17 球体及球面网格划分finish/clear/prep7r0=10sphere,,r0,,90vsbw,allvdele,2,,,1esize,2et,1,solid95mshape,0,3dmshkey,1vsymm,x,all vsymm,y,all vsymm,z,allvglue,allfinish/clear/prep7r0=10sphere,,r0,,90 vsbw,allvdele,2,,,1vdele,allasel,s,loc,x,0asel,a,loc,y,0asel,a,loc,z,0adele,all,,,1allselesize,2et,1,shell63 mshape,0,2d mshkey,1amesh,allarsym,x,allarsym,y,allarsym,z,all nummrg,all!EX3.18 孔板网格划分finish/clear/prep7a0=300b0=800r0=15blc4,,,a0,b0cyl4,a0/4,b0/8,r0 agen,2,2,,,a0/2 agen,2,2,3,1,,b0/8 agen,2,2,5,1,,b0*5/8 asel,s,,,2,9,1cm,a2cm,areaallselwprota,,-90*do,i,1,5wpoff,,,b0/16asbw,all*enddowpoff,,,b0*5/16*do,i,1,5wpoff,,,b0/16asbw,all*enddowprota,,,90*do,i,1,3wpoff,,,a0/4asbw,all*enddowpcsys,-1numcmp,alllsel,s,radius,,r0lesize,all,,,8lsel,invelesize,all,,,4lsel,allet,1,82mshape,0,2dmshkey,1asel,u,loc,y,b0/16,b0*5/16 asel,u,loc,y,b0*11/16,b0*15/16 lsla,slsel,r,tan1,xlesize,all,50,,,,1amesh,allallselamap,21,30,31,54,62 amap,22,31,32,54,64 amap,32,29,30,62,76 amap,33,29,32,64,76amap,42,34,35,76,82 amap,43,35,36,76,84 amap,41,33,34,53,82 amap,44,33,36,53,84amap,19,22,23,49,64 amap,20,23,24,49,60 amap,9,21,22,64,74 amap,31,21,24,60,74amap,38,26,27,74,84 amap,39,27,28,74,80 amap,37,25,26,50,84 amap,40,25,28,50,80amap,17,14,15,43,59 amap,18,15,16,43,63 amap,7,13,14,59,73 amap,8,13,16,63,73amap,34,18,19,73,79 amap,30,17,18,79,44 amap,35,19,20,73,83 amap,36,17,20,44,83amap,15,6,7,39,63amap,16,7,8,39,58amap,4,6,5,63,71amap,5,5,8,58,71amap,27,10,11,71,83 amap,28,11,12,71,78 amap,26,9,10,83,40 amap,29,9,12,40,78!EX3.19 角支架的网格划分finish/clear/prep7blc4,,,150,50blc4,100,,50,-50cyl4,,25,25cyl4,125,-50,25aadd,allnumcmp,allcyl4,,25,10cyl4,125,-50,10 asel,s,,,2,3cm,a1cm,area asel,allasba,1,a1cm lcomb,1,6lfillt,1,2,20 asbl,4,6adele,1,,,1lfillt,3,4,20al,18,19,20 aadd,all numcmp,all!---------------------- wprota,,90 wpoff,,,-25 asbw,all wpoff,,,75 asbw,all wprota,,,90 asbw,all wpoff,,,125 asbw,all wpcsys,-1 wpoff,25 wprota,,,90 asbw,all kwpave,18 asbw,all kwpave,3 wprota,,90 asbw,all kwpave,21 wprota,,-45 asbw,8 wpcsys,-1!-------------------- et,1,plane82 mshape,0,2d mshkey,1 esize,6lesize,33,,,6 lesize,37,,,6 lesize,42,,,6amesh,3,5,2amap,6,9,10,4,23amap,7,9,12,1,23amesh,11,13,2amesh,1,9,8amesh,2,4,2amap,16,3,18,26,28amap,12,14,15,5,28amap,15,13,14,6,28lcomb,19,22lcomb,27,46amesh,10,14,4!========================= !EX3.20 大板小孔的网格划分A finish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1wprota,,90,45asbw,allet,1,plane82mshape,0,2dmshkey,1lsel,s,loc,x,0,a0/20lesize,all,,,8lsel,alllesize,7,,,12,0.1lesize,8,,,12,0.1lesize,9,,,12,10amesh,allcsys,0arsym,x,allarsym,y,allnummrg,all!============================= !EX3.20 大板小孔的网格划分Bfinish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1et,1,plane82mshape,0,2dmshkey,1lesize,5,,,8lesize,1,,,10lesize,4,,,8,5lesize,6,,,8,5amesh,1lesize,7,,,20,0.1lesize,8,,,20,0.1amap,2,50,51,2,4csys,0arsym,x,allarsym,y,allnummrg,all!EX3.21 Z形面网格划分及局部细分finish/clear/prep7blc4,,,15,10blc4,10,6,14,12aadd,allwprota,,-90wpoff,,,6asbw,allwpoff,,,4wprota,,,90wpoff,,,10asbw,allwpoff,,,5asbw,allwpcsys,-1et,1,82esize,2mshape,0,2dmshkey,1amesh,allkrefine,9,10,1,1,,off !不光滑处理(节点位置不变)aclear,allamesh,allkrefine,9,10,1,1,,smooth !光滑处理,节点位置可改变,但既有单元不删除aclear,allamesh,allkrefine,9,10,1,1 !光滑处理,可删除既有单元重细分,节点位置改变!EX3.22 分布小孔环形面的网格划分finish/clear/prep7r1=15r2=25r3=3n=8cyl4,,,r1,,r2,180/ncyl4,0.5*(r1+r2),,r3asba,1,2kl,2larc,6,8,4,0.5*(r1+r2)asbl,all,4lesize,4,,,4lesize,10,,,4lesize,6,,,8lesize,3,,,6lesize,8,,,4lesize,5,,,8lesize,9,,,4lesize,7,,,4et,1,plane82mshape,0,2damap,2,6,7,4,8amap,1,5,6,8,1!lrefine,5,6,1,1,smootharsym,y,allcsys,1agen,n,all,,,,360/nnummrg,all!EX3.23 同心圆轴体的网格划分finish/clear/prep7cylind,50,0,100,150cylind,40,0,200,250cylind,30,0,100,250cylind,20,0,50,250cylind,10,0,0,250vptn,allwprota,,,90vsbw,allwprota,,90vsbw,allet,1,95csys,1lsel,s,length,,50lesize,all,,,5lsel,invelsel,u,radius,,10,50lesize,all,,,4n1=6n2=6n3=8n4=10n5=12lsel,r,loc,x,0,10lesize,all,,,n1,,1lsel,s,radius,,10lesize,all,,,n1lsel,s,radius,,20lsel,s,radius,,30lesize,all,,,n3lsel,s,radius,,40lesize,all,,,n4lsel,s,radius,,50lesize,all,,,n5allselmshape,0,3dmshkey,1vmesh,all!EX3.24A 圆柱与长方体组合模型的网格划分!采用六面体映射网格划分finish/clear/prep7a0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,allaccat,4,6et,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.24B 圆柱与长方体组合模型的网格划分!采用过渡六面体映射网格划分finish/cleara0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,alllsel,s,radius,,r0lsel,a,length,,r0lesize,all,,,6lsel,s,loc,z,0lsel,a,loc,z,h1lesize,all,,,7asel,s,loc,x,a0/2asel,a,loc,y,a0/2accat,allallselet,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.25-任意圈数圆柱形螺旋弹簧的扫略网格划分FINISH/CLEAR/PREP7!1.定义弹簧参数(同上)D=4C=8DZ=C*DT=DZ/6N=4.7*IF,T,LT,D,THEN*ENDIFTKPD=90!2.创建全部螺旋线CSYS,1TDEG=N*360TDEG1=MOD(TDEG,TKPD) N0=(TDEG-TDEG1)/TKPD+1 *AFUN,DEG*DO,I,1,N0CTA=(I-1)*TKPDZ=T/360*CTAK,I,DZ/2,CTA,Z*ENDDO*IF,TDEG1,LT,1.0E-2,THEN *ELSEN0=N0+1CTA=CTA+TDEG1Z=T/360*CTAK,N0,DZ/2,CTA,Z*ENDIF*DO,I,1,N0-1L,I,I+1*ENDDOCM,L1,LINE!3.在螺旋线端部创建簧丝截面CSYS,0WPOFF,DZ/2 WPROTA,,90CYL4,,,D/2wprota,,90asbw,allwprota,,,90asbw,allcm,a1cm,areaVDRAG,a1cm,,,,,,L1et,1,mesh200,7et,2,solid95cmsel,s,a1cmlsla,slesize,all,,,6amesh,allallselesize,d/2vsweep,all!EX3.26 具坑缺陷圆柱的网格划分finish/clear/prep7l0=50r0=10r1=5r2=1h0=5et,1,solid45cyl4,,,r0,,,,l0wpave,0,r0,l0/2wprota,,90con4,,,r1,r2,h0vsbv,1,2wprota,,90vsbw,allwpoff,,,r1+3vsbw,allwpoff,,,-2*(r1+3)vsbw,allwprota,,,90vsbw,allwpoff,,,r1+1vsbw,allwpoff,,,-2*(r1+1)vsbw,allwpcsys,-1wprota,,90vsbw,allwpoff,,,-r1vsbw,allwpoff,,,2*r1vsbw,allwpcsys,-1numcmp,allksel,s,loc,z,l0/2ksel,r,loc,x,0ksel,u,loc,y,0*get,kp1,kp,0,num,min kp2=kpnext(kp1) ksel,s,loc,z,l0/2+r2*get,kp3,kp,0,num,min allselk,1000,,-r0,l0/2+5l,kp3,1000asel,none*get,l1,line,0,num,max arotat,l1,,,,,,kp1,kp2 cm,a1cm,areaallselvsba,all,a1cm!vsel,s,loc,z,17,33!vsel,r,loc,x,-6,6!aslv,s!aplotaccat,128,61accat,160,138 accat,224,202 accat,183,35accat,31,134accat,151,164 accat,215,231 accat,82,195accat,130,5accat,162,145 accat,229,209 accat,192,89accat,12,123accat,117,155 accat,109,219 accat,14,178mshape,0,3d mshkey,1esize,1vmesh,all!EX3.27长方体开圆柱槽的网格划分finish/clear/prep7blc5,,,20,10,30wpoff,,5,5cyl4,,,8,,,,10wpoff,,,10cyl4,,,4,,,,12vsel,s,,,2,3cm,v1cm,voluallselwpcsys,-1vsbv,1,v1cmvoffst,15,5voffst,16,3voffst,17,15vptn,allwprota,,,90vsbw,allwpcsys,-1wpoff,,,5vsbw,allwpoff,,,10vsbw,allwpoff,,,12vsbw,allet,1,mesh200,7esize,2mshape,0,2dmshkey,1asel,s,loc,z,30lccat,6,74lccat,8,68amesh,allasel,s,loc,z,0lccat,3,73lccat,1,67amesh,allallselet,2,solid95mshape,0,3dmshkey,1vsweep,all!ex4.1累加节点自由度约束值finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,1lmesh,alld,1,ux,1e-2,,,,uyd,1,rotzd,2,alldlistnsel,s,d,u,0,0.1dofsel,s,uydscale,2.5dlistdofsel,alldcum,addd,1,ux,2e-2dlistdcum,ignod,2,uy,1.0dlist!ex4.2 对线施加约束并转换finish/clear/prep7et,1,95blc4,,,10,10,10dl,7,,ux,0.1dl,5,,alldl,11,6,symmdl,10,6,asymdl,6,,symmdllistesize,2vmesh,alldtrandlist!ex4.3 施加集中荷载与节点坐标系finish/clear/prep7et,1,beam4k,1k,2,5k,3,10l,1,2l,2,3local,12,0,,,,90nrotat,alldk,1,allfk,2,fy,-1000esize,1lmesh,allnrotat,alllplotfk,3,fy,1000f,6,fx,-1000sbctraneplot!ex4.4a 3d单元sf加载示例finish/clear/prep7et,1,95blc4,,,10,10,20esize,,4vmesh,allasel,s,loc,y,10sf,all,pres,1000asel,s,loc,z,20nsla,s,1sf,all,pres,1000!ex4.4b 2d单元sf加载示例finish/clear/prep7et,1,82blc4,,,100,200blc4,30,60,40,80asba,1,2wprota,,-90wpoff,,,60asbw,allwpoff,,,80asbw,allwprota,,,90wpoff,,,30asbw,allwpoff,,,40asbw,allwpcsys,-1esize,5amesh,all/psf,pres,norm,2sf,all,pres,100sfdele,all,presnsel,s,loc,x,0sf,all,pres,100nsel,s,loc,x,15,20esln,s,1nsel,r,loc,x,15sf,all,pres,110nsel,s,loc,x,40,60nsel,r,loc,y,10,30esln,s,1sf,all,pres,100lsel,s,loc,x,100nsll,s,1esln,snsel,s,loc,x,95sf,all,pres,-100eplot!ex4.5 节点号及其荷载函数finish/clear/prep7et,1,45blc4,,,10,10,20esize,5vmesh,all*dim,mypres,,100*do,i,1,100mypres(i)=i*10.0*enddosffun,pres,mypres(1)nsel,s,loc,y,10sf,all,pres,10sflist*do,i,1,100mypres(i)=i*50.0*enddonsel,s,loc,z,20sf,all,pres,0allselsflist!ex4.6a 2d平面单元plane82 finish/clear/prep7et,1,82blc4,,,10,50esize,2amesh,allsfe,1,4,pres,,100,50sfe,6,4,pres,,200,40nsel,s,loc,x,10esln,ssfe,all,2,pres,,-100/psf,pres,norm,2,0,1eplot!ex4.6b 3d体单元solid95 finish/clear/prep7et,1,95blc4,,,10,10,30esize,5vmesh,all/psf,pres,norm,2,0,1sfe,22,3,pres,,100sfe,22,6,pres,,100,110,120,130 sfelist!ex4.6c 3d壳单元shell63 finish/clear/prep7et,1,63wprota,,90blc4,,,10,10esize,5amesh,all/psf,pres,norm,2,0,1sfe,1,1,pres,,100sfe,1,3,pres,,100sfe,1,6,pres,,100sfe,4,1,pres,,100,110,120,130!ex4.7 在梁单元上施加荷载finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,,10lmesh,all/pnum,elem,1sfbeam,3,1,pres,50,100 sfbeam,5,1,pres,100 sfbeam,7,1,pres,50,100,,,0.2,0.1sfbeam,9,1,pres,100,,,,0.4,-1sfbeam,3,2,pres,50,100,,,0.2,0.1!ex4.8a 在线上施加面荷载finish/clear/prep7et,1,82blc4,,,10,30esize,5amesh,all/psf,pres,norm,2sfl,4,pres,10,60sfl,2,pres,60sftraneplot!ex4.8b 在线上施加面荷载finish/clear/prep7et,1,shell63wprota,,90blc4,,,10,30esize,5amesh,allsfl,3,pres,100sftraneplot!ex4.9 利用荷载梯度在直角坐标系下的施加方法finish/clear/prep7et,1,82blc4,,,10,60esize,2amesh,all/psf,pres,norm,2sfgrad,pres,,y,0,-5nsel,s,loc,x,0nsel,r,loc,y,0,40sf,all,pres,600sfgrad,pres,,y,30,-20nsel,s,loc,x,10sf,all,pres,0allseleplot!ex4.10利用荷载梯度在柱坐标系下的施加方法finish/clear/prep7csys,1k,1,10,-90k,2,10,90k,3,10,90,30l,1,2l,2,3adrag,1,,,,,,2ldele,2,,,1numcmp,allet,1,63esize,2amesh,allsfgrad,pres,1,y,-90,1sfa,all,2,pres,400sftran!可考察下述命令及其施加方法!①!sfgrad,pres,1,y,270,1sfa,all,2,pres,400sftran!②!local,12,1cscir,12,1sfgrad,pres,12,y,270,1sfa,all,2,pres,400sftran!ex4.11 ldread命令的示例finish/clear/filname,ldtest1/prep7et,1,plane82mp,ex,1,2.1e5mp,prxy,1,0.3blc4,,,10,40esize,2amesh,alld,all,allsfl,3,pres,100/solusolve/post1prrsolfinish/filename,ldtest2/solulsclear,allnsel,s,loc,y,0d,all,allallselldread,reac,,,,,ldtest1,rst solve/post1plnsol,s,y!ex4.12 初应力荷载finish/clear/filname,colu1/prep7et,1,plane82mp,ex,1,2e5mp,nuxy,1,0.3blc4,,,1,10esize,2amesh,allnsel,s,loc,y,0d,all,uyd,1,uxnsel,s,loc,y,10sf,all,pres,-10。