User login


You are here

# Finite Element modeling

sp435's picture

High strain rate liver tissue compression

Check out our latest publication on porcine liver biomechanics under high strain rate conditions -

"Mechanical Response of Porcine Liver Tissue under High Strain Rate Compression", Bioengineering, 2019, 6(2), 49;

Emilio Martínez Pañeda's picture

Phase field modelling of crack propagation in functionally graded materials

Dear imechanicians

I hope that this paper is of interest to you. The phase field fracture code developed in FEniCS is available at

Phase field modelling of crack propagation in functionally graded materials

Hirshikesh, S. Natarajana, R.K. Annabattula, E. Martínez-Pañeda

Composites Part B: Engineering 169, pp 239-248 (2019)

When to use equation of state models in modeling high velocity impacts?

Equation of state models (Mie-Gruneisen equation of state etc) have been widely used to study hydrodynamic response of material during explosive deformation, ballistic impacts for a long time. As far as I understand, a hydrodynamic response is required to be incorporated in a model if the pressure of impact leads to an increase volumetric strength. That is the stress wave speed in the material exceeds the speed of sound in the material (v = sqrt(elastic modulus/density)).  Ballistic impacts occur at strain rates close to 10^9 to 10^12/s.

yann.charles's picture

PhD position / La Rochelle - France

A PhD position is available at LaSie (La Rochelle University, France), until oct. 2019.

The work will be dedicated to the caracterisation (numercial and experimental) of the effect of grains boundary on the Hydrogen diffusion in Nickel polycrystals.

Please send your resume and application letter to:

J. Bouhattate (

Postdoctoral Scholar Position in Computational and Experimental Biomechanics at East Carolina University

The Postdoctoral Scholar is responsible for independent and collaborative development, implementation, and analysis of scientific research in the discipline of computational and experimental biomechanics. Responsibilities include management of all aspects of research related to biomechanical characterization and computer modeling of biological tissue and implants in Vahdati Lab. The research projects include investigation of mechanopathologies of soft and calcified tissue and design of personalized medical devices.

MMSLab-CNR's picture

Modal analysis performed by NOSA-ITACA code

The following picture shows a screenshot of the NOSA-ITACA desktop enviroment during a modal analysis. In particular the first and the third mode shape of the San Frediano bell Tower in Lucca (Italy) are shown. The structure is discretized into 45641 eight-node brick elements (element no. 8, NOSA-ITACA element library) with 136923 degrees of freedom. For more details sees: R.M. Azzara, G. De Roeck, E. Reynders, M. Girardi, C. Padovani and D.

A PhD position at the University of Western Ontario, Canada

We are looking for a talented, self-motivated, and enthusiastic student with a BSc or MSc degree in mechanical or materials engineering. Students with strong background in crystal plasticity and computational mechanics are recommended to apply.


Prospective candidates will be assessed based on how well they meet the following criteria:

Solving incompressible finite elasticity without tears

Solving incompressible elasticity has been quite challenging numerically. The conventional approach for handling incompressibility is the so-called penalty method. A volumetric energy term enters into the strain energy and penalizes the volumetric deformation. One straightforward issue is that the penalty parameter goes directly into the tangent matrix. The bigger the penalty parameter, the worse the condition number of the matrix. This is really a manifestation of the ill-posedness of theories based on the Helmholtz free energy, in my opinion [3].

Emilio Martínez Pañeda's picture

Steady-state fracture toughness of elastic-plastic solids: Isotropic versus kinematic hardening

I hope some of you may find this work interesting. We show that kinematic hardening effects play a significant role in monotonic/static fracture.

Steady-state fracture toughness of elastic-plastic solids: Isotropic versus kinematic hardening

K.J.Juul, E.Martínez-Pañeda, K.L.Nielsen, C.F.Niordson

Engineering Fracture Mechanics, 207, pp. 254-268 (2019)

Volume 207, 15 February 2019, Pages 254-268

Kmomeni's picture

PhD position in Multiscale Modeling Materials Under extreme conditions

An immediate fully-funded Ph.D. position is open in Advanced Hierarchical Materials by Design Lab at Louisiana Tech University on multiscale modeling of materials under extreme conditions. The project involves collaboration with the computational research groups at national labs, as well as experimental researchers from MIT and Georgia Tech to test and validate the models and codes.

Emilio Martínez Pañeda's picture

A phase field formulation for hydrogen assisted cracking [code included]

I hope some of you may find this work interesting, the finite element code developed (Abaqus UEL subroutine) can be downloaded from


A phase field formulation for hydrogen assisted cracking


Emilio Martínez-Pañeda, Alireza Golahmar, Christian F. Niordson


Computer Methods in Applied Mechanics and Engineering, 342, pp. 742-761 (2018)


Emilio Martínez Pañeda's picture

Size effects in elastic-plastic functionally graded materials

I hope some of you may find this work interesting:

Size effects in elastic-plastic functionally graded materials

Tittu V. Mathew, Sundararajan Natarajan, Emilio Martínez-Pañeda

Composite Structures, 204, pp. 43-51 (2018)

A post-print is available at

Early Stage Researcher position, Biomedical simulations, Cardiff University

Cardiff University invites applications for an Early Stage Researcher position (Doctoral Candidate) as part of the Rapid Biomechanics and Simulation for Personalised Clinical Design (RAINBOW) MCSA European Training Network. RAINBOW is funded under the European Union’s Horizon 2020 research and innovation programme. The post holder will undertake research on “Meta Modelling for Soft Tissue Contact and Cutting Simulation” leading to a PhD degree award. The post holder will develop numerical methods to simulate the deformations of soft-tissues in the context of computer-aided surgery.

Emilio Martínez Pañeda's picture

Crack Growth Resistance in Metallic Alloys: The Role of Isotropic Versus Kinematic Hardening

We have always modelled crack propagation under monotonic/static loading in metals assuming isotropic hardening. However, we show that anisotropic/kinematic hardening effects play a significant role due to non-proportional straining with crack advance; the isotropic hardening idealization leads to steady state fracture toughness predictions that could be 50% lower. I hope that some of you find this work interesting.

Emilio Martínez-Pañeda, Norman A. Fleck. Crack Growth Resistance in Metallic Alloys: The Role of Isotropic Versus Kinematic Hardening. 

GTN model

I have 2 questions regarding GTN model.

1) What does q1, q2, q3 signify in GTN damage model? What are their physical meanings. I know only their numerical values.
2) How to deal with the problem of GTN in Compression? It is working fine with tension but while in compression it is unable to converge.

