Traslate in matlab equations of rotational displacement
5 次查看(过去 30 天)
显示 更早的评论
Hi everyone. I would like to ask for your help to 'translate' these three equations for seismic motion in Matlab.
These are the equations I want to implement and solve in Matlab to get the values of omega.
I translated it like this and I don't know if this is ok.
equationX = omega_x == 1/2 * (diff(uz, y) - diff(uy, z));
equationY = omega_y == 1/2 * (diff(ux, z) - diff(uz, x));
equationZ = omega_z == 1/2 * (diff(uy, x) - diff(ux, y));
X_solution = solve(equationX, ux);
Y_solution = solve(equationY, uy);
Z_solution = solve(equationZ, uz);
CLARIFY WHAT THE VARIABLES ARE
When we record an earthquake with a seismic station,
we get the values of the displacements (actually oscillations) in the 3 spatial directions: ux uy and uz.
Now I want to calculate the rotational movement around the axes: ωx, ωy and ωz.
0 个评论
回答(1 个)
William Rose
2024-3-27
The equations you have given are for a rigid body. If you are comfortable regarding the earth as a rigid body within a zone of interest*, then you need to measure the x,y,z displacements at three non-collinear points in order to compute (or estimate) the rotation. You can understand this by considering a cube. If I know the x,y,z displacement of one corner of the cube during a movement, I will not be able to determine if the cube rotated, or translated, or did a little of both. But if I know the x,y,z displacements of three corners of the cube, then I will be able to determine how much translation and rotation occurred. If I know the x,y,z displacements of four or more corners, then I can use the "extra" information to estimate whether the cube is in fact moving as a rigid body, or not.
A mathematical way of understanding the requirement for displacements at three points is that if I only have at one point, I will not be able to compute and the other partials. I need at three non-collinear points in order to estimate the different partials.
* In seismology, one often does NOT regard the earth as a rigid body.
2 个评论
William Rose
2024-3-27
You said in your initial post
"we get the values of the displacements (actually oscillations) in the 3 spatial directions: ux uy and uz"
Please note that ux, uy, uz, in the equations for , are the linear velocities, not the linear displacements.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!