Vuhard

"Wendel Melo [email protected] [Abaqus]" <[email protected]> Tue, 12 Mar 2019 18:42:06 +0000 (UTC)
Newsgroups gmane.comp.mathematics.abaqus.user
Message-ID <[email protected]>
Hi everybody, 
My name is Wendel, and I am a master student at University of British Columbia. It is a pleasure to be part of this group; I thank you all for having accepted me. I just started using Vuhard for the plasticity part of my simulation and I have followed Abaqus User Subroutines Manual to write my little code, please see it below:
subroutine vuhard(C Read only -     *     nblock,      *     jElem, kIntPt, kLayer, kSecPt,      *     lAnneal, stepTime, totalTime, dt, cmname,     *     nstatev, nfieldv, nprops,      *     props, tempOld, tempNew, fieldOld, fieldNew,     *     stateOld,     *     eqps, eqpsRate,C Write only -     *     yield, dyieldDtemp, dyieldDeqps,     *     stateNew )C      include 'vaba_param_dp.inc'C      dimension props(nprops), tempOld(nblock), tempNew(nblock),     1   fieldOld(nblock,nfieldv), fieldNew(nblock,nfieldv),     2   stateOld(nblock,nstatev), eqps(nblock), eqpsRate(nblock),     3   yield(nblock), dyieldDtemp(nblock), dyieldDeqps(nblock,2),     4   stateNew(nblock,nstatev), jElem(nblock)CC props(1)=KC Props(2)=nC Props(3)=m      character*80 cmnameC      do km = 1,nblock        yield(km)= props(1)*(eqps(km)**props(2))*(eqpsRate(km)**props(3))        dyieldDeqps(km,1)=props(1)*props(2)*(eqps(km)**(props(2)-1))*(eqpsRate(km)**props(3))     dyieldDeqps(km,2)=props(1)*props(3)*(eqps(km)**props(2))*(eqpsRate(km)**(props(3)-1))   end doC      return      end
As you can see I am using the Ludwick equation for my working hardening behavior and I have set the the parameters K,n and m under the Plastic part (hardening=user). However, I am getting an error like this:

Has anyone else got this error before? Can you give me some advice/ help? 
Thanks in advance,




[Non-text portions of this message have been removed]