24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2005  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

1970-01-01 08:00:00
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

jiaka

铁虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
激光焊接和打孔 热源应该相同的吧  但是如果要模拟打孔的效果的话  估计要用fluent
不屈不挠,一往无前
4楼2014-04-23 14:45:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

小七工作室1

新虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
有焊接的  不清楚和你要的一样不?
2楼2014-04-23 11:00:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

King2012

金虫 (小有名气)

引用回帖:
2楼: Originally posted by 小七工作室1 at 2014-04-23 11:00:01
有焊接的  不清楚和你要的一样不?

谢谢,麻烦你发个。
人生就像做实验总是存在太大的变数,你也不会知道哪个条件的变化对你的最后会产生很大的影响。
3楼2014-04-23 14:05:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaka

铁虫 (初入文坛)

【答案】应助回帖

/FILNAME,Welding residual stress,0
!进入前处理
/PREP7  
ET,1,SOLID226,11 !定义热-结构耦合单元
ET,2,SURF152     !定义结构的表面效应单元
KEYOPT,2,8,1     !设置表面效应单元关键字

U=22          !焊接电压
I=150         !焊接电流
V=0.01        !焊接速度   
YITA=0.7      !焊接热效率
R=0.008       !电弧有效加热半径
Q=U*I*YITA    !电弧热功率
Qm=3/3.1415/R**2*Q !加热斑点中心最大热流密度
!设置材料参数
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000   
MPTEMP,6,1500   
MPTEMP,7,1700   
MPTEMP,8,2500
   
MPDATA,KXX,1,,50
MPDATA,KXX,1,,47
MPDATA,KXX,1,,40
MPDATA,KXX,1,,27
MPDATA,KXX,1,,30
MPDATA,KXX,1,,35
MPDATA,KXX,1,,40
MPDATA,KXX,1,,55
  
MPDATA,DENS,1,,7820
MPDATA,DENS,1,,7700
MPDATA,DENS,1,,7610
MPDATA,DENS,1,,7550
MPDATA,DENS,1,,7490
MPDATA,DENS,1,,7350
MPDATA,DENS,1,,7300
MPDATA,DENS,1,,7090

MPDATA,C,1,,460
MPDATA,C,1,,480
MPDATA,C,1,,530
MPDATA,C,1,,675
MPDATA,C,1,,670
MPDATA,C,1,,660
MPDATA,C,1,,780
MPDATA,C,1,,820   

MPDATA,EX,1,,2.05E11
MPDATA,EX,1,,1.87E11
MPDATA,EX,1,,1.5E11
MPDATA,EX,1,,0.7E11
MPDATA,EX,1,,0.2E11
MPDATA,EX,1,,0.19E2
MPDATA,EX,1,,0.18E2
MPDATA,EX,1,,0.12e2
  
MPDATA,PRXY,1,,0.28
MPDATA,PRXY,1,,0.29
MPDATA,PRXY,1,,0.31
MPDATA,PRXY,1,,0.35
MPDATA,PRXY,1,,0.4  
MPDATA,PRXY,1,,0.45
MPDATA,PRXY,1,,0.48
MPDATA,PRXY,1,,0.5

UIMP,1,REFT,,,20
MPDATA,ALPX,1,,1.1e-5   
MPDATA,ALPX,1,,1.22e-5  
MPDATA,ALPX,1,,1.39e-5  
MPDATA,ALPX,1,,1.48e-5  
MPDATA,ALPX,1,,1.34e-5  
MPDATA,ALPX,1,,1.33e-5  
MPDATA,ALPX,1,,1.32e-5  
MPDATA,ALPX,1,,1.31e-5  


TB,BISO,1,6,2,  
TBTEMP,20   
TBDATA,,225e6,0,,,,
TBTEMP,250  
TBDATA,,185e6,0,,,,
TBTEMP,500  
TBDATA,,85e6,0,,,,  
TBTEMP,750
TBDATA,,45E6,0,,,,   
TBTEMP,1000
TBDATA,,15E6,0,,,,   
TBTEMP,1500
TBDATA,,1E-5,0,,,,   
!建立模型
BLOCK,0,0.18,0,0.005,0,0.15,
BLOCK,0.18,0.23,0,0.005,0,0.15,  
BLOCK,0.23,0.41,0,0.005,0,0.15,
!划分网格
ESIZE,0.005
VMESH,ALL  
!生成表面效应单元
NSEL,ALL
NSEL,S,LOC,Y,0.005   
NPLOT   
TYPE,   2   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
TSHAP,LINE  

ESURF,0

ALLSEL,ALL
!定义焊件之间的接触对
!*  
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,
MAT,1   
MP,EMIS,1,  
R,3
REAL,3  
ET,3,170
ET,4,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,2E6,1.0,,1.0,0.5  
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0
KEYOPT,4,10,2   
KEYOPT,4,11,0   
KEYOPT,4,12,0   
KEYOPT,4,2,0
KEYOPT,3,5,0
KEYOPT,4,1,1
! Generate the target surface   
ASEL,S,,,6  
CM,_TARGET,AREA
TYPE,3  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,11
CM,_CONTACT,AREA
TYPE,4  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,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,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
!*  
!*  
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0   
MAT,1   
R,4
REAL,4  
ET,5,170
ET,6,174
KEYOPT,6,9,0
KEYOPT,6,10,2   
R,4,
RMORE,  
RMORE,,2E6  
RMORE,0
KEYOPT,6,1,1
! Generate the target surface   
ASEL,S,,,12
CM,_TARGET,AREA
TYPE,5  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,17
CM,_CONTACT,AREA
TYPE,6  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,5  
ESEL,A,TYPE,,6  
ESEL,R,REAL,,4  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,5  
ESEL,A,TYPE,,6  
ESEL,R,REAL,,4  
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,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/MREP,EPLOT

)/GOP    ! Resume printing after UNDO process   
)! We suggest a save at this point  

