0 20 40 60 80
Zeit in[s] 03
69 1215
Ain[cm2]
0 20 40 60 80
Zeit in[s] 03
69 1215
scin[cm]
0 20 40 60 80
Zeit in[s] 02
46 108
zcin[cm]
0 20 40 60 80
Zeit in[s] 0
200 400 600 800
Issin[cm4]
0 20 40 60 80
Zeit in[s] 0
500 1000 1500 2000
Izzin[cm4]
0 20 40 60 80
Zeit in[s] 0
200 400 600 800
Izsin[cm4]
Bild B.1:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=1,B=0,Thresholding:pth=11.14kPa}.
78B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.2:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=3,B=0,Thresholding:pth=8.47kPa}.
79
Bild B.3:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=1,B=1,Thresholding:pth=15.35kPa}.
80B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.4:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=3,B=1,Thresholding:pth=15.44kPa}.
81
Bild B.5:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=1,B=2,Thresholding:pth=29.56kPa}.
82B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.6:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=3,B=2,Thresholding:pth=21.01kPa}.
83
Bild B.7:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=1,B=0, Kantendetektion:pe=16.08 Ta x elkPa }.
84B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.8:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=3,B=0, Kantendetektion:pe=12.87 Ta x elkPa }.
85
Bild B.9:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskreten Parameterkombination {M=1,B=1, Kantendetektion:pe=13.74 Ta x elkPa }.
86B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.10:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=3,B=1, Kantendetektion:pe=9.33 Ta x elkPa }.
87
Bild B.11:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=1,B=2, Kantendetektion:pe=8.88 Ta x elkPa }.
88B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.12:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=3,B=2, Kantendetektion:pe=5.45 Ta x elkPa }.
89
Bild B.13:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=1,B=0, Seriell:pe=2.00 Ta x elkPa undpth=12.04kPa}.
90B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.14:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=3,B=0, Seriell:pe=2.00Ta x elkPa undpth=8.93kPa}.
91
Bild B.15:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=1,B=1, Seriell:pe=2 Ta x elkPa undpth=15.43kPa}.
92B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.16:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=3,B=1, Seriell:pe=2.00Ta x elkPa undpth=17.80kPa}.
93
Bild B.17:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=1,B=2, Seriell:pe=2.00 Ta x elkPa undpth=29.04kPa}.
94B Zeitverläufe der optimierten Kontaktgeometrieschätzungen aller diskreten Parameterkombinationen
Bild B.18:Zeitverläufe der geschätzen Kontaktgeometrieparameter sowie deren Referenzverläufe mit der diskre-ten Parameterkombination {M=3,B=2, Seriell:pe=2.00Ta x elkPa undpth=22.56kPa}.
Literatur
[1] Atmel Corporation.SAM4E Explained Pro Design Documentation. Produktdokumentation.
2015.
[2] Atmel Corporation.SAM4E Xplained Pro User Guide. Produktdokumentation. 2015.
[3] Atmel Corporation.SeriesSAM4E - SMART ARM-based Flash MCU. Produktdokumentation.
2015.
[4] Bhattacharjee, T., Rehg, J. M. und Kemp, C. C. „Haptic classification and recognition of objects using a tactile sensing forearm“. In: 2012 IEEE/RSJ International Conference on Intelligent Robots and Systems. IEEE, Okt. 2012, S. 4090–4097.DOI:10.1109/IROS.2012.
6386142.
[5] Bourke, P.Calculating the area and centroid of a polygon. 1988.
http://paulbourke.net/geometry/polygonmesh/(besucht am 21. 06. 2017).
[6] Burger, W. und Burge, M. J.Digitale Bildverarbeitung. 3. Aufl. X.media.press. Berlin, Hei-delberg: Springer Berlin Heidelberg, 2015.DOI:10.1007/978-3-642-04604-9.
[7] Cannata, G., Maggiali, M., Metta, G. und Sandini, G. „An embedded artificial skin for huma-noid robots“. In:2008 IEEE International Conference on Multisensor Fusion and Integration for Intelligent Systems. IEEE, Aug. 2008, S. 434–438. DOI:10.1109/MFI.2008.4648033. [8] Canny, J. „A computational approach to edge detection“. In:IEEE Transactions on Pattern
Anal. Mach. Intell.8.6 (1986), S. 679–698.DOI:http://dx.doi.org/10.1109/TPAMI.
1986.4767851.
[9] Chand, D. R. und Kapur, S. S. „An Algorithm for Convex Polytopes“. In:Journal of the ACM 17.1 (1970), S. 78–86.DOI:10.1145/321556.321564.
[10] Cormen, T. H., Leiserson, C. E., Rivest, R. L. und Stein, C. „Finding the convex hull“. In:
Introduction to Algorithms, Third Edition. 3rd. The MIT Press, 2009, S. 947–957.
[11] Cutkosky, M. R., Howe, R. D. und Provancher, W. R. „Force and Tactile Sensors“. In:Springer Handbook of Robotics. Berlin, Heidelberg: Springer Berlin Heidelberg, 2008, S. 455–476.
DOI:10.1007/978-3-540-30301-5_20.
[12] Dahiya, R. S., Lorenzelli, L., Metta, G. und Valle, M. „POSFET devices based tactile sensing arrays“. In:Proceedings of 2010 IEEE International Symposium on Circuits and Systems. IEEE, Mai 2010, S. 893–896.DOI:10.1109/ISCAS.2010.5537414.
[13] Dahiya, R. S., Metta, G., Valle, M. und Sandini, G. „Tactile sensing-from humans to huma-noids“. In:IEEE Transactions on Robotics26.1 (2010), S. 1–20.DOI:10.1109/TRO.2009.
2033627.
[14] Dargahi, J., Parameswaran, M. und Payandeh, S. „Micromachined piezoelectric tactile sen-sor for an endoscopic grasper - theory, fabrication and experiments“. In:Journal of Micro-electromechanical Systems9.3 (2000), S. 329–335.DOI:10.1109/84.870059.
[15] Eddy, W. F. A New Convex Hull Algorthm for Planar Sets. 1977. DOI:10 . 1145 / 355759 . 355766.
96 Literatur
[16] Elmo Motion Control ltd.SimplIQ - Software Manual. Produktdokumentation. 2011.
[17] Fraden, J.Handbook of Modern Sensors. New York: Springer-Verlag, 2010.
[18] Goger, D., Gorges, N. und Worn, H. „Tactile sensing for an anthropomorphic robotic hand:
Hardware and signal processing“. In:2009 IEEE International Conference on Robotics and Automation. IEEE, Mai 2009, S. 895–901.DOI:10.1109/ROBOT.2009.5152650.
[19] Gonzalez, Rafael C. and Woods, R. E.Digital Image Processing (3rd Edition). Upper Saddle River, NJ, USA: Prentice-Hall, Inc., 2004.
[20] Graham, R. „An efficient algorith for determining the convex hull of a finite planar set“.
In:Information Processing Letters1.4 (1972), S. 132–133.DOI:10.1016/0020-0190(72) 90045-2.
[21] Han, D. „Comparison of Commonly Used Image Interpolation Methods“. In:Proceedings of the 2nd International Conference on Computer Science and Electronics Engineering (ICCSEE 2013). Iccsee. Paris, France: Atlantis Press, 2013, S. 1556–1559.DOI:10.2991/iccsee.
2013.391.
[22] Harder, M.Experimentelle Untersuchung der charakteristischen Eigenschaften einer taktilen Roboter Haut. Bachelorarbeit. Lehrstuhl für Angewandte Mechanik, Fakultät für Maschi-nenwesen, Technische Universität München, 2016.
[23] Hildebrandt, A.-C., Wittmann, R., Wahrmann, D., Ewald, A. und Buschmann, T. „Real-time 3D collision avoidance for biped robots“. In: IEEE/RSJ International Conference on Intelligent Robots and Systems(2014), S. 4184–4190.
[24] Hildebrandt, A. C., Wahrmann, D., Wittmann, R., Rixen, D. und Buschmann, T. „Real-time pattern generation among obstacles for biped robots“. In:IEEE International Conference on Intelligent Robots and Systems(2015), S. 2780–2786.DOI:10.1109/IROS.2015.7353759. [25] Hildebrandt, A.-c., Demmeler, M., Wittmann, R., Wahrmann, D., Sygulla, F., Rixen, D. und Buschmann, T. „Real-Time Predictive Kinematic Evaluation and Optimization for Biped Ro-bots (submitted)“. In:IEEE/RSJ International Conference on Intelligent Robots and Systems (2016).DOI:10.1109/IROS.2016.7759852.
[26] Hooke, R. und Jeeves, T. A. „“ Direct Search” Solution of Numerical and Statistical Pro-blems“. In: Journal of the ACM 8.2 (Apr. 1961), S. 212–229. DOI: 10 . 1145 / 321062 . 321069.
[27] Huang, J. C. „Carbon black filled conducting polymers and polymer blends“. In:Advances in Polymer Technology21.4 (2002), S. 299–313.DOI:10.1002/adv.10025.
[28] Hyung-Kew Lee, Jaehoon Chung, Sun-Il Chang und Euisik Yoon. „Normal and Shear Force Measurement Using a Flexible Polymer Tactile Sensor With Embedded Multiple Capaci-tors“. In: Journal of Microelectromechanical Systems 17.4 (2008), S. 934–942. DOI: 10 . 1109/JMEMS.2008.921727.
[29] Imada CO. ltd.HV-500NII series -Manual Test Stand. Produktdokumentation. 2007.
[30] Jähne, B.Digital Image Processing. Berlin, Heidelberg: Springer Berlin Heidelberg, 2002.
ISBN: 978-3-662-04783-5.DOI:10.1007/978-3-662-04781-1.
[31] Jain, A., Killpack, M. D., Edsinger, A. und Kemp, C. C. „Manipulation in Clutter with Whole-Arm Tactile Sensing“. In:The International Journal . . . 32.4 (2013), S. 458–482.DOI:10.
1177/0278364912471865.
[32] Kajita, S. und Espiau, B. „Legged Robots“. In:Springer Handbook of Robotics. Berlin, Heidel-berg: Springer Berlin Heidelberg, 2008, S. 361–389.DOI: 10.1007/978-3-540-30301-5_17.
Literatur 97
[33] Kalamdani, A., Messom, C. und Siegel, M. „Robots with Sensitive Feet“. In:IEEE Instrumen-tation & Measurement Magazine10.5 (Okt. 2007), S. 46–53. DOI:10.1109/MIM.2007.
4343567.
[34] Kang, W.-X., Yang, Q.-Q. und Liang, R.-P. „The Comparative Research on Image Segmenta-tion Algorithms“. In:2009 First International Workshop on Education Technology and Com-puter Science. IEEE, 2009, S. 703–707.DOI:10.1109/ETCS.2009.417.
[35] Kappassov, Z., Corrales, J. A. und Perdereau, V. „Tactile sensing in dexterous robot hands -Review“. In:Robotics and Autonomous Systems74 (2015), S. 195–220.DOI:10.1016/j.
robot.2015.07.015.
[36] Kheddar, A. und Billard, A. „A tactile matrix for whole-body humanoid haptic sensing and safe interaction“. In:2011 IEEE International Conference on Robotics and Biomimetics, RO-BIO 2011(2011), S. 1433–1438.DOI:10.1109/ROBIO.2011.6181491.
[37] Kinoshita, G., Kimura, T. und Shimojo, M. „Dynamic sensing experiments of reaction force distributions on the sole of a walking humanoid robot“. In:Proceedings 2003 IEEE/RSJ In-ternational Conference on Intelligent Robots and Systems (IROS 2003) (Cat. No.03CH37453) 2.2 (2003), S. 1413–1418.DOI:10.1109/IROS.2003.1248842.
[38] Kinoshita, G., Oota, C., Osumi, H. und Shimojo, M. „Acquisition of reaction force distri-butions for a walking humanoid robot“. In: 2004 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) (IEEE Cat. No.04CH37566). Bd. 4. l. IEEE, 2004, S. 3859–3864.DOI:10.1109/IROS.2004.1390016.
[39] Krishna, G. und Rajanna, K. „Tactile Sensor Based on Piezoelectric Resonance“. In: IEEE Sensors Journal4.5 (Okt. 2004), S. 691–697.DOI:10.1109/JSEN.2004.833505.
[40] Lee, M. und Nicholls, H. „Review Article Tactile sensing for mechatronics—a state of the art survey“. In:Mechatronics9.1 (1999), S. 1–31.DOI:10.1016/S0957-4158(98)00045-2. [41] Lepora, N. F. und Ward-Cherrier, B. „Superresolution with an optical tactile sensor“. In:
IEEE International Conference on Intelligent Robots and Systems. Bd. 2015-Decem. 2015, S. 2686–2691.DOI:10.1109/IROS.2015.7353744.
[42] Liu, H., Greco, J., Song, X., Bimbo, J., Seneviratne, L. und Althoefer, K. „Tactile image based contact shape recognition using neural network“. In:2012 IEEE International Conference on Multisensor Fusion and Integration for Intelligent Systems (MFI). IEEE, Sep. 2012, S. 138–
143.DOI:10.1109/MFI.2012.6343036.
[43] Liu, H., Song, X., Nanayakkara, T., Seneviratne, L. D. und Althoefer, K. „A computationally fast algorithm for local contact shape and pose classification using a tactile array sen-sor“. In:2012 IEEE International Conference on Robotics and Automation. IEEE, Mai 2012, S. 1410–1415.DOI:10.1109/ICRA.2012.6224872.
[44] Lohmeier, S. „Design and Realization of a Humanoid Robot for Fast and Autonomous Bi-pedal Locomotion“. Dissertation. Technische Universität München, 2010.
[45] Maekawa, H., Tanie, K. und Komoriya, K. „A finger-shaped tactile sensor using an optical waveguide“. In:Proceedings of IEEE Systems Man and Cybernetics Conference - SMC. IEEE, 1993, S. 403–408.DOI:10.1109/ICSMC.1993.390885.
[46] Maxon AG.DC Motor RE 35 Typ 118777. Produktdokumentation. 2002.
[47] Maxon AG.Encoder HEDS 5540 Typ 110511. Produktdokumentation. 2017.
[48] Maxon AG.Planetengetriebe GP 32 A Typ 166158. Produktdokumentation. 2017.
98 Literatur
[49] Molchanov, A., Kroemer, O., Su, Z. und Sukhatme, G. S. „Contact localization on grasped objects using tactile sensing“. In: 2016 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS). IEEE, Okt. 2016, S. 216–222. DOI: 10 . 1109 / IROS . 2016 . 7759058.
[50] Nakamoto, H., Fukui, W., Kobayashi, F., Kojima, F., Imamura, N. und Shirasawa, H. „Shape classification based on tactile information by Universal Robot Hand“. In:2009 35th Annual Conference of IEEE Industrial Electronics. IEEE, Nov. 2009, S. 2360–2365. DOI: 10.1109/
IECON.2009.5415419.
[51] Oppermann, K., Arminger, B. R. und Zagar, B. G. „Smart hysteresis compensation of a magneto-elastic force sensor based on Terfenol-D“. In:2010 IEEE Instrumentation & Mea-surement Technology Conference Proceedings. IEEE, 2010, S. 662–667.DOI:10.1109/IMTC.
2010.5488259.
[52] O’Rourke, J. „Convex Hulls in Two Dimensions“. In: Computational Geometry in C. 2nd.
New York, NY, USA: Cambridge University Press, 1998, S. 63–100.ISBN: 0521640105.
[53] Papageorgiou, M., Leibold, M. und Buss, M. „Minimierung einer Funktion mehrerer Varia-blen unter Nebenbedingungen“. In:Optimierung. Berlin, Heidelberg: Springer Berlin Hei-delberg, 2012, S. 67–120.ISBN: 9783540340133.DOI:10.1007/978-3-540-34013-3_5. [54] Preparata, F. P. und Shamos, M. I.Computational Geometry: An Introduction. New York, NY,
USA: Springer-Verlag New York, Inc., 1985.ISBN: 0-387-96131-3.
[55] Schuetz, C., Pfaff, J., Sygulla, F., Rixen, D. und Ulbrich, H. „Motion planning for redundant manipulators in uncertain environments based on tactile feedback“. In: 2015 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS). IEEE, Sep. 2015, S. 6387–
6394.DOI:10.1109/IROS.2015.7354290.
[56] Seminara, L., Capurro, M. und Valle, M. „Tactile data processing method for the recon-struction of contact force distributions“. In:Mechatronics27 (Apr. 2015), S. 28–37. DOI: 10.1016/j.mechatronics.2015.02.001.
[57] Shill, J. J., Collins Jr, E. G., Coyle, E. und Clark, J. „Tactile surface classification for limbed robots using a pressure sensitive robot skin“. In: Bioinspiration & Biomimetics10.1 (Feb.
2015), S. 016012.DOI:10.1088/1748-3190/10/1/016012.
[58] Singh, I. und Neeru, N. „Performance Comparison of Various Image Denoising Filters under Spatial Domain“. In: International Journal of Computer Applications 96.19 (Juni 2014), S. 21–30.DOI:10.5120/16903-6969.
[59] Strohmayr, M. W., Saal, H. P., Potdar, A. H. und Van Der Smagt, P. „The DLR touch sensor I:
A flexible tactile sensor for robotic hands based on a crossed-wire approach“. In:IEEE/RSJ 2010 International Conference on Intelligent Robots and Systems, IROS 2010 - Conference Proceedings. 028056. 2010, S. 897–903.DOI:10.1109/IROS.2010.5650191.
[60] Strohmayr, M. W. und Schneider, D. „The DLR artificial skin step II: Scalability as a prere-quisite for whole-body covers“. In:IEEE International Conference on Intelligent Robots and Systems(2013), S. 4721–4728.DOI:10.1109/IROS.2013.6697036.
[61] Strohmayr, M. W., Worn, H. und Hirzinger, G. „The DLR artificial skin step I: Uniting sen-sitivity and collision tolerance“. In:Proceedings - IEEE International Conference on Robotics and Automation. 2013, S. 1012–1018.DOI:10.1109/ICRA.2013.6630697.
[62] Strohmayr, M. „Artificial Skin in Robotics, A Comprehensive Interface for System-Environment Interaction“. Dissertation. Karlsruher Institut für Technologie, 2012.
[63] Süße, H. und Rodner, E.Bildverarbeitung und Objekterkennung. Wiesbaden: Springer Fach-medien Wiesbaden, 2014. ISBN: 978-3-8348-2605-3.DOI: 10.1007/978-3-8348-2606-0.
Literatur 99
[64] Sygulla, F., Ellensohn, F., Hildebrandt, A.-c., Wahrmann, D. und Rixen, D. „A Flexible and Low-Cost Tactile Sensor for Robotic Applications (PDF Preprint)“. In:IEEE International Conference on Advanced Intelligent Mechatronics (AIM), accepted. 2017.
[65] Takahashi, Y., Nishiwaki, K., Kagami, S., Mizoguchi, H. und Inoue, H. „High-speed pressure sensor grid for humanoid robot foot“. In:2005 IEEE/RSJ International Conference on Intelli-gent Robots and Systems. IEEE, 2005, S. 3909–3914.DOI:10.1109/IROS.2005.1545126. [66] Thrun, S. und Leonard, J. J. „Simultaneous Localization and Mapping“. In:Springer Hand-book of Robotics. 2008, S. 871–889.ISBN: 978-3-540-23957-4.DOI: 10.1007/978-3-540-30301-5_38.
[67] Ulmen, J. und Cutkosky, M. „A robust, low-cost and low-noise artificial skin for human-friendly robots“. In:Proceedings - IEEE International Conference on Robotics and Automation (2010), S. 4836–4841.DOI:10.1109/ROBOT.2010.5509295.
[68] Vukobratovic, M. und Stepanenko, J. „On the stability of anthropomorphic systems“. In:
Mathematical Biosciences 15.1-2 (Okt. 1972), S. 1–37. DOI: 10 . 1016 / 0025 - 5564(72 ) 90061-2.
[69] Wahrmann, D., Hildebrandt, A. C., Wittmann, R., Sygulla, F., Rixen, D. und Buschmann, T. „Fast object approximation for real-time 3D obstacle avoidance with biped robots“. In:
IEEE/ASME International Conference on Advanced Intelligent Mechatronics, AIM2016-Septe (2016), S. 38–45.DOI:10.1109/AIM.2016.7576740.
[70] Wittmann, R., Hildebrandt, A. C., Ewald, A. und Buschmann, T. „An estimation model for footstep modifications of biped robots“. In:IEEE International Conference on Intelligent Robots and Systems(2014), S. 2572–2578. DOI:10.1109/IROS.2014.6942913.
[71] Wittmann, R., Hildebrandt, A.-c., Wahrmann, D., Sygulla, F., Rixen, D. und Buschmann, T.
„Model-Based Predictive Bipedal Walking Stabilization“. In:2016 IEEE-RAS International Conference on Humanoid Robots (Humanoids 2016)(2016).
[72] Yousef, H., Boukallel, M. und Althoefer, K. „Tactile sensing for dexterous in-hand mani-pulation in robotics—A review“. In:Sensors and Actuators A: Physical 167.2 (Juni 2011), S. 171–187.DOI:10.1016/j.sna.2011.02.038.
[73] Yussof, H., Takata, J. und Ohka, M. „Measurement Principles of Optical Three-Axis Tactile Sensor and its Application to Robotic Fingers System“. In: Sensors: Focus on Tactile Force and Stress SensorsDecember (Dez. 2008).DOI:10.5772/6620.
[74] Zwillinger, D.CRC Standard Mathematical Tables and Formulae, 31st Edition. 31. Aufl. CRC Press, 2002.ISBN: 9781420035346.