| 查看: 1154 | 回复: 2 | ||
[求助]
加载高斯热源的一个分析，总是出错，希望大神指点一下
|
|
/clear !清除数据库 /UNITS,SI ！使用国际单位 /PREP7 ！进入前处理器 ET,1,solid70 et,2,surf152 MPTEMP,1,20,250,500,750,1000,1500 MPTEMP,,1700,2500 ！定义材料1 性能的温度表 MPDATA,DENS,1,1,7820,7700,7610,7550,7490,7350 MPDATA,DENS,1,,7300,7090 ！定义各温度下的密度 MPDATA,KXX,1,1,50,47,40,27,30,35 MPDATA,KXX,1,,140,142 !定义各温度下的X上热传导系数 MPDATA,c,1,1,460,480,530,675,670,660 MPDATA,c,1,,780,820 block,0/1000,3/1000,0/1000,3/1000,,1/1000 ！建立一个矩形 block,3/1000,10/1000,0/1000,3/1000,,1/1000 ！建立另一个矩形 vglue,all ！把两个面粘结起来 NUMCMP,ALL ！压缩所有的号码 /pnum,line,1 !打开线编号 vplot !常数定义 PI=3.1415 U=20 ！电压 I=100 ！电流 eff=0.18 ！有效效率 q=eff*U*I ！最大的热流量 Ra=2.1/1000 !电弧半径 QMAX=3*Q/(pi*ra**2) SizeE=0.05/1000 !网格大小 lesize,5,SizeE ！定义1、3线的划分单元的尺寸 lesize,7,SizeE lesize,8,,,40,1/2 ！定义4、2、5 的划分个数和相邻单元的比例 lesize,6,,,40,2 lesize,14,,,40,2 lesize,19,,,15,3 lesize,20,,,15,3 lesize,9,,,25,2 asel,s,loc,z,1/1000 !选中实体上表面 AATT,1, , 2, !指定实体上表面用152号单元 MSHAPE,0,2D !实体单元类型六面体 MSHKEY,1 amesh,all !对上表面划分网格 allsel,all aatt,1,,1 ！ 给相关的元素分配材料属性 MSHAPE,0,3D !实体单元类型六面体 MSHKEY,1 vmesh,all save ！保存 eplot ！ finish ！从当前的处理器中退出 /solu ！进入求解处理器 ANTYPE,TRANS ！进入瞬态分析 tunif,25 ！定义统一温度 TOFFST,273 !定义绝对温度与所选温度的偏差 allsel nsel,s,loc,x,0 dsym,symm,x !施加关于x轴（yz面）的对称边界条件 nsel,s,loc,y,0 dsym,symm,y !表面传热状态 allsel esel,s,type,,2 !选中实体上表面的表面效应单元以方便加对流 sfe,all,,conv,2,10000 allsel asel,s,loc,z,1/1000 !加载热流，选择高斯热源面 nsla,r,1 !选择面上所有节点 esln,r !选择和这些节点有关的单元 eplot !所有单元可视 *do,Xlocatio,0,Ra,SizeE ！定义do循环 *do,Ylocatio,0,Ra,SizeE ！定义do循环 rr=sqrt(xlocatio**2+ylocatio**2) *if,rr,le,ra,then qr=qmax*exp(-3*(rr/ra)**2) ！定义所加的热量？？？？ nsel,s,loc,y,ylocatio ！选择Y坐标为8的节点 nsel,r,loc,x,xlocatio ！再选择X坐标为XLOC的节点 nsel,r,loc,z,1/1000 f,all,heat,qr ！在所选的节点上施加热载荷 *endif *enddo !逐点计算并施加热流 *enddo allsel NSUBST,8,16,5 ！定义载荷步的子步数量 TIME,1e-5 !时间步长(0.001,0.0005) KBC,0 !斜坡 载荷加载方式 AUTOTS,On !自动时间步长 timint,on !Turns on transient effects. OUTRES,ALL !All solution items written to the database OUTPR,ALL,LAST !-- Controls the solution printout TINTP,,,,1 !Defines transient integration parameters TRNOPT,FULL !Specifies transient analysis options EQSLV,ICCG,,2 ！Specifies the type of equation solve save saver LNSRCH,ON ！打开线搜索 solve /POST1 !进入POST1 SET,1,last etable,tempnode,temp,, !将总应变存入ETABLE etable,tempreture,bfe,temp esel,s,etab,tempnode,1350, !选择所有总应变大于或等于0.20的单元 esel,s,etab,tempreture,, FINISH ! /SOLU !重新进入求解器 ANTYPE,，rest !在分析，缺省为new nropt,full !打开牛顿- estif, EKILL,ALL !杀死选择（超过允许值）的单元 ESEL,s,live |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有218人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 一道高斯型的积分问题，求大神指点，想破头了
已经有29人回复
- 谁有ANSYS移动热源模拟T型接头焊接过程的命令流呀？
已经有9人回复
- 求助高斯热源中对流边界的设置
已经有4人回复
- 求大神帮忙，高斯09 “open-new-file”错误
已经有6人回复
- 请教一下，如何利用Ansys在一个长方形中创建50个椭圆。
已经有12人回复
- 施加高斯热源可以用temperature添加吗 还是就得是heat flux函数
已经有11人回复
- 在施加旋转高斯热源时总提示错误，但感觉程序没问题，不知错在哪？
已经有5人回复
- 高斯热源UDF，以及在何处加载？
已经有4人回复
- ANSYS 加载函数化为表格
已经有7人回复
- 新手求助高斯输入文件写法，求大神帮忙
已经有8人回复
- 2种焊接移动高斯热源的程序，有什么区别吗
已经有7人回复
- 打算配置一个单8核的或者两个双4核的机器进行高斯计算，求大神指导如何装备
已经有11人回复
- 能量源项加载问题
已经有19人回复
- 为什么一计算就会出现系统找不到指定文件呢？求高斯大神解答啊....
已经有5人回复
- 【讨论】关于高斯计算不同条件下的热力学函数
已经有8人回复
- 【求助】关于高斯热源峰值功率随时间变化的问题
已经有7人回复
2楼2015-01-23 14:50:13
|
这里有一点点与焊接有关的apdl，共享给大家 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : ansys实例.zip
2016-01-21 16:23:37, 32.58 K
3楼2016-01-21 16:23:46
投票:
30