User login

You are here

mohamedlamine's blog

mohamedlamine'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.

mohamedlamine'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.

mohamedlamine's picture

Transformation Matrix

Hello,

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

mohamedlamine's picture

Publication

Hi all,

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

"CONTRIBUTION TO BRIDGE DAMAGE ANALYSIS" Mohammed Lamine MOUSSAOUI, Abderrahmane KIBBOUA, Mohamed CHABAAT

APPLIED MECHANICS AND MATERIALS Journal,

Volume 704 (2015) pp 435-441

DOI: 10.4028/www.scientific.net/AMM.704.435

 

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

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

and is available at:

mohamedlamine'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.

mohamedlamine's picture

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.

mohamedlamine's picture

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 

 

mohamedlamine'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. 

mohamedlamine's picture

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 

mohamedlamine's picture

academia

Hi all,

 

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

 

                                 www.academia.edu 

Best regards.

mohamedlamine'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 

mohamedlamine's picture

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

mohamedlamine'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 

mohamedlamine's picture

Structural Damping

 

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

mohamedlamine's picture

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

mohamedlamine's picture

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

mohamedlamine's picture

Inverse

It is not recommended to Invert a Matrix in Numerical Analysis. Solving with Factorization methods (L.U, L.D.Lt, ..) is better than Solution with Direct Inverse. This is due to Numerical Errors (roundoff) and/or ill-Conditioned problems. In case of you need to use an Inverse, you are advised to take care about using existing algorithms which can't be exact in all applications.

mohamedlamine's picture

Conditionning

Convergence of Numerical solution may depend on the conditionning of the problem's matrices.

An ill-conditionned matrix produce singulariy cases. This is due to bad shaped elements in

structural mechanics: nearly degenerated elements with zero area (in 2D) or zero volume (in 3D)

The conditionning can be analyzed by numerical algorithms or existing commands in order to

detect these problems.

 

mohamedlamine's picture

Well-Posed Problem in Partial Differentrial Equations

To avoid the Ambiguity appearing in the Formulation Process using the Numerical Methods of Partial Differentrial Equations it is Important to Satisfy the following Conditions given by Fletcher C.A.J. 1989, p18 :

 "The governing Equations and Auxiliary (Initial and Boundary) Conditions are Well-Posed Mathematically if the Following three Condtions are met:

mohamedlamine's picture

Mechanical Vibrations of the Gear Teeth

Recent research work available on request.

 

Subscribe to RSS - mohamedlamine's blog

More comments

Syndicate

Subscribe to Syndicate