激光焊接ANSYS命令流(脉冲热源)

合集下载

焊接模拟ansys实例

焊接模拟ansys实例!下面的命令流进行的是一个简单的二维焊接分析, 利用ANSYS单元生死和热-结构耦合分析功能进!行焊接过程仿真, 计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应力。

finish/clear/filnam,1-2D element birth and death/title,Weld Analysis by "Element Birth and Death"/prep7/unit,si !采用国际单位制!******************************************************et,1,13,4 !13号二维耦合单元, 同时具有温度和位移自由度et,2,13,4!1号材料是钢!2号材料是铝!3号材料是铜!铝是本次分析中的焊料, 它将钢结构部分和铜结构部分焊接起来!下面是在几个温度点下, 各材料的弹性模量mptemp,1,20,500,1000,1500,2000mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11!假设各材料都是双线性随动硬化弹塑性本构关系!下面是各材料在各温度点下的屈服应力和屈服后的弹性模量tb,bkin,1,5tbtemp,20,1tbdata,1,1200e6,0.193e11tbtemp,500,2tbdata,1, 933e6,0.150e11tbtemp,1000,3tbdata,1, 435e6,0.070e11tbtemp,1500,4tbdata,1, 70e6,0.010e11tbtemp,2000,5tbdata,1, 7e6,0.001e11!tb,bkin,2,5tbtemp,20,1tbdata,1,800e6,0.102e11tbtemp,500,2tbdata,1,400e6,0.050e11tbtemp,1000,3tbdata,1, 70e6,0.008e11tbdata,1, 1e6,0.0001e11tbtemp,2000,5tbdata,1,0.1e6,0.00001e11!tb,bkin,3,5tbtemp,20,1tbdata,1,900e6,0.117e11tbtemp,500,2tbdata,1,700e6,0.090e11tbtemp,1000,3tbdata,1,230e6,0.030e11tbtemp,1500,4tbdata,1, 40e6,0.005e11tbtemp,2000,5tbdata,1, 4e6,0.0005e11!!材料密度(假设为常值)mp,dens,1,8030mp,dens,2,4850mp,dens,3,8900! 热膨胀系数(假设为常值)mp,alpx,1,1.78e-5mp,alpx,2,9.36e-6mp,alpx,3,1.66e-5!泊松比(假设为常值)mp,nuxy,1,0.29mp,nuxy,2,0.30mp,nuxy,3,0.30!热传导系数(假设为常值)mp,kxx,1,16.3mp,kxx,2,7.44mp,kxx,3,393!比热(假设为常值)mp,c,1,502mp,c,2,544mp,c,3,385!热膨胀系数(假设为常值)!由于该13号单元还有磁自由度, 此处假设一磁特性, 但并不影响我们所关心的结果mp,murx,1,1mp,murx,2,1mp,murx,3,1!假设焊料(铝)焊上去后的初始温度是1500℃mp,reft,2,1500mp,reft,3,20!****************************************************** !下面建立几何模型csys,0k,1,0,0,0k,2,0.5,0,0k,3,1,0,0 !长1米k,4,0,0.3,0 !厚度0.3米（二维中叫做宽度）k,5,0.35,0.3,0k,6,0.65,0.3,0k,7,1,0.3,0a,1,2,5,4a,2,6,5a,2,3,7,6划分网格esize,0.025type,2mat,2amesh,2!esize,0.05 !网格划分出现问题type,1mat,1amesh,1!mat,3amesh,3eplot!/soluantype,4 ! 瞬态分析trnopt,full在模型的左边界加位移约束nsel,all*get,minx,node,,mnloc,xnsel,s,loc,x,minxd,all,ux,0*get,miny,node,,mnloc,ynsel,r,loc,y,minyd,all,uy,0!*****假设模型的左右边界处温度始终保持在20摄氏度左右*****!其他边界条件如对流和辐射等均可施加，此处因为只是示意而已，故只施加恒温边界条件nsel,all*get,minx,node,,mnloc,xnsel,s,loc,x,minxd,all,temp,20nsel,all*get,maxx,node,,mxloc,xnsel,s,loc,x,maxxd,all,temp,20由于第2个面是焊接所在区域，因此首先将该区域的单元“杀死”nna=2esel,all*get,emax,elem,,num,maxasel,s,area,,nnaesla*get,nse,elem,,count*dim,ne,,nse*dim,nex,,nse*dim,ney,,nse*dim,neorder,,nsemine=0!**********************************************!下面的do循环用于将焊料区的单元按其形心y坐标排序!以便后面模拟焊料由下向上逐步“生长”过程*do,i1,1,nseesel,u,elem,,mine*get,nse1,elem,,countii=0*do,i,1,emax*if,esel(i),eq,1,thenii=ii+1ne(ii)=i*endif*enddo*do,i,1,nse1*get,ney(i),elem,ne(i),cent,y*get,nex(i),elem,ne(i),cent,x*enddominy=1e20minx=1e20*do,i,1,nse1*if,ney(i),lt,miny,thenminy=ney(i)minx=nex(i)mine=ne(i)*else*if,ney(i),eq,miny,then*if,nex(i),lt,minx,thenminy=ney(i)minx=nex(i)mine=ne(i)*endif*endif*endif*enddoneorder(i1)=mine*enddo!************************************************************** max_tem=1500 !按照前面假设，焊料的初始温度为1500℃dt1=1e-3 !用于建立初始条件的一个很小的时间段dt=5 !焊接一个单元所需的时间t=0 !起始时间esel,alleplot/auto,1/replot*do,i,1,nseekill,neorder(i)esel,s,liveeplot*enddoallsel,alloutres,all,allic,all,temp,20kbc,1timint,0,structtimint,1,thermtimint,0,magtintp,0.005,,,1,0.5,0.2!nsub1=2nsub2=40!**************************************************do,i,1,nseealive,neorder(i)esel,s,liveeplotesel,all!******下面的求解用于建立温度的初始条件******t=t+dt1time,tnsubst,1*do,j,1,4d,nelem(neorder(i),j),temp,max_tem*enddosolve!****下面的求解用于保证初始的升温速度为零****t=t+dt1time,tsolve!*********下面的步骤用于求解温度分布***********do,j,1,4ddele,nelem(neorder(i),j),temp*enddot=t+dt-2*dt1time,tnsubst,nsub1solve*enddot=t+50000 !*********下面的步骤用于冷却过程求解***** time,tnsubst,nsub2solvesavefinish后处理过程/post1!**************下面的一系列命令用于生成应力的动画文件******* /seg,dele/cont,1,15,0,1200e6/16,1200e6/dscale,1,1.0avprin,0,0avres,1/seg,multi,stress1,0.1esel,all*do,i,1,nseesel,u,elem,,neorder(i)*enddo*do,i,1,nseesel,a,elem,,neorder(i)set,(i-1)*3+1,1plnsol,s,eqv*do,j,1,nsub1set,(i-1)*3+3,jplnsol,s,eqv*enddo*enddo*do,i,1,nsub2set,(nse-1)*3+4,iplnsol,s,eqv*enddo/seg,off,stress1,0.1/anfile,save,stress1,avi!**********下面的一系列命令用于生成温度的动画文件************ /seg,dele/cont,1,15,0,1500/16,1500/dscale,1,1.0avprin,0,0avres,1/seg,multi,temp1,0.1esel,all*do,i,1,nseesel,u,elem,,neorder(i)*enddo*do,i,1,nseesel,a,elem,,neorder(i)set,(i-1)*3+1,1plnsol,temp*do,j,1,nsub1set,(i-1)*3+3,jplnsol,temp*enddo*enddo*do,i,1,nsub2set,(nse-1)*3+4,iplnsol,temp*enddo/seg,off,temp1,0.1/anfile,save,temp1,avifinish。

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。

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最全命令流解释大全

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小例子

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＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．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的焊接过程模拟分析(含命令流)