R,3,0,0,1, 0.1,0,0  
RMORE,0,0,1e+020, 0,1,0
RMORE,0,2000000,1, 0,1,0.5  
RMORE,,,1, 0,0,
RMORE,,,,,,
RMORE,,,,,1000  
R,4,0,0,0, 0,0,0
RMORE,0,0,0, 0,0,0  
RMORE,0,2000000,0, 0,0,0
RMORE,,,0, 0,0,
RMORE,,,,,,
RMORE,,,,,1000  
!进入求解
/SOL
ANTYPE,4  !设置分析类型为瞬态热-结构分析
TRNOPT,FULL   !使用完全法完成瞬态分析
NLGEOM,1      !激活大变形分析
!设置求解器
EQSLV,SPAR
TUNIF,20,     !设置模型的初始均匀温度
TREF,20,      !设置模型的参考温度
!定义对流边界条件
SFA,4, ,CONV,30,20
SFA,10, ,CONV,30,20
SFA,16, ,CONV,30,20
SFA,2, ,CONV,30,20
SFA,8, ,CONV,30,20
SFA,14, ,CONV,30,20
SFA,5, ,CONV,30,20
SFA,18, ,CONV,30,20
SFA,1, ,CONV,30,20
SFA,7, ,CONV,30,20
SFA,13, ,CONV,30,20
SFA,3, ,CONV,300,20
SFA,9, ,CONV,300,20
SFA,15, ,CONV,300,20
!定义结构约束
DA,5,UX,
DA,5,UY,
DA,5,UZ,
DA,18,UX,
DA,18,UY,
DA,18,UZ,
DA,3,UY,
DA,9,UY,
DA,15,UY,
DA,2,UZ,
DA,8,UZ,
DA,14,UZ,
!定义函数
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS   
*SET,_FNCNAME,'gaosi1'  
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),qm  
*SET,_FNC_C2(1),v   
*SET,_FNC_C3(1),r   
*SET,_FNCCSYS,0
! /INPUT,gaosi1.func,,,1
*DIM,%_FNCNAME%,TABLE,6,23,1,,,,%_FNCCSYS%  
!   
! Begin of equation:  Qm*exp(-3*(({X}-0.18)^2+(0.15-V*{TIME}-{Z})^2)/R^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), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*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, 3, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.18, 0, 0, 2   
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 2, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -3, 17, -1   
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, 18, 3, 1
*SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 0.15, 0, 0, -1
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, -3, 2, -1   
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -5, 2, 4
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, -1, 17, -3  
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, 19, 17, -1  
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 4, -2   
*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:  Qm*exp(-3*(({X}-0.18)^2+(0.15-V*{TIME}-{Z})^2)/R^2)
!-->
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS   
*SET,_FNCNAME,'GAOSI2'  
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),QM  
*SET,_FNC_C2(1),V   
*SET,_FNC_C3(1),R   
*SET,_FNCCSYS,0
! /INPUT,GAOSI2.func,,,1
*DIM,%_FNCNAME%,TABLE,6,23,1,,,,%_FNCCSYS%  
!   
! Begin of equation: Qm*exp(-3*(({X}-0.23)^2+(V*({TIME}-15)-{Z})^2)/R^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), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*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, 3, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.23, 0, 0, 2   
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 2, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -3, 17, -1   
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 15, 0, 0, 1
*SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, 1, 2, -1
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, 18, 3, -3   
*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, -1, 2, 4
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, -3, 17, -1  
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, 19, 17, -1  
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 4, -2   
*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: Qm*exp(-3*(({X}-0.23)^2+(V*({TIME}-15)-{Z})^2)/R^2)  
!-->
!定义第1个载荷步文件
ESEL,ALL
ESEL,S,TYPE,,2  
EPLOT   
SFE,all,1,HFLUX, , %GAOSI1%
ALLSEL,ALL
OUTRES,ALL,ALL,
TIME,15
AUTOTS,1
NSUBST,75,150,50,1
KBC,0   
LSWRITE,1
!定义第2个载荷步文件
ESEL,ALL
ESEL,S,TYPE,,2  
EPLOT  
SFEDELE,ALL,1,HFLUX
SFE,all,1,HFLUX, , %GAOSI2%
ALLSEL,ALL
TIME,30
AUTOTS,1
NSUBST,75,150,50,1
KBC,0  
LSWRITE,2
!定义第3个载荷步文件
ESEL,ALL
ESEL,S,TYPE,,2  
EPLOT  
SFEDELE,ALL,1,HFLUX
ALLSEL,ALL
TIME,200
AUTOTS,1
NSUBST,600,,50,1
KBC,0  
LSWRITE,3
!定义第4个载荷步文件
TIME,3600
AUTOTS,1
NSUBST,100,150,50,1
KBC,0  
LSWRITE,4
!使用载荷步文件求解
LSSOLVE,1,4,1
不屈不挠,一往无前
5楼2014-04-23 14:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见