Quantcast
Channel: iMechanica - User subroutines in ABAQUS
Viewing all 63 articles
Browse latest View live

About Dflux of abaqus

$
0
0

I would like to use Dflux to give a Temperature and dT/dt related heat flux. Can Dflux be used for this purpose? I am not sure how to get dT/dt ? Or other subroutine can fulfill this purpose?

 

Thank you 


A Very silly abaqus question... applying a wind load

$
0
0

Dear All,

Sorry to ask such a silly question but I am quite new to Abaqus... I am trying to model a tower exposed to wind. I know how the wind pressure varies according to height, but what I dont know is how to model it in abaqus... Also, if i apply a pressure load, then the load will act normally at all points on the surface, right? But in actuality wind doesnt act that way... Could anyone help me out on this one?

 

Thanks and Regards,

Lekshmi 

Abaqus "system error code 884"

$
0
0

Hello again

 

Im running a simulation on Abaqus and i got the next text message:

 

Abaqus Error: The executable C:\SIMULIA\Abaqus\6.11-PR3\exec\standard.exe 

aborted with system error code 884. 

Please check the .dat, .msg, and .sta files for error messages if the files 

exist.  If there are no error messages and you cannot resolve the problem, 

please run the command "abaqus job=support information=support" to report and 

save your system information

 

 do you know what is this error? ive been looking on ABAQUS documentation and in the web and i cant find information

 

Thanks for your help

 

Juan Carlos Martinez

 

How to use USDFLD request variables from other elements in ABAQUS

$
0
0

Hi all,

 I want to use USDFLD to have diffusion and mechanical coupling analysis. The model is building a cohesive zone between two parts. Before the calculation is  executed in each increment, I need to get the stress information of last increment of one element and the other two element next to it in the same cohesive zone. After geting this, I can update the parameters of cohesive zone and a following mechanical analysis will be executed. How can I save the variables of desired elements for each increment? Or how can I request other elements' information in the USDFLD subroutine ?

Non-uniform expansion coefficient alpha in UEXPAN (ABAQUS)

material points will change while using subroutine in ABAQUS

$
0
0

This subroutine using VUSDFLD is to
determine position of damage and provide evolution of damage surface, when some
elements failed. The “outer” damage status will transfer to “inner” damage
status.

Now the problem is that when I read the
coordinates of material points, using coordMp
k,1),coordMpk,1),coordMpk,1), it returns wrong locations. Then I
tried a very simple loop:

Do
k=1,nblock
        stateNew(k,7)=nblock
end Do

I found the number of material points
(namely, nblock) varied! First time
I ran this loop, I got the contour of state variable 7 (SDV7) shown below. The
SDV7 were various. But from the subroutine they were supposed to be same.

 

Then I coarsen the mesh, and get a
different result, SDV7 of every element were same (equal to the material point
number=125)! Shown below:

 

Then I substituted coordMp(1,1) for nblock
in the above-mentioned loop, and found similar phenomenon.

 

I wonder if the number of material points
and their coordinates will vary in a loop? Those results are beyond my expectation.
Is that Fortran grammar problem or misunderstanding on subroutine in Abaqus?

 

I attached the subroutine "position.for" and its "Job-1.inp" below. But I switch their extensions into .txt for uploading.

Thank you all.

Calling a File inside of the UEL

$
0
0

Hi everybody,

I'm writing a User Difined Element (UEL) for my research.

Would you please let me know if it's possible to call a file (txt, dat, etc) inside of the UEL or not?

Example:

    IF (LFLAGS(3)=1) THEN
        OPEN (UNIT = 10, FILE = "AAA.dat", STATUS = "OLD")
        READ (10, *), StiffnessTable
        AMATRX=StiffnessTable
    END IF

Thanks in advance,

Ashkan

Sample code for Abaqus subroutines


Properties of Polycarbonate in Abaqus

UEL visualization and strain energy for UEL

$
0
0

I am only seeing a "x" for a UEL at its first node. And I dont see any other node of UEL. Is it correct or some problem is there?
I also want to calculate the strain energy for whole model(including UEL). But in history output I am only seeing strain energy of Abaqus defined elements under "whole model strain energy". Can some one please help me in calculating strain energy of whole model including UELs 

