能否帮我看一下有无语法问题？

饭远长 · 发表于 2014-10-15 10:41:01

我是学机械的，Fortran真的是零基础学起。。最近要用ABAQUS进行二次开发，用Fortran改编了一个子程序，如下。但是ABAQUS调用子程序时提示错误：
Error in job 10_15sub: Problem during linking - Single Precision Abaqus/Explicit User Subroutines. This error may be due to a mismatch in the Abaqus user subroutine arguments. These arguments sometimes change from release to release, so user subroutines used with a previous release of Abaqus may need to be adjusted.
Job 10_15sub aborted due to errors.

我觉得可能是语法有错误，但是不知道哪里有问题，实在很头疼，敬请各位版主和朋友帮我解决一下，多谢！！

subroutine vusdfld(
c Read only -
   * nblock, nstatev, nfieldv, nprops, ndir, nshr,
   * jElemUid, kIntPt, kLayer, kSecPt,
   * stepTime, totalTime, dt, cmname,
   * coordMp, direct, T, charLength, props,
   * stateOld,
c Write only -
   * stateNew, field )
c
   include 'vaba_param.inc'
c
   dimension props(nprops),
   *       jElemUid(nblock), coordMp(nblock, *),
   *       direct(nblock, 3, 3), T(nblock,3,3),
   *       stateOld(nblock, nstatev),
   *       stateNew(nblock, nstatev),
   *       field(nblock, nfieldv)
   character*80 cmname
c Properties array
c    props(1) -> Transverse tensile strength, Yt
c    props(2) -> Matreix compressive strength, Yc
c    props(3) -> Ply shear strength, Sc
c    props(4) -> Fiber buckling strength, Xc
c    props(5) -> Initial shear modulus, G12
c    props(6) -> Nonlinear shear factor, alpha
c    props(7) -> stress max
c
   character*3 cData(maxblk*6)
   dimension jData(maxblk*6)
   dimension stress(maxblk*6),strain(maxblk*6)
c Read properties
   yt = props(1)
   yc = props(2)
   sc = props(3)
   xc = props(4)
   g12 = props(5)
   alpha = props(6)
   delet = props(7)
c Get stresses and strains from previous increment
   jStatus = 1
   call vgetvrm( 'S', stress, jData, cData, jStatus )
   jStatus = 1
   call vgetvrm( 'LE', strain, jData, cData, jStatus )
c
   call evaluateDamage( nblock, nstatev,
   *    nfieldv, ndir, nshr,
   *    yt, yc, sc, xc, g12, alpha,
   *    delet,
   *    stress, strain,
   *    stateOld,
   *    stateNew, field )
c
   return
   end
c
   subroutine evaluateDamage ( nblock, nstatev,
   *    nfieldv, ndir, nshr,
   *    yt, yc, sc, xc, g12, alpha,
   *    delet,
   *    stress, strain,
   *    stateOld,
   *    stateNew, field )
c
   include 'vaba_param.inc'
c
   dimension stress(nblock,ndir+nshr),
   *    strain(nblock,ndir+nshr),
   *    stateOld(nblock,nstatev),
   *    stateNew(nblock,nstatev),
   *    field(nblock,nfieldv)
   parameter ( zero = 0.d0 )
   delet = props(7)
c
c initialize failure flags from statev.
   do k = 1, nblock
      stateNew(k,1) = stateOld(k,1)
      stateNew(k,2) = stateOld(k,2)
      stateNew(k,3) = stateOld(k,3)
c
      em    = stateOld(k,1)
      efs = stateOld(k,2)
      damage = stateOld(k,3)
c
      s11 = stress(k,1)
      s22 = stress(k,2)
      s12 = stress(k,4)
*
      e12 = 2.0*strain(k,4) !e12 is engineering strain
c
c damage index: = 0 if no strain to prevent divide by zero
c
      damage = 0.d0
      if (stress(k,2).gt.delet) then
      stateNew(k,2) = zero
   else
      if (e12.ne.0) then
   *       damage = (3.d0*alpha*g12*s12**2-2.d0*alpha*(s12**3)/e12)/
   *                (1.d0+3.d0*alpha*g12*s12**2)
c
      f1 = s12**2/(2.d0*g12) + 0.75d0*alpha*s12**4
      f2 = sc**2 /(2.d0*g12) + 0.75d0*alpha*sc**4

c
c matrix tensile/compressive failure
      if (em .lt. 1.d0) then
         if (s22 .lt. 0.d0) then
            em = sqrt((s22/yc)**2 + f1/f2)
         else
            em = sqrt((s22/yt)**2 + f1/f2)
         endif
         stateNew(k,1) = em
      endif
c
c fiber-matrix shear failure
      if (efs .lt. 1.d0) then
         if (s11 .lt. 0.d0) then
            efs = sqrt((s11/xc)**2 + f1/f2)
         else
            efs = 0.d0
         endif
         stateNew(k,2) = efs
      endif
c
c state transition diagram
c
c fv1: matrix compr/tens failure
c fv2: fiber/matrix shear failure
c fv3: material damage (shear nonlinearity)
c                         fv1  fv2  fv3    e1  e2  nu12  g12
c
c (0) no failure             0 0 0  -->  e1  e2  nu12  g12
c (1) matrix (compr/tens)    1 0 0  -->  e1 0 0 g12
c (2) fib/mtx shear          0 1 0  -->  e1  e2 0    0
c (3) matrix & f/m shear    1 1 0  -->  e1 0 0    0
c (4) pure damage          0 0 1  -->  e1  e2  nu12 0
c (5) mtrx & damage          1 0 1  -->  e1 0 0    0
c (6) f/m shear & damage    0 1 1  -->  e1  e2 0    0
c (7) mtrx, f/m shr & damage  1 1 1  -->  e1 0 0    0
c
c    update field variables
c
      field(k,1) = 0.d0
      field(k,2) = 0.d0
      if (em .ge. 1.d0) then
         field(k,1) = 1.d0
      end if
      if (efs .ge. 1.d0) then
         field(k,2) = 1.d0
      end if
      field(k,3) = damage
      stateNew(k,3) = field(k,3)
      end if
   end do
c
   return
   end

楚香饭 · 发表于 2014-10-15 11:33:00

1.你没有给出各include 代码，我无法帮你检查。
2.编译器会帮你检查。我认为你的错误不是语法错误，而是编译器与 abaqus 的设置问题。
3.我不懂 abaqus

[求助] 能否帮我看一下有无语法问题？

新人勋章

爱心勋章

水王勋章

元老勋章

热心勋章