User login


You are here

mohammedlamine's blog

mohammedlamine's picture

Efficient Constant Strain Tetrahedron Shape Functions

Dear researchers,

Hi all,

I am pleased to post my following publication:

Study on Dynamic Detection of Reinforced Concrete Bridge Damage by Finite Element Model Updating


mohammedlamine's picture

My Scientific New Paper

Hi All,

My New Scientific Paper is :

Structural Probabilistic Health Monitoring of a Potentially Damaged Bridge

Mohammed Lamine MOUSSAOUI , Mohamed CHABAAT , Abderrahmane KIBBOUA

International Journal of Structural and Civil Engineering Research

Volume 10, No. 3, pp. 119-122, 2021


mohammedlamine's picture

Scientific Paper

Hi All,

Please Find in this Post my Scientific Paper :

Numerical Analysis of Damage Zones in a Bridge

Mohammed Lamine MOUSSAOUI , Mohamed CHABAAT

International Journal of Structural Integrity, Emerald Publishing,

Volume 11, No. 1, pp. 1-12, 2020



Best Regards

mohammedlamine's picture

An efficient convergence test for the fixed point method

The fixed point method consists to find the solution of F(X)=X.

One can not get fixed with the convergence condition |F'(X)|<1 because if the function has an optimum then |F'(X)|=0 even if the solution is not yet reached.


We introduce an efficient convergence test with the condition:

|Xn+1 - Xn| ≤ epsilon1 And |F(Xn+1)-Xn+1| ≤ epsilon2

mohammedlamine's picture

Concrete stress-strain curve

The concrete stress-strain curve corresponds to the compression case with negative values plotted on positive axis. If one needs to plot the tensile curve for this concrete material this will be almost similar with {0.1*compression} stresses with positive values. Furthermore if one needs to plot on the same axis the tensile steel and concrete stresses for a composite mixture he needs to consider this last case. Fortunately the tensile and compression steel stresses are quite similar with opposite values.

mohammedlamine's picture

Failure Envelope

A finite element software computes the stresses and/or strains sets values at each dynamic iteration step. These values can be plotted on a stress-space giving several failure envelopes. For an isotropic or an homogeneous mixture the Von Mises criterion represents an ideal envelope which may not coincide with the available envelopes but the stresses values can be checked with the material strength.

mohammedlamine's picture

Transformation Matrix


The attached file describes the development of the transformation matrix of a beam element:

mohammedlamine's picture


Hi all,

I am pleased to post the references of my recent paper entitled :



Volume 704 (2015) pp 435-441

DOI: 10.4028/


N.B: The equation (7) has been reviewed with:

[P2](nxg)=[[I](nxn)| [0](nxf)]

and is available at:

mohammedlamine's picture

Reinforced Concrete Strength

The composite strength of Reinforced Concrete which is a mixture of two main materials can be calculated by the formula given by Tsai : X=(vf+vm*Em/Ef)*Xf  where vf is the fiber (steel) fraction and vm is the concrete (cement+sands+gravels) fraction. Xf is the yield strength of Steel and Xm is the compressive strength of the Concrete matrix. For the application case: vf=0.01911, vm=0.98089, Xf=270 MPa, Xm=40 MPa, Ef=207.000 GPa, Em=30.798 GPa we obtain X=44.563 MPa.

mohammedlamine's picture

Neutralization of errors


The problem in numerics is the neutralization of the errors for example plus (+)error and minus (-)error. For that reason it is difficult to locate a convergent process. Some numerical algorithms avoid these omissions like leat squares methods which handle the square of the errors in the minimization function. In the other applications like solving processes truncating errors and roundoff errors can be analyzed successfully.

mohammedlamine's picture

Shape Functions of a Constant Strain Tetrahedron


I have developped a New Method to Compute Numerically the Shape Functions of a Constant Strain Tetrahedron from the Geometrical Information of its Vertices using the Planes' Equations Defining the Barycentric Coordinates versus the Cartesian Coordinates. The Algorithm for Solving the Linear Systems of Equations has been Implemented Efficiently with the Partial Pivoting Gauss LU Method for any arbitrary Shaped and Non-Degenerated Element.

References: FEOFS2013, CanCNSM2013 and FDM2013.


Best Regards.


Mohammed Lamine Moussaoui 


mohammedlamine's picture

Technical Publication

Hi all,

Here are the references of the following technical publication: 

Title: " Dynamic Detection of Reinforced Concrete Bridge Damages by Finite Element Model Updating ".

Author(s) : Mohammed Lamine MOUSSAOUI, Mohamed CHABAAT, Abderrahmane KIBBOUA.

Area: Non-Smooth Vibrations in Structural Dynamics II.

Conference: 4th Canadian Conference on Nonlinear Solid Mechanics (CanCNSM 2013)

Location: McGill University, Montreal, Canada

Date: July 23-26, 2013

Status: Accepted / Done


Best regards. 

mohammedlamine's picture


Hi all,

There are several formulae allowing the computation of tetrahedral volume elements. The weakness of some of them is that they give negative volumes. This is due to the vector dot product and its direction. You must take care when using them.

You have to check your formula with some simple examples like the vertices of the following one: (1,0,0) , (0,1,0) , (0,0,1) , (0,0,0). 


Mohammed lamine Moussaoui 

mohammedlamine's picture


Hi all,


If you have published works before the era of internet or want to share research, find below a useful link:



Best regards.

mohammedlamine's picture

Convergence of a Solution


The Solution of a Differential Equation or a Set of Differential Equations Converges vers the Exact if it is Consistant and it is Stable : Lax's Theorem. Since the Exact is not always known it is convenient to apply this Theorem. Numerical instabilities are a result of Roundoff errors and Truncation errors. The Domain of Stability can be obtained from a Von Newman analysis in the Complex domain. This implies a condition (relation) between the variation steps. 


Mohammed Lamine 

mohammedlamine's picture


" Damage Detection in Bridges using a Mathematical Model by an Updating Method "

Mohammed Lamine MOUSSAOUI , Mohamed CHABAAT , Abderrahmane KIBBOUA

The 9th International Conference on Fracture and Strength of Solids, FEOFS 2013, June 2013, Jeju, South Korea

Status: Accepted

mohammedlamine's picture

free-free vibrations v.s F.E. softwares

Hi All,

Do you know why some finite element softwares can not solve free-free vibrations of many applications unless you fix at least one degree of freedom ?

In this case I do not see how they are able to analyze free-free structures as planes.


Mohammed Lamine 

mohammedlamine's picture

Structural Damping


Which best Describes Structural Damping : Dissipative Energy, Linearization from Stiffness and Mass Matrices, or Strain Energy  

mohammedlamine's picture

Are Finite Element Softwares Reliable ?


Some Finite Element Softwares use the Formulation of an Element and Jump to another one

by Merging Nodes : Quadrilateral to Triangular and Hexahedral to Tetrahedral without modifying

the Shape Functions.

This Method is Mathematically and Numericaly Forbidden because they use the Same Shape Functions

when Merging the Nodes for the Defined First Element (see Element Definition and after Mesh Process).

mohammedlamine's picture


 The Simplex K with 4 Vertices in Dimension 3 is a Tetrahedron and :

 K = determinant((x2-x1) , (x3-x2) , (x4-x3) ; (y2-y1) , (y3-y2) , (y4-y3) ; (z2-z1) , (z3-z2) , (z4-z3))

 The Volume of this element is :

 V(K) = abs(K) / 3! 

 If V(K) ~= 0  then K is Nearly Degenerated




Subscribe to RSS - mohammedlamine's blog

Recent comments

More comments


Subscribe to Syndicate