deformation gradeint in VUMAT Abaqus Explicit subroutine

$
0
0

Hello everyone,

I'm modelling VUMAT for a composite. So, in my model I have inotropic and anisotropic behaviour. I have modelled isotropic behaviour but for anisotropic behaviour I need to use deformation gradient, which is causing errors. CAN ANYONE HELP with DEFORMATION GRADIENT. How to use that, do we need to rotate stress back in corotational form as abaqus manual says if yes then How?

How can we get Rotation matrix in VUMAT, if we can multiply it with stretch I think we have the result. May be an alternative.....

Please help, stuck badly :(

 

Thanks in advance

Troubles with ABAQUS subroutine compiler after update of Visual Studio 2015

$
0
0

I always use "Microsoft Visual Studio 11.0" and "Intel\Composer XE 2011 SP1" to apply all the user subroutines.

Recently, the university administrator updated the Visual Studio to 2105 (14.0) on my lab computer. 

When I tried to run an OLD simulation which WORKED perfectly before, it returned me the error message below:

 

Begin Linking Abaqus/Standard User Subroutines

   Creating library standardU.lib and object standardU.exp

standardU_static.lib(SMAStsUsubs_IF.obj) : error LNK2005: DllMain already defined in MSVCRT.LIB(dllmain.obj)

standardU.dll : fatal error LNK1169: one or more multiply defined symbols found

Abaqus Error: Problem during linking - Abaqus/Standard 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.

Abaqus/Analysis exited with errors

 

Although Visual Studio has been updated, I still called the "Microsoft Visual Studio 11.0" and "Intel\Composer XE 2011 SP1".

Can anyone help me on this issue?

Thank you very much!

-Ninggang

 

Abaqus User defined Routines

$
0
0

I have been using an Abaqus dll created using a user-defined routine. Things work well on a serial case. The issue comes when I try to use it in parallel. I am trying to do this on windows 7. The error I get is

forrtl: severe (59): list-directed I/O syntax error, unit 12,

in Fortran I/O but I do not get this error when I run the simulation serially. Any ideas will be helpful.

Thanks

 

some errors still exists for me during importing subroutine file in abaqus/cae

$
0
0
AttachmentSize
Image icon2.jpg131.2 KB
Image icon3.jpg117.88 KB
Image icon11.jpg96.42 KB

Thank you asim sir,

I am sorry to say that i am still getting the error message.I have attached the picture of  procedure and error  messages i gets.please correct my procedure.

 

I use

Normal
0

false
false
false

EN-IN
X-NONE
X-NONE

MicrosoftInternetExplorer4

/* Style Definitions */
table.MsoNormalTable
{mso-style-name:"Table Normal";
mso-tstyle-rowband-size:0;
mso-tstyle-colband-size:0;
mso-style-noshow:yes;
mso-style-priority:99;
mso-style-qformat:yes;
mso-style-parent:"";
mso-padding-alt:0cm 5.4pt 0cm 5.4pt;
mso-para-margin-top:0cm;
mso-para-margin-right:0cm;
mso-para-margin-bottom:10.0pt;
mso-para-margin-left:0cm;
line-height:115%;
mso-pagination:widow-orphan;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-ascii-font-family:Calibri;
mso-ascii-theme-font:minor-latin;
mso-fareast-font-family:"Times New Roman";
mso-fareast-theme-font:minor-fareast;
mso-hansi-font-family:Calibri;
mso-hansi-theme-font:minor-latin;
mso-bidi-font-family:"Times New Roman";
mso-bidi-theme-font:minor-bidi;}

C:\SIMULIA\Documentation\docs\v6.8\books\eif\ufricxxy.f

as the inpu file

Normal
0

false
false
false

EN-IN
X-NONE
X-NONE

MicrosoftInternetExplorer4

/* Style Definitions */
table.MsoNormalTable
{mso-style-name:"Table Normal";
mso-tstyle-rowband-size:0;
mso-tstyle-colband-size:0;
mso-style-noshow:yes;
mso-style-priority:99;
mso-style-qformat:yes;
mso-style-parent:"";
mso-padding-alt:0cm 5.4pt 0cm 5.4pt;
mso-para-margin-top:0cm;
mso-para-margin-right:0cm;
mso-para-margin-bottom:10.0pt;
mso-para-margin-left:0cm;
line-height:115%;
mso-pagination:widow-orphan;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-ascii-font-family:Calibri;
mso-ascii-theme-font:minor-latin;
mso-fareast-font-family:"Times New Roman";
mso-fareast-theme-font:minor-fareast;
mso-hansi-font-family:Calibri;
mso-hansi-theme-font:minor-latin;
mso-bidi-font-family:"Times New Roman";
mso-bidi-theme-font:minor-bidi;}

C:\SIMULIA\Documentation\docs\v6.8\books\eif\ufricxxy.f

 as subroutine file

 

 

UEL ABAQUS Subroutine

$
0
0

Hi,

I'm stucking with my subroutine. May ask anybody to have a look the the file. 

I run a simple model that and try to run it,but an error message was appeared as below;   ***ERROR: USER SUBROUTINE UEL MISSING

 

I don't know what is the problem. I did locate the input (.inp) and user (.f) in the same folder and run it (abaqus job=.... user= ...ff int )

 

Hope to get help from anybody...

 

Aida


ABAQUS VUHARD subroutine problem

$
0
0

Hello

I am working on the simulation of ring rolling. I have written the subroutine for the plasticity part in the subroutine VUHARD and have set the PLASTIC part in material definition to USER.
I dont have HARDENING PROPERTIES so i have set it BLANK. i have set the path of the subroutine in the job defintion part too. But abaqus crashes whenever i try to do DATA CHECK
It says ABAQUS KERNEL EXITED WITH AN ERROR.I dont know why this problem arises?. I am stuck.

The curve i want to define is yieldstress=K*plasticstrain^m*plasticstrainrate^n

I have defined the yield,Dyieldeqps(k,1)Dyieldeps(k,2) as abaqus has told.
Abaqus has told to define Statenew too.I dont know how to define it.

This is my subroutine :

c-----subroutine for plasticity
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.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)

