友情提示:如果本网页打开太慢或显示不完整,请尝试鼠标右键“刷新”本网页!阅读过程发现任何错误请告诉我们,谢谢!! 报告错误
一世书城 返回本书目录 我的书架 我的书签 TXT全本下载 进入书吧 加入书签

高斯赋值-第2章

按键盘上方向键 ← 或 → 可快速上下翻页,按键盘上的 Enter 键可回到本书目录页,按键盘上方向键 ↑ 可回到本页顶部!
————未阅读完?加入书签已便下次继续阅读!



autots;on
kbc;1
deltim;0。01;0。01;0。2          !时间步长                   
!***********************************以下为施加热流密度载荷  
*do;i;nemin;nemax;1
  *if;esel(i);eq;1;then       !如果该单元在上面的集合里,则
   xsy=centrx(i)               !读取该单元的中心坐标
   ysy=centry(i)
   zsy=centrz(i)
   rr=abs(sqrt(xsy*xsy+(ysy…b)*(ysy…b)+(zsy…c)*(zsy…c)))    !该单元中心距离热源中心
   *if;rr;le;radius;then                  !该单元中心在加热半径范围内
    qr=qmax*exp(…k*rr*rr)           !该单元中心处的热流大小
    sfe;i;nmface(i);hflux;;qr           !在每个单元所指定的面上施加热源载荷
   *endif
  *endif
*enddo
  
allsel                                        !全选求解
outres;all;all
solve
   
esel;s;type;;1                   
esel;r;live                          
vsel;s;;;v1                              !选择焊缝1上表面单元    
eslv;r
nsle 
asel;s;;;a1 
nsla;r;1
esln;r;0  
    
!!!!!!!!!!!!!!输入先暂时停止一下
!插入〃/post1plnsol;temp/psf;hflux;;1;1;on/replot〃的内容进行查看
!问题:热流载荷的分布和温度场的分布不一致(该完全一致才对),有点偏移?
!若要继续,读入〃finish/solu〃
  
*do;i;1;6
  sfedele;all;i;hflux                       !删除每个单元六个面上的热载荷
*enddo 
       
*enddo

alls
esel;s;live
eplot                                            !查看温度场                          
save;'1。11thermal_analysis_alive1';'db'
!=
!生成节点温度分布云图动画文件的程序段
!=
/post1
/seg;dele                             !允许图形数据存储到局部终端存储器
/dscale;1;1。0            !显示位移时,设置放大系数
avprin;0;0                !规定主项量和总向量如何被计算
avres;1                    !规定当powergraphies可用时,结果数据如何被平均
/seg;multi;1。11temp…10;0。2    !存动画;分节存储后来的显示;动画名字;每个画面间延迟系数

esel;s;live
!*************************************************读取初始稳态温度场结果
set;;;1;;dt0    !从结果文件读;载荷步;子步;比例因子实部;比例因子虚部;时间
plnsol;temp   
!*****************************************读取第1道焊缝过程中的温度场结果
*do;im;0;lnum;1 
tm=dt1+im*tinc 
set;;;1;;tm   
plnsol;temp                           
*enddo 

/seg;off;1。11temp…10;0。2
anim;1;1          !*********************************生成动画完毕

!****************************************存储第1道焊缝过程中的温度场结果图片
esel;s;live
*do;im;0;lnum;1 

*if;im;eq;0;then
   tm=dt1
   set;;;1;;tm   
   plnsol;temp
   /image;save;1。11temp0;jpeg 
*endif
  
tm=dt1+im*tinc 
set;;;1;;tm   
plnsol;temp
  
*if;tm;eq;dt1+lnum/2*tinc;then   
   set;;;1;;tm   
   plnsol;temp6。5;jpeg 
*endif   
                     
*enddo
/image;save;1。11temp12。5;jpeg   
!温度场完毕以第一道为例
返回目录 上一页 回到顶部 0 0
未阅读完?加入书签已便下次继续阅读!
温馨提示: 温看小说的同时发表评论,说出自己的看法和其它小伙伴们分享也不错哦!发表书评还可以获得积分和经验奖励,认真写原创书评 被采纳为精评可以获得大量金币、积分和经验奖励哦!