基于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命令流查询与解释》【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激光焊接命令流

/REP,FAST
/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，这里面有很多的命令流是废话，本人懒，就不做剪切了，其中重要的地方我都做了解释。

将命令流复制粘贴到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平台焊接模拟中不同焊接热源的比较

图C 面热源稳态温度场分布图/ 面热源加载形式
图?
点热源加载形式
!"#
电弧有效半径
焊接电弧有 01232 中进行焊接热过程模拟时，
图E 体热源加载形式
效半径的大小对模拟结果的精确性有着直接的影响。采用点热源加载形式的电弧有效半径可以认为近似为零，在这种情况下所得到结果通常不是很精确，在热源附近会出现局部高温现象，故而，一般采
热源加载形式
利用 01232 软件进行焊接模拟时，热源的加载
点热源、面热源和体热源，它们的形式一般有 / 种：载荷形式分别由热流率、热流密度和热生成率表现4/5。热流率作为节点集中载荷，主要用于线单元模型中，如果输入的值为正，代表热流流入节点，即单元获取热量。一般在焊接模拟中，采用点热源所形成的温度场分布易出现局部温度梯度过大等现象，故而不提倡使用此种模型；热流密度是一种面载，当通过单位面积的热流率已知或利用 01232 的 6789:901;6< 计算求得时，可以在模型相应的外表面施加热流密度，此种加载形式仅适用于实体和壳单元；第三种载荷形式是以热生成率作为体载荷施加于单元上 = 即为单位体积的热流率。以下是利用 / 种不同的加载形式在 01232 中模拟平板堆焊的情况>见图 ?@图 AB，
・ /F ・
用面热源或者体热源，或将两者结合后的组合热源。焊接模拟时电弧有效半径变小，热量相对集中，焊接熔池的熔深变大，熔宽变小；反之亦然。
!"#$%&’$ (#")’*+ ,-$.’*#
研究与设计
梁晓燕等：基于 /0121 平台焊接模拟中不同焊接热源的比较
第9期

ANSYS热分析教程及命令流算例

＝＝＝＝＝【热力耦合分析单元简介】＝＝＝＝＝＝SOLID5－三维耦合场实体具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。

本单元由8个节点定义，每个节点有6个自由度。

在静态磁场分析中，可以使用标量势公式（对于简化的RSP，微分的DSP，通用的GSP）。

在结构和压电分析中，具有大变形的应力钢化功能。

与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。

INFIN9－二维无限边界用于模拟一个二维无界问题的开放边界。

具有两个节点，每个节点上带有磁向量势或温度自由度。

所依附的单元类型可以为PLANE13和PLANE53磁单元，或PLANE55和PLANE77和PLANE35热单元。

使用磁自由度（ＡＺ）时，分析可以是线性的也可以是非线性的，静态的或动态的。

使用热自由度时，只能进行线性稳态分析。

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激光焊接命令流

/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
/REP,FAST
SET,LIST,999
SET,NEXT
/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
solve
*enddo
FINISH
/POST1
SET,FIRST
SET,FIRST
SET,LIST,2
PLDISP,0
!*
/EFACET,1
PLNSOL, TEMP,, 0
/EFACET,1
PLNSOL, TF,SUM, 0
/EFACET,1
PLNSOL, TF,SUM, 0
SET,LIST
SET,FIRST
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,19,1), 0.0, 99, 0, 1, -2, 0, 0

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。