Vuhard

classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|

Vuhard

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

Reply | Threaded
Open this post in threaded view
|

Re: Vuhard

Wendel Melo
 The text I have written about the error I am getting when using V-uhard was not displayed in the way I would like to. So, I have attached a copy of my little script (.f), the inp. file and the error file that I am getting. I would appreciate if someone could give me some help or suggestions. Thanks...  
    On Wednesday, March 13, 2019, 3:59:02 AM PDT, Wendel Melo [hidden email] [Abaqus] <[hidden email]> wrote:

  
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]


  #yiv1166828722 #yiv1166828722 -- #yiv1166828722ygrp-mkp {border:1px solid #d8d8d8;font-family:Arial;margin:10px 0;padding:0 10px;}#yiv1166828722 #yiv1166828722ygrp-mkp hr {border:1px solid #d8d8d8;}#yiv1166828722 #yiv1166828722ygrp-mkp #yiv1166828722hd {color:#628c2a;font-size:85%;font-weight:700;line-height:122%;margin:10px 0;}#yiv1166828722 #yiv1166828722ygrp-mkp #yiv1166828722ads {margin-bottom:10px;}#yiv1166828722 #yiv1166828722ygrp-mkp .yiv1166828722ad {padding:0 0;}#yiv1166828722 #yiv1166828722ygrp-mkp .yiv1166828722ad p {margin:0;}#yiv1166828722 #yiv1166828722ygrp-mkp .yiv1166828722ad a {color:#0000ff;text-decoration:none;}#yiv1166828722 #yiv1166828722ygrp-sponsor #yiv1166828722ygrp-lc {font-family:Arial;}#yiv1166828722 #yiv1166828722ygrp-sponsor #yiv1166828722ygrp-lc #yiv1166828722hd {margin:10px 0px;font-weight:700;font-size:78%;line-height:122%;}#yiv1166828722 #yiv1166828722ygrp-sponsor #yiv1166828722ygrp-lc .yiv1166828722ad {margin-bottom:10px;padding:0 0;}#yiv1166828722 #yiv1166828722actions {font-family:Verdana;font-size:11px;padding:10px 0;}#yiv1166828722 #yiv1166828722activity {background-color:#e0ecee;float:left;font-family:Verdana;font-size:10px;padding:10px;}#yiv1166828722 #yiv1166828722activity span {font-weight:700;}#yiv1166828722 #yiv1166828722activity span:first-child {text-transform:uppercase;}#yiv1166828722 #yiv1166828722activity span a {color:#5085b6;text-decoration:none;}#yiv1166828722 #yiv1166828722activity span span {color:#ff7900;}#yiv1166828722 #yiv1166828722activity span .yiv1166828722underline {text-decoration:underline;}#yiv1166828722 .yiv1166828722attach {clear:both;display:table;font-family:Arial;font-size:12px;padding:10px 0;width:400px;}#yiv1166828722 .yiv1166828722attach div a {text-decoration:none;}#yiv1166828722 .yiv1166828722attach img {border:none;padding-right:5px;}#yiv1166828722 .yiv1166828722attach label {display:block;margin-bottom:5px;}#yiv1166828722 .yiv1166828722attach label a {text-decoration:none;}#yiv1166828722 blockquote {margin:0 0 0 4px;}#yiv1166828722 .yiv1166828722bold {font-family:Arial;font-size:13px;font-weight:700;}#yiv1166828722 .yiv1166828722bold a {text-decoration:none;}#yiv1166828722 dd.yiv1166828722last p a {font-family:Verdana;font-weight:700;}#yiv1166828722 dd.yiv1166828722last p span {margin-right:10px;font-family:Verdana;font-weight:700;}#yiv1166828722 dd.yiv1166828722last p span.yiv1166828722yshortcuts {margin-right:0;}#yiv1166828722 div.yiv1166828722attach-table div div a {text-decoration:none;}#yiv1166828722 div.yiv1166828722attach-table {width:400px;}#yiv1166828722 div.yiv1166828722file-title a, #yiv1166828722 div.yiv1166828722file-title a:active, #yiv1166828722 div.yiv1166828722file-title a:hover, #yiv1166828722 div.yiv1166828722file-title a:visited {text-decoration:none;}#yiv1166828722 div.yiv1166828722photo-title a, #yiv1166828722 div.yiv1166828722photo-title a:active, #yiv1166828722 div.yiv1166828722photo-title a:hover, #yiv1166828722 div.yiv1166828722photo-title a:visited {text-decoration:none;}#yiv1166828722 div#yiv1166828722ygrp-mlmsg #yiv1166828722ygrp-msg p a span.yiv1166828722yshortcuts {font-family:Verdana;font-size:10px;font-weight:normal;}#yiv1166828722 .yiv1166828722green {color:#628c2a;}#yiv1166828722 .yiv1166828722MsoNormal {margin:0 0 0 0;}#yiv1166828722 o {font-size:0;}#yiv1166828722 #yiv1166828722photos div {float:left;width:72px;}#yiv1166828722 #yiv1166828722photos div div {border:1px solid #666666;min-height:62px;overflow:hidden;width:62px;}#yiv1166828722 #yiv1166828722photos div label {color:#666666;font-size:10px;overflow:hidden;text-align:center;white-space:nowrap;width:64px;}#yiv1166828722 #yiv1166828722reco-category {font-size:77%;}#yiv1166828722 #yiv1166828722reco-desc {font-size:77%;}#yiv1166828722 .yiv1166828722replbq {margin:4px;}#yiv1166828722 #yiv1166828722ygrp-actbar div a:first-child {margin-right:2px;padding-right:5px;}#yiv1166828722 #yiv1166828722ygrp-mlmsg {font-size:13px;font-family:Arial, helvetica, clean, sans-serif;}#yiv1166828722 #yiv1166828722ygrp-mlmsg table {font-size:inherit;font:100%;}#yiv1166828722 #yiv1166828722ygrp-mlmsg select, #yiv1166828722 input, #yiv1166828722 textarea {font:99% Arial, Helvetica, clean, sans-serif;}#yiv1166828722 #yiv1166828722ygrp-mlmsg pre, #yiv1166828722 code {font:115% monospace;}#yiv1166828722 #yiv1166828722ygrp-mlmsg * {line-height:1.22em;}#yiv1166828722 #yiv1166828722ygrp-mlmsg #yiv1166828722logo {padding-bottom:10px;}#yiv1166828722 #yiv1166828722ygrp-msg p a {font-family:Verdana;}#yiv1166828722 #yiv1166828722ygrp-msg p#yiv1166828722attach-count span {color:#1E66AE;font-weight:700;}#yiv1166828722 #yiv1166828722ygrp-reco #yiv1166828722reco-head {color:#ff7900;font-weight:700;}#yiv1166828722 #yiv1166828722ygrp-reco {margin-bottom:20px;padding:0px;}#yiv1166828722 #yiv1166828722ygrp-sponsor #yiv1166828722ov li a {font-size:130%;text-decoration:none;}#yiv1166828722 #yiv1166828722ygrp-sponsor #yiv1166828722ov li {font-size:77%;list-style-type:square;padding:6px 0;}#yiv1166828722 #yiv1166828722ygrp-sponsor #yiv1166828722ov ul {margin:0;padding:0 0 0 8px;}#yiv1166828722 #yiv1166828722ygrp-text {font-family:Georgia;}#yiv1166828722 #yiv1166828722ygrp-text p {margin:0 0 1em 0;}#yiv1166828722 #yiv1166828722ygrp-text tt {font-size:120%;}#yiv1166828722 #yiv1166828722ygrp-vital ul li:last-child {border-right:none !important;}#yiv1166828722  
  ----------

         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 )

      include 'vaba_param.inc'

      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)
