CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > System Analysis

Deformable mesh in 1D elastic body

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 23, 2013, 08:40
Smile Deformable mesh in 1D elastic body
  #1
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

Reply

Tags
deformable mesh elastic


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 Off
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
[snappyHexMesh] Layers:problem with curvature giulio.topazio OpenFOAM Meshing & Mesh Conversion 10 August 22, 2012 10:03
[Gmsh] 2D Mesh Generation Tutorial for GMSH aeroslacker OpenFOAM Meshing & Mesh Conversion 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 19:10


All times are GMT -4. The time now is 09:22.