CFD Online Logo CFD Online URL
Home > Forums > System Analysis

Deformable mesh in 1D elastic body

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   November 23, 2013, 08:40
Smile Deformable mesh in 1D elastic body
New Member
Alexander Kazakov
Join Date: Nov 2013
Posts: 1
Rep Power: 0
alexanderkazak is on a distinguished road
Good afternoon!
I'm solving the pde's of 1D deformable solid body (just an elastic case without viscosity)

du/dt + A*du/dx = 0,

where u is a 2D vector u = (velocity, tension) = (v, t)
As I understand, there are two methods: Euler's and Lagrange's
In Euler's coordinates, "x", which are not moving in the space, the matrix A is

( v -1/d
-E v ), d is density, E is Young's modulus

In Lagrange's coordinates, "y", which are moving with points of body

A = ( 0 -1/d
-E 0 )
Here equation is the same, du/dt + A * du/dx = 0, (not du/dt + A * du/dy = 0 !), but
but the coordinate of the particle is changed: dx = v*dt

1) Is it correct?

I use a grid-characteristic method to solve it.
In Euler's method, I diagonalize A and transfer the values of two invariants along the corresponding characteristics.
In Lagrange's method, I do the same with Lagrange's matrix (and formulas for recalculation invariants from (v, t) is the same!) and after that I do dx = v*dt for each node of my mesh.
This methods are equal, to my mind.

2) Is it correct?

3) What's the exact solution? It's not like a "f(x - at) + g(x + at)", where a = sqrt(E/d), because of "dx = v*dt" in Lagrange's method and diagonal "v" in Euler's method!

I'm looking forward for any answers!
alexanderkazak is offline   Reply With Quote


deformable mesh elastic

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Mesh motion with Translation & Rotation Doginal CFX 2 January 12, 2014 07:21
Layers:problem with curvature giulio.topazio OpenFOAM Native Meshers: snappyHexMesh and Others 10 August 22, 2012 09:03
2D Mesh Generation Tutorial for GMSH aeroslacker Open Source Meshers: Gmsh, Netgen, CGNS, ... 12 January 19, 2012 04:52
[ICEM] Problem making structural mesh on a surface froztbear ANSYS Meshing & Geometry 1 November 10, 2011 09:52
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 18:10

All times are GMT -4. The time now is 07:34.