24小时热门版块排行榜    

查看: 1140  |  回复: 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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

中兵大将

新虫 (初入文坛)

这个就是为了加载高斯热源,将高于熔点的单元杀死,不知道哪出错了,算不出来,希望大神给指点指点啊
2楼2015-01-23 14:50:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

中兵大将

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 中兵大将 at 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 中兵大将 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 关于代码变化问题,想知道的进来 +9 且听虎啸 2026-08-07 12/600 2026-08-07 21:53 by 理想三旬2025
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见