ilinca's picture

Post-doctoral position (computational mechanics) at Rice University

Post-doctoral position available in computational mechanics at Rice University. The announcement and instructions to apply are posted at

simulate inclusion-matrix debonding under tensile load with abaqus cohesive element

Dear all:

  I want to simulate the  inclusion-matrix debonding under tensile load with abaqus cohesive element.But when I added the cohesive element around the inclusion, the stress filed was totally changed. I think it is because the stiffness of the cohesive element (K/G1/G2)was not right. Do you know how to choose it so it will not affect the stress field of the oringinal structrure? or can you give me an example? Thanks very much.

Modelling Delamination in composite material using cohesive elements

Hi All,

Basic problem in ABAQUS/CAE

I am simulating a cube of sides 10 units. I have applied the
load as pressure of 16 units at top. The cube is meshed by usign a
global seed size of 5. It means, there are 9 elements and 27 nodes
(using C3D8R). I have used the material property of steel for instance.
The base of the cube is fixed and the sides are allowed to move
vertically.  Now i want to calculate the nodal force at the centre of
the cube i.e. at the centroid of the cube (not in the centre of any of
the exterior sides). 


The question about ABAQUS element: CPS4R

The element is added displacement load in one direction, at the steptime=0, the nodal displacement is zero, and so as the integrated points' strain and stress. However, the internal force NFORC is nonzero. When we calculate internal force by ourselves, it is obtained by the integral for BT*stress. But the stress is zero, how Abaqus calculates the nonzero initial internal force?

The details can be found  in the attachment.

Abaqus Eigen Value Analysis Error


I am having issues running an Eigen value simulation on a bridge model that I have developed. The model consists of a bridge super-structure supported by pinned connections on the ends of the bridge length and a roller support in the middle.

The model runs without errors when I do not request an Eigen value analysis, but aborts when I request an Eigen analysis. The message file sites the cause as a factorization error. Warnings also appear stating a zero pivot when releasing the 1st dof of the nodes associated with the roller connection.

I have attache the .inp file of my model. Do you know what may be causing this error?




Mechanics journal impact factors 2012

Some funding agencies require that articles be published only in journals with impact factors greater than 1. Here's a list from 2012 that tells you which to choose (if you care about impact factors).
-- Biswajit


simply supported boundary condition in Ansys Autodyn

How could i apply simply supported boundary condition in Ansys Autodyn?

SIMULIA's picture

Interested in learning how to automate multi-step simulation workflows?

Join us on Thursday for this FREE eSeminar and learn how Learn how Isight could be used to:

•Automate a series of functions to create a Sim-flow
•Add components to a Sim-flow
•Execute a Sim-flow
•Visualize Sim-flow results
•Assess the impact of the various design variables on the objectives
and identify significant interactions among the design variables
•Perform material model calibration

Register Today!

Abaqus:warning message,Tangent directions used to enforced contact constraints for contact pair

Well I tried to search for a related thread and could not find one, so creating a

Announcing Wear of Materials 2015 - Website has just gone live

20th International Conference on Wear of Materials
12-16 April 2015
The Sheraton Centre, Toronto, Canada

The 20th International Conference on Wear of Materials will focus on both the fundamental and applied aspects of wear and friction of materials at the macro-, micro-, and nano-scale. It will address the understanding of tribological phenomena; particularly the progress in recent decades and a special session will concentrate on modeling of wear.

Paper: On finding cohesive law parameters for foam-metal interfaces

A detailed report of our work on trying to find cohesive laws for interfacial fracture is attached.  I will add the experimental data after a shorter version of this report has been published.

The reason for posting this report on iMechanica is so that people who want to replicate the tests or perform similar tests know what is involved.  As we have found out, many details and potential problems faced by an experimentalist are hard to determine from the concise papers that are typically published as journal articles.

Please let me know if you have any questions.

-- Biswajit


 On finding cohesive law parameters for foam-metal interfaces


Faculty Position - Computational Materials Science

