激光焊接ANSYS命令流(脉冲热源)
ansys命令流整理
ansys命令流整理
前处理
节点复制:
NGEN,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE
是一个节点复制命令,
它是将一组节点在现有坐标系统下复制到其它位置。
ITIME: 复制的次数,包含自己本身。
INC: 每次复制节点时节点号码的增加量。
NODE1,NODE2,NINC: 选取要复制的节点,即要对哪些节点进行复制。
DX,DY,DZ: 每次复制时在现有坐标系统下,几何位置的改变量。
SPACE:间距比,是最后一个尺寸和第一个尺寸的比值
单元复制
EGEN, ITIME, NINC, IEL1, IEL2, IEINC, MINC, TINC, RINC, CINC, SINC, DX, DY, DZ
循环:
*do,I(循环控制变量),1(变量初值),100(循环变量结束变量值),1(步长)
循环内部命令流
*enddo(循环命令结束)
节点选择
NSEL, type,item,comp,vmin,vmax,vinc,kabs
例子:nsel,s,loc,x,-0.1,0.1
S---开始选择节点
U---在已选择节点内要删除部分节点
Loc---利用坐标选择节点
定义截面面积
SECTYPE, SECID, Type, Subtype, Name, REFINEKEY
SECID---ID(编号)
节点自动填充
FILL, NODE1, NODE2, NFILL, NSTRT, NINC, ITIME, INC, SPACE。
ANSYS14.0常用命令流
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,,ALLSFL,3,PRES,1000SFL,2,PRES,1000......SOLVE !求解标识FINISH !求解模块结束标识/POST1 !进入通用后处理器标识....../POST26 !进入时间历程后处理器……/EXIT,SAVE !退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE !指定绕轴旋转视图/DIST !说明对视图进行缩放/DEVICE !设置图例的显示,如:风格,字体等/REPLOT !重新显示当前图例/RESET !恢复缺省的图形设置/VIEW !设置观察方向/ZOOM !对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例:K,1,0,0,02.在激活坐标系生成直线LSTR,关键点P1,关键点P2例STR,1,23.在两个关键点之间连线L,关键点P1,关键点P2例,1,2注:此命令会随当前的激活坐标系不同而生成直线或弧线4.由三个关键点生成弧线LARC,关键点P1,关键点P2,关键点PC,半径RAD例ARC,1,3,2,0.05注:关键点PC是用来控制弧线的凹向5.通过圆心半径生成圆弧CIRCLE,关键点圆心,半径RAD,,,,圆弧段数NSEG例:CIRCLE,1,0.05,,,,46.通过关键点生成样条线BSPLIN,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6例:BSPLIN,1,2,3,4,5,67.生成倒角线LFILLT,线NL1,线NL2,倒角半径RAD例FILLT,1,2,0.0058.通过关键点生成面A,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6,P7,P8...例:A,1,2,3,49.通过线生成面AL,线L1,线L2,线L3,线L4,线L5,线L6,线L7,线L8,线L9,线L10例:AL,5,6,7,810.通过线的滑移生成面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,02.通过矩形中心定位点生成面BLC5,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH注:与上条命令的不同就在于矩形的定位点不一样例:BLC5,2.5,1.5,5,3,03.通过在工作平面定义矩形X.Y坐标生成面RECTNG,矩形左边界X坐标X1,矩形右边界X坐标X2,矩形下边界Y坐标Y1,矩形上边界Y坐标Y2 例:RECTNG,0,5,0,3生成圆面4.通过中心定位点生成实心圆面CYL4,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,圆面的内半径RAD1,内圆面旋转角度THETA1,圆面的外半径RAD2,外圆面旋转角度THETA2,圆面的深度DEPTH注:如要实心的圆面则不用RAD2,THETA2,DEPTH例:CYL4,0,0,5,3605.生成扇形圆面命令介绍如上例1实心扇形:CYL4,0,0,5,60例2扇形圆环:CYL4,0,0,5,60,10,60例3整的圆环:CYL4,0,0,5,360,10,360注:同时可通过定义圆面的深度以生成柱体6.通过在工作平面定义起始点生成圆面CYL5,开始点X坐标XEDGE1,开始点Y坐标YEDGE1,结束点X坐标XEDGE2,结束点Y坐标YEDGE2,圆面深度DEPTH例:CYL5,0,0,2,2,7.通过在工作平面定义内外半径和起始角度来生成圆面PCIRC,内半径RAD1,外半径RAD2,起始角度THETA1,结束角度THETA2例CIRC,2,5,30,1808.生成面与面的倒角AFILLT,面1的编号NA1,面2的编号NA2,倒角半径RAD例:AFILLT,2,5,2第五天目标:掌握多边形面和体的生成1.生成多边形面命令:RPR4,多边形的边数NSIDES,中心定位点X坐标XCENTER,中心定位点Y坐标YCENTER,中心定位点距各边顶点的距离RADIUS,多边形旋转角度THETA例:RPR4,4,0,0,0.15,30注:这条命令可通过定义不同的NSIDES生成三边形,四边形,...,八边形2.生成多边形体命令:RPR4,多边形的边数NSIDES,中心定位点X坐标XCENTER,中心定位点Y坐标YCENTER,中心定位点距各边顶点的距离RADIUS,多边形旋转角度THETA,多边形的深度DEPTH例:RPR4,4,0,0,0.15,30,0.1注:多边形体和面命令唯一的不同就在于深度DEPTH的定义到此,关键点,线,面的生成讲解已结束第六天目标:掌握体的生成命令1.通过关键点生成体命令:V,关键点P1,关键点P2, P3, P4, P5, P6, P7, P8例:V,4,5,6,7,15,24,252.通过面生成体命令:VA,面A1,面A2, A3, A4, A5, A6, A7, A8, A9, A10例:VA,3,4,5,8,103.通过长方形角上定位点生成体命令:BLC4该命令前面在讲生成面的时候已作介绍,唯一的不同在于深度DEPTH的定义.4.通过长方形中心定位点生成面命令:BLC55.通过定义长方体起始位置生成体命令:BLOCK,开始点X坐标X1,结束点X坐标X2, Y1, Y2, Z1, Z2例:BLOCK,2,5,0,2,1,36.生成圆柱体基本命令通生成圆形面,不同在于DEPTH的定义基本命令:CYL4基本命令:CYL5基本命令:CYLIND7.生成棱柱基本命令通生成多边形,不同在于DEPTH的定义基本命令:RPR48.通过球心半径生成球体命令:SPH4,球心X坐标XCENTER,球心Y坐标YCENTER,半径RAD1,半径RAD2例:SPH4,1,1,2,59.通过直径上起始点坐标生成球体命令:SPH5,起点X坐标XEDGE1,起点Y坐标YEDGE1,结束点X坐标XEDGE2,结束点Y坐标YEDGE2 例:SPH5,2,5,7,610.在工作平面起点通过半径和转动角度生成球体命令:SPHERE,半径RAD1,半径RAD2,转动角度THETA1,转动角度THETA2例:SPHERE,2,5,0,6011.生成圆锥体命令:CONE,底面半径RBOT,顶面半径RTOP,底面高Z1,顶面高Z2,转动角度THETA1,转动角度THETA2 例:CONE,10,20,0,50,0,180第七天目标:掌握常用的布尔操作命令1.沿法向延伸面生成体命令:VOFFST,面的编号NAREA,面拉伸的长度DIST,关键点增量KINC例:VOFFST,1,2,,2.通过坐标的增量延伸面生成体命令:VEXT,面1的编号NA1,面2的编号NA2,增量NINC,X方向的增量DX,Y方向的增量DY,Z方向的增量DZ, RX, RY, RZ例:VEXT,1,5,1,1,2,2,3.面绕轴旋转生成体命令:VROTAT,面1的编号NA1,面2的编号NA2,NA3, NA4, NA5, NA6,定位轴关键点1编号PAX1,定位轴关键点2编号PAX2,旋转角度ARC,生成体的段数NSEG例:VROTAT,1,2,,,,,4,5,360,44.沿线延伸面生成体命令:VDRAG,面1的编号NA1,面2的编号NA2, NA3, NA4, NA5, NA6,导引线1的编号NLP1,导引线2的编号NLP2, NLP3, NLP4, NLP5, NLP6例:VDRAG,2,3,,,,,8,5.线绕轴旋转生成面命令:AROTAT,线1的编号NL1, NL2, NL3, NL4, NL5, NL6,定位轴关键点1的编号PAX1,定位轴关键点2的编号PAX2,旋转角度ARC,生成面的段数NSEG例:AROTAT,3,4,,,,,6,8,360,46.沿线延伸线生成面命令:ADRAG,线1的编号NL1,NL2, NL3, NL4, NL5, NL6,导引线1的编号NLP1, NLP2, NLP3, NLP4, NLP5, NLP6例:ADRAG,3,,,,,,87.同理可以延伸关键点,相应的命令如下:LROTAT, NK1, NK2, NK3, NK4, NK5, NK6, PAX1, PAX2, ARC, NSEGLDRAG, NK1, NK2, NK3, NK4, NK5, NK6, NL1, NL2, NL3, NL4, NL5, NL6各选项的含义雷同于上.8.延伸一条线命令EXTND,线的编号NL1,定位关键点编号NK1,延伸的距离DIST,原有线是否保留控制项KEEP例EXTND,5,2,1.5,09.布尔操作:加命令COMB,线编号NL1,线编号NL2,是否修改控制项KEEP例COMB,2,5注:对面和体的相应为:VADD,AADD.选项的含义都类似10.布尔操作:粘接和搭接搭接的核心关键字为:OVLAP,随实体的不同略有不同,如:对体为VOVLAP对面为AOVLAP对线为LOVLAP粘接的核心关键字为:GLUE,随实体的不同略有不同,如:对体为VGLUE对面为AGLUE对线为LGLUE但其他的选项的含义是类似的,这里就不再累述.目标:掌握体素的移动,复制,删除,映射一.移动关键点命令:KMODIF,关键点编号NPT,移动后的坐标X,移动后的坐标Y,移动后的坐标Z例:KMODIF,5,0,0,2二.移动复制关键点命令:KGEN,复制次数选项ITIME,起始关键点编号NP1,结束关键点编号NP2,增量NINC,偏移DX,偏移DY,偏移DZ,关键点编号增量KINC,生成节点单元控制项NOELEM,原关键点是否被修改选项IMOVE例:KGEN,2,1,10,1,2,2,2,,,,注:IMOVE选项说明,设置为0时,不修改原关键点,即为复制,设置为1时,修改原关键点,即为移动,从而通过控制IMOVE选项实现移动或复制.三.移动复制线命GEN,ITIME,NL1,NL2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上四.移动复制面命:AGEN,ITIME,NA1,NA2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上五.移动复制体命令:VGEN,ITIME,NV1,NV2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上六.修改面的法向方向命令:ANORM,面的编号ANUM,单元的法向方向是否修改选项NOEFLIP例:ANORM,2七.体素的删除基本的命令为:*DELE组合不同的关键字形成不同的命令如:KDELE,LDELE,ADELE,VDELE基本的命令格式为:*DELE,起始体素编号N*1,结束体素编号N*2,增量NINC,是否删除体素下层的元素选项KSWP如DELE,2,5,1,1八.体素的映射基本的命令为:*SYMM组合不同的关键字形成不同的命令如:KSYMM,LSYMM,ARSYM,VSYMM基本的命令格式为:*SYMM,映射轴选项NCOMP,起始体素编号N*1,结束体素编号N*2,增量NINC,关键点编号增量KINC,NOELEM, IMOVE如:VSYMM,X,1,10,1,,,,ANSYS 的单位(讲得十分透)ANSYS 软件并没有为分析指定系统单位,在结构分析中,可以使用任何一套自封闭的单位制(所谓自封闭是指这些单位量纲之间可以互相推导得出),只要保证输入的所有数据的单位都是正在使用的同一套单位制里的单位即可。
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杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK180 3D有限应变杆EPCDFGB 另可考虑粘弹塑性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命令流(入门必备)
ANSYS命令集/EXIT,Slab,Fname,Ext,Dir Slab=ALL 保存所有资料Slab=NOSA VE所有更改资料不保存Slab=MODEL保存实体模型,有限元模型,负载的资料(系统默认)例:/EXIT,ALL-------------------------------------------------------- /FILNAM,Fname Fname=工作文件名称,不要扩展名例:/FILNAM,Sanpangzi--------------------------------------------------------/SA VE,Fname,Ext,Dir 保存目前所有的Datebase资料,即更新Jobname.db--------------------------------------------------------/RESUME,Fname,Ext,Dir,NOPAR 回到最后SA VE时的Datebase状态--------------------------------------------------------/CLEAR 清除所有Datebase资料-------------------------------------------------------- LOCAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2 定义区域坐标系统KCN 区域坐标系统代号,大于10的任何号码KCS=0,1,20=笛卡儿坐标1=圆柱坐标2=球面坐标XC,YC,ZC 该区域坐标原点与整体坐标原点的关系THXY,THYZ,THZX 该区域坐标与整体坐标XYZ轴的关系例:LOCAL,11,1,1,1,0-------------------------------------------------------- CSYS,0,1,2声明当前坐标系统例:CSYS,0-------------------------------------------------------- /UNITS,LABEL 声明系统分析时所用的单位LABEL=SI (米,千克,秒)LABEL=CGS (厘米,克,秒)LABEL=BFT (英尺)LABEL=BIN (英寸)例:LABEL,SI-------------------------------------------------------- /PREP7进入通用前处理器-------------------------------------------------------- N,NODE,X,Y,Z,THXY,THYZ,THZX 定义节点NODE 节点号码X,Y,Z 节点在当前坐标系中位置例:N,1,2,3,4-------------------------------------------------------- NDELE,NODE1,NODE2,NINC 删除已建立的节点NODE1,NODE2 删除从NODE1到NODE2的节点,如1到100 NINC 间隔号码,1为1到100全删,2为1,3,5 (99)例:NDELE,1,100,2--------------------------------------------------------NPLOT,KNUM 将节点显示在图形窗口中KNUM=0不显示节点号码KNUM=1显示节点号码--------------------------------------------------------NLIST,NODE1,NODE2,NINC将节点资料列在窗口中例:NLIST--------------------------------------------------------NGEN,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE 复制节点ITIME 复制次数,包括本身INC复制时节点号码增量NODE1,NODE2,NINC 要复制的节点DX,DY,DZ 复制出的节点的位置改变量例:NGEN,4,5,1,5,1,1,2,3 将节点1到5复制4次,每次复制X,Y,Z方向分别移动1,2,3单位长度-------------------------------------------------------- FILL,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE 填充节点(默认为均分填充)例:FILL,1,100 在节点1到100之间填充2,3 (99)-------------------------------------------------------- ET,ITYPE,Ename,KOPT1……KOPT6,INOPR 定义元素ITYPE 元素类型编号Ename 所使用元素名称KOPT1-KOPT6 元素特性编码例:ET,1,LINK1 第1类元素为LINK1单元-------------------------------------------------------- MP,Lab,MAT,C0,C1,C2,C3,C4 定义材料特性材料特性为固定值,其值为C0材料特性随温度变化,由C1-C4控制Lab 材料特性类别MAT 对应ET所定义的元素类型编号ITYPELab=EX,EY,EZ 杨氏系数Lab=DENS 密度Lab=PRXY,PRYZ,PRZX 泊松比Lab=GXY,GYZ,GZX 剪力模数例:MP,EX,1,207E9 第一类元素的杨氏系数为207E9-------------------------------------------------------- R,NSET,R1……R6 定义元素类型几何特性NSET 属性组别号码(系统默认值1)R1-R6 所定义元素类型几何特性值例:R,1,1E-4,2.09E-10 ,0.005-------------------------------------------------------- E,I,J,K,L,M,N,O,P定义元素连接方式I-P 定义元素节点顺序的号码例:E,1,2,5,7 四节点元素的节点顺序为1,2,5,7-------------------------------------------------------- EPLOT,KNUM 将元素显示在图形窗口中ENUM=0 不显示元素ENUM=1 显示元素-------------------------------------------------------- ELIST 将元素资料列在窗口中-------------------------------------------------------- EDELE,IEL1,IEL2,INC 删除已建立的元素IEL1,IEL2,INC 欲删除元素的范围例:EDELE,1,10,1-------------------------------------------------------- EGEN,ITIME,NINC,IEL1,IEL2,IEINC,MINC,IINC,RINC,CINC,DX,DY,DZ 复制元素ITIME 复制次数,包括本身NINC复制时节点号码增量IEL1,IEL2,IEINC 欲复制的元素范围DX,DY,DZ 复制出的元素的位置改变量例:EGEN,6,12,1,4,1 将元素1到4复制6次-------------------------------------------------------- /PNUM,Label,KEY 在图形中显示号码Label=NODE,ELEM,KP,LINE,AREA,VOLUKEY=0 不显示号码KEY=1 显示号码例:/PNUM,ELEM,1-------------------------------------------------------- /SOLU 进入解题处理器--------------------------------------------------------ANTYPE,Antype,Status 声明分析类型Antype=STATIC or 0 静态分析(系统默认)Antype=BUCKLE or 1 屈曲分析Antype=MODAL or 2 振动模态分析Antype=HARMIC or 3 调和外力动力系统分析Antype=TRANS or 4 瞬时动力系统分析例:ANTYPE,STATIC-------------------------------------------------------- F,NODE,Lab,V ALUE,V ALUE2,NEND,NINC定义节点上的集中力NODE 节点号码Lab 外力形式Lab=FX,FY,FZ,MX,MY,MZ 结构力学Lab=HEAT 热学的热流量Lab=AMP,CHRG 电学的电流,电荷Lab=FLUX 磁学的磁通量V ALUE 外力大小NODE,NEND,NINC 施力节点范围例:F,1,FY,,-200,5,1 =F,ALL,FY,-200 节点1-5 的Y方向定义集中力-200(注意FY,,表明V ALUE2默认)-------------------------------------------------------- FDELE,NODE,Lab,NEND,NINC 删除节点集中力例:FDELE,1,FY,5,1 = FDELE,ALL-------------------------------------------------------- D,NODE,Lab,V ALUE,V ALUE2,NEND,NINC,Lab2,Lab3,Lab4,Lab5,Lab6 定义节点自由度的限制NODE,NEND,NINC 选取自由度约束节点的范围Lab 相对元素每一个节点受自由度约束的形式结构力学Lab=UX,UY,UZ(直线位移)Lab=ROTX,ROTY,ROTZ(旋转位移)例:D,1,UX,,,5,1 节点1到5 X方向约束D,1,UX,,,5,1,UY 节点1到5 X Y方向约束D,1,ALL,,,5,1 节点1到5 全部自由度约束注意:使用命令前要先定义元素--------------------------------------------------------DDELE,NODE,Lab,NEND,NINC 取消节点自由度约束例:DDELE,ALL-------------------------------------------------------- DLIST,NODE1,NODE2,NINC 列出节点自由度约束--------------------------------------------------------DL,LINE,AREA,Lab,V ALUE1,V ALUE2定义线自由度限制LINE 线号AREA 线所属面积例:DL,8,3,ALL 定义面积3上面线8的约束注意:同时有DLLIST,DLDELE命令--------------------------------------------------------SFBEAM,ELEM,LKEY,Lab,V ALI,V ALJ,V AL2I,V AL2J,IOFFST,JOFFST定义分布力作用于梁元素的方式及大小ELEM 分布力所作用的元素编号LKEY 梁元素的4个面中分布力所在面号码Lab=PRES (表示分布压力)V ALI,V ALJ 在I,J点分布力的值例:SFBEAM,1,1,PRES,60,30 元素1上1号面作用分布力SFBEAM,1,1,PRES,-30,60--------------------------------------------------------SFE,ELEM,LKEY,Lab,KV AL,V AL1,V AL2,V AL3,V AL4定义分布力作用于元素上的方式及大小ELEM 分布力所作用的元素编号LKEY 分布力作用边,面的号码Lab=PRES (表示分布压力)V AL1-V AL4 分布力在元素边,面上节点的值例:SFE,4,2,PRES,,20,60元素4的第2边,面作用分布力--------------------------------------------------------SF,NLIST,Lab,V ALUE,V ALUE2 定义节点间分布力NLIST 分布力作用边或面上所有节点Lab=PRES例:SF,ALL,PRES,10注意:SFE适用于非均匀分布力,作用在元素的边或面上SF适用于均匀分布力,作用在节点之间-------------------------------------------------------- SFLIST,NODE,Lab 显示分布力-------------------------------------------------------- NSEL,Type,Item,COMP,VMIN,VMAX,VINC,KABS节点选择命令Type 选择方式Type=S 选择一组节点为ACTIVE点Type=R 在现有ACTIVE点中选出一部分作为ACTIVE点Type=A 在选择一部分节点,加入ACTIVE点中Type=U 在现有ACTIVE点中,排除某些节点Type=ALL 选择全部节点作为ACTIVE点Item =NODE 用节点号码选取┇Item =LOC 用节点坐标选取COMP=无┇COMP=X 以节点X坐标为准VMIN,VMAX,VINC 节点选取范围例:NSEL,S,NODE,,1,13,1选1-13中奇数点为ACTIVE点NSEL,A,NODE,,14,20,1 选14-20加入ACTIVE点--------------------------------------------------------OUTPR,Item,FREQ,Cname 分析结果是否显示于输出窗口Item=ALL 所有结果Item=NSOL 节点自由度结果FREQ 负载的次数FREQ=ALL 最后负载例:OUTPR,ALL,ALL注意:仅用于小程序--------------------------------------------------------SOLVE 开始解题--------------------------------------------------------/POST1 进入后处理器--------------------------------------------------------PLDISP,KUND 显示结构变形结构KUND=0 显示变形后结构形状KUND=1 显示变形前后结构形状KUND=0 显示变形前后结构形状,但仅显示外观--------------------------------------------------------PLESOL,Item,Comp显示元素的解答Item(何种解答)Comp(Item分量)S X,Y,Z,XY,YZ,ZX 应力S 1,2,3 主应力S EQV,INT 等效应力F X,Y,Z 结构力M X,Y,Z 结构力矩例:PLESOL,S,X,Y,Z--------------------------------------------------------PLNSOL,Item,Comp显示节点的解答Item(何种解答)Comp(Item分量)S X,Y,Z,XY,YZ,ZX 应力S 1,2,3 主应力S EQV,INT 等效应力F X,Y,Z 结构力M X,Y,Z 结构力矩U X,Y,Z,SUM 位移ROT X,Y,Z,SUM 旋转位移例:PLNSOL,S,Y PLNSOL,U,X--------------------------------------------------------PRESOL,Item,Comp 打印元素解答Item(何种解答)Comp(Item分量)S X,Y,Z,XY,YZ,ZX 应力F X,Y,Z 结构力M X,Y,Z 结构力矩例:PRESOL,S,X--------------------------------------------------------PRNSOL,Item,Comp 打印节点解答Item(何种解答)Comp(Item分量)U X,Y,Z 位移U X,Y,Z方向及总向量方向位移S COMP 应力S PRIN 主应立,等效应力例:PRNSOL,U PRNSOL,S注意:查看结果通常使用PLDISP,1 PLNSOL,U,Y……PRNSOL,S其中PLNSOL中Comp不能省略,PRNSOL中可以省略--------------------------------------------------------TYPE,Itype 声明建立元素时,元素形式号码(对应ET的Itype)例:ET,1,LINK1 ET,2,PLANE42TYPE,1建立LINK1元素TYPE,2建立PLANE42元素--------------------------------------------------------REAL,NSET 声明建立元素时,元素几何参数属性编号(对应R,NSET)--------------------------------------------------------MAT,MAT 声明建立元素时,元素材料特性属性编号(对应MP,MAT)--------------------------------------------------------LSWRITE,LSNUM 多重负载资料保存至文件Jobname.S0i 例:LSWRITE 自动编号,不用输入--------------------------------------------------------LSSOLVE,SLMIN,LSMX,LSINC读取所定义的多重负载并解答SLMIN,LSMX,LSINC 读取负载的范围例:LSSOLVE,1,4 获得1-4负载的解答--------------------------------------------------------DDELE,NODE,Lab,NEND,NINC删除定义的约束条件NODE,NEND,NINC 删除约束的节点范围Lab删除约束的节点方向例:DDELE,1,UY,5,2 删除1,3,5节点的Y向约束--------------------------------------------------------FDELE,NODE,Lab,NEND,NINC删除定义的集中力NODE,NEND,NINC 删除集中力的范围Lab删除外力的方向例:FDELE,1,FY,5,2 删除1,3,5节点的Y向集中力--------------------------------------------------------SFDELE,Nlist,Lab 将已定义的面载荷删除Nlist 面负载所含的节点(由NSEL选择,设Nlist=ALL)Lab=PRES (结构力学)--------------------------------------------------------SFEDELE,ELEM,LKEY,Lab 将负载从元素上删除ELEM 元素编号LKEY 负载作用边,面的号码Lab=PRES (结构力学)--------------------------------------------------------SET,Lstep,SBSTEP,FACT,KIMG,TIME,NGLE,NSET检查负载结果Lstep=负载编号例:SET,2 检查第二负载的结果PLDISP,1……注意:此命令要在/POST1中使用-------------------------------------------------------- FILE,Fname,Ext,Dir 读取分析后的结果文件例:FILE,TEST,RST--------------------------------------------------------ANTYPE,Antype,Status 声明分析类型例:ANTYPE=MODAL or 2 模态分析-------------------------------------------------------- MODOPT,Method,NMODE,FREQB,FREQE,PRMODE,NUMKEY 选择模态分析方法Method=REDUC 降阶法Method=SUBSP 次空间法NMODE 欲求模态个数(降阶法小于主自由度一半)FREQB,FREQE 欲探讨振动频率范围(默认全部)PRMODE 分析后,模态结果保存到结果文件的个数例:MODOPT,SUBSP,5-------------------------------------------------------- M,NODE,Lab1,NEND,NINC,Lab2。
激光焊接ansys小例子
nsel,s,loc,x,x_center-0.002,x_center+0.002
nsel,r,loc,y,y_center,y_center+0.002
nsel,r,loc,z,0
sf,all,hflux,%laser% !在选择的节点上面施加定义好的热流密度函数
allsel,all
nsel,r,loc,z,-0.01,0
ine,all,3 !对激光光斑作用的区域进行网格细化
save
x_center=x0+rb !设置光斑的初始位置
y_center=y0
然后通过gui方式调用函数编辑器定义要加载的热流密度函数,函数的形式为
QC*exp(-2*(({X}-x_center)^2+({Y}-y_center)^2)/(rb^2)),保存后退出,利用函数加载器
nsubst,5
solve
*enddo
不好选择,所以采用圆形区域的外接正方形区域来近似.
再输入第二部分命令流
/sol
antype,4,new !定义分析类型
trnopt,full
toffst,0 !设置采用的温度单位为K
tunif,300 !设置初始温度为300K
!---------------------------
*set,tim,0
QC=Q/(PI* rb*rb)
et,1,solid70 ! 定义单元类型
mp,dens,1,7833 !定义密度
mp,c,1,465 !定义比热
mp,kxx,1,54 !定义热传导系数
block,0,0.052,0,0.025,-0.025,0 !建长方体模型
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命令流(脉冲热源)
/CONFIG,NRES,10000/FILNAME,Thermal3,1/TITLE,Laser weldding/UNITS,SI/prep7LENGTH=0.050 !焊缝长度WIDTH=0.050 !焊接宽度的一半Rlaser=0.0003 !激光最小半径Vlaser=0.3/60 !激光焊接速度LSIZE=0.00025 !划分最小单元尺寸TT=LENGTH/Vlaser !总的焊接时间F=30 !激光脉冲频率TINC=1/F !激光作用周期L_ZONE=0.001 !网格细化分宽度QMAX=35.38e8 !激光最大输出能量et,1,shell57 !采用shell57单元R,1,0.0005 !厚度!定义材料属性!导热系数MPTEMP,,,,,,,,MPTEMP,1,293MPTEMP,2,373MPTEMP,3,473MPTEMP,4,573MPTEMP,5,673MPTEMP,6,773MPTEMP,7,873MPTEMP,8,973MPTEMP,9,1073MPTEMP,10,1173MPTEMP,11,1273MPTEMP,12,1373MPTEMP,13,1500MPTEMP,14,1600MPTEMP,15,1700MPTEMP,16,1800MPTEMP,17,1900MPTEMP,18,2000MPTEMP,19,2200MPTEMP,20,2400MPTEMP,22,2600 MPTEMP,23,2700 MPTEMP,24,2800 MPTEMP,25,2900 MPTEMP,26,3000 MPDATA,KXX,1,,10.8 MPDATA,KXX,1,,11.9 MPDATA,KXX,1,,12.9 MPDATA,KXX,1,,14.2 MPDATA,KXX,1,,15.8 MPDATA,KXX,1,,17.6 MPDATA,KXX,1,,19.5 MPDATA,KXX,1,,21.3 MPDATA,KXX,1,,23.2 MPDATA,KXX,1,,24.8 MPDATA,KXX,1,,26.3 MPDATA,KXX,1,,27.4 MPDATA,KXX,1,,28.2 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1!定义密度MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,8890!定义比热容MPTEMP,,,,,,,, MPTEMP,1,293 MPTEMP,2,373 MPTEMP,3,473MPTEMP,5,673 MPTEMP,6,773 MPTEMP,7,873 MPTEMP,8,973 MPTEMP,9,1073 MPTEMP,10,1173 MPTEMP,11,1273 MPTEMP,12,1373 MPTEMP,13,1500 MPTEMP,14,1600 MPTEMP,15,1700 MPTEMP,16,1800 MPTEMP,17,1900 MPTEMP,18,2000 MPTEMP,19,2200 MPTEMP,20,2400 MPTEMP,21,2500 MPTEMP,22,2600 MPTEMP,23,2700 MPTEMP,24,2800 MPTEMP,25,2900 MPTEMP,26,3000MPDATA,C,1,,423 MPDATA,C,1,,454 MPDATA,C,1,,472 MPDATA,C,1,,489 MPDATA,C,1,,505 MPDATA,C,1,,522 MPDATA,C,1,,538 MPDATA,C,1,,553 MPDATA,C,1,,569 MPDATA,C,1,,584 MPDATA,C,1,,598 MPDATA,C,1,,613 MPDATA,C,1,,830 MPDATA,C,1,,1650 MPDATA,C,1,,660 MPDATA,C,1,,670 MPDATA,C,1,,680 MPDATA,C,1,,688 MPDATA,C,1,,695 MPDATA,C,1,,699MPDATA,C,1,,705MPDATA,C,1,,705MPDATA,C,1,,705MPDATA,C,1,,705MPDATA,C,1,,705!建立模型S_DIS=0.0005RECTNG,S_DIS,0.050+S_DIS,0,0.0008 RECTNG,S_DIS,0.050+S_DIS,0.0008,0.0018 RECTNG,S_DIS,0.050+S_DIS,0.0018,0.005 RECTNG,S_DIS,0.050+S_DIS,0.005,0.050 aglue,allnumcmp,all!划分网格!------------------------ 面4LESIZE,12,,,20LESIZE,13,,,20LESIZE,7,,,20MSHKEY,1MSHAPE,0,2DAMESH,4!------------------------- 面3LESIZE,10,,,2LESIZE,11,,,2MSHKEY,1MSHAPE,0,2DAMESH,3LREFINE,5,,,1,1,OFF!------------------------- 面2LESIZE,8,,,2LESIZE,9,,,2MSHKEY,1MSHAPE,0,2DAMESH,2LREFINE,3,,,1,1,OFF!------------------------ 面1LESIZE,2,,,4,LESIZE,4,,,4,MSHKEY,1MSHAPE,0,2DAMESH,1ALLSEL,ALLNUMCMP,ALLALLSEL,ALL/REPLOTFINISH/SOLUANTYPE,TRANSTRNOPT,FULLTUNIF,293nlgeom,offnropt,fullKK=3N_ONCE=TT/TINCMAX_TIME=N_ONCE+1MAX_X=(LENGTH+2*S_DIS)/LSIZE+1MAX_Y=2*L_ZONE/LSIZE+1H_MAX_Y=L_ZONE/LSIZE+1!------------------------定义载荷并求解-------------------------------SFA,2,,CONV,400,293SFA,3,,CONV,250,293SFA,4,,CONV,50,293*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME*DO,K,1,MAX_TIME,1*DO,I,1,MAX_X,1FLUX2(I,0,K)=(I-1)*LSIZE*ENDDO*DO,J,1,MAX_Y,1FLUX2(0,J,K)=-(H_MAX_Y-J)*LSIZE*ENDDO*DO,I,1,MAX_X,1*DO,J,1,MAX_Y,1XCENTER=(K-1)*Vlaser/FYCENTER=0XCOORD=(I-1)*LSIZEYCOORD=-(H_MAX_Y-J)*LSIZEDISTANCE=SQRT((XCOORD-XCENTER)**2+(YCOORD-YCENTER)**2)*IF,DISTANCE,LE,Rlaser,THENFLUX2(I,J,K)=QMAX/EXP(KK*(DISTANCE**2)/(Rlaser**2))*ELSEFLUX2(I,J,K)=0*ENDIF*ENDDO*ENDDOFLUX2(0,0,K)=(K-1)*TINCOUTRES,BASIC,ALLKBC,1TSRES,ERASETIMINT,ONLNSRCH,ONTM1=0.005+(K-1)*TINCTIME,TM1SFA,1,,HFLUX,%FLUX2%AUTOTS,-1NSUBST,4/GST,ON/REP,FAST/PSF,HFLUX,,1 !表面上热流加载用带颜色填充曲面表示/REPLOTSOLVESFADELE,1,,HFLUXTM2=K*TINCTIME,TM2SFA,1,,CONV,400,293AUTOTS,-1NSUBST,4/GST,ONSOLVESFADELE,1,,CONVSAVE*ENDDOSAVE!***********************************!冷却阶段!************************************DO,k,1,30,1TIME,TT+5*k !载荷步结束时间NSUBST,25KBC,1SOLVESAVE*ENDDOoutres,all,allfinish/CONFIG,NRES,10000/FILNAME,stress3,1/TITLE,Laser weldding/UNITS,SI/PREP7et,1,solid45 !把热单元转化成结构单元R,1,0.0005!------------------------设置结构分析的材料属性----------------------- MPTEMP,,,,,,,,MPTEMP,1,293MPTEMP,2,373MPTEMP,3,473MPTEMP,4,513MPTEMP,5,573MPTEMP,6,673MPTEMP,7,773MPTEMP,8,873MPTEMP,9,973MPTEMP,10,1073MPTEMP,11,1173MPTEMP,12,1273MPTEMP,13,1373MPTEMP,14,1423MPTEMP,15,1473MPTEMP,16,1596MPTEMP,17,1673MPTEMP,18,1873MPTEMP,19,2073MPTEMP,20,2273MPTEMP,21,2473MPTEMP,22,2673MPTEMP,23,2873MPTEMP,24,3073MPTEMP,25,3273MPDATA,ALPX,1,,0.184E-5MPDATA,ALPX,1,,1.17E-5MPDATA,ALPX,1,,1.21E-5MPDATA,ALPX,1,,1.28E-5MPDATA,ALPX,1,,1.31E-5MPDATA,ALPX,1,,1.35E-5MPDATA,ALPX,1,,1.40E-5MPDATA,ALPX,1,,1.47E-5MPDATA,ALPX,1,,1.55E-5 MPDATA,ALPX,1,,1.60E-5 MPDATA,ALPX,1,,1.65E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPTEMP,,,,,,,, MPTEMP,1,293 MPTEMP,2,373 MPTEMP,3,473 MPTEMP,4,573 MPTEMP,5,673 MPTEMP,6,773 MPTEMP,7,873 MPTEMP,8,973 MPTEMP,9,1073 MPTEMP,10,1173 MPTEMP,11,1273 MPTEMP,12,1373 MPTEMP,13,1473 MPTEMP,14,1596 MPTEMP,15,1673 MPTEMP,16,1873 MPTEMP,17,2073 MPTEMP,18,2273 MPTEMP,19,2473 MPTEMP,20,2673 MPTEMP,21,2873 MPTEMP,22,3073 MPTEMP,23,3273 MPDATA,EX,1,,2.05E11 MPDATA,EX,1,,2.03E11 MPDATA,EX,1,,1.98E11MPDATA,EX,1,,1.92E11 MPDATA,EX,1,,1.86E11 MPDATA,EX,1,,1.80E11 MPDATA,EX,1,,1.78E11 MPDATA,EX,1,,1.67E11 MPDATA,EX,1,,1.59E11 MPDATA,EX,1,,1.50E11 MPDATA,EX,1,,1.41E11 MPDATA,EX,1,,1.2E11 MPDATA,EX,1,,0.8E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPTEMP,,,,,,,,TB,BISO,1,6,2,TBTEMP,298TBDATA,,354e6,5.89e9,,,,TBTEMP,873TBDATA,,212e6,1.7e9,,,,TBTEMP,1073TBDATA,,199e6,1.4e9,,,,TBTEMP,1473TBDATA,,100e6,1e8,,,,TBTEMP,1596TBDATA,,50e6,1e7,,,,TBTEMP,2173TBDATA,,5e6,1e7,,,,savefinish/SOLUantype,transtrnopt,full !指定为完全瞬态分析nropt,full,,on !完全牛顿拉夫逊法自适应下降nlgeom,1 !大变形autots,on !激活时间步自动阶跃timint,on,struct !打开结构瞬态分析tintp,0.005,,,1,0.5,0.2kbc,0CNVTOL,F,,0.005,2,1.0e-6CNVTOL,U,,0.005,2,1.0e-6OUTRES,ERASEOUTRES,NSOL,LASTnsel,s,loc,y,0.015,0.05d,all,all,0allsel,allnsel,s,loc,y,0dsym,symm,y,,allsel,all*SET,ii,1*DO,T,1E-6,1,TINC/2TM1=TINC*(ii-1)/2+1E-6TIME,TM1ldread,temp,,,TM1,,Thermal3,rth nsubst,20solve*SET,ii,ii+1*enddo*SET,aa,1*DO,T,1,30.5,TINC/2TM2=TINC*(aa-1)/2+1TIME,TM2ldread,temp,,,TM2,,Thermal3,rth nsubst,6solve*SET,aa,aa+1*enddo*SET,k,1*DO,T,30.51,55.51,0.1TM3=(k-1)*0.1+30.51TIME,TM3ldread,temp,,,TM3,,Thermal3,rth nsubst,6solve*SET,k,k+1*enddonsel,s,loc,y,0.03d,all,ux,0d,all,uz,0allsel,allnsel,s,loc,y,0dsym,symm,y,,allsel,all*SET,bb,1*DO,T,55.52,57.52,0.02TM4=(bb-1)*0.02+55.52 TIME,TM4ldread,temp,,,TM4,,Thermal3,rth nsubst,30solve*SET,bb,bb+1*enddo*SET,jj,1*DO,T,57.53,180,2TM5=(jj-1)*2+57.53TIME,TM5ldread,temp,,,TM5,,Thermal3,rth nsubst,6solve*SET,jj,jj+1*enddosavefinish。
焊接ANSYS命令流
/CLEAR !清除数据库,新建文件/FILNAME, EXAMPLE30 !改变工作名!开始热分析!建立有限元模型/PREP7 !进入前处理器ET,1,PLANE77 !平面热单元ET,2,SOLID90 !空间热单元V=0.005 !焊接线速度LENGTH=0.05 !钢板沿焊缝方向长度L_NUMBER=LENGTH/V !沿焊缝方向单元段数 MPTEMP,1,0,300,600,1000,1300,2000 !定义材料模型,温度表 MPDATA,DENS,1,1,7859,7770,7659,7600,6750,6750 !密度MPDATA,DENS,2,1,7859,7770,7659,7600,6750,6750MPDATA,C,1,1,450,514,967,1370,2800,2000 !比热MPDATA,C,2,1,450,514,967,1370,2800,2000MPDATA,KXX,1,1,67,53,39,31, 20,20 !传导系数MPDATA,KXX,2,1,67,53,39,31,20,20RECT,0,0.10,0,0.01 !因为对称,创建 1/2实体模型 K,5,0.002K,6,0.002,0.002K,7,0.00662,0.01K,8,0.00339,0.01112K,9,0,0.0115LARC,7,9,8A,1,5,6,7,9AGEN,2,2ASBA,1,2AGLUE, ALLK,20,0,0.00385CIRCLE,20,0.008,,,90CIRCLE,20,0.012,,,90ASBL,3,4ASBL,2,5ESIZE, 0.00075SMRTSIZE,7MSHAPE,1AMESH,3,5,1ESIZE, 0.003AMESH,1 !划分单元EXTOPT,ESIZE,L_NUMBEREXTOPT,ACLEAR,1MA T,1VEXT,3,5,1,,,LENGTH !挤出MA T,2VEXT,1,,,,,LENGTHNUMMRG,NODE !合并两次挤出产生的位置重合的节 点ALLSFINISH!热分析求解!作稳态热分析施加初始温度/SOLU !进入求解器AUTOTS,ON !打开自动时间步长OUTPR,NSOL,ALL !输出控制OUTRES,NSOL,ALLKBC,1 !阶跃载荷ANTYPE,TRANS !瞬态分析TIMINT,OFF !关闭时间积分NSUBST,4 !子步数目D,ALL,TEMP,20 !施加初始温度TIME,0.01 !时间SOLVE !解答!杀死焊缝单元ESEL,S,MAT,,1 !选择焊缝单元EKILL,ALL !杀死单元ESEL,ALLTIME,0.02SOLVE!开始激活焊缝单元,进行瞬态热分析TIMINT,ON,THER !打开时间积分,进行瞬态分析 DDELE,ALL,TEMP !删除前面用 D 命令施加的温度载 荷TIME=0.02 !变量 TIME 的初始值*DIM,QQ,ARRAY,8 !定义数组,存储单元节点热流密 度*DIM,E_NODE,ARRAY,20 !定义数组,存储单元所属节点编 号*GET,ELEM_MAX, ELEM, 0, NUM,MAX !获取最大单元号*GET,ELEM_MIN, ELEM, 0,NUM,MIN !获取最小单元号*GET,NODE_MAX, NODE, 0, NUM,MAX !获取最大节点号*GET,NODE_MIN, NODE, 0,NUM,MIN !获取最小节点号*DIM,NODE_A VE,ARRAY,NODE_MAX !定义数组,存储节点相邻单元数 量*DO,J,NODE_MIN,NODE_MAX !计算节点相邻单元数量,循环开 始ALLS*IF,NSEL(J),EQ,1,THENNSEL,S,,,J !选择节点JESLN,S !选择节点J 相邻单元*GET,NODE_A VE(J), ELEM, 0,COUNT !获取单元数量*ENDIF*ENDDO !计算节点相邻单元数量,循环结 束ALLS*DO,KK,1,3 !计算三层焊缝温度场,循环开始 R=0.02 !电弧有效加热半径*IF,KK,EQ,1,THEN !第一层焊缝的基本参数Q=3200*0.7 !输入热量QMAX=3*Q/3.1415926/R/R !变量XC= 0.11190E02 !第一层焊缝质心X 坐标YC= 0.21451E02 !第一层焊缝质心Y 坐标V_NUM=2 !第一层焊缝的体号TIME_INC=LENGTH/L_NUMBER/V !第一层焊缝时间增量*ENDIF*IF,KK,EQ,2,THEN !第二层焊缝的基本参数Q=4500*0.7QMAX=3*Q/3.1415926/R/RXC= 0.20373E02YC= 0.60454E02V_NUM=1TIME_INC=LENGTH/L_NUMBER/1.5/V*ENDIF*IF,KK,EQ,3,THEN !第三层焊缝的基本参数Q=5800*0.7QMAX=3*Q/3.1415926/R/RXC= 0.29574E02YC= 0.94733E02V_NUM=3TIME_INC=LENGTH/L_NUMBER/2/V*ENDIF*DO,J,1,L_NUMBER !从头至尾激活第 KK 层焊缝,开 始SFEDELE,ALL,ALL,HFLUX !删除以前施加的热流密度VSEL,S,,,V_NUM !选择第 KK 层焊缝所对应的体ESLV,S,1 !选择体上单元即第 KK 层焊缝单 元NSEL,R,LOC,Z,(J1)*LENGTH/L_NUMBER+0.001,J*LENGTH/L_NUMBER0.001!选择第 J 行单元的中间节点 ESLN,R !选择第 KK 层第J 行单元EALIVE,ALL !激活单元ALLSZC=(J0.5)*LENGTH/L_NUMBER !第 KK 层第J 行单元质心 Z坐标*DO,I,ELEM_MIN,ELEM_MAX !计算施加单元热流密度,循环开 始*IF,ESEL(I),EQ,1,THENESEL,S,,,I !选择第 I个单元*DO,K,1,8 !计算第I单元前8个节点热流密度*GET,E_NODE(K),ELEM,I,NODE,K !获取第 I 个单元第 K 个节点的编 号SS=E_NODE(K)DISP=SQRT((NX(SS)XC)*(NX(SS)XC)+(NY(SS)YC)*(NY(SS)YC)+(NZ(SS)ZC)*(NZ(SS) ZC)) !计算节点到热中心点(质心)的距 离EEE=3*DISP*DISP/R/R !变量*IF,EEE,LT,25,THENQQ(K)=QMAX*EXP(EEE)/NODE_A VE(E_NODE(K)) !第 K 个节点的热流密 度*ELSEQQ(K)=0*ENDIF*ENDDO*IF,QQ(1)+QQ(2)+QQ(3)+QQ(4),NE,0,THENSFE,I,1,HFLUX,,QQ(1), QQ(2), QQ(3),QQ(4) !在第 I 个单元第一个面上加热流密 度*ENDIF*IF,QQ(1)+QQ(2)+QQ(5)+QQ(6),NE,0,THENSFE,I,2,HFLUX,,QQ(1), QQ(2), QQ(6),QQ(5) !在第 I 个单元第二个面上加热流密 度*ENDIF*IF,QQ(6)+QQ(2)+QQ(3)+QQ(7),NE,0,THENSFE,I,3,HFLUX,,QQ(2), QQ(3), QQ(7),QQ(6) !在第 I 个单元第三个面上加热流密 度*ENDIF*IF,QQ(1)+QQ(4)+QQ(5)+QQ(7),NE,0,THENSFE,I,5,HFLUX,,QQ(1), QQ(4), QQ(7),QQ(5) !在第 I 个单元第五个面上加热流密 度*ENDIF*IF,QQ(5)+QQ(6)+QQ(7)+QQ(8),NE,0,THENSFE,I,6,HFLUX,,QQ(5), QQ(6), QQ(7),QQ(8) !在第 I 个单元第六个面上加热流密 度*ENDIFALLS*ENDIF*ENDDO !计算施加单元热流密度,循环结 束ALLSASEL,S,LOC,Z,0 !选择边界面ASEL,A,LOC,Z, LENGTHASEL,A,,,20,22,1ASEL,A,,,17SFA,ALL,1,CONV,30,20 !施加对流边界条件ALLSTIME=TIME+TIME_INC !时间变量NSUBST,5,10,3 !步长LNSRCH,ON !打开线性搜索TIME,TIME !设定分析时间SOLVE*ENDDO !从头至尾激活第 KK 层焊缝,结 束TIME=TIME+3600 !空冷 3600sTIME,TIMESFEDELE,ALL,ALL,HFLUX !空冷SOLVE !解答*ENDDO !计算三层焊缝温度场,循环结束 SA VE !保存FINI!开始结构分析!前处理/PREP7ETCHG,TTS !将热单元转换为结构单元ET,3,COMBIN14 !弹簧单元R,1,1000 !弹簧刚度MPTEMP !以下定义材料模型MPTEMP,1, 0,300,600,1000,1300,1400,2000 !温度表MPDATA,EX,1,1,2E11,1.86E11,1.35E11,2E10,2E7, 2E7,2E7 !弹性模量MPDATA,EX,2,1,2E11,1.86E11,1.35E11,2E10,2E7, 2E7,2E7MPDATA,PRXY,1,1,0.33,0.33,0.35,0.36,0.4,0.4 !泊松比MPDATA,PRXY,2,1,0.33,0.33,0.35,0.36,0.4,0.4MPDATA,ALPX,1,1,12E6,12.6E6,13.0E6,13.4E6,13.7E6,14E6,13.7E6 !线膨胀系数 MPDATA,ALPX,1,1,12E6,12.6E6,13.0E6,13.4E6,13.7E6,14E6,13.7E6MP,REFT,1,20 !参考温度MP,REFT,2,20TB,BKIN,1,5TBTEMP,15TBDATA,1,280E6,2E9 !屈服极限,切向模量TBTEMP,800TBDATA,1,10E6,1E7TBTEMP,2500TBDATA,1,10E6,1E6TBCOPY,BKIN,1,2 !复制材料 BKIN 数据!以下创建弹簧单元,目的是防止结构产生较大的刚性位移,导致计算不收敛,由于弹 簧!刚度极小,弹簧力不大,对分析不会产生影响K,50,0.11,0.01K,51,0.11,0.01,LENGTHL,50,3L,51,25TYPE,3REAL,1ESIZE,,1LMESH,49,50,1FINISH!结构分析的求解/SOLUDK,10,ALL !在关键点上施加全约束,不妨碍热变 形DA,8, UX !在对称面上施加垂直该面方向约 束DA,14, UXDA,18, UXANTYPE,TRANS !瞬态分析OUTPR,ALL,ALL !输出控制OUTRES,ALL,ALLKBC,0AUTOTS,ON!杀死焊缝单元TIME,0.02ESEL,S,MAT,,1EKILL,ALL !杀死单元ESEL,ALLNSUBST,1SOLVETIME=0.02!以下顺序激活焊缝单元*DO,KK,1,3*IF,KK,EQ,1,THENV_NUM=2TIME_INC=LENGTH/L_NUMBER/V*ENDIF*IF,KK,EQ,2,THENV_NUM=1TIME_INC=LENGTH/L_NUMBER/1.5/V*ENDIF*IF,KK,EQ,3,THENV_NUM=3TIME_INC=LENGTH/L_NUMBER/2/V*ENDIF*DO,J,1,L_NUMBERVSEL,S,,,V_NUMESLV,S,1NSEL,R,LOC,Z,(J1)*LENGTH/L_NUMBER+0.001,J*LENGTH/L_NUMBER0.001 ESLN,REALIVE,ALLALLSTIME=TIME+TIME_INCBFDELE,ALL,TEMP !删除以前施加的温度载荷 LDREAD,TEMP,,,TIME,,,RTH !从热分析的结果文件读温度载荷 NSUBST,10,50,3LNSRCH,ONTIME,TIMENLGEOM,ON !打开大变形选项SOLVE*ENDDOTIME=TIME+3600BFDELE,ALL,TEMPLDREAD,TEMP,,,TIME,,,RTHTIME,TIMESOLVESA VE*ENDDOFINI!查看结果/POST1 !进入普通后处理器SET, , , , , 1.02 !从结果文件读 1.02s 的结果 PLNSOL, BFE,TEMP, 0,1.0 !温度云图PLNSOL, S,EQV, 0,1.0 !应力云图FINI!0.02s 至 10.02s 为焊接过程,10.02s 至 3610.02s 为冷却过程。
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命令流
令查看。
建立元素前必须先行定义使用者欲选择的元素型号、元素材料特性、元素几何特性等,为了程序的协调性一般在/PREP7后,就定义元素型号及相关资料,只要在建立元素前说明使用哪种元素即可。
相关命令ET,ITYPE,Ename,KOPT1,KOPT2,KOPT3,KOPT4,KOPT5,KOPT6,INOPR 元素类型(Element Type)为机械结构系统的含的元素类型种类,例如桌子可由桌面平面单元各桌脚梁单元构成,故有两个元素类型。
ET命令是由ANSYS 元素库中选择某个元素并定义该结构分析所使用的元素类型号码。
ITYPE:元素类型的号码Ename:ANSYS元素库的名称,即使用者所选择的元素。
KOPT1~KOPT6:元素特性编码。
Menu Paths:Main Menu>Preprocessor Element Type>Add/Edit/DeleteMP,Lab,MAT,C0,C1,C2,C3,C4定义材料的属性(Material Property),材料属性为固定值时,其值为C0,当随温度变化时,由后四个参数控制。
MAT:对应ET所定义的号码(ITYPE),表示该组属性属于ITYPE。
Lab:材料属性类别,任何元素具备何种属性在元素属性表中均有说明。
例如杨氏系数(Lab=EX,EY,EZ),密度(Lab=DENS),泊松比(Lab=NUXY,NUXYZ,NUZX),剪切模数(Lab=GXY,GYZ,GXZ),热膨胀系数(Lab=ALPX,ALPY,ALPZ)等。
Menu paths:Main Menu>Preprocessor>Matial Props>IsotropicR,NSET,R1,R2,R3,R4,R5,R6定义”实常数”,即某一单元的补充几何特征,如梁单元的面积,壳单元的厚度。
所带的的参数必须与元素表的顺序一致。
Menu paths:Main Menu>Preprocessor>Real ConstantsE,I,J,K,L,M,N,O,P定义元素的连接方式,元素表已对该元素连接顺序作出了说明,通常2-D平面元素节点顺序采用顺时针逆时针均可以,但结构中的所有元素并不一定全采用顺时针或逆时针顺序。
(完整版)史上最全的ANSYS命令流大全
《史上最全的ANSYS命令流查询与解释》【1】*************************************************************************************对ansys主要命令的解释1,/PREP7 ! 加载前处理模块2,/CLEAR,NOSTART ! 清除已有的数据, 不读入启动文件的设置(不加载初始化文件)初始化文件是用于记录用户和系统选项设置的文本文件/CLEAR, START !清除系统中的所有数据,读入启动文件的设置/FILENAME, EX10.5 ! 定义工程文件名称/TITLE, EX10.5 SOLID MODEL OF AN AXIAL BEARING ! 指定标题4,F,2,FY,-1000 ! 在2号节点上施加沿着-Y方向大小为1000N的集中力6,FINISH ! 退出模块命令7,/POST1 ! 加载后处理模块8,PLDISP,2 ! 显示结构变形图,参数“2”表示用虚线绘制出原来结构的轮廓9,ETABLE,STRS,LS,1 ! 用轴向应力SAXL的编号”LS,1”定义单元表STRSETABLE, MFORX,SMISC,1 ! 以杆单元的轴力为内容, 建立单元表MFORXETABLE, SAXL, LS, 1 ! 以杆单元的轴向应力为内容, 建立单元表SAXLETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容, 建立单元表EPELAXLETABLE,STRS_ST,LS,1 !以杆件的轴向应力“LS,1”为内容定义单元表STRS_STETABLE, STRS_CO, LS,1 !以杆件的轴向应力“LS,1”定义单元表STRS_COETABLE,STRSX,S,X ! 定义X方向的应力为单元表STRSXETABLE,STRSY,S,Y ! 定义Y方向的应力为单元表STRSY*GET,STRSS_ST,ELEM,STEEL_E, ETAB, STRS_ST !从单元表STRS_ST中提取STEEL_E单元的应力结果,存入变量STRSS_ST;*GET, STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO”从单元表STRS_CO中提取COPPER_E单元的应力结果,存入变量STRSS_CO10 FINISH !退出以前的模块11, /CLEAR, START ! 清除系统中的所有数据,读入启动文件的设置12 /UNITS, SI !申明采用国际单位制14 /NUMBER, 2 !只显示编号, 不使用彩色/NUMBER, 0 ! 显示编号, 并使用彩色15 /SOLU ! 进入求解模块:定义力和位移边界条件,并求解ANTYPE, STATIC ! 申明分析类型是静力分析(STA TIC或者0)OUTPR, BASIC, ALL ! 在输出结果中, 列出所有荷载步的基本计算结果OUTPR,BASIC,ALL !指定输出所有节点的基本数据OUTPR,BASIC,LAST ! 选择基本输出选项,直到最后一个荷载步OUTPR,,1 ! 输出第1个荷载步的基本计算结果OUTPR,BASIC,1 ! 选择第1荷载步的基本输出项目OUTPR,NLOAD,1 ! 指定输出第1荷载步的内容OUTRES,ALL,0 !设置将所有数据不记录到数据库。
ansys激光焊接命令流
/DIST,1,0.924021086472,1
/REP,FAST
SET,LIST,999
SET,NEXT
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
SET,LAST
SET,LAST
/DIST,1,1.08222638492,1
/REP,FAST
!*
/EFACET,1
PLNSOL, TEMP,, 0
SET,LAST
/DIST,1,1.08222638492,1
/REP,FAST
/BATCH
/COM,ANSYS RELEASE 12.1 UP20091102 17:19:45 02/24/2012
/input,menust,tmp,'',,,,,,,,,,,,,,,,1
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1
/REPLOT,RESIZE
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
基于ANSYS的焊接过程模拟分析(含命令流)
基于ANSYS的某焊件两焊缝在顺序焊接过程中分析(含命令流)目录一、前言 (3)二、基本方法 (3)三、条件 (3)四、物理性能表 (3)五、基本数据 (3)六、分析结果 (3)七、附录: (4)一、前言本文是对一焊件两个焊缝的凝固过程的温度场进行分析,焊缝及两钢板的材料为钢。
二、基本方法使用有限元分析方法并且使用ANSYS分析软件。
采用三维八节点热分析SOLID70单元,利用生死单元技术,对两个焊缝连续凝固的过程进行分析,本分析分6步进行,首先杀死一个焊缝的所有单元,进行稳态分析,得到温度的初始条件;进行瞬态分析,分析右焊缝的液固相变的转换过程;进行瞬态分析,分析右焊缝的凝固过程;激活焊缝的所有单元,进行短时间的瞬态分析,得到温度初始条件;进行瞬态分析,分析左焊缝的液固相变的转换过程;进行瞬态分析,分析左焊缝的凝固过程。
分析时,采用英制单位。
三、条件初始条件:焊接件的温度为70℉,焊缝的温度为3000℉。
对流边界条件:对流系数0.00005 Btu/(s·in2·℉),空气温度70℉。
四、物理性能表五、基本数据底板尺寸:2in×1.2in×0.34in上板尺寸:1in×1.2in×0.34in焊角尺寸:R=0.17in六、分析结果根据以上的有限元模型,完成对焊件的稳态分析,焊缝相变分析,焊缝凝固分析后,得到温度分布图,见附录。
七、附录:分析模型:网格化分:温度分布图(2s):温度分布图(1000s):温度分布图(1100s):焊接过程中温度随时间变化曲线图:命令流:FINISH/FILNAME,QMH/PREP7ET,1,SOLID70MP,KXX,1,.5e-3MP,C,1,.2MP,DENS,1,.2833MPTEMP,1,0,2643,2750,2875,3000 MPDATA,ENTH,1,1,0,128.1,163.8,174.2,184.6 MP,KXX,2,.5e-3MP,C,2,.2MP,DENS,2,.2833MP,KXX,3,.5e-3MP,C,3,.2MP,DENS,3,.2833MPTEMP,1,0,2643,2750,2875,3000 MPDATA,ENTH,3,1,0,128.1,163.8,174.2,184.6 BLOCK,-0.17,0.17,0,0.34,0,1.2BLOCK,0.17,0.34,0,0.34,0,1.2BLOCK,0.34,1,0,0.34,0,1.2BLOCK,-0.17,0.17,0.34,0.51,0,1.2 BLOCK,-0.17,0.17,0.51,1.34,0,1.2WPA VE,0.17,0.34,0CYLIND,0.17,0,0,1.2,0,90WPA VE,0,0,0CSYS,0FLST,3,3,6,ORDE,3FITEM,3,2FITEM,3,-3FITEM,3,6VSYMM,X,P51X,,,,0,0VGLUE,ALLVPLOT,ALLV ATT,1,1,1VSEL,S,,,1VSEL,A,,,12,17,1V ATT,2,1,1VSEL,S,,,11V ATT,3,1,1ALLSEL,ALLESIZE,0.05VSWEEP,ALLESEL,S,MA T,,3 TOFFST,460EKILL,ALL ALLSEL,ALL/SOLUANTYPE,TRANS TIMINT,OFFESEL,S,MA T,,1NSLED,ALL,TEMP,3000 NSEL,INVED,ALL,TEMP,70 TIME,1KBC,0ALLSEL,ALLSOLVEDDELE,ALL,TEMP TIMINT,ON TINTP,,,,1TIME,100DELTIME,1,0.5,10 AUTOTS,ONKBC,1OUTRES,ERASE OUTRES,ALL,ALL ASEL,S,EXTASEL,S,EXTASEL,U,LOC,Y,0 SFA,ALL,,CONV,5E-5,70 ALLSEL,ALLSOLVETIME,1000 DELTIME,50,10,100 AUTOTS,ONSOLVEEALIVE,ALL ALLSEL,ALLESEL,S,MA T,,3NSLED,ALL,TEMP,3000 TIME,1001 DELTIME,1,1,1 ALLSEL,ALLSOLVEDDELE,ALL,TEMPTIME,1100DELTIME,1,0.5,10SOLVETIME,2000DELTIME,100,10,200 SOLVE/POST1 SET,,,,,,,1PLNSOL,TEMP,,0 SET,,,,,,,25PLNSOL,TEMP,,0 SET,,,,,,,35PLNSOL,TEMP,,0SET,NEXTPLNSOL,TEMP,,0 SET,,,,,,,58PLNSOL,TEMP,,0 SET,,,,,,,64PLNSOL,TEMP,,0/POST26NSOL,2,4727,TEMP,,TEMP_2 STORE,MERGENSOL,3,4752,TEMP,,TEMP_3 STORE,MERGENSOL,4,4808,TEMP,,TEMP_4 STORE,MERGENSOL,5,4833,TEMP,,TEMP_5 STORE,MERGENSOL,6,4883,TEMP,,TEMP_6 STORE,MERGENSOL,7,4908,TEMP,,TEMP_7 STORE,MERGENSOL,8,5088,TEMP,,TEMP_8 STORE,MERGENSOL,9,5308,TEMP,,TEMP_9 STORE,MERGE/AXLAB,X,TIME/AXLAB,Y,TEMPERATURE /XRANGE,0,2000PLV AR,2,3,4,5,6,7,8,9/EXIT,NOSA V。
ANSYS点焊构件有限元模拟命令流
以下是本人做的一个点焊有限元分析的例子的命令流,注意本人用的是ANSYS,这里面有很多的命令流是废话,本人懒,就不做剪切了,其中重要的地方我都做了解释。
将命令流复制粘贴到ANSYS的命令窗口按键盘上的回车键即可运行。
这个案例用的是刚性短梁模拟焊点,希望对某些同志有一些参考意义。
还要啰嗦一句,分析结果没有用实验验证过,正确与否,本人未知。
Finish 这三句是必须粘贴进去的/clear/prep7ET,1,SHELL181 这两行是单元类型!*ET,2,BEAM188!*R,1,0.0017, , , , , , 钢板厚度RMORE, , , , , , ,!*!*MPTEMP,,,,,,,, 这段是材料参数MPTEMP,1,0MPDATA,EX,1,,2.07E11MPDATA,PRXY,1,,0.3TB,BISO,1,1,2,TBTEMP,0TBDA TA,,300E6,490E6,,,,MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,2.0E11MPDATA,PRXY,2,,0.3TB,BISO,2,1,2,TBTEMP,0TBDA TA,,400E6,1000E6,,,,SECTYPE, 1, BEAM, CSOLID, , 0SECOFFSET, CENTSECDATA,0.001,0,0,0,0,0,0,0,0,0 用于模拟焊点的短梁的半径,焊点的大小可在这里定义,这里是1mmRECTNG,0,0.1,0,0.16, 这段是模型的几何建模,本人是先画钢板,后画/DIST,1,1.08222638492,1 线的,这短线代表了焊点/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/REP,FAST/DIST,1,0.924021086472,1/REP,FASTFLST,3,1,5,ORDE,1 FITEM,3,1AGEN,2,P51X, , ,0,0,0.0015, ,0 /DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/VIEW,1,-1/ANG,1/REP,FAST/VIEW,1,,,1/ANG,1/REP,FASTFLST,3,1,5,ORDE,1 FITEM,3,2AGEN, ,P51X, , ,0.06,0,0, , ,1/REP,FAST/VIEW,1,,1/ANG,1/REP,FAST/VIEW,1,,,1/ANG,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FASTK,9,0.08,0.08,0, 画关键点K,10,0.08,0.08,0.0015,LSTR, 9, 10 连线AESIZE,ALL,0.005, 这段是划分网格MSHKEY,0FLST,5,2,5,ORDE,2FITEM,5,1FITEM,5,-2CM,_Y,AREAASEL, , , ,P51XCM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_Y!*AMESH,_Y1!*CMDELE,_YCMDELE,_Y1CMDELE,_Y2!*CM,_Y,LINELSEL, , , , 9CM,_Y1,LINECMSEL,S,_Y!*!*CMSEL,S,_Y1LATT,2,1,2, , , ,1CMSEL,S,_YCMDELE,_YCMDELE,_Y1!*FLST,5,1,4,ORDE,1FITEM,5,9CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_Y!*LESIZE,_Y1, , ,3, , , , ,1!*MSHKEY,0LMESH, 9/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FASTNUMCMP,NODE 这行是压缩节点/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST!*!*/COM, CONTACT PAIR CREATION - START 这段是创建接触对,用modeling>create>contact pair就可以创建了,注意要把里面的算法换成MPC算法,就可以创建一个MPC接触对了CM,_NODECM,NODECM,_ELEMCM,ELEMCM,_KPCM,KPCM,_LINECM,LINECM,_AREACM,AREACM,_VOLUCM,VOLU/GSA V,cwz,gsav,,tempMP,MU,1,MA T,1MP,EMIS,1,7.88860905221e-031R,3REAL,3ET,3,170R,3,,,1.0,0.1,0,RMORE,,,1.0E20,0.0,1.0,RMORE,0.0,0,1.0,,1.0,0.5RMORE,0,1.0,1.0,0.0,,1.0RMORE,10.0KEYOPT,4,4,0KEYOPT,4,5,0KEYOPT,4,7,0KEYOPT,4,8,0KEYOPT,4,9,0KEYOPT,4,10,2KEYOPT,4,11,0KEYOPT,4,12,5KEYOPT,4,2,2KEYOPT,3,5,0! Generate the target surfaceASEL,S,,,1ASEL,A,,,2CM,_TARGET,AREATYPE,3NSLA,S,1ESLN,S,0ESLL,UESEL,U,ENAME,,188,189NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) ESURFCMSEL,S,_ELEMCM! Generate the contact surfaceLSEL,S,,,9CM,_CONTACT,LINETYPE,4NSLL,S,1ESLN,S,0ESURF*SET,_REALID,3ALLSELESEL,ALLESEL,S,TYPE,,3ESEL,A,TYPE,,4ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1! Reverse target normals ESEL,NONEESEL,A,TYPE,,3 ESEL,R,REAL,,3 ESURF,,REVERSE ESEL,ALLESEL,S,TYPE,,3 ESEL,A,TYPE,,4 ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOT!! Reverse contact normals ESEL,NONEESEL,A,TYPE,,4 ESEL,R,REAL,,3 ESURF,,REVERSE ESEL,ALLESEL,S,TYPE,,3 ESEL,A,TYPE,,4 ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOT!ESEL,ALLESEL,S,TYPE,,3 ESEL,A,TYPE,,4 ESEL,R,REAL,,3 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM CMSEL,S,_VOLUCMCMDEL,_VOLUCM/GRES,cwz,gsavCMDEL,_TARGETCMDEL,_CONTACT/COM, CONTACT PAIR CREATION - END /MREP,EPLOT/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FASTFINISH/SOL!*ANTYPE,0LSEL,S, , , 4NSLL,S,1FLST,2,33,1,ORDE,4FITEM,2,1FITEM,2,54FITEM,2,74FITEM,2,-104!*/GO 这段是加载及约束边界D,P51X, ,0, , , ,ALL, , , , , 这行是对左边界全约束LSEL,S, , , 6NSLL,S,1FLST,2,33,1,ORDE,3FITEM,2,695FITEM,2,715FITEM,2,-746!*/GOD,P51X, ,0, , , ,UZ, , , , , 这行是对最右边缘节点的Z方向约束FLST,2,33,1,ORDE,3FITEM,2,695FITEM,2,715FITEM,2,-746!*/GOD,P51X, ,0.001, , , ,UX, , , , , 这行是对最右边缘的节点的X方向加载一个1mm的位移ALLSEL,ALLSA VE,'5m1','db','D:\MYDOCU~1\ANSYS\' 这行是文件的保存名字/STA TUS,SOLUSOLVE 这行是求解。
ansys命令流操作大全.
ansys——ANSYS命令流(Ⅰ)1. A,P1,P2,…,P17,P18(以点定义面)2. AADD,NA1,NA2,…NA8,NA9(面相加)3. AATT,MAT,REAL,TYPE,ESYS,SECN(指定面的单元属性)【注】ESYS为坐标系统号、SECN为截面类型号。
4. *ABBR,Abbr,String(定义一个缩略词)5. ABBRES,Lab,Fname,Ext(从文件中读取缩略词)6. ABBSAVE,Lab,Fname,Ext(将当前定义的缩略词写入文件)7. ABS,IR,IA,--,--,Name,--,--,FACTA(取绝对值)【注】*************8. ACCAT,NA1,NA2(连接面)9. ACEL,ACEX,ACEY,ACEZ(定义结构的线性加速度)10. ACLEAR,NA1,NA2,NINC(清除面单元网格)11. ADAMS,NMODES,KSTRESS,KSHELL【注】*************12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC【注】*************13. ADD,IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC(变量加运算)14. ADELE,NA1,NA2,NINC,KSWP(删除面)【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。
15. ADRAG,NL1,NL2,…,NL6,NLP1,NLP2,…,NLP6(将既有线沿一定路径拖拉成面)16. AESIZE,ANUM,SIZE(指定面上划分单元大小)17. AFILLT,NA1,NA1,RAD(两面之间生成倒角面)18. AFSURF,SAREA,TLINE(在既有面单元上生成重叠的表面单元)19. *AFUN, Lab(指定参数表达式中角度单位)20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE(复制面)21. AGLUE,NA1,NA2,…,NA8,NA9(面间相互粘接)22. AINA,NA1,NA2,…,NA8,NA9(被选面的交集)23. AINP,NA1,NA2,…,NA8,NA9(面集两两相交)24. AINV,NA,NV(面体相交)25. AL,L1,L2,…,L9,L10(以线定义面)26. ALIST,NA1,NA2,NINC,Lab(列表显示面的信息)【注】Lab=HPT时,显示面上硬点信息,默认为空。
ANSYSMechanical在焊接仿真中的应用
ANSYSMechanical在焊接仿真中的应用ANSYSMechanical在焊接仿真中的应用1前言焊接作为现代制造业必不可少的工艺,在材料加工领域一直占有重要地位。
焊接是一个涉及到电弧物理、传热、冶金和力学等各学科的复杂过程,其涉及到的传热过程、金属的融化和凝固、冷却时的相变、焊接应力和变形等是企业制造部门和设计人员关心的重点问题。
焊接过程中产生的焊接应力和变形,不仅影响焊接结构的制造过程,而且还影响焊接结构的使用性能。
这些缺陷的产生主要是焊接时不合理的热过程引起的。
由于高能量的集中的瞬时热输入,在焊接过程中和焊后将产生相当大的残余应力和变形,影响结构的加工精度和尺寸的稳定性。
因此对于焊接温度场合应力场的定量分析、预测有重要意义。
传统的焊接温度场和应力测试依赖于设计人员的经验或基于统计基础的半经验公式,但此类方法带有明显的局限性,对于新工艺无法做到前瞻性的预测,从而导致实验成本急剧增加,因此针对焊接采用数值模拟的方式体现出了巨大优势。
ANSYS作为世界知名的通用结构分析软件,提供了完整的分析功能,完备的材料本构关系,为焊接仿真提供了技术保障。
文中以ANSYS为平台,阐述了焊接温度场仿真和热变形、应力仿真的基本理论和仿真流程,为企业设计人员提供了一定的参考。
2焊接数值模拟理论基础焊接问题中的温度场和应力变形等最终可以归结为求解微分方程组,对于该类方程求解的方式通常为两大类:解析法和数值法。
由于只有在做了大量简化假设,并且问题较为简单的情况下,才可能用解析法得到方程解,因此对于焊接问题的模拟通常采用数值方法。
在焊接分析中,常用的数值方法包括:差分法、有限元法、数值积分法、蒙特卡洛法。
差分法:差分法通过把微分方程转换为差分方程来进行求解。
对于规则的几何特性和均匀的材料特性问题,编程简单,收敛性好。
但该方法往往仅局限于规则的差分网格(正方形、矩形、三角形等),同时差分法只考虑节点的作用,而不考虑节点间单元的贡献,常常用来进行焊接热传导、氢扩散等问题的研究。
ANSYS命令流解释大全
ANSYS命令流解释大全A N S Y S命令流解释大全Document serial number【LGGKGB-LGG98YT-LGGT8CB-LGUT-一、定义材料号及特性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=,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,TB,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热分析教程及命令流算例
=====【热力耦合分析单元简介】======SOLID5-三维耦合场实体具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。
本单元由8个节点定义,每个节点有6个自由度。
在静态磁场分析中,可以使用标量势公式(对于简化的RSP,微分的DSP,通用的GSP)。
在结构和压电分析中,具有大变形的应力钢化功能。
与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。
INFIN9-二维无限边界用于模拟一个二维无界问题的开放边界。
具有两个节点,每个节点上带有磁向量势或温度自由度。
所依附的单元类型可以为PLANE13和PLANE53磁单元,或PLANE55和PLANE77和PLANE35热单元。
使用磁自由度(AZ)时,分析可以是线性的也可以是非线性的,静态的或动态的。
使用热自由度时,只能进行线性稳态分析。
PLANE13-二维耦合场实体具有二维磁场、温度场、电场和结构场之间有限耦合的功能。
由4个节点定义,每个节点可达到4个自由度。
具有非线性磁场功能,可用于模拟B-H曲线和永久磁铁去磁曲线。
具有大变形和应力钢化功能。
当用于纯结构分析时,具有大变形功能,相似的耦合场单元有SOLID5、SOLID98和SOLID62。
LINK31-辐射线单元用于模拟空间两点间辐射热流率的单轴单元。
每个节点有一个自由度。
可用于二维(平面或轴对称)或三维的、稳态的或瞬态的热分析问题。
允许形状因子和面积分别乘以温度的经验公式是有效的。
发射率可与温度相关。
如果包含热辐射单元的模型还需要进行结构分析,辐射单元应当被一个等效的或(空)结构单元所代替。
LINK32-二维传导杆用于两节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
可用于二维(平面或轴对称)稳态或瞬态的热分析问题。
如果包含热传导杆单元的模型还需进行结构分析,该单元可被一个等效的结构单元所代替。
LINK33-三维传导杆用于节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
ansys命令流入门教程
ansys命令流入门教程在 ANSYS 中,荷载包括边界条件和作用力,对结构分析可以是以下内容:位移、力、压力、温度、重力一般可将荷载分为六类,如表 4-1 所示。
★ 荷载即可施加在几何模型(关键点、硬点、线、面、体)上,也可施加在有限元模型(节点、单元)上,或者二者混合使用。
★ 施加在几何模型上的荷载独立于有限元网格,不必为修改网格而重新加载;★ 施加在有限元模型上且要修改网格,则必须先删除荷载再修改网格,然后重新施加荷载。
★ 不管施加到何种模型上,在求解时荷载全部转换(自动或人工)到有限元模型上。
在结构分析中自由度共有7 个,自由度的方向均依从节点坐标系。
约束可施加在节点、关键点、线和面上。
一、施加自由度约束1. 节点自由度约束及相关命令(1) 对节点施加自由度约束命令:D, NODE, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6NODE - 拟施加约束的节点号,其值可取 ALL、组件名。
Lab - 自由度标识符,如UX、ROTZ等。
如为ALL,则为所有适宜的自由度。
VALUE - 自由度约束位移值或表式边界条件的表格名称。
VALUE2 - 约束位移值的第二个数,如为复数输入时,VALUE 为实部,而 VALUE2 为虚部。
NEND,NINC - 节点编号范围和编号增量,缺省时 NEND=NODE,NINC=1。
Lab2,Lab3,Lab4,Lab5,Lab6 - 其它自由度标识符,VALUE 对这些自由度也有效。
各自由度的方向用节点坐标系确定,转角约束用弧度输入例如:D,ALL,ALL ! 对所选节点的全部自由度施加约束D,18,UX,,,,,UY,UZ ! 对节点 18 的 3 个平动自由度全部施加约束D,20,UX,1.0e-4 ! 对节点20 的UX 施加约束,且约束位移值为1.0e-4D,22,UX,0.1,,25,,UY,ROTY ! 对节点 22~25 的 UX,UY,ROTY 施加约束,且位移值均为 0.1(2) 在节点上施加对称和反对称约束命令:DSYM, Lab, Normal, KCNLab - 对称标识,如为 SYMM 则生成对称约束,如为 ASYM 则生成反对称约束。
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=,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,TB,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命令得到更详细的解释。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
/CONFIG,NRES,10000/FILNAME,Thermal3,1/TITLE,Laser weldding/UNITS,SI/prep7LENGTH=0.050 !焊缝长度WIDTH=0.050 !焊接宽度的一半Rlaser=0.0003 !激光最小半径Vlaser=0.3/60 !激光焊接速度LSIZE=0.00025 !划分最小单元尺寸TT=LENGTH/Vlaser !总的焊接时间F=30 !激光脉冲频率TINC=1/F !激光作用周期L_ZONE=0.001 !网格细化分宽度QMAX=35.38e8 !激光最大输出能量et,1,shell57 !采用shell57单元R,1,0.0005 !厚度!定义材料属性!导热系数MPTEMP,,,,,,,,MPTEMP,1,293MPTEMP,2,373MPTEMP,3,473MPTEMP,4,573MPTEMP,5,673MPTEMP,6,773MPTEMP,7,873MPTEMP,8,973MPTEMP,9,1073MPTEMP,10,1173MPTEMP,11,1273MPTEMP,12,1373MPTEMP,13,1500MPTEMP,14,1600MPTEMP,15,1700MPTEMP,16,1800MPTEMP,17,1900MPTEMP,18,2000MPTEMP,19,2200MPTEMP,20,2400MPTEMP,22,2600 MPTEMP,23,2700 MPTEMP,24,2800 MPTEMP,25,2900 MPTEMP,26,3000 MPDATA,KXX,1,,10.8 MPDATA,KXX,1,,11.9 MPDATA,KXX,1,,12.9 MPDATA,KXX,1,,14.2 MPDATA,KXX,1,,15.8 MPDATA,KXX,1,,17.6 MPDATA,KXX,1,,19.5 MPDATA,KXX,1,,21.3 MPDATA,KXX,1,,23.2 MPDATA,KXX,1,,24.8 MPDATA,KXX,1,,26.3 MPDATA,KXX,1,,27.4 MPDATA,KXX,1,,28.2 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1 MPDATA,KXX,1,,41.1!定义密度MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,8890!定义比热容MPTEMP,,,,,,,, MPTEMP,1,293 MPTEMP,2,373 MPTEMP,3,473MPTEMP,5,673 MPTEMP,6,773 MPTEMP,7,873 MPTEMP,8,973 MPTEMP,9,1073 MPTEMP,10,1173 MPTEMP,11,1273 MPTEMP,12,1373 MPTEMP,13,1500 MPTEMP,14,1600 MPTEMP,15,1700 MPTEMP,16,1800 MPTEMP,17,1900 MPTEMP,18,2000 MPTEMP,19,2200 MPTEMP,20,2400 MPTEMP,21,2500 MPTEMP,22,2600 MPTEMP,23,2700 MPTEMP,24,2800 MPTEMP,25,2900 MPTEMP,26,3000MPDATA,C,1,,423 MPDATA,C,1,,454 MPDATA,C,1,,472 MPDATA,C,1,,489 MPDATA,C,1,,505 MPDATA,C,1,,522 MPDATA,C,1,,538 MPDATA,C,1,,553 MPDATA,C,1,,569 MPDATA,C,1,,584 MPDATA,C,1,,598 MPDATA,C,1,,613 MPDATA,C,1,,830 MPDATA,C,1,,1650 MPDATA,C,1,,660 MPDATA,C,1,,670 MPDATA,C,1,,680 MPDATA,C,1,,688 MPDATA,C,1,,695 MPDATA,C,1,,699MPDATA,C,1,,705MPDATA,C,1,,705MPDATA,C,1,,705MPDATA,C,1,,705MPDATA,C,1,,705!建立模型S_DIS=0.0005RECTNG,S_DIS,0.050+S_DIS,0,0.0008 RECTNG,S_DIS,0.050+S_DIS,0.0008,0.0018 RECTNG,S_DIS,0.050+S_DIS,0.0018,0.005 RECTNG,S_DIS,0.050+S_DIS,0.005,0.050 aglue,allnumcmp,all!划分网格!------------------------ 面4LESIZE,12,,,20LESIZE,13,,,20LESIZE,7,,,20MSHKEY,1MSHAPE,0,2DAMESH,4!------------------------- 面3LESIZE,10,,,2LESIZE,11,,,2MSHKEY,1MSHAPE,0,2DAMESH,3LREFINE,5,,,1,1,OFF!------------------------- 面2LESIZE,8,,,2LESIZE,9,,,2MSHKEY,1MSHAPE,0,2DAMESH,2LREFINE,3,,,1,1,OFF!------------------------ 面1LESIZE,2,,,4,LESIZE,4,,,4,MSHKEY,1MSHAPE,0,2DAMESH,1ALLSEL,ALLNUMCMP,ALLALLSEL,ALL/REPLOTFINISH/SOLUANTYPE,TRANSTRNOPT,FULLTUNIF,293nlgeom,offnropt,fullKK=3N_ONCE=TT/TINCMAX_TIME=N_ONCE+1MAX_X=(LENGTH+2*S_DIS)/LSIZE+1MAX_Y=2*L_ZONE/LSIZE+1H_MAX_Y=L_ZONE/LSIZE+1!------------------------定义载荷并求解-------------------------------SFA,2,,CONV,400,293SFA,3,,CONV,250,293SFA,4,,CONV,50,293*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME*DO,K,1,MAX_TIME,1*DO,I,1,MAX_X,1FLUX2(I,0,K)=(I-1)*LSIZE*ENDDO*DO,J,1,MAX_Y,1FLUX2(0,J,K)=-(H_MAX_Y-J)*LSIZE*ENDDO*DO,I,1,MAX_X,1*DO,J,1,MAX_Y,1XCENTER=(K-1)*Vlaser/FYCENTER=0XCOORD=(I-1)*LSIZEYCOORD=-(H_MAX_Y-J)*LSIZEDISTANCE=SQRT((XCOORD-XCENTER)**2+(YCOORD-YCENTER)**2)*IF,DISTANCE,LE,Rlaser,THENFLUX2(I,J,K)=QMAX/EXP(KK*(DISTANCE**2)/(Rlaser**2))*ELSEFLUX2(I,J,K)=0*ENDIF*ENDDO*ENDDOFLUX2(0,0,K)=(K-1)*TINCOUTRES,BASIC,ALLKBC,1TSRES,ERASETIMINT,ONLNSRCH,ONTM1=0.005+(K-1)*TINCTIME,TM1SFA,1,,HFLUX,%FLUX2%AUTOTS,-1NSUBST,4/GST,ON/REP,FAST/PSF,HFLUX,,1 !表面上热流加载用带颜色填充曲面表示/REPLOTSOLVESFADELE,1,,HFLUXTM2=K*TINCTIME,TM2SFA,1,,CONV,400,293AUTOTS,-1NSUBST,4/GST,ONSOLVESFADELE,1,,CONVSAVE*ENDDOSAVE!***********************************!冷却阶段!************************************DO,k,1,30,1TIME,TT+5*k !载荷步结束时间NSUBST,25KBC,1SOLVESAVE*ENDDOoutres,all,allfinish/CONFIG,NRES,10000/FILNAME,stress3,1/TITLE,Laser weldding/UNITS,SI/PREP7et,1,solid45 !把热单元转化成结构单元R,1,0.0005!------------------------设置结构分析的材料属性----------------------- MPTEMP,,,,,,,,MPTEMP,1,293MPTEMP,2,373MPTEMP,3,473MPTEMP,4,513MPTEMP,5,573MPTEMP,6,673MPTEMP,7,773MPTEMP,8,873MPTEMP,9,973MPTEMP,10,1073MPTEMP,11,1173MPTEMP,12,1273MPTEMP,13,1373MPTEMP,14,1423MPTEMP,15,1473MPTEMP,16,1596MPTEMP,17,1673MPTEMP,18,1873MPTEMP,19,2073MPTEMP,20,2273MPTEMP,21,2473MPTEMP,22,2673MPTEMP,23,2873MPTEMP,24,3073MPTEMP,25,3273MPDATA,ALPX,1,,0.184E-5MPDATA,ALPX,1,,1.17E-5MPDATA,ALPX,1,,1.21E-5MPDATA,ALPX,1,,1.28E-5MPDATA,ALPX,1,,1.31E-5MPDATA,ALPX,1,,1.35E-5MPDATA,ALPX,1,,1.40E-5MPDATA,ALPX,1,,1.47E-5MPDATA,ALPX,1,,1.55E-5 MPDATA,ALPX,1,,1.60E-5 MPDATA,ALPX,1,,1.65E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPDATA,ALPX,1,,1.8E-5 MPTEMP,,,,,,,, MPTEMP,1,293 MPTEMP,2,373 MPTEMP,3,473 MPTEMP,4,573 MPTEMP,5,673 MPTEMP,6,773 MPTEMP,7,873 MPTEMP,8,973 MPTEMP,9,1073 MPTEMP,10,1173 MPTEMP,11,1273 MPTEMP,12,1373 MPTEMP,13,1473 MPTEMP,14,1596 MPTEMP,15,1673 MPTEMP,16,1873 MPTEMP,17,2073 MPTEMP,18,2273 MPTEMP,19,2473 MPTEMP,20,2673 MPTEMP,21,2873 MPTEMP,22,3073 MPTEMP,23,3273 MPDATA,EX,1,,2.05E11 MPDATA,EX,1,,2.03E11 MPDATA,EX,1,,1.98E11MPDATA,EX,1,,1.92E11 MPDATA,EX,1,,1.86E11 MPDATA,EX,1,,1.80E11 MPDATA,EX,1,,1.78E11 MPDATA,EX,1,,1.67E11 MPDATA,EX,1,,1.59E11 MPDATA,EX,1,,1.50E11 MPDATA,EX,1,,1.41E11 MPDATA,EX,1,,1.2E11 MPDATA,EX,1,,0.8E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,EX,1,,0.08E11 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.33 MPTEMP,,,,,,,,TB,BISO,1,6,2,TBTEMP,298TBDATA,,354e6,5.89e9,,,,TBTEMP,873TBDATA,,212e6,1.7e9,,,,TBTEMP,1073TBDATA,,199e6,1.4e9,,,,TBTEMP,1473TBDATA,,100e6,1e8,,,,TBTEMP,1596TBDATA,,50e6,1e7,,,,TBTEMP,2173TBDATA,,5e6,1e7,,,,savefinish/SOLUantype,transtrnopt,full !指定为完全瞬态分析nropt,full,,on !完全牛顿拉夫逊法自适应下降nlgeom,1 !大变形autots,on !激活时间步自动阶跃timint,on,struct !打开结构瞬态分析tintp,0.005,,,1,0.5,0.2kbc,0CNVTOL,F,,0.005,2,1.0e-6CNVTOL,U,,0.005,2,1.0e-6OUTRES,ERASEOUTRES,NSOL,LASTnsel,s,loc,y,0.015,0.05d,all,all,0allsel,allnsel,s,loc,y,0dsym,symm,y,,allsel,all*SET,ii,1*DO,T,1E-6,1,TINC/2TM1=TINC*(ii-1)/2+1E-6TIME,TM1ldread,temp,,,TM1,,Thermal3,rth nsubst,20solve*SET,ii,ii+1*enddo*SET,aa,1*DO,T,1,30.5,TINC/2TM2=TINC*(aa-1)/2+1TIME,TM2ldread,temp,,,TM2,,Thermal3,rth nsubst,6solve*SET,aa,aa+1*enddo*SET,k,1*DO,T,30.51,55.51,0.1TM3=(k-1)*0.1+30.51TIME,TM3ldread,temp,,,TM3,,Thermal3,rth nsubst,6solve*SET,k,k+1*enddonsel,s,loc,y,0.03d,all,ux,0d,all,uz,0allsel,allnsel,s,loc,y,0dsym,symm,y,,allsel,all*SET,bb,1*DO,T,55.52,57.52,0.02TM4=(bb-1)*0.02+55.52 TIME,TM4ldread,temp,,,TM4,,Thermal3,rth nsubst,30solve*SET,bb,bb+1*enddo*SET,jj,1*DO,T,57.53,180,2TM5=(jj-1)*2+57.53TIME,TM5ldread,temp,,,TM5,,Thermal3,rth nsubst,6solve*SET,jj,jj+1*enddosavefinish。