| 查看: 1592 | 回复: 4 | |||
旗木卡卡西a铁虫 (初入文坛)
|
[求助]
ANSYS怎么样用数组在圆周周围加载变载荷 已有1人参与
|
|
二维平面圆周上加载激光，用数组来编写，思路是用数组来编写，用角度控制选择不同的节点，然后加上对应的热流密度。但是我对于数组用的很不通，最后计算出来的温度根本不对，求指导。 /COM,Laser !commence 开头的意思 /UNITS,SI !si国际单位制（m，kg，s，K） /TITLE,Laser ultrasound /FILNAME,thermal,1 !文件名 !******************前处理******************** /PREP7 ET,1,PLANE55 !定义材料属性 MPTEMP,,,,,,,, !设置材料参数 MPTEMP,1,300 !设置温度表格 MPTEMP,2,730 MPTEMP,3,2000 MPDATA,DENS,1,,2703 !密度 MPDATA,DENS,1,,2542 MPDATA,DENS,1,,2329 MPDATA,C,1,,926.67 !比热容 MPDATA,C,1,,1136.51 MPDATA,C,1,,1756.27 MPDATA,KXX,1,,225.5 !热导率 MPDATA,KXX,1,,191.1 MPDATA,KXX,1,,170.5 MPDATA,EX,1,,7.02E10 !基体弹性模量 MPDATA,EX,1,,7.02E10 MPDATA,EX,1,,7.02E10 MPDATA,PRXY,1,,0.34 !定义基体泊松比 MPDATA,PRXY,1,,0.34 MPDATA,PRXY,1,,0.34 UIMP,1,REFT,,,300 !定义基体线膨胀系数参考温度 MPDATA,ALPX,1,,2.35E-5 !定义基体线膨胀系数 MPDATA,ALPX,1,,2.35E-5 MPDATA,ALPX,1,,2.35E-5 !***************建立模型********************* CSYS,1 *AFUN,DEG RO=3 RI=2.5 PCIRC,RO,RI,0,8 PCIRC,RO,RI,8,180 PCIRC,RO,RI,180,360 AGLUE,ALL NUMCMP,ALL !*****************划分网格******************* ALLS LSEL,S,,,1 LSEL,A,,,3 a=0.05 LESIZE,All,,a !按照角度划分网格。密网格0.05° LSEL,A,,,6,9,1 b=0.5 LESIZE,All,,b !疏网格0.5° ALLS LSEL,S,,,2 LSEL,A,,,4 LSEL,A,,,5 LESIZE,ALL,,,50 ALLS MSHAPE,0,2D MSHKEY,1 MAT,1 !选择材料号 TYPE,1 !选择单元类型号 AMESH,ALL !****************建立激光参数**************** A=0.08 !吸收率 rlaser=300e-6 !光斑半径 elaser=13.5e-3 !单脉冲能量 TUP=10e-9 !脉冲上升时间 n=4 !单脉冲作用时间与上升时间的关系 tmax=(n*tup) !脉冲作用时间，假设去上升时间的3倍 P=(elaser/(3.14*rlaser*rlaser*tup)) !脉冲的峰值功率密度 PI=3.1415926 dm=RO*1000*PI*a/180 dt=4e-9 !步长 !**************建立激光源数组模型************ *SET,TNUM,10 !定义热流密度加载的时间的段数40ns/4ns *SET,YNUM,120 !定义热流密度加载的位置的段数0.001/20e-6 *DIM,FLUX2,TABLE,TNUM+1,YNUM+1,,TIME,Y !定义一个二维的数组 *DO,I,1,YNUM+1,1 !基于角度的循环次数 *SET,FLUX2(0,I),(I-1)*a !写入角度Y的坐标值 *ENDDO *DO,I,1,TNUM+1,1 !循环tt次，一个脉冲分为tt段加载（10次） *SET,FLUX2(I,0),(I-1)*dt !写入时间TIME的坐标值 *ENDDO *DO,I,1,TNUM+1,1 *DO,J,1,YNUM+1,1 *SET,DISTANCE,(J-1)*dm !I时刻，激光x坐标值 *SET,TTT,(I-1)*dt !I时刻，时间值 *IF,DISTANCE,LE,rlaser*10/3,THEN SPACEE=EXP(-(DISTANCE**2)/(rlaser**2)) TIMEE=(TTT/TUP)*EXP(-TTT/TUP) *SET,FLUX2(I,J),A*P*SPACEE*TIMEE *ELSE *SET,FLUX2(I,J),0 *ENDIF *ENDDO *ENDDO !**************非线性求解设置**************** /solu !进入求解器 alls antype,trans !瞬态分析(等价于antype,4 ) trnopt,full !完全瞬态积分法 solcontrol,on !自动控制 nropt,full,,on !定义完全牛顿-拉普森方法，激活自使用下降 pred,on !打开预测校正 lnsrch,on !线性搜索 tintp,,,,1,, ! 一阶瞬态积分参数 timint,on !时间积分器打开，瞬态效果打开 tref,300 !参考温度 neqit,50 !规定每个子步中最大迭代次数,默认为25 outres,all !输出结果：除SVAR和LOCI外的所有解决方案。 !*****************稳态分析******************* /solu !进入求解器 OUTRES,ERASE !重新设定结果数据的保存 timint,off !关闭瞬态效应 alls !全选 ic,all,temp,300 !initial condition，选中所有节点加温度20摄氏度 time,10e-10 !稳态分析的时间 deltim,1e-10,1e-10,1e-10 !指定用于此加载步骤的时间步大小 kbc,0 !阶跃方式加载 allsel,all solve !求解 /PBF,TEMP,,1 /REPLOT !*******************瞬态分析***************** TSRES,ERASE !删除之前时间子步的设置 T=4E-9 TIMINT,ON KBC,0 !阶跃方式加载 *DO,K,1,TNUM,1 TIME,T DELTIM,DT,DT,DT !用于指定该载荷步的时间步大小 NSUBST,10,,,1 !设定子载荷步数 outres,ALL,LAST !保存最后一子步的所有数据 SFL,1,hflux,%FLUX2% !对线6施加激光作用等效热流 ALLSEL SOLVE T=T+dt *ENDDO |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
2楼2018-03-06 09:28:11
旗木卡卡西a
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 13.9
- 散金: 70
- 帖子: 47
- 在线: 32.5小时
- 虫号: 7482999
- 注册: 2017-11-16
- 专业: 交叉学科中的光学问题
3楼2018-03-06 16:47:00
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
4楼2018-03-12 09:35:04
旗木卡卡西a
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 13.9
- 散金: 70
- 帖子: 47
- 在线: 32.5小时
- 虫号: 7482999
- 注册: 2017-11-16
- 专业: 交叉学科中的光学问题
5楼2018-05-18 10:06:29
30