Hi everyone,

I'm using UEL in ABAQUS. My model just has one 4-node element (defined in UEL) fixed in one end. When I'm running the UEL, it's working without any problm. But when I' trying to print out the global stiffness matrix I'm getting this error: "aborted with system error code 1073740940".

This is the way that I'm trying to print it out:

*STEP, NAME=STEP-1
*MATRIX GENERATE, STIFFNESS
*MATRIX OUTPUT, STIFFNESS
*END STEP

Would you please let me know what's the problem and how I can fix it?

Regards,

Hi everyone,

I'm working on UEL and studying wave propagation in beam. I have a special purpose beam element defined in UEL which has four nodes. When the beam is horizontal and a vertical load is being applied, it's working perfectly:

*Node
1, 0., 0., 0.
2, 0, 0.25, 0.
3, 0., 0., 0.
4, 0., 0.25, 0.

*USER ELEMENT,NODES=4,TYPE=U201,PROPERTIES=15,COORDINATES=3,VARIABLES=100, UNSYMM
1,2,6
*ELEMENT,TYPE=U201,ELSET=Set-1
1, 1, 2, 3, 4

*Cload, Amplitude=TB40KR
FR, 2, 1.

But, when the beam has a rotation (let's assume 45 degree wrt horizontal line), and the same load (vertical) is being applied, something is going wrong and I'm getting exactly the same answer with the horizontal one (they shouldn't be the same. It has been tested by using B32 elements from ABAQUS library).

*Node
1, 0., 0., 0.
2, 0.176776695, 0.176776695, 0.
3, 0., 0., 0.
4, 0.176776695, 0.176776695, 0.

*USER ELEMENT,NODES=4,TYPE=U201,PROPERTIES=15,COORDINATES=3,VARIABLES=100, UNSYMM
1,2,6
*ELEMENT,TYPE=U201,ELSET=Set-1
1, 1, 2, 3, 4

*Cload, Amplitude=TB40KR
FR, 2, 1.

I'm pretty much sure it will raise lots of questions in your mind about why I'm not using beam elements from ABAQUS library and etc. I should say that, I HAVE TO define and use my own elements. :(

It would be very kind of you to let me know what's wrong here.

Thanks in advance.

Best Regards,

Hi everyone,

I am working in UEL. I can get the stiffness matrix from fortran to be printed in .msg file.How can I print the global stiffness matrix of the model by using ABAQUS. My model is not big or complicated, so the global stiffness matrix shouldn't be too big!
Thanks!

Normal
0

false
false
false

EN-US
X-NONE
AR-SA

/* 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-parent:"";
mso-padding-alt:0in 5.4pt 0in 5.4pt;
mso-para-margin-top:0in;
mso-para-margin-right:0in;
mso-para-margin-bottom:10.0pt;
mso-para-margin-left:0in;
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-hansi-font-family:Calibri;
mso-hansi-theme-font:minor-latin;}

Hi All,

Recently I've got a problem with coordinate transformation matrix in ABAQUS.

Let's say you have a beam element with 2 nodes and 3 DOFs per node (imagine a Timoshenko beam). When you are rotating this beam to have some angle with horizontal line, based on DOFs and rotation angle, ABAQUS is applying a coordinate transformation matrix to the beam to give you the results in global coordinate system (transform between global and local coordinate system).

Is there any way that I can change/edit that coordinate transformation matrix that ABAQUS is applying or not?

I would appreciate any kind of help for this.

Thanks,

Hi everybody,

I'm an ABAQUS UEL user. I have Intel Parallel Studio 2013 installed on my computer with Visual Studio 2013 and ABAQUS 6-12.3

I'm almost a begginer in this area. I realized that Intel Composer XE 2013 (which has been intalled on my machine) contains MKL which has BLAS and LAPACK libraries in it.

I tried to use some of their functions (let's say eigenvalue solver and SVD solver). When I'm
trying to run it (abaqus job=test user=test1 int) I'm getting this error
message:

LNK2019: unresolved external symbol dgemm referenced in function uel

LNK1120: unresolved externals

To me, it seems that the problem is in linking between ABAQUS, MKL and
LAPACK. If I'm right, would you please let me know how I can resolve it?

I really appreciate if you give me a step by step solution. I'm under pressure due to an upcomin deadlline.

Best Regards,

Hello everyone,

I want to simulate contact problem by using UEL.I define the surface in the inp file.But abaqus don not

recognize it. How can I solve this problem?

PS. Sorry for my English.I hope it is clever enough to understand the problem.

Regards.

J.H. Zhai

Hi everybody,

I'm working on UEL. The stiffness matrix (AMATRX) which I have
defined is consisting of complex numbers (real and imaginary part).

When I'm running it I'm getting this error message:

Abaqus Error: The executable standard.exe aborted with system
error code 1073741819. Please check the .dat, .msg and .sta files for
error messages if the files exist. ...

I checked the files and no error message was found.

When stiffness matrix is just real numbers, it's working perfectly.

Would you please let me know what should I do?

Thanks in advance.

Hi everybody,

I have a UEL for truss element (T3D2). I'm calling and reading some .dat files inside of the UEL. I'm just reading these files inside of the UEL without doing anything to them or with them (It's just a test case). Here is an example:

OPEN(101, file="C:\Users\ak1070\Desktop\UMAT-UEL-UELMAT\UEL TEST\V1.dat", status='old',form='formatted',&
&action='read',position='rewind')
READ(101,*) AEVectors(1,:)
CLOSE(101)

.

.

.

OPEN(164, file="C:\Users\ak1070\Desktop\UMAT-UEL-UELMAT\UEL TEST\V64.dat", status='old',form='formatted',&
&action='read',position='rewind')
READ(164,*) AEVectors(64,:)
CLOSE(164)

I'm supposed to call and read 64 files. But it's just working for 38 files. When I'm calling 38 files it's working perfectly. I Mean:

OPEN(101, file="C:\Users\ak1070\Desktop\UMAT-UEL-UELMAT\UEL TEST\V1.dat", status='old',form='formatted',&
&action='read',position='rewind')
READ(101,*) AEVectors(1,:)
CLOSE(101)

.

.

.

OPEN(138, file="C:\Users\ak1070\Desktop\UMAT-UEL-UELMAT\UEL TEST\V38.dat", status='old',form='formatted',&
&action='read',position='rewind')
READ(138,*) AEVectors(38,:)
CLOSE(138)

But, when I'm calling 39 files the results change to sth weird. I thought maybe sth is wrong with data so, I called one of the previous files twice to make the total number of files 39 and still the results change to sth weird.

OPEN(101, file="C:\Users\ak1070\Desktop\UMAT-UEL-UELMAT\UEL TEST\V1.dat", status='old',form='formatted',&
&action='read',position='rewind')
READ(101,*) AEVectors(1,:)
CLOSE(101)

.

.

.

OPEN(138, file="C:\Users\ak1070\Desktop\UMAT-UEL-UELMAT\UEL TEST\V38.dat", status='old',form='formatted',&
&action='read',position='rewind')
READ(138,*) AEVectors(38,:)
CLOSE(138)

OPEN(138, file="C:\Users\ak1070\Desktop\UMAT-UEL-UELMAT\UEL TEST\V38.dat", status='old',form='formatted',&
&action='read',position='rewind')
READ(138,*) AEVectors(38,:)
CLOSE(138)

It would be very kind of you to let me have your advice.

Thanks in advance.

Best Regards,

Hi every body,

I'm a new user of UEL and have a questions. It would be very kind of you to help me with it.

I have a MATLAB code which gives me some big matrices. Is it
possible to save them as a file (.txt, .dat, etc) and call them inside
of the UEL? (It's almost impossible to change all of those codes to
FORTRAN or copy-paste them to UEL.)

Thanks in advance,

Best Regards,

-------------

Ashkan

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

Hi there,

Recently I had lots of problem with linking ABAQUS and Fortran compiler together. I did all of the things which has been described in forums but none of them worked for me. Finally I found the solution.

I'm using ABAQUS 6.12-3, Visual Studio 2012 Professional and Intel Composer XE Suites 2013. For linking them simply follow these steps:

1- find the file ifortvars.bat and its path.

2- startmenu --> Programas --> right click on ABAQUS Command --> Shortcut --> add the path of ifortvars.bat file inside of quotation mark to the very end of the target --> make an space and add this: intel64 vs2012 -->save and exit.

in my case it is: "C:\Program FIles (x86)\Intel\Composer XE 2013 SP1\bin\ifortvars.bat" intel64 vs2012

I hope this method works for you as well.

Hi all,

I have been trying to model the surface interaction between two blocks of material using surface-based cohesive elements (*cohesive behavior). This works well in ABAQUS CAE (standard) when I set the material properties as "Mechanical, Elastic" so no problem with that.

However, I would now like to use the material properties from a fortran code. To do this, I was planning on defining the elements using *UEL with an overlay to record all the values. Is this compatible with cohesive zone elements?

Many thanks for your help,

Lucy

Hello Experts,

I am trying to do modal analysis on user-defined piezoelectric shell element. My eigen value probelm of free vibration is given in the attachment.

My issue is Mass matrix [M] is not positive definite due to electrical degrees of freedom having no mass associated with them and the stiffness matrix [K] given is non-positive definite due to negative dielectric terms. I know , it is possible to condense out the above equations and find the real eigen values and eigen vectors but how to do this in UEL subroutine using Lanczos solver. I introduced arbitrary/artificial mass terms in Mass matrix to get the eigen values and eigen vectors. Is there any other appropriate method? / ideas?

I need to resolve this issue asap. I would really appreciate any help regarding this.

Thanking you in advance,

Somu

Hello Experts,

I have written a uel subroutine where I have totally 12 degree of freedom. Among these 12 degree of freedom, 6 are standard abaqus degree of freedom (3 Displacements and 3 Rotations) and another 6 are User defined degree of freedom.

I assign these these degree of freedom in following way.

*USER ELEMENT,TYPE =U9, NODES=9,PROPERTIES=24,coordinates=3
1,2,3,4,5,6
9,1,2,3,4,5,6,7,8,9,10,11,12

Although above mentioned degree of freedom like 9 is potential and 11 is temperature. I think I can overide them in UEL. Correct me if I am wrong.

I can get the output at the nodal degree of freedom 1,2,3,4,5,6,9 by using following input statement.

*node print
u
epot

But how to get the output at the other user-defined degree of freedom in .dat file i.e 7,8,10,11,12

I tried with nt, it didn't print the output for dof 11 as per abaqus.

I didn't find in abaqus manual any information on this?

Please help me!

Thanks and Regards,

Somu

Hello Experts,

I am trying to extract the eigenvalues of the User-Defined Element. I am getting following error when I try to extract frequencies using Lanczos method

*Frequency, eigensolver=Lanczos, acoustic coupling=off, normalization=displacement
3,,,, , ,

Please help me to resolve this error. Please find the message file in attachment.

SHIFT IS DETERMINED IN ORDER TO CLOSE GAP BETWEEN TRUST INTERVALS
NEW SHIFT 1.00000e+000

NUMBER OF NEGATIVE EIGENVALUES 11
NUMBER OF POSITIVE EIGENVALUES 197
NUMBER OF NEGATIVE EIGENVALUES NOT ASSOCIATED WITH LAGRANGE MULTIPLIERS 11
NUMBER OF NEGATIVE EIGENVALUES ASSOCIATED WITH LAGRANGE MULTIPLIERS 0

***ERROR: NUMBER OF EIGENVALUES FOUND ON THE INTERVAL EXCEEDS THE EXPECTED
NUMBER. THIS MAY HAPPEN IF A SYSTEM IS OVERCONSTRAINED. RUN A STATIC
PERTURBATION ANALYSIS TO PERFORM THE OVERCONSTRAINT CHECK. ALSO THIS
MAY HAPPEN IF A SYSTEM HAS RIGID BODY MODES AND NO SHIFT IS GIVEN

Hello Experts,

I have written a Uel subroutine and trying to implement on a plate constrained at all the edges.

My UEL Subroutine executes without any errors when the total number of elements are just 15 to 20. But when I increase the number of elements and mesh it more fine, My Uel subroutine aborts with ERROR Code 5.

If I print the JELEM, UEL Subroutine is aborting always at 52th Element, i.e when JELEM=52. Even when I remesh the structure with different size, Uel subroutine consistantly aborts for the 52ed element.

Does any one have any idea what could be the reason? My Uel subroutine works fine on simple structures like beam fixed at one end.

Morever, there aren't any error messages inside .dat file and .msg file.

Abaqus Error: The executable C:\Programme\SIMULIA\Abaqus\6.9-EF2\exec\standard.exe
aborted with system error code 5.
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. Use the same command to run Abaqus that you
used when the problem occurred. Please contact your local Abaqus support
office and send them the input file, the file support.log which you just
created, the executable name, and the error code.
Abaqus/Analysis exited with errors

Thanks and Regards,

Somu

Hello All,

Does anyone have any idea how to automate meshing using USER DEFINED ELEMENTS (for Example UEL) in Abaqus using python. Actually I need to mesh an arbitrary shaped structure with 9-Node Quadrilateral Element. Abaqus doesn't support meshing of 9 Node Quad Element (S9R) either.

One idea would be to generate mesh with 8 Node Quad Element (S8R) and then interpolate the 9th node. Do file handling operations in Abaqus python to modify the input file for 9 Node Quad element. Is there any other better workaround for this? I am sure many would have come accross this situation of meshing of UEL element.

Any advise/help is greatly appreciated !

Thanks and Regards,

Somu

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

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

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

Hi

I have made a user element, UEL of a cohesive zone element. which is a 8 noded brick element with zero thickness. After solving models using this element I would like postprocess the model in abaqus CAE 6.10 in the way that I am able to:

1. Plot element and node labels
2. Plot undeformed user elements together with the rest of the elements model
3. Plot deformed user elements together with the rest of the elements model
4. Plot any of the SVARS values element wise

I have searched for an answer to this on existing forum discussions and in the abaqus manual but I did not find any solution.

Thanks in advance for you help.

Best regards

Brian Bak

Dear Friends,

I have written a UEL for Strain Gradient Plasticity and is working well when i give Boundary condition as Concentrated force, But when i apply displacement boundary condition the it leads to 'system error 144'.

Iknow that system error 144 corresponds to division by zero or infinity, but then why it is coming in on displacement boundary condition is i am worried about.

I will be extremely helpful to anybody for any type of help on above problem.

Thanking you in advance

Soham

Dear Friends,

I have written an UEL for Strain Gradient Plasticity but it shows following error :

FE SOLVER AND MESH HAVE INCOSISTENT NUMBER OF DOFS

I will be very thankful if someone can throw some light on this error.

Thanking in advance.

Soham