# mohammedlamine's blog

## 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

## 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

## 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

## 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

## 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.

## Transformation Matrix

Hello,

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

## 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.

## Neutralization of errors

Hello,

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.

## Shape Functions of a Constant Strain Tetrahedron

Hello,

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

## 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)

Date: July 23-26, 2013

Status: Accepted / Done

Best regards.

## Weakness

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

Hi all,

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

Best regards.

## 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

## Publication

" 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

## 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

## Structural Damping

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

## Are Finite Element Softwares Reliable ?

Hello,

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).

## Tetrahedron

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

Moussaoui

## Pages 