24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

旗木卡卡西a

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 13.9
散金: 70
帖子: 47
在线: 32.5小时
虫号: 7482999
注册: 2017-11-16
专业: 交叉学科中的光学问题

[求助] 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

回复此楼

» 猜你喜欢

导师想让我从独立一作变成了共一第一已经有8人回复
博士读完未来一定会好吗已经有23人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
申请2026年博士已经有6人回复

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

1楼2018-03-05 15:52:26

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

旗木卡卡西a

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 13.9
散金: 70
帖子: 47
在线: 32.5小时
虫号: 7482999
注册: 2017-11-16
专业: 交叉学科中的光学问题

引用回帖:

2楼: Originally posted by shgao20 at 2018-03-06 09:28:11
似乎你的表数组定义有问题。
按照 ANSYS 的规定，对于瞬态分析，使用表数组加载时，该表数组应该有两列，行数为时间点数。
表数组的第 0 列 (表数组的行、列编号都是从 0 开始) 为时间点值，1 列为载荷值。而 ...

您能帮我改一下么，或者加个联系方式指导下我修改。对于数组不怎么会用