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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 有没有哪位仁兄有ansys激光打孔的教程或者是APDL语言？
51/1返回列表
查看: 2289  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

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的回帖
查看全部 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

铁虫 (初入文坛)

【答案】应助回帖

感谢参与，应助指数 +1
激光焊接和打孔 热源应该相同的吧  但是如果要模拟打孔的效果的话  估计要用fluent
赞一下 回复此楼
不屈不挠，一往无前
4楼2014-04-23 14:45:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