C
C props(1)=K
C props(2)=n
C props(3)=m
C props(4)=Sigma Yield
C props(5)=Elastic strain
      character*80 cmname
C
      do km = 1,nblock
        yield(km)= props(4)+props(1)*((props(5)+eqps(km))**props(2))*(eqpsRate(km)**props(3))
        dyieldDeqps(km,1)=props(1)*props(2)*((props(5)+eqps(km))**(props(2)-1))*(eqpsRate(km)**props(3))
            dyieldDeqps(km,2)=props(1)*props(3)*((props(5)+eqps(km))**props(2))*(eqpsRate(km)**(props(3)-1))
          end do
C
      return
      end

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

Reply | Threaded
Open this post in threaded view
|

Re: Vuhard

francisco.avs
Hi Wendell,

 Are you brazilian?
 

 please, try this and let me know if it'll be okay:
 

 

 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 )

include 'vaba_param.inc'

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)
C

 C
 real*8 :: val_K
real*8 :: val_n
real*8 :: val_m
real*8 :: Syield
real*8 :: Eel

 C
 C
 val_K = props(1)
val_n = props(2)
val_m = props(3)
Syield = props(4)
Eel = props(5)
 
character*80 cmname
C
 

 parameter ( one = 1.0d0)


 do km = 1,nblock


 yield(km)= Syield + val_K*((Eel+eqps(km))**val_n)*(eqpsRate(km)**val_m)


 dyieldDeqps(km,1)=val_K*val_n*((Eel + eqps(km))**(val_n - one)) * (eqpsRate(km)**val_m)


 dyieldDeqps(km,2)=val_K*val_m*((Eel+eqps(km))**val_n)*(eqpsRate(km)**(val_m - one))


 end do
C
return
end