Notes from April 30 – Thursday Many engineering problems can be solved using a matrix approach. Although the Cartesian Coordinate System is the most common for the simpler problems, there are many more suitable systems for more advanced problems. Most unusually, many physical problems have a “Natural Coordinate System.” A complicated formulation in the Cartesian system can be much more simplified using the Natural Coordinate System. For example, 200 coupled simultaneous equations in the Cartesian System can be simplified into 200 independent equations in the Natural system. In terms of matrix algebra, a fully coupled full matrix can be simplified into a diagonal matrix. In the loose set of examples which follow, some important applications are touched upon. (1) The concept of principle stresses. (2) Orthogonal Transformation made possible by an eigenvalue problem. (3) Transforming coupled ordinary differential equations of a structural analysis into independent single degree of freedom differential equations. Other applications: (4) Space curves, the tangent-normal- binormal coordinate system. (5) How to deal with slowly convergent series. (6) Finite Difference formulation for Partial Differential Equations.
15
Embed
Notes from April 30 – Thursday Many engineering problems ...
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Notes from April 30 – Thursday
Many engineering problems can be solved using a matrix approach. Although the Cartesian Coordinate System is the most common for the simpler problems, there are many more suitable systems for more advanced problems.
Most unusually, many physical problems have a “Natural Coordinate System.” A complicated formulation in the Cartesian system can be much more simplified using the Natural Coordinate System. For example, 200 coupled simultaneous equations in the Cartesian System can be simplified into 200 independent equations in the Natural system. In terms of matrix algebra, a fully coupled full matrix can be simplified into a diagonal matrix.
In the loose set of examples which follow, some important applications are touched upon. (1) The concept of principle stresses. (2) Orthogonal Transformation made possible by an eigenvalue problem. (3) Transforming coupled ordinary differential equations of a structural analysis into independent single degree of freedom differential equations. Other applications: (4) Space curves, the tangent-normal-binormal coordinate system. (5) How to deal with slowly convergent series. (6) Finite Difference formulation for Partial Differential Equations.
We started with a symmetric matrix A and obtained its eigenvectors as the columns of matrix Q. It can be shown Q*Q’ or (Q’)*Q is an identity matrix. That means Q is an orthogonal transformation. The transformation (Q’)*A*Q would transform A in to a diagonal matrix D which contains its eigenvalues.
Here is a strategy often used to solve complicated problems. Consider this simple matrix example first: Ax=b. “x” is in the xyz system, so the matrix A is full. Transform x to a new system using Qx=y. Then equation system becomes, Ax=AQy=b. Now pre-multiply the entire system by Q’, then (Q’)*A*Qy=(Q’)y, or Dy=(Q’)b. Normally, a Greek letter, like eta, is used for y. In this transformed coordinate system, the matrix equation is diagonal. It would not be worth the effort to use this method for such a simple problem, but for complicated structural response calculation, this is the most common method. Having a 3000 degrees of freedom system, we can obtain a small sample of eigenvalues and eigenvectors, perhaps 100 or less, the results would be good enough. There is no need to do direct inversion of the 3000 by 3000 matrix.
These are the methods used by engineers worldwide. Create “stick” models of complicated structures, such as a nuclear reactor containment building, and use modal (eigenvalue) analysis to calculate the responses.
There are many cases in transportation engineering which require a special coordinate system, the tnb system. It has three unit vectors in the tangent, normal and binormal directions. These unit vectors change as a vehicle traverses the path. Examples are (1) Roller Coasters. The new ones are totally complicated and the vehicle path goes up and down with many curves. (2) Fighter Jet airplanes travel in erratic paths. (3) On and off ramps from highway or freeways. Very exotic designs of highway bridges, clover leave highway interchanges.
The vectors have a very simple form, the velocity vector has only one component, the tangent component. It means the path is design for the vehicle to pass through it in the tangent direction. The acceleration has two components, the tangential change of speed and the centripetal acceleration calculated as v^2/rho, in which rho is the radius of curvature at that location. The binormal unit vector forms the complete set for force related calculations.
We are lucky in ce108 to have introduced the variable spacing interpolation formulas. We would need those delta functions. One parameter describes the path and the distance “s” marks the location of a point along the path. That point is prescribed with x(s), y(s) and z(s). The derivatives are not simple but we can calculate all of them using divided differences. We have all the tools, we can say with a Jamaican accent, “no problem, mon.”
Oh, have you notice about the “superelevation” design put on the race tracks or on the ramps to fast freeways? To counteract the centripetal acceleration using the weight of the vehicle, a slight slope (or large slope for Olympic Bicycle Racing Tracks) were put in to make driving more pleasant and to save some tire wear. The weight of the vehicle in the negative z direction, but the normal force holding the weight of the vehicle is in the binormal direction, we need orthogonal transformation matrices along the path too. With computer programming, all of these things can be accomplished routinely. It is fun and exciting to code a computer. But seriously, DO NOT WAIT FOR Computer Science graduates to do this. They don’t know any Physics or math higher than Calculus III.
Numerical Determination of Unit Vectors
Thex, y, andz coordinates of a space curve can be defined as a function of parameter,ξ, as
x = x(ξ) , y = y(ξ) , z = z(ξ) .
The tangential unit vector,et, can be obtained as
et =d~r
ds=
d~r
dξ
dξ
ds=
dξ
ds
dx
dξ
dy
dξ
dz
dξ
,
in whichds =
√(dx)2 + (dy)2 + (dz)2 ,
and
ds
dξ=
√(dx
dξ
)2
+(
dy
dξ
)2
+(
dz
dξ
)2
=√
∆ ordξ
ds=
1√∆
The above formulation guarantees that|et| = 1. To find the normal unit vector, perform anotherderivative with respect tos,
the sum after 100 terms is 1.56761, and after 1000 terms is 1.57048, compared to the exact value of1.5708. Needless to say, the series is very slowly convergent. Perform now the forward differencesas shown below: