| 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]