求教--Dload如何实现在不同载荷步不同加载平面的加载？

virgilee

各位高手，我想在两个载荷步里分别对模型的不同部分施加动态载荷，请问如何实现，我现在是不知道如何把所需的载荷步和加载面的信息传递到Dload里边，请各位指教。下边的子程序可以编译运行，但是没有应力结果，感觉应该是if里边参数调用的问题，请赐教。


SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,
     1 COORDS,JLTYP,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION TIME(2), COORDS (3)
      CHARACTER*80 SNAME
C
      PARAMETER(P0=5e4,Dm0=250,Dmiu=0.002,speed=20.0,DL=6.0,width=0.2)
      Pi=3.1415926535897932384626433832795
      omega=2.0*Pi*speed/DL !角速度 speed为移动速度
      Disp=TIME(1)*speed !移动距离
      xMin=Disp+0.002 !0.002 设置误差范围，积分点坐标和节点坐标不一样
      xMax=Disp+width-0.002 ! width 加载宽度
      IF(SNAME.EQ.'ASSEMBLY_SURF_Y'.AND.KSTEP.EQ.1) THEN ！想实现在surf_y面上在载荷步1中施加载荷
                   if(COORDS(1).ge.xMin.and.COORDS(1).lt.xMax) then
                        F=P0+Dm0*Dmiu*omega*omega*sin(omega*TIME(1)) !定义载荷大小
                    end if
      ELSE
         F=0
      END IF
      IF(SNAME.EQ.'ASSEMBLY_SURF_z'.AND.KSTEP.EQ.2) THEN ！想实现在surf_z面上在载荷步2中施加载荷
                 if(COORDS(1).ge.xMin.and.COORDS(1).lt.xMax) then
                     F=P0+Dm0*Dmiu*omega*omega*sin(omega*TIME(1)) !定义载荷大小
                 end if
       ELSE
          F=0
      END IF
        
      RETURN
      END

tengyuan

请问你的问题解决了么,我也遇到了同样的问题，而且我以后主做Abaqus二次开发，可以交流一下，微信：lty1002845535

tengyuan

答主你好，请问子程序Dload加载不同平面不同分析步的问题，你解决了么？想要请教

kb24lal

dload和vdload有什么区别？