yann.charles's picture

2 years PostDoctoral Position - Universities Lille & Paris 13 , Paris, France

A 2 years post Doctoral position is avalaible in France, on material science numerical modeling fields. This is a join project between Univ. Lille (UMET lab), Paris 13 (LSPM lab), CEA (IRFM lab) and Marseille univ. (PIIM lab).

The present postdoc subject is focused on the modeling of the interactions between impurities (H, He...) and tungsten. 

Emilio Martínez Pañeda's picture

Non-local plasticity effects on notch fracture mechanics (code included)

I hope some of you find this work interesting, the code with the control algorithm (Segurado and Llorca, 2004) to avoid convergence problems in cohesive zone model analyses of crack propagation can be downloaded from (documentation and Abaqus input files are provided).

Non-local plasticity effects on notch fracture mechanics

Emilio Martínez-Pañeda, Susana del Busto, Covadonga Betegón 

Theoretical and Applied Fracture Mechanics (2017)

claudiosaul.lopes's picture

PhD position: multiscale computational mechanics of advanced composites

The Composites Materials Group ( at IMDEA Materials, led by Dr. Carlos González and Dr. Cláudio S. Lopes, invites applications for a PhD student position in the field of multiscale computational mechanics of advanced composites. Candidates with interest and knowledge on computational damage mechanics, model development, numerical simulation and optimization methods for advanced composite materials and structures are strongly encouraged to apply.

PhD positions in the area of materials modelling and characterization: Western University, Canada

PhD positions are available:

I am looking for talented, self-motivated, and enthusiastic students with MSc degree in mechanical or materials engineering. Students with a background in any of the following areas are recommended to apply:

Flavio Stochino's picture

A simplified model for railway catenary wire dynamics

In this paper, a simplified analytic model for the dynamic behaviour of railway catenary wire is presented. The model is discussed and validated with the help of numerical results obtained by a finite element code, for the case of a typical Italian railway installation. The simplified model is convenient from the computational point of view and is useful for sensitivity analysis. Some parametric studies have been developed by considering as free parameters the velocity and the distance of the train pantographs and looking at their effect on catenary dynamics.

Noel ODowd's picture

Fellowships (3 year) available in computational mechanics of composite materials

Expressions of interest are invited for applicants interested in applying for a 3 year post-doctoral fellowship in the area of computational mechanics of composite materials.

rajan_prithivi's picture

Local mesh refinement using Gmsh

Hello All,

I am using Gmsh for creating a finite element mesh. After I create the volume mesh, Is there a way to do a local mesh refinement around a particular point (say x,y,z) and within a sphere of given radius.

Kindly help !



Michael A. McCarthy's picture

Marie Curie Early Stage Researcher (PhD student), Improving the crashworthiness of composite transportation structures

Applications are invited for a PhD position at the University of Limerick, Ireland, funded through the EU Marie Curie Early Stage Researcher Programme.

SALARY: €37,903 p.a. (plus family allowance of €4,589 p.a. where appropriate in accordance with Marie Skłodowska-Curie programme rules).

Topic: Development of discrete fastener-less joining technologies for impact protection devices


Subscribe to RSS - # Finite Element modeling

Recent comments

More comments


Subscribe to Syndicate