Skip to main content

Andreas Burger's blog

Abaqus implicit, rigid body-deformable body contact, negative dot product

Submitted by Andreas Burger on

Hello,



I am running a simulation where I have a rigid body (shell elements) pushing into a deformable part.  Unfortunately, the simulation aborts with the following error:



The dot product -0.19758 cannot be negative for element 30050 instance xzy with record ID 844 at step 2 inc 1 kitgen 5.



This element (C3D4H) belongs to the deformable part and is within the body, so no not close to the contact at all.  At the point where the simulation aborts, this element is not deformed and the stress is basically 0 (E-6).

Bergstroem Boyce Model - inelastic part

Submitted by Andreas Burger on

Hi, 

 I was wondering if anyone has experience in implementing the Bergstroem-Boyce Model (Constitutive modeling of the time-dependent and cyclic loading of elastimers and application to soft biological tissue, 2001)? I am struggling with the implementation of the viscoplastic part (Network B) of the model:

NTVDM CPU has encoutnered an illegal instruction

Submitted by Andreas Burger on

Hi,

I'v been using the Abaqus command window to run my jobs without any problems.

Today, however, I suddenly got the following message:

Intel Visual Fortran Compiler 9.1 Build Environment for 32-bit applications - abaqus job=eightchain user=eightchain_Jaumann

The NTVDM CPU has encountered an illegal instruction.

CS:117a IP:0204 OP:63 65 6c 27 3a Choose 'Close' to terminate the application

The header of the message says: 16 bit MS-DOS Subsystem

 Has anyone had a similar problem before???

 Thanks,

writing data to dat file and plotting values during job execution

Submitted by Andreas Burger on

Hi,

I am in the process of debugging a UMAT.  In my UMAT code I am currently using write statements to output certain variables to the dat file.  I was wondering if there are any other ways to output data during the execution of a job.  Is it possible to print/plot variables, for example using Matlab (maybe even realtime)?

Also, has anyone experience in accessing the fill files (e.g converting them to ASCII and then use Matlab to look at data?).

Thanks,

Andreas

Pressure over closure function for Polyurethane or (nearly) incompressible materials in general

Submitted by Andreas Burger on

Hi,

Is anyone aware of a pressure-overcosure function for (nearly) incompressible materials.

 I am aware of functions for PE and articular cartilage but they do not work with v=0.5.   I have been trying to find similar functions for elastomers, but can't seem to find anything.

 I am grateful for any input.

 Andreas

UMAT: Calling a another subroutine in UMAT, Errors during job execution

Submitted by Andreas Burger on

Hi,



I am in process of implementing the Bergstroem-Boyce Model. In a first
step I am just considering the hyperlastic part of it (Network A:
8-chain). I have implemented this part in a UMAT using a numerical
approach for the elasticity tensor (Forward Difference Approximation
with perturbation of F).



I have two problems:



1.

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

In order to keep my UMAT as easy to read as possible I have tried to keep the FDA in a separate file using



call KFDA_Jaumann(dfgrd1,......)

Perturbation technique based on forward difference approximation

Submitted by Andreas Burger on

Hi,

 I am trying to implement (in Abaqus) the elasticity tensor using a numerical approach.  I have come across a paper by C.Miehe (Numerical computation of algorithmic (consistent) tangent moduli in large-strain computational inelasticity, 1996).

There are a few things I don't quite understand.  

1 PK, 2PK vs true, nominal stress

Submitted by Andreas Burger on

Hi,

 if I have a material model defined in cauchy stress then I use test data with true stress (calculated from nominal stress: sigma_true=sigma_nom(1+epsilon_nom).

If I wanted to use nominal stress data, which stress tensor would I have to use to define my material model? First or second Piola Kirchhoff?

 Andreas