You are not logged in.
Pages: 1
Hi Bob,
Many thanks for your assistance, this has been great revision for me. I am happy with the results using
tanE = cosBsinC / (cosCcosD - sinBsinCsinD).
The problem is solved just fine and is "close enough for government work".
I can now go ahead and finish my project.
Thanks again, John.
I have been through the other examples and all check out when set up on the mill.
A B C D E Measured C
20 -60 -30 10.3141 -17.8894 ~ -30.0
15 -20 -30 14.1327 -30.4869 ~ -30.0
12 -20 -10 11.2955 -9.7052 ~ -10.0
18 -25 -33 16.4085 -33.7231 ~ -33.0
32 -41 -20 31.228 -23.5713 ~ -20.0
Now just need to clear up the sign/sense problems in Rot(-B)x.
John.
The large difference in A was bothering me because this was correct before. I have reworked the math referring to diagram Johns+box+v3. I have also decided that we are rotating the object, not the axes, I suspect that I have my signs/notations mixed somewhere, see below.
Rotate PQ onto OY around X by angle -B
Rot(-B)x Q R T Q' R' T'
[ 1 0 0 ] [ 0 cosA cosC ] [ 0 cosA cosC ]
[ 0 cosB sinB ] x [ cosB 0 sinC ] = [ 1 sinAsinB cosBsinC ]
[ 0 -sinB cosB ] [ sinB sinA 0 ] [ 0 cosBsinA -sinBsinC ]
In order to make Q' lie on the Y axis I seem to have to reverse the X rotation matrix. I don't understand why and the sign of D is also reversed althought the magnitude is correct.
Vice angle D = angle between plane Q'R'T' and OX in plane XOZ
tanD = R'z / R'x = sinAcosB / cosA = tanAcosB
Next make P'Q'R' parallel to XOY by rotating around Y by angle D. I think the mistake before was using angle C whose projection onto plane XOZ has been reduced to D by the first rotation B. This gave rise to the error in the horizontal vice rotation about Z.
Rot(D)y T' T''
[ cosD 0 sinD ] [ cosC ] [ cosCcosD - sinDsinBsinC ]
[ 0 1 0 ] x [ cosBsinC ] = [ cosBsinC ]
[ -sinD 0 cosD ] [ -sinBsinC ] [ -sinDcosC - cosDsinBsinC ]
Rotate vice about Z by angle E = angle between plane ZST'' and plane XOZ
tanE = T''y / T''x = cosBsinC / (cosCcosD - sinBsinCsinD)
I made another model with these angles:
A B C D E measured C
20 -60 -30 10.3141 -17.8894 ~ -30.0
Best I can tell this seems to work apart from the problems above.
When time permits I will work through the other models to check.
John.
Hi Bob,
I believe I did have X and Y axes mixed up so you are right. I also had not realised we are rotating the axes and not the object if I understand the conventions correctly.
I made another model according to your last diagram (Johns+box+v3).
This time the angles are A = -15, B = 20, C = 30. Photos attached. I set A -ve because if I have it correct +ve angles are X->Y, Y->Z and Z->X. I calculated D = -19.3701 and E = 29.3200. Using these numbers I set up a test block and milled a fin (inverse of a groove). Still the results don't seem to quite agree, particularly in XZ plane this time. PR should have been 15 to X but turns out to be ~20.4. ST should have been at 30 to X in the XY plane but turns out to be ~27.2. The differences can be seen in the photos.
I really don't get what is going on.
John.
Bob I have edited my last post, my fuzzy thinking got it wrong. I hope you were not confused by it. May be it has something to do with having emigrated from Australia to Canada 43 years ago and still not having my world turned around!
Your diagram above is correct. C is the angle between a line from S and plane XZ on the top surface of the prism. A side wall of the groove is the plane normal to surface PQR and containing the line from S.
We were close the first go around but C turned out to be in PQR and not in the top surface of the prism.
John.
Bob I checked out your cube and Excel sheet, thanks. No doubt the math works.
I made another test model with the first values for A, B and C (photo included).
You can see angle C (10°) lies on the XY'' plane. I machined a small area of this plane so the angle can be seen. This is all just as it should be.
After making this it became apparent that I have been missing the fact that what is desired is the line of intersection of one vertical wall of the slot (plane XZ''') and the top surface of the cube (plane XY) be at 10° to plane XZ measured in plane XY. This turns out to be close to 15° in the current test, I would like it to be the angle C (10°).
I have spent most of the day trying to figure out how to achieve this without any result yet.
John.
I must add that I did not define the problem correctly in the beginning. It should now be clearer with an image. The last edit probably muddied the waters some more. What is desired is a skewed groove in the rectangular block oriented such that the bottom surface of the groove extends onto two adjacent sides at angles A and B and that a groove side wall intersects the top surface at angle C (measured in the top surface) to the side in which B is, or its complement (to the side containing B).
John.
Bob,
Here are a couple of photos of the milling setup. Shown is a trial cut in a piece of wood.
1. The block of wood is twisted by angle A around an axis parallel to the lenght of the vice.
2. The vice body is tilted by angle D to the XY plane. This is done by placing a spacer of thickness 4" x sine E under one end of the vice. The centres of the cylindrical surfaces are 4" apart. This gets our A and B angled lines parallel to the XY plane (vice base).
3. The vice is rotated around Z by angle E from the transverse axis of the macine. This is parallel to the riser blocks you see under the vice.
This whole assembly can be translated in / out (X) and left / right (Y). The quill carrying the end mill cutter above the block of wood can be moved up and down (Z).
Hope this is clear.
John.
Hi Bob,
Thanks for checking my arithmetic. I was going to check that our two formulae were equivalent but it slipped my mind. Too bad there wasn't some error that explains the descrepancies, I fear this is not a trivial problem.
Careful setup, milling and measurenment still seems to show the theory is not fitting the practice.
I have been doing a little reading and it appears we may not treating the problem properly. Consider what happens when we do the two rotations to make plane ROQ cincident with the XY plane. These two totations can be combined into one whose axis lies in the XY plane so determining D works ok. The next step is to rotate OP around Z so that it lies in the YZ plane. The three discreet rotations can be combined in to one whose axis in general will not lie in the XY plane. This means that the angles used for the first two rotations will change.
Check this link
http://www.euclideanspace.com/maths/algebra/matrix/orthogonal/rotation/index.htm
If this is too complicated I'll just set my object up by eye and call it good enough.
John.
The matrix revision was great, it turned out to be pretty simple really given the 3 rotation matrices. I would not thought of the way of evaluating angles D and E by setting the appropriate components to 0. I worked through it all and I believe there is a mistake in evaluating tan E.
cosEcosAsinC - sinEcosDcosc - sinEsinDsinAsinC = 0 is correct but this
=> tan E = cosAtanC / (cosD + sinAtanCsinD)
I recalculated E with the following results:
A B C D (Calc) E (Meas) E (Calc) E (Recalc)
12 20 10 -19.60 ~ 5.4 10.51 11.04
18 25 33 -23.92 ~24.9 36.56 36.56 (why this is the same I am not sure)
32 14 20 -11.94 ~15.2 17.13 18.21
The correction has not made a huge difference. I set up two of these examples in the mill and made some cuts. I cannot measure any significant differences from angles A, B and C so I beleive all is well. I guess I have to conclude that my ability to measure some of the angles is quite poor. Makes sense when trying to measure the orientation of a line of intersection of two planes at small angles to each other.
I quite enjoyed the exercise and many thanks for your assistance.
John.
Bob many thanks for all the info. and the Wiki link, it is going to take me some time to digest all this. I'll let you know.
John.
Bob I tried your function for E and made some more measurements as follows:
A B C D (Calc) E (Meas) E (Calc)
12 20 10 -19.60 ~ 5.4 10.51
18 25 33 -23.92 ~24.9 36.56
32 14 20 -11.94 ~15.2 17.13
I am not sure what to make of this. I am pretty certain that my statement "Intuitively if A=B and C=45° then Φ=45°." is wrong.
I have been assigning angles as positive as measured from their respective axis. If I understand the usual convention YOQ is the only positive angle measured from 0, maybe this is the problem. I am interpreting your diagram above as having X coming forwards and the corner formed by the three planes as concave.
Could you post the notation you are using, or a refence to it, defining the sense in which angles are positive and the matrices plus the multiplication reminder.
Thanks, John.
Bob I did suspect it would be much easier using matrices and vectors. I did some introductory courses in linear algebra about 25 years ago but have never had to use it, so its gone now. I am now retired from working in the airborne geophysical business searching mainly for metal deposits and mapping geology. One of my hobbies is metal machining by traditional hand operated methods. I am sure this problem comes up often and is no doubt solved in software in CAD/CAM and CNC machining packages.
You may not have a mistake with the sign because I did not pay any particular attention to signs of angles, I am really interested in magnitude because when viewing the actual part one can see from layout lines how it should be oriented.
I have made carefull measurement of C in one case only. The error would be about +/- 0.25 degrees.
A=12°, B=20°, C=10°. θ=atan(tanBcosA)=19.5966°. Rotation about Z, Φ=~5.4°(measured)
Intuitively if A=B and C=45° then Φ=45°. If A=B=0° or 90° then something most likely will be undefined. I suspect there may be trouble if C>45° also.
Take your time, I've been messing with this for about 3 weeks now.
John.
Yes Bob that is the gist of it. I should have made a diagram but am only capable of drawing on paper with a pencil and then scanning it.
Angle A would be ROX, B YOQ and C POY.
The verical jaw face of the vice is the XZ plane. The prism would be twisted about the Y axis such that OR coincides with OX.
Next the vice is tilted around axis X so that OQ lies in plane XY. This angle is given by atan(tanBcosA). This is the common compound angle problem found in miter cuts for crown moldings and the like.
The third rotation about the Z axis is to bring line OP into the YZ plane.
Rotation of the prism about OY and OX will shorten the projection of OP on the XY plane. So it seems to me that the third angle of rotation about OZ must be some function of A, B and C.
With my rather rusty high school trig I don't seem to be able to figure this angle out.
The problem relates to machine tool setup.
Suppose we have cartesian axes X, Y & Z and a rectangular right angled prism with adjacent sides a, b & c, has lines at angles A, B, C on each surface respectively all intersecting at the common corner. The requirement is to orient the prism such that the plane containing lines at angles A and B is parallel to X-Y.
I have been able to derive this. Twist the prism such that line at angle A is parallel to X-Y, then using this line as an axis tilt the prism by angle atan(tanBcosA) X-Y (compound angle problem).
Next the block must be rotated about Z by some angle so that a vertical plane containing the line at angle C is parallel or perpendicular to X - Z or Y - Z.
How to find this angle eludes me.
Help or a link that shows how will be much appreciated.
Thanks, John.
Pages: 1