| 查看: 1066 | 回复: 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 |
回复此楼» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有200人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 一道高斯型的积分问题，求大神指点，想破头了
已经有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
精华III:
20