parameter( one=1.0d0,small=1.0d-8)
C
character*80 cmname

double precision a,b,c

k=100.0d0
m=0.30d0
n=0.15d0
C
C
do km = 1,nblock

plstrain=eqps(km)
plstrnrt=eqpsRate(km)

yield(km)=k*(plstrain**m)*(plstrnrt**n)
dyieldDeqps(km,1)=k*m*(plstrain**(m-one))*(plstrnrt**n)
dyieldDeqps(km,2)=k*n*(plstrain**m)*(plstrnrt**(n-one))

end do
C
return 
end

 

I would be grateful if someone could help me figure out the error. 

Umats with Abaqus

Storing Abaqus sensor values

$
0
0

Hello

I have a sensor in VUAMP subroutine to measure the diameter of a ring. I want to compare the current value with the sensor value 10 seconds ago.

Is there a way to store the sensor values and retrieve them in abaqus?

 

 

Thanks

ravikumar. 

How to update solution-dependat variable (SVARS) during iterations

$
0
0

Hello everyone,

 I am implementing a uel user subroutine in ABAQUS. This uel is for modeling contact in X-FEM framework. Thus, it's a nonlinear problem and I need to store some solution dependent variables during each iteration. I know I can store SDVs in SVARS. But this vector will update at the end of each step, and not during iterations. Anybody knows about this problem?

 I really appreciate if anybody can help me with this.

    Thanks

    Ehsan 

Modelling moving discontinuity using XFEM

$
0
0

Greeting to all iMechanicians,

I am trying to model motion of discontinuity (grain boundary) using XFEM as it saves me from the computational expense required for mesh adaptation. The material behaviour on the both sides of the discontinuity is different. After any step the discontinuity might cut through an element. Then to compute the stiffness of the element I will have to integrate seperately for the splitted regions and then add them. This is easy when I have a 1-D case but seems difficult for 2-D and 3-D case. Since the discontinuity may take any arbitrary shape, I will have to find out the points in the element that lie on either side of it for which I can use signed distance function. But how do I implement this in ABAQUS. I have a vague idea that I might have to use UEL or UMAT. If somebody has faced a similar situation, please help me out.

Best regards,

Luv

 

Viewing all 63 articles
Browse latest View live