RE: How to update YoungÆs modulus according to change in plastic strain

I have checked the USDFLD subroutine in the PDF file. Becouse I am new to write user subroutine in ABAQUS, I couldn't change that for this special problem. Would you please put the subroutine here and let me know how to use that?

The other thing that i should mention is that the analysis that I want to use the subroutine is in ABAQUS STANDARD, please let me know if there is any differences between ABAQUS STANDARD and Explicit to use this subroutine.

Sorry again because of this request. I have tried a lot but I couldn't use that.

C,A

RE: How to update YoungÆs modulus according to change in plastic strain

I made for You very simple example + USDFLD subroutine used by the model.Please check attached *.zip file.

A material change Young module respect to elastic strain.The idea is to use field variable, since meaning of field variable is set by a user in the example field variable are elastic strains.The subroutine read EE11 output for each integration point and save it as field variable, next abaqus base on the field variable value choose Young module.The example is with truss elements since there is only one direction for stress and strain so it is easy to check results.

Of course the same approach you can use with plastic strain, You need to change only subroutine to read plastic strain.Just add new output in the model ('PE') and in the subroutine read 'PE' output instead 'EE'.For plastic strain outputs will be save in 'array' table in following order:array(1) = PE11, array(2) = PE22, array(3) = PE33array(4) = PE12, array(5) = PE13, array(6) = PE23array(7) = PEEQ, array(8) = PEMAG

>>> please let me know if there is any differences between ABAQUS STANDARD and Explicit to use this subroutine.USDFLD is for Standard and VUSDFLD is for Explicit.Usually we can find small difference in subroutines for explicit and standard. I do not know is any differences in this case. You need to check it in documentation.

RE: How to update YoungÆs modulus according to change in plastic strain

I have performed the analysis which is a nonlinear plane stress problem and Young's Modulus (E) is changed with Equivalent Plastic Strain (PEEQ). To validate the model I want to plot PEEQ versus E. I couldn't find a relation between them, so I have made a new post about that. I made a decision to write this question here too, to ask you about that. Please let me know if you know how to get E in a nonlinear plane stress problem to plot PEEQ versus E.

Cheers,A

RE: How to update YoungÆs modulus according to change in plastic strain

RE: How to update YoungÆs modulus according to change in plastic strain

RE: How to update YoungÆs modulus according to change in plastic strain

I want to define a material in which the Young's modulus (E) can change trough thickness. I wrote a USDFLD subroutine for that to use a field variable. Would you please help me solve the problems that I have as below:

• As you can see in the attached file, I want to change the E in thickness direction (Y). So I defined COORD (2) as a field variable. I think this way is not correct because the amount of coordinate (2) will change by deformation of sheet, but I want to change the E trough thickness. How can I define the field variable for this problem?

• Instead of writing exact numbers for E, I want to have a relation between E and field variable (1) like:E=69600(1+Field (1)/2.5)How can I use this equation?

I appreciate if you please let me know if there is a better way to change E with thickness and this way is not correct.

RE: How to update YoungÆs modulus according to change in plastic strain

I've prepared a pack for you with my solution. The pack includes *.inp and *.for files.This is your model with small modification (new step, new definition for field variable in *ELASTIC section).Please take a look for it since I put some explanations in comments.

c do it for all others inc. and steps elsec assign initial values calculated in first inc. field(1) = statev(1)

end if

The main idea is to:1. define field variable as E instead thickness as in your case.2. point 1 allow us to calculate E respect to your equation.3. thickness (y-coordinate) I read for material points not for nodes. It easier since coordinates for material points are pass into subroutinewith dummy arguments lists (coord array) and we do not have to use GETVRM routine to read node coordinates.4. E is calculated only for first increment when all elements are in initial positions, for all others increment I assign initial values.So I have no influence of element deformation for my material data during solution.It is true only if in first increment of my analysis I have no loads which could deform geometry.5. remember that USDFLD subroutine is call at the end of an increment, it means E values assign for elements are not correct for first increment.To solve it I see two possibilities:a) add one artificial step with one increment where we do not generate any stress/strains at the beginning of the analysis.b) set initial values for field variable.Take a look for abaqus documentation 30.2.1 Initial conditions in Abaqus/Standard and Abaqus/Explicit (Defining initial values of predefined field variables).I think the best is to use external odb file to read initial field variables values. In this case You do not need USDFLD subroutine for other runs with the same part.You need to only run one analysis with subroutine to set field variables.

Value of E assigned for each element you can plot in Viewer with FV1 field output.

RE: How to update YoungÆs modulus according to change in plastic strain

You put wonderful comments for me. Regarding to you comments would you please answer my questions as below:

• What do you mean by material point coordinate? Is it the coordinate of integration point?• When I wanted to plot FV1 to see E for each element the minimum and maximum were 92800 and 185600 respectively. In the equation for E, the minimum and maximum should be 69600 and 208800 (69600*3). I thought maybe we don't have enough elements trough thickness, so I increased the number of element to 50 trough thickness!! But the minimum and maximum are 70990 and 207400 respectively which are not exact yet. Would you please let me know what is the reason that we can not get the exact amount for E in this analysis?

Regards,A

RE: How to update YoungÆs modulus according to change in plastic strain

Quote:

What do you mean by material point coordinate? Is it the coordinate of integration point?

Yes, you are right, this is the same points.

Quote:

When I wanted to plot FV1 to see E for each element the minimum and maximum were 92800 and 185600 respectively.

The values are correct. The equation for E is calculate respect to y-coordinate of material point.For your mesh there are 3 material points in y-direction (in middle of element height):1. For bottom elements it is 0.83333 -> E=69600*(1+0.83333/2.5) = 928002. For middle elements it is 2.49999 -> E=69600*(1+2.49999/2.5) = 139200 3. For top elements it is 4.16666 -> E=69600*(1+4.16666/2.5) = 185600

Quote:

In the equation for E, the minimum and maximum should be 69600 and 208800

You can get these values only for node position. It does not make sense, E value you can assign only for elements not for nodes.

One more point from my side, I forgot about this in previous post.The subroutine works ok for reduced elements (with one integration points) if you want to use fully integrated elements then you need to modify the subroutine. I believe the best solution is to calculate E value for all material points in an element and assign one average value.

Regards,Bartosz

RE: How to update YoungÆs modulus according to change in plastic strain

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Resources

ENGINEERING.com surveyed 364 product development professionals to glean insights into how their teams are performing today, and what technologies are most useful to support their processes. Download Now

3D printing has quickly expanded beyond the realm of prototyping and into the world of end part manufacturing. Desktop 3D printers need not be used for making parts themselves, but as a secondary process for producing tools. Download Now

Creating tooling for injection molding is an expensive and complex proposition. Until recently, the only way to avoid tooling defects causing air traps, voids, shrinkage and the like was to employ a time-consuming and tedious physical prototyping process. Download Now