The Department of Mechanical Engineering
at the University of Wyoming invites applications for a tenure-track faculty
position at the Assistant or Associate Professor level in the area of computational
materials science.  The successful
applicant will be expected to establish an externally funded research program, develop
collaboration within the Department and College, mentor graduate students, and teach
at both the undergraduate and graduate levels. 
Minimum qualifications include an earned doctorate in mechanical

Calling a File inside of the UEL

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?


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

Thanks in advance,


niuwenkai's picture

material points will change while using subroutine in ABAQUS

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

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:

CISM Advanced Course "Ferroic Functional Materials"

We would like to draw your attention to the advanced course at the International Centre for Mechanical Sciences (CISM) to be held in Udine (Italy) from September 8 through 12, 2014 on the topic of: "FERROIC FUNCTIONAL MATERIALS: EXPERIMENT, MODELING AND SIMULATION". The advanced course will comprise of distinguished lectures from the following experts:

- Kaushik Bhattacharya (Caltech, USA)
- Manfred Fiebig (ETH Zurich, Switzerland)
- John Huber (University of Oxford, UK)
- Doru C. Lupascu (University of Duisburg-Essen, Germany)
- Christopher S. Lynch (UCLA, USA)
- Jörg Schröder (University of Duisburg-Essen, Germany)

For more information please visit:

Job gets aborted due to cpset and nset history output definition?

Dear all,

 I am modelling an interactin between deformable and analytical bodies by utilisign Abaqus/Explicit. In order to define the contact between those objects I have selected surface-to-surface technique type and applied kinematic contact method as formulation method. As well I have selected individual output for a specially created set of my primary interest. Nevertheless, when I submit a job for analysis it gets aborted and the following warning errors appear:

The cpset or the surface parameter must be specified for history contact output

The nset parameter can only be used for bondstat or bondload output

 Can anybody suggest me what should be done to overcome the problem? Thanks in advance!

Pengfei Liu's picture

When is strong and when is weak by using extended finite element method?

As we know, the extended finite element method has been well used to simulate the localization band propagation and crack propagation without remeshing. It can deal with both strong discontinuity problems (displacement discontinuity) and weak discontinuity problems (strain or acceleration discontinuity). Now, I have always a question: when is strong and when is weak discontinuity using XFEM? In other words, what is the difference for the interpolation function and level set function or other function between strong and weak discontinuity? Thanks.

Liu pengfei

Nuwan Dewapriya's picture

Modelling fracture of graphene using Griffith’s criterion and quantized fracture mechanics

In armchair graphene sheets, crack propagates perpendicular to the applied strain, whereas crack propagation in zigzag sheets occurs at an angle to the straining direction. This occurs due to different bond structure along armchair and zigzag directions as shown in Fig. 1. Videos 1 and 2 show the fracture of armchair and zigzag sheets, respectively.


Fig. 1: Armchair and zigzag directions of graohene

Ajit R. Jadhav's picture

MWR for the first- and third-order differential equations

Hi all,

In engineering sciences, we usually end up using either the second- or the fourth-order differential equations, and the MWR (the method of weighted residuals) works pretty well for them.

The question is: how about the first- and the third-order differential equations? Why don't we see any applications of MWR for these odd-ordered differential equations? What gives?

Error in finding connector element result please help

I am getting displacement and force both zero. I am attaching .cae file please have a look where i am making mistakes.

crack simulation in a body

i want to do crack anlysis in composite.i am thinking to use freefem++ or matlab for this purpose.please suggest me which one should i use.if any other,please recommend.

genepak's picture

Senior Researcher/Post-Doc Position in Solid Mechanics at AICT, Seoul National University, Korea

Applications are invited for the position of senior researcher/postdoctoral fellow in the following areas with special focus on multiscale defect mechanics of electronic materials and devices:

mohsenzaeem's picture

Call for Abstracts in Mini-Symposium on "Computational Modeling of Microstructural Evolution II”- ASME Congress 2014

Dear Colleague,

We would like to invite you to participate in a mini-symposium entitled “Computational Modeling of Microstructural Evolution II” in ASME 2014 International Mechanical Engineering Congress & Exposition (
ASME Congress and Exposition will be held in Montreal's downtown convention center, 14-20 November 2014.

siegmund's picture

Post Doc Position

Syndicate content