应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 焊接仿真计算
101/1返回列表
查看: 8631  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tuzuoqi

银虫 (小有名气)

[求助] 焊接仿真计算 已有1人参与

/PREP7
et,1,solid70
pi=3.1415926
LSIZE=0.005
V=0.007
Tinc= LSIZE/V
e=0.8
Q=160*7
af=0.002
ar=0.004
b=0.002
c=0.002
ff=0.667
fr=1.333
Qmax11=6*e*3**(1/2)*ff*Q/(pi*pi**(1/2)*af*b*c)
Qmax12=6*e*3**(1/2)*fr*Q/(pi*pi**(1/2)*ar*b*c)
MPTEMP,,,,,,,,  
MPTEMP,1,297,693,1093,1293,1493,1693
MPTEMP,7,1743,1793
MPDATA,C,1,1,450,610,670,670,700,1490
MPDATA,C,1,7,9260,810
MPDATA,C,2,1,450,590,640,660,710,3160
MPDATA,C,2,7,7390,810
MPDATA,DENS,1,1,7770,7640,7480,7400,7310,7190
MPDATA,DENS,1,7,7010,6690
MPDATA,DENS,2,1,7800,7660,7500,7410,7310,7170
MPDATA,DENS,2,7,7010,6980
MPDATA,KXX,1,1,16.71,21.42,26.4,28.94,31.36,33.49
MPDATA,KXX,1,7,30.79,31.48
MPDATA,KXX,2,1,16.09,21.02,26.4,28.93,31.51,32.94
MPDATA,KXX,2,7,30.71,31.43
Block,0,0.1,0,0.1,0,0.004
Wpoff,0.05,0,0.004
Block,-0.002,0.002,0,0.1,0,0.0005
Wprota,0,0,-90
Wpoff,0,0,0.002
vsbw,all,,delete
Wpoff,0,0,-0.004
vsbw,all,,delete
Wprota,0,0,90
Wpoff,-0.002,0,0.001
K,,0.05,0,0.0055
Larc,13,14,25
FLST,2,2,4  
FITEM,2,17  
FITEM,2,33  
AL,P51X
voffst,18,-0.1
FLST,2,2,6,ORDE,2   
FITEM,2,2   
FITEM,2,4   
VADD,P51X   
FLST,2,2,5,ORDE,2   
FITEM,2,10  
FITEM,2,27  
AADD,P51X
FLST,2,2,5,ORDE,2   
FITEM,2,9   
FITEM,2,18  
AADD,P51X   
vglue,all
lsel,s,line,,20
lesize,all,,,40
lsel,s,line,,34
lsel,a,line,,17
lesize,all,,,4
lsel,s,line,,43
lesize,all,,,8
mat,2
vmesh,2
alls
lsel,s,line,,7
lsel,a,line,,37
lsel,a,line,,4
lsel,a,line,,2
lesize,all,,,30,0.1
lsel,s,line,,10
lesize,all,,,5
mat,1
vmesh,1
alls
lsel,s,line,,28,29
lsel,a,line,,25,26
lesize,all,,,30,0.1
lsel,s,line,,12
lesize,all,,,5
mat,1
vmesh,3,5
alls
finish
/solu
j=0
*do,i,1,21,1
j=j+tinc
XX=(i-1)*LSIZE
antype,4
trnopt,full
tunif,293
outres,all,all
time,j
nsubst,3,,,1
autots,1
kbc,1
lnsrch,1
tsres,erase
timint,1,therm
nsel,s,ext
sf,all,conv,200,293
allsel,all
local,11,2,0.05,XX,0.0055,,,,2,1
NSEL,S,LOC,X,,0.002                    
NSEL,R,LOC,Y,,180                    
NSEL,R,LOC,Z,-90,0                    
CM,NODEF,NODE   
local,12,0,0.05,XX,0.0055,,,,1,1
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNCCSYS   
*SET,_FNCNAME,'front'   
*DIM,_FNC_C1,,1
*SET,_FNC_C1(1),qmax11  
*SET,_FNCCSYS,12
! /INPUT,.\Desktop\front1.func,,,1  
*DIM,%_FNCNAME%,TABLE,6,23,1,,,,%_FNCCSYS%  
!   
! Begin of equation: Qmax11*exp(-3(({x}/0.002)^2+({y}/0.002)^2+({z}/0.002)^2))  
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*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
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 0.002, 0, 0, 2  
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, 2, 4, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -2, 17, -1   
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 0.002, 0, 0, 3  
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, 3, 4, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -2, 17, -1  
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 0.002, 0, 0, 4
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 4, 4, -2
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, -4, 17, -2  
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -1, 1, -5   
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 3, -4   
*SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, 17, 3, -1   
*SET,%_FNCNAME%(0,23,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qmax11*exp(-3(({x}/0.002)^2+({y}/0.002)^2+({z}/0.002)^2))
!-->
BF,nodef,HGEN,%front%
CMDELE,nodef   
ALLSEL,ALL
local,13,2,0.05,XX,0.0055,,,,2,1
NSEL,S,LOC,X,,0.002                    
NSEL,R,LOC,Y,-180,0                    
NSEL,R,LOC,Z,-90,0                    
CM,NODEB,NODE   
local,14,0,0.05,XX,0.0055,,,,1,1
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNCCSYS   
*SET,_FNCNAME,'back'
*DIM,_FNC_C1,,1
*SET,_FNC_C1(1),qmax12  
*SET,_FNCCSYS,14
! /INPUT,.\Desktop\back.func,,,1
*DIM,%_FNCNAME%,TABLE,6,23,1,,,,%_FNCCSYS%  
!   
! Begin of equation:  Qmax12*exp(-3(({x}/0.002)^2+({y}/0.004)^2+({z}/0.002)^2))
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*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
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 0.002, 0, 0, 2  
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, 2, 4, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -2, 17, -1   
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 0.004, 0, 0, 3  
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, 3, 4, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -2, 17, -1  
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 0.002, 0, 0, 4
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 4, 4, -2
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, -4, 17, -2  
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -1, 1, -5   
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 3, -4   
*SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, 17, 3, -1   
*SET,%_FNCNAME%(0,23,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation:  Qmax12*exp(-3(({x}/0.002)^2+({y}/0.004)^2+({z}/0.002)^2))   
!-->
BF,nodeb,HGEN,%back%
CMDELE,nodeb  
ALLSEL,ALL
SOLVE
BFDELE,all,HGEN
*ENDDO


焊接仿真的命令流，出现了图中的问题，那位大神解释一下怎么回事

焊接仿真计算
12.png
回复此楼

» 猜你喜欢
1楼 2017-12-19 15:38:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
tuzuoqi: 金币+10, ★有帮助 2017-12-20 08:43:39
在你的循环：
    *do,i,1,21,1
  中，定义了一个函数数组，名为 FRONT，由于循环的缘故，所以每个循环都要重新定义该表数组，导致上述错误。
  可以把定义表数组的命令：
    *DIM,%_FNCNAME%,TABLE,6,23,1,,,,%_FNCCSYS%  
  移到循环之前，并改为：
    *DIM,FRONT,TABLE,6,23,1,,,,%_FNCCSYS%
赞一下(1人) 回复此楼
2楼2017-12-20 07:41:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuzuoqi

银虫 (小有名气)
引用回帖:
2楼: Originally posted by shgao20 at 2017-12-20 07:41:57
在你的循环：
    *do,i,1,21,1
  中，定义了一个函数数组，名为 FRONT，由于循环的缘故，所以每个循环都要重新定义该表数组，导致上述错误。
  可以把定义表数组的命令：
    *DIM,%_FNCNAME%,TABLE,6,23,1,,, ...

那个图中左下角为什么会出现红色的一块，能给解释一下吗
赞一下 回复此楼
3楼2017-12-20 08:45:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)
抱歉，不知道是什么。不过，应该不影响计算的。
赞一下 回复此楼
4楼2017-12-21 07:48:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuzuoqi

银虫 (小有名气)
引用回帖:
4楼: Originally posted by shgao20 at 2017-12-21 07:48:35
抱歉，不知道是什么。不过，应该不影响计算的。

能在请教您一个问题吗，就是对流的施加只能在外表面节点施加，但是对于多层的堆焊，堆焊完第一层时，需要在第一层表面施加冷却过程的对流，可是在焊接的有限元模型中，堆焊的几层是一块建模的，后面层会覆盖前面层，这个怎么施加对流
赞一下 回复此楼
5楼2017-12-22 23:22:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)
没有这样做过。不过，由于堆焊的材料，一般是预先建一个最终的完整模型，并杀死所有需要堆焊的材料对应的单元，然后一层一层的激活堆焊材料的单元。可以试试如下的两种方法：
    方法 1：  在建模时最好将每一层堆焊的材料都建为一个实体。相应的堆焊层就会具有对应的表面 (areas)，然后，每激活一层堆焊单元，就在相应的表面上施加对流载荷。添加新的堆焊单元后注意删除前面的对流载荷。
    方法 2：  如果建模时没有为各层堆焊材料创建相应的实体和表面，那么可以在每激活一层堆焊单元后，选择该层单元表面的所有节点，然后对这些节点施加 (表面) 对流载荷 (SF 命令)，软件会自动将施加到这些节点上的表面载荷转换到与这些节点相连的单元表面上。同样，添加新的堆焊单元后注意删除前面的对流载荷。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

