|
[Sponsors] |
November 23, 2013, 08:40 |
Deformable mesh in 1D elastic body
|
#1 |
New Member
Alexander Kazakov
Join Date: Nov 2013
Posts: 1
Rep Power: 0 |
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! |
|
Tags |
deformable mesh elastic |
|
|
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 |