求助：怎样移动加载啊

吴鑫

下面是我的程序，do循环了，但是坐标系不跟着移动，热源也就不动。期待大家的建议：
/solu
antype,trans
tunif,20                  !对所有节点指定一个均布温度        
kbc,1                     !指定载荷为阶跃载荷
timint,on                 !打开瞬态效应
autots,on                 !自动时间步长跟踪
!选择单元加载对流条件

。。。。。。。。。。。。。
*dim,laserbodysource,table,11,5,1,z,x   !定义激光体热源数组,11行5列1
*tread,laserbodysource,t1,txt,,0           !从t1读取数据laserbodysource

v=10/60
j=0   
   
*do,t,1e-6,0.0002*50,0.0002
        *if,j,gt,0,then
  !卸载
    cmsel,,lasersource              !用来选中之前创建的集合,和asel等类似
    bf,all,hgen,0                   !BFE或BF来加载热源，hgen生热率
    cmdel,lasersource
      

    local,13,1,0,y+distance,0,,,,1,1  !distance
    nsel,s,loc,x,0,Rarc               !Rarc
    nsel,r,loc,z,0
    sfdele,all,hflux                  !将定义的面负载删除，hflux热流heat flux
    sf,all,conv,%convection%,20       !定义节点间对流载荷
    allsel
*else
allsel
    *endif
j=j+1
            
        time,100000*t                !指定荷载步结束时间为t
! fcum,add                     !Specifies that force loads are to be accumulated
    nsubst,50,,,14        !指定载荷步中需要的子步数，此荷载步的子步数为50            
    autots,0              !是否使用自动时间步长,0是什么意思了
    deltim,,,,2           !在本载荷步中指定时间步长大小
    kbc,0                 !线性插值
    outres,all,all        !规定写入数据库的求解信息，全写入

    y=v*t                     !不同时间时刻对应不同y，不同y处建立局部坐标系和加载热源
   
    local,11,1,0,y,0,,,,1,1,  !定义局部坐标,坐标系代号11，
    nsel,s,loc,x,0,Re+remain  !选取从光斑中心到Re+remain的表面节Re为光斑半径，remain=re/20
    nsel,r,loc,z,0,0.0025     !选取表面到z=0.0025处节点
    cm,lasersource,node       !定义组元（集合）lasersource，将几何元素分组形成组元，node为组元类型
    nplot
    bf,all,hgen,%laserbodysource%  !BFE或BF来加载热源，在节点上施加载荷

    allsel         
    solve
!  save
!  finish  

        
       /solu
       antype,trans,rest                  !读入当前最高载荷步的重启动文件

        /post1
        set,last             !读取求解结果
       etable,tmp,temp       !或者ttot
       esel,s,etab,tmp,150000  !melt=1000
       finish

       /solu
       resume
       antype,trans,rest                  !读入当
       *get,elemnum,elem,0,count          !Number of elements in the selected set
       *if,elemnum,gt,0,then
       ekill,all
       esel,s,live
!      nsle,s,1
!      nsel,invert
!      d,all,all,0
       *else
       allsel
       *endif
      
       time,100000*t+10
      
       allsel
       solve


!    /solu
!     resume
!     antype,trans,rest                  !读入当前最高载荷步的重启动文件

*enddo

吴鑫

大家帮帮忙啊~~先谢谢啦，急急急

守护蓝天

利用函数加载器怎么样？ 2# 吴鑫

吴鑫

您好，能不能具体说说啊？新手不太明白。不胜感激~ 3# 守护蓝天