tuzuoqi
6楼2017-12-23 09:01:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuzuoqi

银虫 (小有名气)
送红花一朵
引用回帖:
6楼: Originally posted by shgao20 at 2017-12-23 09:01:49
没有这样做过。不过，由于堆焊的材料，一般是预先建一个最终的完整模型，并杀死所有需要堆焊的材料对应的单元，然后一层一层的激活堆焊材料的单元。可以试试如下的两种方法：
    方法 1：  在建模时最好将每一层堆 ...

谢谢。
回复此楼
7楼2017-12-23 09:54:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuzuoqi

银虫 (小有名气)
引用回帖:
6楼: Originally posted by shgao20 at 2017-12-23 09:01:49
没有这样做过。不过，由于堆焊的材料，一般是预先建一个最终的完整模型，并杀死所有需要堆焊的材料对应的单元，然后一层一层的激活堆焊材料的单元。可以试试如下的两种方法：
    方法 1：  在建模时最好将每一层堆 ...

再请教您一下，多道焊接时，焊接第一道时，作为死单元的第二道也有温度，而且还比较高，这个正常吗？
赞一下 回复此楼
8楼2017-12-23 19:07:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)
由于死单元实际上也参与计算，只是它的材料属性被改为很小的值，因而死单元也会有温度，但是不应该太高。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

tuzuoqi
9楼2017-12-24 07:07:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuzuoqi

银虫 (小有名气)
送红花一朵
引用回帖:
9楼: Originally posted by shgao20 at 2017-12-24 07:07:53
由于死单元实际上也参与计算，只是它的材料属性被改为很小的值，因而死单元也会有温度，但是不应该太高。

那死单元温度偏高的原因一般是什么呢？
赞一下 回复此楼
10楼2017-12-24 10:41:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tuzuoqi 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于代码变化问题，想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了，今天 +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
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