<?xml version="1.0" encoding="UTF-8"?>
<Worksheet><Version major="6" minor="1"/><View-Properties><Hide name="Section Range"/><Hide name="Group Range"/><Zoom percentage="100"/></View-Properties><Styles><Layout alignment="left" bullet="none" firstindent="0.0" leftmargin="0.0" linebreak="space" linespacing="0.0" name="Heading 3" rightmargin="0.0" spaceabove="0.0" spacebelow="0.0"/><Layout alignment="centred" bullet="none" linespacing="0.0" name="Author" spaceabove="8.0" spacebelow="8.0"/><Layout alignment="left" bullet="none" firstindent="0.0" leftmargin="0.0" linebreak="space" linespacing="0.0" name="Heading 2" rightmargin="0.0" spaceabove="7.9992003" spacebelow="2.0016"/><Layout alignment="left" bullet="none" firstindent="0.0" leftmargin="0.0" linebreak="space" linespacing="0.0" name="Heading 1" rightmargin="0.0" spaceabove="7.9992003" spacebelow="4.0032"/><Layout alignment="left" firstindent="0.0" leftmargin="0.0" linebreak="space" linespacing="0.0" name="Normal" rightmargin="0.0" spaceabove="0.0" spacebelow="0.0"/><Layout alignment="centred" bullet="none" firstindent="0.0" leftmargin="0.0" linebreak="space" linespacing="0.0" name="Title" rightmargin="0.0" spaceabove="12.0024" spacebelow="12.0024"/><Layout alignment="left" bullet="dot" linespacing="0.0" name="Bullet Item" spaceabove="3.0" spacebelow="3.0"/><Font background="[0,0,0]" bold="true" executable="true" family="Monospaced" foreground="[255,0,0]" name="Maple Input" opaque="false" size="12"/><Font background="[0,0,0]" family="Serif" foreground="[0,128,128]" hyperlink="true" name="Hyperlink" opaque="false" size="12" underline="true"/><Font background="[0,0,0]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" name="Text" opaque="false" size="12" subscript="false" superscript="false" underline="false"/><Font background="[0,0,0]" bold="false" family="Times New Roman" foreground="[0,0,0]" italic="false" name="Bullet Item" opaque="false" size="12" underline="false"/><Font background="[0,0,0]" bold="true" family="Arial" foreground="[0,0,0]" italic="true" name="Heading 3" opaque="false" size="14" subscript="false" superscript="false" underline="false"/><Font background="[0,0,0]" bold="true" family="Arial" foreground="[0,0,0]" italic="false" name="Heading 2" opaque="false" size="16" subscript="false" superscript="false" underline="false"/><Font background="[0,0,0]" bold="true" family="Arial" foreground="[0,0,0]" italic="false" name="Heading 1" opaque="false" size="18" subscript="false" superscript="false" underline="false"/><Font background="[0,0,0]" bold="true" family="Arial" foreground="[0,0,0]" italic="false" name="Title" opaque="false" size="36" subscript="false" superscript="false" underline="false"/><Font background="[0,0,0]" bold="false" family="Times New Roman" foreground="[0,0,0]" italic="false" name="Normal" opaque="false" size="12" underline="false"/></Styles><Group><Input><Text-field layout="Title" style="Title"><Image height="78" width="800">TUZOV3RLVWI8b2I8Uj1NRExDZE5WWlpKOnROPkg6eFhWRXJwczo7Qk5TRE9FVGxNWGxnd2dpVzttRFtVVVVXVXNLaXRVZl1XZnZfaXZtaXhvWUtFVmNzSXl1eXZheXZVSXZfaW9peG9PV2tneHdpeXdPdmVDSHdnSXhpSXhteXFBWXNdSXdnWXRVaXVJWHBDSUZpU0lhQkFBc2E7R2JZeXZjaXhxeXhlWXdleXVZeXVXZE1XVHVVWXV5eXl5QTs6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6OjpaamlmRHFFdGtdYE5cXEBOZFxcUWdxeEhgandoU1dEUVZ5UHhQTEFJWFVgd3l5eVNVdW5gcltETlpdV21ValB1SlpdWV9sU0xxcVdpb1N4d3d0TEVRbEBVTkdpT0NAWHlRalhMWWJJdk48eHdhTG5BdD11T1pkUW5BdEU8U0lkUW5RSkxZUklkcTpgeEpZcnlxSkJoeU5Gdkw/Xl5Zb09BW3lZZWxvZmlHYnQ/d1t3W1BoZEs/Z1NPXkRHcExZZUpwXXQ/ZmpIb1xcSV86eW87SF1cXGBcXDpHb0RGXWBocUVodD13W0ZfYWxTPXdVVG9UdE9IUHdDYm9yWVt3Oj1FcFlkUllyWU1DaEtkRT9CRG1pZEtHPVFzQ19ZUm1IblFCTFlyP1FlRV9YX2tyaWdlOltpQlljZl9ERGFHZVNzXFxlVFBPYl93WXJ3c1hpcmRJdmlHYk53R107VFllVEttZ3l3dkpHQnN5Q3ldVmxtRmV5RVF3Y1g9amp5eDpgc1FNUF5cXFlQaG9fVGs+eE1zbXRzSVBNaEttWUx3TVh3SVdYcU14cVBJVWtFUVQ/bW9EaHRIRW9fbFlAbVFIUXBaRHlMVXJZSHBuPHlSdXRuSFV2PGxweEtZUFd3SVhSX3BgSWBwWGZXT3l5PmVNeV9KV3U9cWFSPnBwVnhPXmZ1bnI/R2BIdl5RaWFddnV1b2NKcHdVUWRUZ2RgX21leF1UdmZcXHhmcmhkYlh2cEllX0hzW0lpSD5uVW9udkBiS3BpWkh0WGBpYmhmS09gSkZkUFBrSXF2eV5xPD9tQHZ1dkFba2BmRGhia1lkTnF4al9jPl9mT2Z2X3dkeF5fRT91WVh5UUBvbEZxWUlmO19lXUl5UFZxbm9zZlB5SkFePWFzdXFbamBaUj9rRV55akh0SFFnT0h4U25cXHdZb0loYFRZXFxGZ2BSeGBJcVtWd3E6QF1UeXliUXh2QF1rPmtpdmRhWVxcdWlcXGRXaXJuW1BxclRncFBZYnhedHZGZmtXWmJpaGxZYT5eYktAd1R3c1Fodk95YkA/XWdxaHdvbW5nXz5vZz0+d3BHYXJBY11oaWJBeVhAZUxvZ1FuaGx5a0Q/czxfY1xcPmJAUXV2QV5reG1ecHBBWHZqVlpzRl5BRm9eX25WVmZsaXhyaWZoYXFpP2JISVxcSmZfXU9dc15gbHlzc0FzcF9iPUlaXWFrZFBtSm5pQXZeUG5hTkF3OkdpPlZxbWZ2Ukl1eUZfW05tcFFqYz9wSXFeUFdqaUZkVVlyYz5nbFBxaFBbQj9qTE5xS0F3eXhuVmhxXFxhallRXlpGVlF4az9lO19mQFViSVNzPz9UPGFCdz1mS11VeVl5W29SQU15Uj1Id2l3RVVIZm1SUFN0eV1Uc1N0YkFIeFN1WU1zXnlHS1V1PUlCPVF4ZW1VQT1ycndJO2FJWD1CSj9iXnNzW19UWEVZVENlRWt1R2dDTmdlRUtZOnl4RUtCTFdibXVCSGt2ak9ndmFjSV9XPV9kR2t0UmVnWXdyXVdGUT95VEtCQlV3SVtIVFlyQnlHanlGXFxXYndndnddU3hhd2FhV3M7eUF3VENBU155eGQ/WGQ9c0JneVJhRFc9RGpzVDo9aFxcS2dtTUdbYXZcXEtkXXNUSkVjdltkVjtmdmNoO3dTOl9Ea1l1XVF3T0NkTztzZz15b2V5dFNHYGtJbXNGeW9nXj94RU9CTENGVmlESWdJQEd5XW90XmlyUDtISz9oWk9zamdTXFxvSD9FVVN1REdNVUF1RkpJSGlfRktTV1V3UlRbaG89U3VjY147SXNfVlRVRT1JQ29TSXN3Q1dxUlpRRzxfaVVhY3JDZWhPY2FJUld1c3BxUmZZVEBjY2ZNdWhzeUNXclltSVBLSWJRaGRDZWhxeFxcc3Q/XURHXUVxTUlGWWZXXXJnVUNicXZJR1Nnb2ZMV2dgYUhKS2RsdUVxRWV1PWl4a3dRU3RyU1d0V2djZ3dKU0lHa3Veb3hnS1Z5UVdaRXReZ0JlS0daS3hjZWQ9SWRUT2hKRWZSW3hyTUJrS2debUdKXUhjW3RyT1Q6X1I/ZUZkX0ZWQ1haQ0Q/UUNxU1hdWWV0R0Y8RXVRZVVmY0NMTWhqR3ZWS3NfU1RrVXdeXUNFVUVsW2Y8P2hORXdkb0g/TVdmXUZiZXNQS1Vda2dIXWJTRVM7UVZWP2hIcWRUPWNlP2JwX2g9R0d1cUdEW3lAU1U9SUZLRXhFZVVXQWhOTVhed2RSWUZJTWV2S2VIWVdTc0NsW0hHYXVbQUVaaVJfaVRKVURTXVlja1hzb1Y+XUdCcWJAO1ZNPURsdVZIZ1Z1UWVEcXhMVUVdXVdTQVJfb0I/b3hMZ3I9PXZxa1I/TWNQQUVHXVdCS1ZQW0hWT0k+SXJFdUJrVWNxU2NrQ3dwc0ZvX1JjP2VCW2hoQ1hZU3JGQ2hWQVN0W1VVV1dzXWNlWUJoeUQ+YVVUTVdaO3ZEb1I8TWlnUUR1X1R0Q1V1UWVUcVhMT0k+UVZfQ2lJXXdfQ3J1RUhvc1J3b0ZmP0VjUWlKP2JoPHJUdVhbWG0+UU4/WXROZHBQUU1TeFVNPExxPXFASU5CQUtFVEpCaHhTdExzRXE6XFxWbVlNY0VKdkxNYFxcam9BV0tsdkxgb1RFeGJxUmBgdVJxSzs9UFg8TEF1c0NoTz9AbU5FamVhUF1JU1docEB5V2xcXFdjPXk8UWxQWEpRdVN3bFc9eHRZeXZKSE9UdEs7VFc+bE9JRE9EVEpaeU5vVVBSTHhIbHdQZWxLeFQ7dG9SRXZcXEFsY11rYnBwZmB5b2x5dlB2T01reERLPl11XFxFVkNdTkFhbkFZYz1GXktfdWRnZFtRXlZpXkRyPlt0Uj1IXFxhR18/R1RbcnRTcnVbWEJ1R0RzVUthZz9RVUVHRUtDaWdjR01lWW9HQjxVUkJJYltlYnZZRkFLYkd5R0s9Q01RQ1FdQ15bVWtVVEZjWFZFaD1dZzxbVkRvQkFJZ095WENnc1FzZGBDRmM9dWpRSEtdWWM7eE9PaUBZeGxPRlhFYnhPR2VDczxraGFJUlZJZ09tcz1lVE9JeVB5cmZ5QnRxdFZ1eVJFeTpvclBjZVxcSWdxa2JWTVVaQVg+c0hzVXVPa1lxQWdDPXN5b1lBc3ZgS0NoQVheV1JdX3h2Y0Zfa1JnQWNeSWNQW1NJX0RbVWY9TWVyb2ZjUUdvWUJmY2NhW1RpRVR2YWU7PUhHY3RhcXZ1V0hkO0liQWlKRVlkeXRHP2hCb3JkVVRYV0M7ZWJjaXNMW1V4WURDW2drUXJITWdIcWVibXJpa3ZqP0hyU2lQeXJja3hrQ3dRQURDb0llRXZiVWJHQWJvUWhYRWhbO2RcXEtISG1iX09GdFd1X3liX1VjUk90bndiUVVIakV1TD1VcF9SYj1VWUFnVU1FPmdDQWdDaXlTRW9zRVFVR3FkV01XcT1jP0VyS01XSUdGd09DZUd3Xz9jQFlCTT1zYFFlcz9VYEd2RElHdV1EaF9VXFxhRUNRQ2tpZ2BLWV5JZDxVckZTR2RpZENRZD93dmpzZ2pvY2Bhdj9BQlVjQ3FrRGJnVVFtWWRXVHlVSElFSV4/dk89eHJJWUBfSVhLeHlleTpXeV1ZUnJ1eERpR2lTdj91aUhHYlZRbDpES1tIbXJwUEhQcV5QbEVcXGtBTWt2bUxMeWxGQW9rbGpjZXY9bHFpPFlXdFJMZXdxSVFQXW51VGpTcXZvXXhndHI6YFREZW9zYHFzWG9VdHNePFFWQUtkPWxIRXdSUU9WWXFFeVRvYFk6YVlOUEtoXFxWPkFzTlF4O1R4cmRXXllKXnRqYVxcdkhkbmxVa1Jla29ZSnZYT1Zlc09xbFVNTkBtUG5VUG9YbVRAanRtVWRwS29IeGhtdUQ9UUhld2tcXG5CbE9odXFXWG93eXM8XFxWUGRrWkFKZ0VSb2BSQGV2W2V2cFRxYGFTeDxOVUF2eVVyYV1udnRSaUh1QlRRSVRzXFxxVlxceUxsXXJhWHRcXEBQQ0hTQHRyO1xcWG1EUztYbUZwVlJ5WHVIaklNS0JgbVppdlc9TkhUU0FcXHNycHRnTW1JQU5xZXVZQHFKTU9GaHJ4RUxpXXZvbXJQQGtnXVRFTVNORVhyZWxteXJveGtVXFxZeU1sbT1LYEF2dmFYV21LUXFtQTxRVElVXklRaG13XklZSFFxXlxcc1dsbEVdc2tscz1RWXdBUnRwVVBIVldpbj5US3llcWBETGtZUEQ9Vll4T1VpdTxRVG89dT1QVWNZWEx5a0VNbUJIWXd1T1Nkc211dV9kUm1dV2xwTElcXHhLbHF5QEteQU86SUpcXGFvO1lzZEh4UkhwT0B5REBMP0lwTHhyZFVwX0h2Y3B2Z111RVFWS1h3dmRucGBWTnFWVkB0W2xMX2lvO3FPSVVOd0xTZmlKOm10PnlWVHFOZ01Wb2FvUF1STmlWRFFPYEBWRGlzZEh5d3RZO0BWZnRMcVlsc3RyRTx2aG1yQmltVU1yPkVrSkF1R3h2WWlZSm1XeHhZcWRqR3hLbF1UQFFQVllPWWBMSmBtYGFqTlxcTUJNVlFteXNMTkRZc3E8b3BZeVNEbWBBdlBAcUJIbFBpT1xcQXg8cW9cXEBZZVhyU0hQUkBWZVlWR0FTcnhRWllQR3R0c1BrXWVMRWhXbzxQR0FQP1F4WkxYWDx1Y01TXFxsSnlkU1NNd0dAa1FMdmpBTVdUeVV0b3hVTGtVUFhUdTxQUVxcaHNhUGtkUEtOaHVIZGtBdHVDUVBaUU1LRFN2UVZQWXlwTFJUeHlNVFBWTVVVaHFzbW1EcG5jWWxYPU5xbHFreFJwZFBPZWtSeE1wQGtTbFVdSFc7eHQ/PVNfTG09QXRuYExxVVFFRVZCQVdCVW5uPXRCTFh4cHRGYE5TQVRkVU5HSFRFPFdOSU5QeFdOSVJHZXdKVE53SHU9TVlWQHVFQEs8TU1fZVNHRWtfRFBAYVZAbWxAPUxcXEV1Q1B2Y3l3U3BrYUB1W3RRaERwPGVTQEF2bUBVXFxNdj1BUVpEVztNTXdVa1E9bT1hVE1BWV1Mb3V1UzpNTlxceVBzXFxRWG1QVkl3dnFKb1BNVElwckFzXlFSdmxzU110TmRLQ0VsP3hUd21uW3Byal1XTXhLT0ROSUlMXkBzbjxZZkVrWEhRTmRRTXRXTGFQTVFMcVVUP015WnRXVGFSQ1Vsa0BYO0BtV2RLP21uVkRsRkB4dnRMVlFRc0ltXlRXczxvWD1LYUV1RWhZRUx0YF1RckB5VEhYUnR4Qlh1dkRyWmR0Xk1uVkhYSWFSeHFMS0xTR0FTTUh3XUBqZHlyTlRNPT1zZV1yYF1vRzxLXFw9VlBgWUtEalh1VHVJakVAd0NRU3hkTV5Ad1BQUz1Nc2I9az5MT1tcXG87dHNtXVVDSUxkRVZqPFM7UFRFaU5VTVZNbW9QdUpHTFRIVE5HcFhLUEtnRHBKbFVUSEtzdW88UGNKbl9jeHBARndhZ1pOWV9XcGVNYHFBV2ddaFxcZklzQWBiWl5hdGl0dXdgPkFpYXloW1ByRVFpZ3BiTU93eXZhSnZ4OkhnYllnO1htPE9yTW9nZFB3XFw+Xj8/a05GYVZYcUdQXmR5WndGcldHeEtuW2tmZ0w+YEdZblBZa2R3YktxYllYcHBoaE9Hcz55W1tGc3VWXUF2Okd1VktHWD9ydG1iVT1VQXlCWFFWSU93RHFFS0lFc29lOmFkXWtYSmF2UkdkZFtCd01jdUVZXFxlRURleF1bZHhPZV1BdVJJZEJLdlNbRD5DWGdBVkhdQndVVUdzQllpeEJ5ZlZ3dnJrU2FdQm1HV2djZnFfaWw7aGdpZzxBUnVhckh1aE5RZHFrSGtLV3FBZHBFY0dvR0NNSV5pd2FXY1d5RlNtU2xxc0lfV21nR2NxZVZpc21xYm9XaFlXR2RXU0ZHQ253dmtPSD9rRlBDc2pvaGFhRG95Z0I7Y2VzUlI9SFBBY0FbaXZhck5nckdReWY7Y2VbY1NhYj5vRjxlWEpdYkE/cmdJVUpLZXRtc2NxQ3FFVUxtVTxNZmVRd3J5dEtzWHRzdmd3ZlI9c3NlZV1lVnI9aVhnVHVzUkVLSXJhd1ZhZXNvRFtdeXFBZ1Q9VkU7WHdnVFg7Q3NpY21zaFpJeWhddGY7RD5hZ1lrU1o7Y0VlZU9VYnBnQktdQ15jY2g9WFtXY3c/VEtdRF9rd1pjcl5hSU1haXZdZ0M/dVJbaUFlR249Z1dJUkRZeGFBZ0x5aVVbY29zZlRNc2BXSEVDZVdpaFpTaVtTU0g/dkRdQ0NTeUpTZUxhd3A9ZEJfeE9PeEtxZ0BjeV4/YmJfd0pTQ2BVWVRbYk5FR2RdSFppSU9BYlZPckplRWNLSWo/Rm5hV2dBZWxTQ0FfV0tddFxcQVVBcVNNTURoU1ZSZ2VwP0JGTXVda1JEa3ZPS3dgXWl5XVNzQXhuVXQ+U2ZpO0lXdXV1R2V2PXJoV2RAcWdnWUZIdWNyPWZmb3ZMUUJfW2U+PWllQXdqQWRcXGloRFV0T1NIRmVXXndpQWVXbk1DY3NYc21oXW94Qz9nXz9mWWdybnNTdmd1dU1iW0lmP0NmYG1SUEV2WXVzbDtYUT1TcW1EUkFzdWN0UEd2c2d4SGdCS3VpXFw7RnZdR0xBRF5bQ29NRT4/U09XVWlnVEtVWVRZRFRjSVZDR3Q7YEZPXVVma1dWXkd3a2FXYHVma0tQWnQ/XFxEd15vT19cXF5pd0Bibl9tdk5eXT92S3FjYHZgdFZzPEB5Zk5fZWBcXER2c0ZnYGlZY0xIdHdfYkt3WmteaEtWbmhWZm1nXFxzSXdRUGZBR2BQbmh1aWtqd3c9YXRYUGJkd2VOQGpOV3VHYGFnaWJeVmlZZmVhUGlyP2c7X1s9T3I6eFxcdnFpPWlxWmlaVHZrQGFmWHZwZHFsU1dcXGlweWRhcllHXFxKRm09WV94Z2U/X1pCcWNESF9Od3JOWHhrUG49V1tsWW49YV89SFpbaW4/UG5BYXlUP3l4dmlpeGJoWHI7aWdkSGlvb1FdZURXZmJpRkBTeD1jdHZ3dVtRV0lpSV5zRV1XclBbRWJNdWdXUm5zQmRbeVJbZD9ZVT1NeWJjaGFkcVI9bG9NdndEU1dZUFd4UlxcPUxPZFdqPExJbXVMXXc8WE5BRVJPaHdLSEx4bUtmdFN4WVVmPEtFTXh4dExObFlLeHI+dVF0cFRWWVV0XFxtdWRUYVxcdF1lblpEUXBgWWtxU0lAb11xT3BgbG9dUEN4eHhBa3J4eF9Ec1RJbVFIU2F5VGhhUkVcXEpMdHB0cVJBdVF5WE9OZFV5VU47YXhbTVdzeHhOPXQ/PUpfeVU/UEs/UXg6SG5fWU5rPFNMSXlcXFVVRFlPaGh2Ok11Zm1Na0lPclVKcVhzbERyVUB3WXFRRWhMS0FrW1hua3F2X114Y1VLbD1UYUVuVUVOUG12SGlwVnBTSFl0Y1BRS3hUSVV3YFhYZUFxW2lYeWhxWGRXXXF1U2xrdVVOVFF3W2h5V2FsUkVTUWhLdWRNUHVSTkhRVlRwWjxWQUhuc0hsWXh2Xm1tVkhKdERZb1R3XWF0XFxYWF5ZeTxMakpga1ZAb110Uz1YSnFNUV1FcHFgc1ZJSl5wV3RQS2NMbWtkTklodVo9V1l5dmhUTnJwVFBpTXdUUWFxcVpFUm5YamU9bklgVXg9cFB1dHRdWF5Ad0dVS0x4als9cUpRblldVk1MVUZNcnFgd0lEc1hoUGNAUVxcQXlydFR1SE49SWtmZFVlcU9bbExxUXBIdHY9cUo+cU5ObVJVVUxeXXdSVExEQHFdYHNXaHhrVEs8aE5eZHJ4UW1rQVJGUXJWQXBKWFlybExXaFZLcXNgcHhSTFR3dVFqZHFmPHdSPWxPbXVfXU5cXGFYQ2lzZ21KdERYc2RVdFVMPFBrYFhwXjxVUlB1RVxcVHlAVUdES1hNS2xRTTpwdDxgbkh5dmc8a0V5bl5sVkZJUT5xUG51d0JleXJ1VG1IWW1YSl91ckRwS3FJUnBQTEx4UlZdSnRMU2t1anhtb2tFbHhNdXhBWE5ZV2NoUD1oUlh4VWpwdnFxbkdFbnY9WWpRclpUUF5lcFRxamFcXFBGRFludWNTUGRpZmxWZ1tqd1tId1xcal5fb3d1VlBkZmdfQ1ZnZFhuSGhoa1F3TVZzaGdaVHhvbFliaGBvakhxd19gZVhaXFw+d1hPbmVgbT9nb0xfd09uYEJ3YGFfdmZ5eVhHdUpHdWdmc29gbWdpdnRIbVheY1NwbVFhZlxcXl1ueWg9b1p4X3dQWG5PaXRyaWI6WHdZT2hwV3ldcWRsV3Z1YGNYPl1KZ3NtPnRkcXNzbl9GP2FuZnlOaFpLZ2dcXE5nRHlwO0FoOl9saEFzW3Z0REZcXE13aFxcZ3dCQWxbeWJNWF4/V2VdWWRFblp3aHk6UXc+YXV0QF9sT2w6PmhmZ2FveHVGUWJLbmJZSHBIUW9id15DX25XQHFEbnBjUXFFR2F3VlxcYEBybnBjbGhjaz5eWEdkTkBxZEF1W0ZmVUledV9cXDpgcWZ2cT9fc29HXFxVZ3VBQVxcQW5da1BsRk5kQkBzS1ZwZE50SF5nQWZvaXBkYUdkRUdsUHdiSlB0W09zUW5eVU5gbUZaWnZsbm9iPnlnTF53V1ltXFxWaGVWZU1HalBockpIZW5JYnBAeFxcd2VdWG9jYGBocGVgeHA6dnVYd2VNWWdbUHFUcG5pSGBvb1tKZ110P3NpQGBwdm9mSXRzblxcXklkYG92VmFnQXFsYUl4VkBdalZdZHZhUUZhbF9oYm93QU94RGBfYVlqSmhsb3FrV1lsSl5mQWZiaT5sTVBgUU5mW2dyWD5yQF9uSF9qX2FeVE52b3hpSlZyc15ldVBjb0BcXFFPW09fcEU+Z1lQbUBfbW9QXlVRX0JwZkVOY0hgak1uWmlZdG14YFZPZ3h2XFxmT3Fob2RAeW9XQW9ITmteV2JDWWRzT2hyeWdKbmRLdnFWWGJSQF1pPmpBSHlXXl1oP11meGdDSWNObl1Jb15sTndIRmY+QGdZQWtRVmNEQGlCP1xcVUdyVFZeaGZqRGlmZ155dEF5SXZcXFFgbXlWeGB2X0RRWlxcSHh0YF5RcV9zUW1AaGRDbnRUQGM9eGZnQGBVWW9cXFl4eGZwZ1lqSElgZGdnWW9fcWB0aEleV2BhPUdyQmhlVVZvUHdraHh5ZEdaU15ucD95Rl5tR2hodmhdVElePHFod3FeSEZcXHNRcFZHdG9vW0dhYklWXFxmQGZCeXdDP2pPd29HRlxcY0Z5cW5tbU5oZXduOndrZnhvYU9pcGhvOndfXnddR1hpQF54aVFlcUZpT25fZ0Feb1ZwVVluPE54RWdsP0lpZ2lgWlFobEd1V292QV94bmFcXFhOczp5Yl9QcHJYXkdpdjtGaHF4ZzxJdGVeZERGYWpIZlN2b1FZaT9XeFpQZGNJX05HbT1pWl5JdmA+ZFk/cD1xaG1QcD0+XU9gYklRd05nZWxRZD9WYllAaV9PYFxcSWJESWVaZnJtcGJsdmxaZlp5XnN2bnNuSWhtTmhbYXBqVmJtVmZVZlpbQXQ9YGZCZ3ZLZmdXeGtiP2Nmb2pkR3ZyaGlMZnZgWV9DX2RpcGdYd29DWHRzSGxcXG5dTlBabU9geVdcXGVeaFRfeERGbGhgW1BJXlpucFdwbURnWj1fY0dmY2NWdlpubkpZa1ZvZmdeaGxXdz5wYVxcYGxNcGZIUGpDUGo+R2ZuT2BUPmljdl1bZmpAdmt0UHJvbnltUFRnZGJQcD55TmRUbnBkb3JRbVRheTpEb0N4cmBpUD5RTGNoTkBEWERUcnlYeUldakc8dUVoTEBUdUE9bGV5T2ZcXFhyUHRacE9pbVlpcVFnUXJxWHZybFdxdXF0TFB2UXlwXXI/UVlOdXJgdU5XZVJgeHFeSGpGaXBnVVlYRHNBWXhOVHJtQE9KcVBJWW46ZW9GTVhZRXRjTFBScUxHSHdLbG5hVU1wSG9jTXdOPXlaVmZFXl5fSXFcXEZuQ19jVEhobldzVz9vTkBuYlBdXWhydkdicz9vcW5tQl9hW3h2bj5mY0BfRU9pPlhoTmZwdVZhQHhoTkljXm9ybUlxZmZvRl9tZkhjZ3lkZ05dX192OkhbcEBfUl5eS2BdZWFvTWlqUldcXFpPeWBoaEFwZVVwbUJoYEBAbVFgbltmaE1xYlRfb09maEVwaHRBbGA/bG5pZGJAdmh3YnFeeFZHY21vYG5odUphdmtOWkleYDxnZ2hIb0Q/YGhfbUlReUlIcUZWZ1txX1Rfbl1YYl5IXFxTZl5dbnRgd2ZXb19WbnZMSGZTbmJzWXl1aWBBUWdjcV1EPnNLQGZxZltDaGlTX1p2ZmZbSHRCUmxnRmA/c09TU3JxdmRBSUg/eHRtVV11V11HR0RhUmlTRj9NVHBLdm9raGJZRj5RaGhLVU4/RHk9VmZ3dj1FVkdXaXZvZUtfdUxhZ0RdclZDdE0/ZD89cz1HUlFvWWZJWTs9VXFvZXJNRHNnSUxvUz5FWDxtQ2tfR0ZBWVxcaUV2SURnZ0VBcWVocWdkV2hTYVNPVWRBVUJxZXlNTWNvO2dmXXdPX3Nsc0dFYUlsVXRvZWdnP0VhSVRrU3VgUXJ3b3Z4S3lea2VMXXJdR3ZqW1dKc2hPX2lnQURCd2ddYHZTbUxMSXNSSG9DZExsSE1GUHhmcGtMUHJgXWptUE9GUXh2SWtkTG4+eFZReXNEdHNFTXB0eHZLYU5gYHRobFJsYHFJYFhpaW1TZUo6XFxTaUxVPG1SXUhyYl14QW11ZFRrV2VQUWl1UTxsb3Vydl10TXRXVWFqaXlvOnBuVz1QRT1vTERRXnlMPmBvR2BqWHRMcGx4c2BtamB0W0hPRFBua0RPdGV0dEhxVkhLVT1UZ2xOYWRSXnhVS1lSO1huOzxZR2BQdFFZWFRPbFB0U0RRTkVsaU13XmRzdklwZnRQVWlrSGFsPkxWcXhrWHhSU2hVQ2BUaXB2YUFOO3Bzc1BPRVFwdGxMXm1ySz1NS3lUQ2B1Yz1tRWRwUmB1PUFxZXhqckBSa1VQcVhVR0FMbl1sPkh2Y1hLREB5Y0l2O1FXS3RTVW1VX1VSQmBrUmBNQkRsWFRuYkxucT1ZWFVuVlB0dXR5T2BReGBKZm1VYm1RR3hLbFFtVFFsTFVNaT1SWlhtRmVOYGxYVTx1UWF0eERZWHFRTGFWVVhyZHlVS015QWhrTVFxVERqdXFTVHhwSlRLQnF2XlF5bmxLXWRYbEBzWHBvZ3hSXnFPWGR2R3htWVlXSU1uZlxcT1l1dWVMTT9pUm5lUnd5d2NRUnJQbkZMUEc9VlJQT15pS108UnFwTVlscXVBeFdwWXdxS211ZnlmaXE/XmNAX0RObm5Zd0BvXWs+W15pdUJmaFFJeXJOdUBpbmxwW19XZnVeX1lwXUVBYWxfeV52ZV9gYlBeYWdnYT1BZVJ3YUhBalFvZURPeVZeZV1vcVVHYHk/Y2h3XT1OeHhPXFx3VndabmRrP2JWPnBQVmpQWWpEbmdAeGM9cXBRX2NPSFxcUT5cXF9mYEdmbV5vZG5RYD5YWmRXZVxcR21GRnZzUHVhcHA9YGxGVnRoaVtCeWE9TnFxRnhZWWNtcVxccGZzYUd4O0FcXERRbGRQd0JQcU1mXFxVb25zVnVoSGhYXmZjTm9Rb2NlT2toSXN0aV47cXRaR3RPbnhYP152cWRZUHBqb11oUGw8cWNAVmtUaGVHT3dFb2tuQFtOYWl5Xm1CUWpcXElyY09tT25rPG9rXXF0PXFyPVZvSVZfZF9qWHZ1cWlvd1h3W25jeEhfTFlddGdwdW5kUE5lR3Bpaz52W2l2UkZ0X0d1S25xRkl3XFxuW0RwaHRgdmdJa1phdl9RXU5xdmJBa21IcUFxZ1VPcV9IXlpYakU+Y0BPYklYanRub2ZHbz9xcV5nZWs+dXRAaXVgdHA/bFxcSXBPaHdDYXNcXGFvbmZpSXZtP1BkcVFnTU5oVlFqRmdlREFkU0ZmXFxWZURBZ2hxbFJHXFxVRnFNWWVwcF54dmxQUHZeR15lQGRGP15MQG1LSXVyUV5eeXhLRmlLbm1YR2o6SXZSeW12PmxnT2xyP3E8aHRKSV1sd1tecXVMQF9cXEdmdWdwPD94R3dbbF93P2llZ19pam5cXEReXFxlSW5Cd1s8Z1puSFxcb1FfRkFiVm5hRmFhUmBfPkdgckhuYV5za1BbQEZqaGloPmBiQnBiZ3Z2ckFhQnhwPUFxb3FiQHdfal5xO3FhXFxWdUZxW0B4W1ZAWz9BZ0peW2t4ZXJQZW5ocWZXYWZXeUBxZm1xdmN3WmtuZ1thd3FAb25hXkVHXUtncmJZXV9gXU54XmZHWnJPeGVBYFJGaDpAd0ZOamZoWl1OY2dpcT1ZZGRXXFxiXm1AQGVyRmNncVtMUW1SQXBlQVtxQXhBPnJUR3JCSG1KcHhveGF5cXFpYWdWT2BjSW5bPlteWHJASWhLPmI6WXdtXmtvTnhxZ1tYYGc8UXlDX2dNYFt4X2Z5PnRcXEBmYnhiO0FgaF5oZFFpS1h1ZkZuRkhdTGhxY2hqXmZ5T2VDW2Zjb3dZc3ZScUk7a1dbO3RlR1ZeP1RJVXdXW2V0R0RPbWhsVWZbdWJmY3lpXUlNb2dpd0heVWg6b3dBd2ZXYVJwbXRDQ3U9SUZkQWVeT2hzS3J1ZWVaRVNaVUJeY2luTUlmd0VfbWZza3hKc0VuP2ZPYWNAP0dKPXVYa0RJaUI6cVlPV0JVS1I/PXU/cWdOYXRoP3ZkR0ZQP3VgVXl3XWZDY2NVT3R2VUdZY2k/Q3U9Z1hwXXhJY0l5T3JBbXl4WXhJY1VOS0ZwUURZWVI7Q3dCPXNQcWlGP2VGP0M+VVRZZWY+PVJKTVRcXEVCQ0dVRGloaEtVUl12Z0FlblV2WmlDZWliV2tDY1VCb3Fkb29jd3lDQ19yQGlmO1NzcWtzd19EVm9CdXlZO3lyZHFmdVNjaEFZQVVncVN5U2tnWnVXbWtEPVFjdll3WVtia3NySnNFYUtXaVd1TT11X2djOkVnX1lmTVl0bnFIb2dSaXNFaT14cHNYcV1ncklocldGaUVDeVNVOm1SZV94Z0VCc1tDXFx5Um8/YzpzaWZDZ25bWEhlRFVnV0tvaWljRVBZSXFxSHhrV3dBZ1ljaHFnRGpnSTthRnBLRFdRdl9PQlpZdD5Xd2RBVVl1Ukw7Y295cU51VmhdVXFwdGRpaztwWGR5a3hUUUdsVFF0dlhgWVY9cTpAVm1AcmZUcklUbFR5all0eG9Aam5VU1NIU3NdcV91eXFdc0NlV05NdUhtSz5cXGpXPWx3ZE5aZG9CYFdzcG9QRExvSW1CVUphYVhYbVhQYUxTZW9uPVV1YW1RcW1wRVVXSHNMVHl2RVhUZU9paVl5aXF5PXFxXVhIQVJIYE1fSVFjPFlBZW5IeVFUSW55TWxSWUtLeGtmRVh5QW1ZeVFUPXB0aHlbZHV4QW1DQW1BeVl3VW9zZHlBYFBrUWtVdE95eHJgSVJedG1NTHZlUEw7RFlaZW5yeXZ1PFlHdU5yYHhrUHFQeE1ybW1DPVNAXFx2YXlKa014bnV0Y2F2PkBqVkR2YGxxZG1MO2l3ZlRrQmVMcWxVV1FQZlVYUVh5WnVMZ21zT3lRPUhWQHVOeGVRdXVtckBRT015a0RNY0xtcz1PSEBzUUlVPUhRbV1yeHVqbHV5R3FweFl1TnFVRF1ZeGFyO2F2Q2BtSFxcVTt0UDtsV2t1UXNJSm14TUdNTHZNTmFZUFNgbXJZT2hxWUpAU3JlUXhdVGJIT3JIbGNJbT9FeXBJam5obGtFd2dZbkZ0WG9oSlNVc2xMcnc8a3FgWGNYdUhtUzpJUFBNV15cXG5lTE5yRVlURUxraHlXRFc6cHZAUFJOQVNBcW9xYFlMeFV5dFhaTXZyYFRqbG1rVXE9aFFQUVJ2ZHRTRVhUWVV1VUxxQFl4QVQ/eVNMbVVWZE5vUU5WXFxvS1F3b1VYPW15SkB2XnlyZklzSmRPUj1YPj1XWEBVUkhqeUlWSXVYSUxMU2Bua1xcTHNkTEp4VlhwbWpQcnB5SjpZcWdNU1V1ckZYcWVhV0R4cFdgWTo9UnBQTExhbz9UUU1ETmxgUGRZcDxtWD55WFhFTEVgd1FtVD5dUWV0eVpFVUV4a3I9Uj95cmFJWFRJWWp1dWZpbnJpS1xcQHlTcUtFeFheXXhdYHJYPVI6TFZocFBCbFBlWVlcXGRMYlVOVEBQWGVTXmF3Wml3QXVralhXPDxtZWRqd1BgcXdlWE5ya29oX1lvPWhia1ddd2B3RUBlTldibHdraFhcXGheck5vYk1Rcnl3aUtZa1l5YHlXdF1pZmh5bHReaHlBeFxceFxceUhrVFlrQnhza195akFoPWh2W19eeWZ1YV9kXFw/a3R4YXI+am93Yll5WmxRYVBWbWpeaXdQcT9gW0FIXl5RYHBBeWB3ZFpgY01weXZpXjtGcm14cXl3dFBwd1dOYHFxYVtRZGxYeEpfaHdnWkxfckh4a2lHYVpXcVxceHFNSHdgeHZOQWhJX3Bqd3NMQHBNV3J4Xm9ocWlaP29oeXdeeHg8T3RRV1tjcWxyRmNab2FcXE5dU1diS3dxP2h5a3dqPXldYF9mZU5zZl9aW0BpXnh2eFhgaXZed15mWFdkUXZgdD9icVl2ZlZecXBoa2dubV9oT0lsZm5cXGZ5bm1weFA/eWZ2cHVpZXBuaHlmcnV5aT8+cVlvXnJ4dXhBdTtoeGVpYV93eVlmW0l2XWBIaXE/bkhXeER5XklQWnlobT9RYVZvaldhcFBnbW5JX3hgbl95aW1AYGp5cj9PeXlXXFxgeHFBQHVYTmFIV1xcP3BkQXlrdXZlSUBjO0FxbUl5Q3dcXGVBZj9Pb1dhcFBHeElIcHh5alJmanlha05gZ1VGeGFPY1N2eVlBbD9xZWZudTtmdklgeHRvd0hIcnBXaF5pY1lPaXRxeTxwZHlPbVM+eWluYFxceWE9aG50dm9bZ3RtcHlZVmpreWlPP2J0RmZ1YXhlX2V0eWhsWFp5eVo/eHBBeHlcXD5rWUZ2O1BbdG9gYkl2bEBuX2lmTUlhXVBgSXl5YFBjd154anlxRW5rc19iOm9qZElpPGZ4cFZ5cVF0S1lyeXlhYkhzS3ZcXEU/bFlndmpvdXU+ZHdeZXRIeWpYYHM/eVRxWlVZdVJeXlhgWj95a3hvaXQ/cW4/c0xgXztobGo/XV54ZHNoZ3N5cW1Bb053W3deeXlOW0J2XzxBZ09mdGloYFNJYT1oX2ZwYnhAdUlBZEl2bEhWXmY/YmRZZkBod1NfXzp2aF5Jd19Jb3lQb1ZZYkl2XFw9Vl9KPnA7Rm1oWWVKPnhhbl1iWW94b2l0Q1ZbYkJLZXFxcm9lR0JDST1NSGxReWNJd1tRdjs/Y0FheEo9UnhheV9FZlhLWXJ5eTpveHZjZHI9VFZZQ0N1Q3c/WDpJWDtDSXJJaUFvRWhTRXRpV2txRXQ9dz90S3dcXHg8dmpYbml1XnlBdl1BWWNOaWVkUGdqRDo7al5QTmFMTlFFTmpENUI=</Image></Text-field></Input></Group><Group><Input><Text-field layout="Title" style="Title">Analysis of a Dual-Input, Multi-Loop Servo</Text-field><Text-field layout="Author"><Font background="[0,0,0]" encoding="ISO8859-1" family="Times New Roman">\251 2004 J.M. Redwood</Font></Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">Introduction</Text-field></Title><Group><Input><Text-field layout="Normal" style="Text">This worksheet demonstrates the use of Maple to develop transfer functions for a naval gunmounting.</Text-field><Text-field layout="Normal" style="Text">The following Maple techniques are highlighted:</Text-field><Text-field firstindent="0.0" layout="Bullet Item" leftmargin="0.0" linebreak="space" rightmargin="0.0" style="Bullet Item"><Font executable="false" subscript="false" superscript="false">Developing the model</Font></Text-field><Text-field firstindent="0.0" layout="Bullet Item" leftmargin="0.0" linebreak="space" rightmargin="0.0" style="Bullet Item"><Font style="Text">Plotting results from the transfer functions</Font></Text-field><Text-field firstindent="0.0" layout="Bullet Item" leftmargin="0.0" linebreak="space" rightmargin="0.0" style="Bullet Item"><Font style="Text">Creating a Nyquist diagram</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">Part 1</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Synopsis</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">Transfer functions are developed in this worksheet for a dual-input, multi-loop servo representing the training motion of a naval gunmounting and are analyzed by classical methods using Bode attenuation and phase diagrams.  The closed loop responses are analyzed in Part II.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Introduction</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The purpose of the servo analyzed in this worksheet was to train a naval gunmounting so that its twin guns were aimed continuously at the target as directed by the ship's gun fire control system.  (Training is the angular motion in a plane parallel to the ship's deck.)  When a gun fires, the propellant gases acting on the gun's breech produce a force about the gunmounting's training axis, so causing a load torque tending to throw the guns off aim.  Thus, the training servo had two inputs - a training signal from the gun fire control system, and a load torque.  The training signal was transmitted by a synchro transmitter in the gun fire control system and the gunmounting's actual training was measured by a synchro control transformer located on the gunmounting.  The difference between the two provided the error signal input to the servo.</Text-field><Text-field layout="Normal" style="Normal">  The servo had three feedback loops - the gunmounting's angular position, angular velocity and angular acceleration.  The block diagram shown below illustrates the inputs, loops and main components of the servo.</Text-field><Text-field layout="Normal" style="Normal">  The ship is assumed to be in flat calm water, so that it does not suffer from sea-induced linear, or angular, motions.  (This avoids the complications that arise when such motions have to be taken into account.)</Text-field><Text-field layout="Normal" style="Normal">  The details of the servo have been adapted from those given in the reference, which are believed to have been used on a gunmounting similar to the one discussed in this worksheet.</Text-field><Text-field layout="Normal" style="Normal">  The servo's transfer functions are derived in this worksheet, and the open loop performance is analyzed using Bode attenuation and phase diagrams.  Closed loop performance is analyzed separately in Part II.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field bookmark="Block Diagram" layout="Heading 2" style="Heading 2">Block Diagram</Text-field></Title><Group><Input><Text-field alignment="centred"><Image height="346" width="573">TUZOV3RLVWI8b2I8Uj1NRExDZE5OWndqOmRvOkg6eFZddHFYQjpHQk5TRE9FPHl5eXlZY214YXVgQG5fQ19iO19kck91Z2FtO3dnVWZiTj5qWVF1aFdwR0BcXFhfbkJfblFOcHdScmdkRD9TS2VES2dibFVSZE9GSD1jTUdkcmN4Tm91PFVjakFVREtGUVtSeT1CZGtpTkNYOj1UcmN2SltmZz9kRHNiRD15OklGQT9HV2FUc1t4SmNGVTtma3FmU3V3aldIbW9jZ2l4VHFVSXFTS1dkeEdJZnFXQHF4cXdXWGFpc1VIaWlzaF1IREtTRElWSFtET1tCWztDOz1DPD1EPD9Ebk1uRGRxeUx0XFxYeHZobUtcXFRgTHZmdVZiTU1vZExjQVNlZGtUdW9LTVhkZUtMcGtYQUxrZUw7UEpPTUx3RFRfRXU/dVRfQXdJVHNwWW9pVFFTRFJoSHJIWEx5RHdmYVJCYFlxWWtVPHR5PHlDYXNadVJdYXlSdU5WSU1nXFxmTEFjVUhdTl5ub190QlZvZWZGeVZoPXZYWUlWX3FnRVd3bGpPVU5JYVR2cHl0PXdyaFdlSXhePU9WbVNoPHhlWEtzZXJnbVRtPHdCeUxZUGtwQVFTPFJWYWtaVVN2SVJsVHFLWG9tTWtAQHN5PHFCQFJKTGpYeHJfSVhfUG5iRVlGRVRTQUs8PUw8Pl9PR2hXZmJaZnVTZnBgaFxcZmdiWkhic2ZqblFcXGNQZ1RwalBoa1peWls+O19WUjtmQE1GdFN1W1l0ZUN0PkVIX2FYckV3SD9mYENZXTt2YnFUQTtXeEtkeUtlTVVEbW15Q1lFO0Nybk13YGdJbGtCZUFXc2NDYmlkeEdpTHVzWj1mSW9SSHFjPm9XcUlTYD1UOllkYmVFY21ScT11U1N1Wz1mQUt2VUdVXkt0czt3a01FX2V4eWdjYFVoWVl3PHFmbD9HYl1jQkNScHlUdHN3Q3VVTEdUcGNVQE92THV2anF0cGFWOkVyVVF2XFx3Z1A7aUZZQjpFeXVheVxcP3NJV3U/cVhYeVRRT0Vtb2ZyYVhjZ1RzRWlrU1VHbWhmVWlfZ3dDd2VTSXRJX1U6eUJ5b3R5b2dvX3U7P1d3PVVzVWdQU3dFdVM6X0Naa1ZePVNOV1JvU3ZpbXVRd1JAd3lAc0hNa1hXa1RVXXVxQXd4YVhNO1ZwX3JbPWZjPWZvP3dlZVhoeVdsYWY7U0RbSVltS3lYYVJuX0lgaUdUYWZnU3NpbVdAW2VLO2Zqd3ZuPURTd3ZCUWVGaVk8cXI+d1JXVXJPU3Y9P1lQXVZrbVVaUUVIXVdHYUVEQ1Q6S2VdZXl3bVNyd1NYd0NZUWk6U1RVZURiU2hVc2d2PUVfYXlFd3RzSUdmX3haZVZPcURPd3JfTWNhYXl5aVlbO1lUa3VmT1dUPWh4UUVtdUJgZ1Rnb3VyU1NCRXdVZ1dqXUdVUVl0bUhmZVNHUWJUaVJCU0dpdUlEY1JwR2l4P3VgPVJSeUg/S1VzYWRvY1c9W2hKZXdOSUlQZ2hbdURDcUdlXWdpS2lOZ1VlSXY/UUY6VUhxcUdEVVhtO3ZCXUNLc0RuV2NLPWhLV2N0SUlJUVlOW1hJXXR0TXJ1O2NoY3R1W2lBW0hHdWJ3UXRcXGl5O0lFVUViSkl2clN3U01pdldXdltYXFxvVVhTY3ZNRFdtdW5rSFc/ZVdDd1dtdV1Rd0ZVeV5TSG1pZXBVaE5TaFpBVmB1UmVTdT1TYlhvY15zWXh5Uzp3dV15SVplZHlnZVpndEJddUNdSGZDU0pRVFZFRHBrR3dBRG5NaDtld113aWN5Yl1ndmtNVkh5aWFxUjpdaElVV0FxQlF3VVJLUlJreF9dWTp1RDs9QzxcXEp3eXFIeVFHRXdzQWxeQXBacG1JRFNkcHRLRFQ/TFhcXElVTURTSnV0WGVSWlRLbFxcTENddFdkVEVNajxkVUtQcEVEWEpxTEtdVlNxdHU9UXR5UjtcXFh4WXg9aVJhVE5CeFZsbHd2bXBhbVBXQGtpREt0WHNIdEp4eHBKYUpnQXM7QFlESEpQYXFQQU1ZZWxqQWxBYHNySVk8YU5qbXZpSW1LUW12WVZYdFd0WXI6UHhOaFFdRE10PVJBaWpcXHlPS01zVFFwSWhQb1xcU01ReHZkVD9dTXJJcFNVVE9gVHB0U1pFbXFJbktNV19QUF5AUnZEal9ddTpIcGRVVWRYU0BcXHJ2QU9PaVlNcXB1WHFkPVdFcVF0cGxWSFNYbFg8YGprWHVqZFJIRXVieXheeFRdVE5vcVhAWVBPSXRzTG5iZXJQdXRrTVFUXFxUbmVwRVl2WEhxS0hWO11MdWVvTXVNaW1RV3lYZFBLOk11dVlxSXR2V2BvS0lsUUx1R3FrYXFZdXFPWWV4YEF2WmlSQ0VrUFlyYGhPd3VzcHlMUlFLPjxuW1R1bmxwb0V2dmhWWUVwUnlNVVRqcWVTdV1QZ0hYQFxccGw8dWV1VGB0WUJQbDtYVF9ET0pBVUpFTj5Qal9oWTx5eERZb188TU5pWXd0VXZEVkdNcFRxTndRTktgdEtUbWFpc2JVc0FIc1d4bE9teG49dWNxSkZZS3NsclV4T1NBYz9wZm5hdmlIc1ZAdUpvc18/YEJBcE52dmBxY2VBd2lhX3hualNXWkV5XFxQUG5zZ3dgR11PaG1fX2Q8SW94V2FdYXhheFtLXmpaUmQ/UmF5Z1lrR1RLQ2xVRmp1dmxxZXM/ZUtNRlldZ3NnaVBZVUthUkhbcm1nYmVtV0V5WHBbRjtdSVU9d0g9SVFLZW13WFdNUkFFdkxrZV9VQ0VBdl5BcmJHUnVZZk9Vcl93Q0FBR0lxRExrZEZhaVVhWVxcaUh0XWhkPXdgR3g/eWlTcVRbU2VmV3ZLU2lLW2hwWVZuQ1NKc1lqa1NxZ3NvV1ZnQVNwUUVNPXdLaWhBS1U9cXdwV0ZTX2dZW1lzYVV0R0RcXE9DXkFHPENUT01ES1VWTmFkd2F3cEt3ak1iRkFIaUdWdEFEOz1DPE5iSk5qSl50WV5sTnhaaXFpUGBfPz9kbVhsSndoSmlcXFBvbjpmXnd5d2Zoa3JnYnFfY1VWaXV5WlpOaVRpXk9HZ0lIZ0VpZUNRZ0VxY0JmYkVHcERGXFxNcXhpT2VvRm5oR15jXm5WXlpMWXRmWFxcTWFyblhcXENmaFdXXUtHZ1xcT3U8T19dWF5yPmdLXmpaUkJLWmR5ZHhZWmFodUJxXnNPYlZnOzxKOmBuXFx0TlxcdFRbPFA7Mjs=</Image></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Filter" style="Hyperlink">Return to Filter</Hyperlink>      <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Transfer Function - Input" style="Hyperlink">Return to H1</Hyperlink>     <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#H2" style="Hyperlink">Return to H2</Hyperlink></Text-field></Input></Group><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Definitions</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  R</Font> is the reference input (angular position in training required by the gun fire control system).</Text-field><Text-field layout="Normal" style="Normal">  <Font family="SymbolPi">e</Font> is the error between the controlled output and the reference input.</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">L</Font> is an external disturbing torque applied to the motor shaft that drives the load (gunmounting).</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">C</Font> is the controlled output (angular position of the gunmounting in training).</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">C/R</Font> is the closed loop response of the controlled output to the reference input.</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">C/L</Font> is the closed loop response of the controlled output to an external disturbing torque applied to the motor shaft that drives the load.</Text-field><Text-field layout="Normal" style="Normal">  </Text-field><Text-field layout="Normal" style="Normal">N.B.  The servo is assumed to be linear throughout the operating region and the principle of super-position therefore applies.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Transfer Functions shown in the Block Diagram</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The "Bode2" package will be required for plotting the Bode attenuation and phase margin diagrams.  It is therefore necessary to copy bode2.lib, bode2.ind and bode2.hdb to Maple's lib directory, for example, "C:\Program Files\Maple 9.5\LIB".  It is also necessary to call the "plots" package and to change the imaginary unit to <Font bold="true" italic="true">j</Font>, instead of <Font italic="true">I</Font>.  A reduction to 4 decimal places shown makes the display easier to read, so "displayprecision" is set to 4, (computation, however, is at the precision set by "Digits").</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Maple Input">restart; 
interface(imaginaryunit=j,displayprecision=4): 
with(plots):
with(Bode2);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">a.  <Font underline="true">Pre-amplifier</Font> (including gains of synchros and error-phase-discriminating circuits).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K1;</Font><Font italic="true"> </Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  K1</Font> has units of volts/radian.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">b.  <Font underline="true">Power Amplifier and Amplidyne</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K2G2 := K2/((Tf*s + 1)*(Tq*s + 1));</Font><Font italic="true"> </Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  K2</Font> is gain of the power amplifier and amplidyne generator (volt/volt).</Text-field><Text-field layout="Normal" style="Normal"><Font italic="true">  Tf</Font> and <Font italic="true">Tq</Font> are the time constants (sec) of the amplidyne generator field and armature quadrature currents respectively.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">c.  <Font underline="true">Combined Motor and Load</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K3G3 := K3/(s*(Tm*s + 1));</Font>  </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"><Font italic="true">  K3</Font> is the combined motor and load velocity constant (radian/sec/volt).</Text-field><Text-field layout="Normal" style="Normal"><Font italic="true">  Tm</Font> is the motor time constant = the combined motor and load inertias divided by the combined friction and armature current losses (sec).</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">d.  <Font underline="true">Tacho-generator</Font> (gunmounting's angular velocity feedback).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K4G4 := K4*s;</Font> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  K4</Font> is the tacho-generator constant (volt/radian/sec).</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">e.  <Font underline="true">Accelerometer</Font> (gunmounting's angular acceleration feedback).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K5G5 := K5*s^2;</Font> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">  <Font italic="true">K5</Font> is the accelerometer constant (volt/radian/sec squared).</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field bookmark="Filter" layout="Normal" style="Normal">f.  <Font underline="true">Filter</Font>.   <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Block Diagram" style="Hyperlink">Block Diagram</Hyperlink></Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K6G6 := K6*s^2*(Tc*s + 1)/((Ta*s + 1)*(Tb*s + 1)*(Td*s + 1));</Font> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  K6</Font> is the filter constant (sec squared).</Text-field><Text-field layout="Normal" style="Normal"><Font italic="true">  Ta</Font>, <Font italic="true">Tb</Font>, <Font italic="true">Tc</Font> and <Font italic="true">Td</Font> are time constants (sec) of the filter network.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">g.  <Font underline="true">Applied load torque to motor voltage</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K7;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  <Font italic="true">K7</Font> is the reciprocal of the DC motor torque-per-volt constant, and its units are volt/Nm.  It transforms the external load torque applied at the motor shaft into the equivalent voltage that would be required at the terminals of the motor's armature to produce that torque.  (Note that it is voltage that is summed at the input to the DC motor.)</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">h.  <Font italic="true" underline="true">G7</Font>.</Text-field><Text-field layout="Normal" style="Normal">  Let <Font italic="true">G7</Font> be the transfer function of the forward loop in the block diagram from the input to the power amplifier and amplidyne to the output.  <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Block Diagram" style="Hyperlink">Block Diagram</Hyperlink></Text-field><Text-field layout="Normal" style="Normal"> <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Observations" style="Hyperlink">Return to Observations</Hyperlink>  </Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">G7 := K2G2*K3G3;</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">i. <Font italic="true" underline="true">K54G54</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Let <Font italic="true">K54G54</Font> be the ratio of accelerometer transfer function to that of the tacho-generator (sec).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K54G54 := K5G5/K4G4;</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field bookmark="Transfer Function - Input" layout="Normal" style="Normal">j.  <Font italic="true" underline="true">H1</Font></Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Let <Font italic="true">H1</Font> be the transfer function of the inner feedback loop (comprising the tacho-generator, accelerometer and filter).  <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Block Diagram" style="Hyperlink">Block Diagram</Hyperlink>  </Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Maple Input">H1 := K4G4*K6G6*(1 + K54G54);</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Component Parameters</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The servo components had the following parameters (expressed in the SI units mentioned above). </Text-field><Text-field layout="Normal" style="Normal"> </Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">Params := {K1 = 14.8, K2 = 2300, K3 = 0.526, K4 = 0.98, K5 = 0.01368, K6 = 0.193,
 Tf = 1/31.5, Tq = 1/17.6, Tm = 1/9.1, Ta = 1/0.986, Tb = 1/5.26, Tc = 0.016, Td = 0.0016,
 K7 = 0.63062};</Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Loop Transfer Functions</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Transfer Function between Reference Input &amp; Controlled Output</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  When determining the transfer function from the reference input, <Font italic="true">R</Font>, to the controlled output, <Font italic="true">C</Font>, the externally applied torque is set to zero.</Text-field><Text-field layout="Normal" style="Normal">  Let the loop transfer function from the pre-amplifier output, to the controlled output, <Font italic="true">C</Font>, be <Font italic="true">G8</Font>.  Then</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">G8 := G7/(1 + G7*H1);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Define <Font italic="true">K1G8</Font> as the loop transfer function from the pre-amplifier input, <Font family="SymbolPi">e</Font>, to the controlled output, <Font italic="true">C</Font>.  Thus</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K1G8 := K1*G8;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Hence, the closed loop transfer function from the reference input, <Font italic="true">R</Font>, to the controlled output, <Font italic="true">C</Font>, is</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/R` := simplify(K1*G8/(1+K1*G8)):</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Transfer Function between Disturbing Torque and Controlled Output</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  When determining the transfer function from the externally applied disturbing torque, <Font italic="true">L</Font>, to the controlled output, <Font italic="true">C</Font>, the reference input <Font italic="true">R</Font> is set to zero.</Text-field><Text-field layout="Normal" style="Normal">The loop is then opened at the output from <Font italic="true">K2 G2</Font>, the power amplifier and amplidyne.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field bookmark="H2" layout="Normal" style="Normal">   Let <Font italic="true">H2</Font> be the transfer function of <Font italic="true">K1</Font> in parallel with the inner feedback loop, <Font italic="true">H1</Font>, determined above.    <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Block Diagram" style="Hyperlink">Block Diagram</Hyperlink></Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">H2 := K1+H1;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"> Let <Font italic="true">H3 </Font>be the transfer function of the feedback loop, thus</Text-field><Text-field bookmark="H3" layout="Normal" style="Normal">     <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Observations" style="Hyperlink">Return to Observations</Hyperlink> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">H3 := K2G2*H2;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"> Define <Font italic="true">G9</Font> as the open loop transfer function from the external torque, <Font italic="true">L</Font>, to the controlled output, <Font italic="true">C</Font>, thus</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">G9 := K7*K3G3;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"> Hence, the closed loop transfer function between the externally applied torque, <Font italic="true">L</Font>, and the controlled output, <Font italic="true">C</Font>, is</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/L` := G9/(1+K3G3*H3);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The dimensions of <Font italic="true">C/L</Font> = <Equation input-equation="[K7]*[K[3]*G[3]]/(1+[K[3]*G[3]]*[K2*G[2]]*[H[2]]);" style="2D Comment">NiMqKDcjJSNLN0ciIiI3IyomJiUiS0c2IyIiJEYmJiUiR0dGK0YmRiYsJkYmRiYqKEYnRiY3IyomJSNLMkdGJiZGLjYjIiIjRiZGJjcjJiUiSEdGNUYmRiYhIiI=</Equation> are <Equation input-equation="volt/(N*m);" style="2D Comment">NiMqJiUldm9sdEciIiIqJiUiTkdGJSUibUdGJSEiIg==</Equation><Font encoding="ISO8859-1"> \327 </Font><Equation input-equation="radian/volt" style="2D Comment">NiMqJiUncmFkaWFuRyIiIiUldm9sdEchIiI=</Equation>.  This equation is normalized by dividing both sides by 1 <Equation input-equation="radian/(N*m);" style="2D Comment">NiMqJiUncmFkaWFuRyIiIiomJSJOR0YlJSJtR0YlISIi</Equation>. </Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Bode Diagrams</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Attenuation Diagram for Inner Loop</Text-field></Title><Group><Input><Text-field bookmark="G7H1" layout="Normal" style="Normal">   <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Observations" style="Hyperlink">Return to Observations</Hyperlink></Text-field><Text-field layout="Normal" style="Normal">  The transfer function, <Font italic="true">G7H1</Font>, for the inner loop from the power amplifier's input is</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">G7H1 := G7*H1;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Inserting the parameter values in <Font italic="true">G7H1</Font> to yield <Font italic="true">G7H1_atten</Font></Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">G7H1_atten := eval(G7H1,Params);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  <Font italic="true">G7H1_atten</Font> is now plotted on the Bode attenuation diagram together with the straight-line asymptotic approximations thus,</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">k := 228.8215720:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">num := [.016,.01395918367]:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">den := [.03174603175,.05681818182,.1098901099,1.014198783,.1901140684,.0016]:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">a := linelog(k,2,num,den,.01,5000,top):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">b := attenlog(G7H1_atten,.01,5000):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(a,b);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Since the normalizing time constant, <Font italic="true">T</Font> = 1 second, the values on the abscissa may be read as frequencies in 1/sec (= radians/second).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The gain of <Font italic="true">G7H1</Font> increases at 12 dB/octave from very low frequencies, breaks twice before reaching a maximum of 57 dB at a frequency of <Font family="SymbolPi">w</Font> = 6/sec.  The gain then declines at an average rate of about 11.5 dB/octave becoming 12 dB/octave when the frequency is greater than <Font family="SymbolPi">w</Font> = 625/sec.  The break points are located at <Font family="SymbolPi">w</Font> = 0.986/sec, 5.26/sec, 9.1/sec, 17.6/sec, 31.5/sec, 62.5/sec, 71.6/sec and 625/sec.  These correspond to the reciprocals of the time constants <Font italic="true">Ta</Font>, <Font italic="true">Tb</Font>, <Font italic="true">Tm</Font>, <Font italic="true">Tq</Font>, <Font italic="true">Tf</Font>, <Font italic="true">Tc</Font>, <Font italic="true">K54</Font>, and <Font italic="true">Td</Font>, where <Font italic="true">K54</Font> is <Font italic="true">K5/K4</Font>.</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">Ta</Font> and <Font italic="true">Tb</Font> are used to reduce the gain rate at frequencies between 1/sec and 10/sec, while <Font italic="true">Tc</Font> and <Font italic="true">K54</Font> are used to reduce the attenuation rate at frequencies of 62.5/sec and 71.6/sec until 625/sec, when the attenuation rate again increases to 12 dB/octave.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Since <Equation input-equation="K1G8 = K1/(1+K1G7H1);" style="2D Comment">NiMvJSVLMUc4RyomJSNLMUciIiIsJkYnRiclJ0sxRzdIMUdGJyEiIg==</Equation>, this is approximately <Equation input-equation="1/G7H1;" style="2D Comment">NiMqJiIiIkYkJSVHN0gxRyEiIg==</Equation> when <Font italic="true">K1G7H1</Font> is large.  Thus the shape of the main forward loop, <Font italic="true">K1G8</Font>,<Font italic="true"> </Font>is the same as that of the reciprocal of <Font italic="true">G7H1</Font>.  This is very convenient, since the latter is easier to use when choosing time constants to produce good performance.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Loop Attenuation Diagram for Reference Input</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  Inserting the parameters in <Font italic="true">K1G8</Font>, the open loop gain for the reference input can be plotted.  However, to plot the straight-line approximations, <Font italic="true">K1G8</Font> must be factored as follows.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K1G8_atten := eval(K1G8,Params):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K1G8_atten := factor(K1G8_atten);
</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">  If a line segment of the plot starts and ends at the same frequency, <Font family="SymbolPi">w</Font> = 1/<Font italic="true">T</Font>, the "linelog" function in the Bode package produces an error message "plot ranges must be non-empty".  That would occur at each of the double break points, unless they are separated very slightly as follows: </Text-field><Text-field layout="Normal" style="Normal">  a.  "Digits" is increased, and</Text-field><Text-field layout="Normal" style="Normal">  b.  The second of each approximate factor of the quadratic's approximate factors in the denominator is decreased very slightly from that of the first.  That is the second approximate <Font family="SymbolPi">w</Font> = 1/<Font italic="true">T</Font> for each quadratic is reduced slightly.</Text-field><Text-field layout="Normal" style="Normal">  As a result, "linelog" plots an extremely short line segment between the two very slightly separated frequencies.  To the eye, this appears as a double break point. </Text-field><Text-field layout="Normal" style="Normal"> "Digits" is restored to 10 after plotting the Bode attenuation diagram.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">Digits := 13:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">k := 90331642.97:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">num := [625.0000001, 5.259999996, .9860000007]:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">den := [79.07244591, 59.53034908, sqrt(796354.8919), sqrt(796354.89189),
         sqrt(.4362516033e-2), sqrt(.43625160329e-2)]:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">a := linelog2(k,1,num,den,.001,5000):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">b := attenlog(K1G8_atten,.001,5000):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(a,b);</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">Digits := 10:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Using the "apparent gain" of 90331642.97, "linelog2" calculated the true gain as 85 dB (17905.04) by multiplying the apparent gain of <Font italic="true">K1G8_atten</Font> by all its numerator time constants and dividing it by all its denominator time constants.  This is the gain of <Font italic="true">K1G8</Font> when expressed in canonical form, and matches the gain indicated by extrapolating the low frequency asymptote of <Font italic="true">K1G8</Font> to the gain axis.  It is known as the velocity error coefficient, <Font italic="true">Kv</Font>, and has units of 1/sec.  <Font italic="true">Kv</Font> determines the closed loop system's error when following a constant velocity input, which is an important operational requirement for a gunmounting.  If the constant velocity input is <Font family="SymbolPi">W</Font> radian/sec, the error after all the transients have decayed, is <Equation input-equation="Omega/K[v] = Omega/17905.04;" style="2D Comment">NiMvKiYlJk9tZWdhRyIiIiYlIktHNiMlInZHISIiKiZGJUYmLSUmRmxvYXRHNiQiKC8weiIhIiNGKw==</Equation> radians = .5585e-4 <Font family="SymbolPi">W</Font> radians.  This will be examined further in Part II. </Text-field><Text-field layout="Normal" style="Normal">Since the normalizing time constant, <Font italic="true">T</Font>, is 1 second, the values on the abscissa may be read as frequencies in 1/sec (= radians/second).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The break points of <Font italic="true">K1G8_atten</Font> occur at <Font family="SymbolPi">w</Font> = 0.986/sec, 5.26/sec, 59/sec, 79/sec, 625/sec.</Text-field><Text-field layout="Normal" style="Normal">  The double break points at <Equation input-equation="omega = sqrt(.4362516033e-2);" style="2D Comment">NiMvJSZvbWVnYUctJSVzcXJ0RzYjLSUmRmxvYXRHNiQiK0xnXmlWISM3</Equation> = 0.066/sec and <Equation input-equation="sqrt(796354.8919);" style="2D Comment">NiMtJSVzcXJ0RzYjLSUmRmxvYXRHNiQiKz4qW04neiEiJQ==</Equation>= 892/sec.  Break point frequencies that can be adjusted are 1/<Font italic="true">Ta</Font> = 0.986/sec and 1/<Font italic="true">Tb</Font> = 5.26/sec; the others are affected by the values of 1/<Font italic="true">Tc</Font> = 62.5/sec, 1/<Font italic="true">K54</Font> = 71.6/sec and <Font italic="true">K6</Font> = 0.193/<Equation input-equation="sec^2;" style="2D Comment">NiMqJCUkc2VjRyIiIw==</Equation>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The very high gain at low frequencies is reduced rapidly by the -18 dB/octave segment between <Font family="SymbolPi">w</Font> = 0.066/sec and 0.98/sec, and by the -12 dB/octave segment that follows it from <Font family="SymbolPi">w</Font> = 1/<Font italic="true">Ta</Font> = 0.98/sec to <Font family="SymbolPi">w</Font> = 1/<Font italic="true">Tb</Font> = 5.26/sec.  Together with the overall gain setting of 85 dB, <Font italic="true">Ta</Font> and <Font italic="true">Tb</Font> produce gain crossover over (also known as feedback cutoff) about 1/3 down the 3.5 dB long, -6 dB/octave segment from <Font family="SymbolPi">w</Font> = 5.26/sec to 59/sec.  This indicates a very stable system that should follow the control signal well, while rejecting much of the noise in it.  After a -18 dB/octave section from <Font family="SymbolPi">w</Font> = 59/sec to 625/sec, the loop gain then declines at -24 dB/octave, ensuring that unwanted high frequency signals are well attenuated.  There is a second very short, -18 dB/octave section between <Font family="SymbolPi">w</Font> = 625/sec and 892/sec, before the loop gain decreases at the final rate of -24 dB/octave.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  It will be seen that at frequencies less than 59/sec, the attenuation is essentially the same as the reciprocal of the inner loop's transfer function, <Font italic="true">G7H1</Font>, plotted previously.  As noted above this can be helpful when choosing the time constants <Font italic="true">Ta</Font>,<Font italic="true"> Tb</Font> and <Font italic="true">Tc</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The frequency at which gain crossover over occurs is obtained by substituting <Font italic="true">j</Font><Font family="SymbolPi">a</Font> = <Font italic="true">s</Font>, where <Font family="SymbolPi">a</Font> = <Font family="SymbolPi">w</Font><Font italic="true">T</Font>, setting <Font italic="true">K1G8</Font> = 1, and solving for <Font family="SymbolPi">a.</Font></Text-field><Text-field layout="Normal" style="Normal">  </Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">omega*T = abs(fsolve(subs(s=j*alpha,K1G8_atten)=1,alpha));</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Alternatively, the cursor may be used to read the frequency at gain crossover over, thus log10(<Font family="SymbolPi">w</Font><Font italic="true">T</Font>) = 1.16, whence <Font family="SymbolPi">w</Font><Font italic="true">T</Font> = 14.4, which is sufficiently close for all practical purposes.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Loop Phase Margin Diagram for Reference Input</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The loop phase margin is shown by</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">phaselog(K1G8_atten,.001,5000);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The sudden loss of phase margin seen at about log10(<Font family="SymbolPi">w</Font><Font italic="true">T</Font>) = -1.4 occurs when the open loop attenuation changes slope from -6dB/octave to -18 dB/octave.  At gain crossover, log10(<Font family="SymbolPi">w</Font><Font italic="true">T</Font>) = 1.17, the phase margin is 43 degrees.  Although, this indicates a very stable system, it appears to be only conditionally stable. If the gain were to reduce, and gain crossover-over occurred at about log10(<Font family="SymbolPi">w</Font><Font italic="true">T</Font>) = 0.4, i.e. <Font family="SymbolPi">w</Font> = 2.4/sec, the phase margin would fall to zero.  For this to happen, the loop gain would need to reduce by about 28 dB i.e. by a factor of 25.  This is unlikely since it is greater than the pre-amplifier gain (<Font italic="true">K1</Font> = 14.8) and because the power amplifier and amplidyne gains are most unlikely to halve at the same time as the pre-amplifier loses all of its gain.</Text-field><Text-field layout="Normal" style="Normal">  Nevertheless, a reduction of phase margin to, say, 12 degrees caused by the pre-amplifier's gain halving, would cause sluggish, oscillatory performance.  This is illustrated in Part II, when the transient responses are examined.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Loop Attenuation Diagram for Disturbing Torque Input</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  When considering the feedback loop from the (load) output to the motor, the reference input is set to zero.  Before examining this loop, it is worth looking briefly at <Font italic="true">K3G3_atten</Font>, the motor and load transfer function, which is shown below.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K3G3_atten := eval(K3G3,Params);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">a := linelog(0.526,1,[ ], [.1098901099],.01,100):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">b := attenlog(K3G3_atten,.01,100):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(a,b);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  <Font italic="true">K3G3_atten</Font> demonstrates the usual characteristics of a motor and load - small gain at low frequencies, and attenuation at 12 dB/octave at frequencies above the single break-point at <Font family="SymbolPi">w</Font> = 1/<Font italic="true">Tm</Font> = 9.1/sec (again, the normalizing time constant, <Font italic="true">T</Font> = 1 second).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The transfer function for the complete open loop feedback from the (load) output to the motor is <Font italic="true">H3</Font>,</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">H3;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Inserting the component parameters, this becomes</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">H3_atten := factor(eval(H3,Params));</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The attenuation diagram for the open loop feedback transfer function from the (load) output is plotted below.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">Digits := 13:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">k := 174605.1797:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">n := [91.70572349,15.84188672,.9721856834,sqrt(248.0537289),sqrt(248.05372889)]:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">d := [624.9999999,31.49999999,17.60000001,5.260000000,.9859999999]:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">a := linelog2(k,0,n,d,0.01,5000):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">b := attenlog(H3_atten,0.01,5000):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(a,b);</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">Digits := 10:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The very large low frequency gain of <Font italic="true">H3_atten</Font> remains virtually constant and never dips below 81 dB, thus ensuring that even a small movement of the load shaft produces a very large opposing torque at all frequencies.  This is caused by the feedbacks from the tacho-generator and accelerometer offsetting the small low frequency gain and attenuation of <Font italic="true">K3G3</Font>, the motor and load transfer function.</Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Nyquist Diagram for Reference Input</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  Before using the Nyquist stability criterion for a multi-loop system such as this, it is necessary to determine the number of poles, <Font italic="true">P</Font>, of the open loop transfer function that lie in the right half of the <Font italic="true">s</Font> plane.  This is given by,</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">solve(K1G8,{s});</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Since <Font italic="true">Ta</Font>, <Font italic="true">Tb</Font> and <Font italic="true">Td</Font> are all positive real constants, <Font italic="true">K1G8</Font> has no poles in the right-half plane, i.e. <Font italic="true">P</Font> = 0.  Moreover, this result is independent of all gain settings.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  (Maple may not be able to find the roots of the open loop transfer function symbolically.  However, it should be able to find numerical solutions after the servo's parameters have been inserted in <Font italic="true">K1G8</Font>.  In such cases, it is necessary to check that gain changes do not affect the number of poles, or their location in the left, or right-half planes.)</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Nyquist's stability criterion derives from the theorem that <Font italic="true">N = Z - P</Font>, where</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"><Font italic="true">    Z</Font> is the number of zeros in the right-half plane of the denominator of the closed loop transfer function (1 + <Font italic="true">K1G8</Font>),</Text-field><Text-field layout="Normal" style="Normal">    <Font italic="true">P</Font> is the number of poles in the right half plane of the open loop transfer function (<Font italic="true">K1G8</Font>), and</Text-field><Text-field layout="Normal" style="Normal"><Font italic="true">    N</Font> is the number of net clockwise encirclements of the [-1, j0] point described by a representative point as it travels from <Font family="SymbolPi">w = </Font>-<Font encoding="ISO8859-1" family="SymbolPi">\245</Font> to <Font encoding="ISO8859-1" family="SymbolPi">+\245</Font> around the locus of the open loop transfer function in the complex plane.  (A negative value of <Font italic="true">N</Font> signifies the number of net anti-clockwise encirclements of the [-1, j0] point.)</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  For stability, <Font italic="true">Z</Font> must be identically zero.  Otherwise, the closed loop transfer function would contain poles in the right-half plane, and the controlled output, <Font italic="true">C</Font>, would increase without bound.  Thus, the system is stable if, and only if, <Font italic="true">N</Font> <Font encoding="ISO8859-1" family="SymbolPi">\272</Font> <Font italic="true">-P</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  In this case <Font italic="true">P</Font> = 0, so <Font italic="true">N</Font> must be zero for stability.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  It is now convenient to obtain a function <Font italic="true">f</Font> (<Font italic="true">j</Font><Font family="SymbolPi">w</Font>) for the loop transfer function in the complex plane thus,</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K1G8_jw := eval(K1G8,Params):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K1G8_jw := subs(s=j*omega,K1G8_jw):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">f := unapply(K1G8_jw,omega):</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The Nyquist diagram is plotted below with some representative points shown on the locus of <Font italic="true">K1G8_jw</Font>.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">a := complexplot([seq(evalc(f(i)),i={12,15,20,30,65,100,-12,-15,-20,-30,-65,-100})],x=-1..0,
style=point,symbol=circle,colour=black):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">b := complexplot(f,10..1000,colour=blue):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">c := complexplot(f,-10..-1000,colour=red,linestyle=3):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">d := plot([[-1,0]],style=POINT,symbol=cross,symbolsize=10,color=black):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">e := coordplot(polar,[0..1,Pi/2..3*Pi/2],grid=[4,13],color=[black,black],linestyle=[2,2]):
f := textplot([-1,-0.05,"[-1, j0]"],font=[TIMES,ROMAN,11],colour=black):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(a,b,c,d,e,f,scaling=constrained,title="Nyquist Diagram\nPositive frequencies (blue), Negative frequencies (red)",titlefont=[TIMES,ROMAN,12]);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The Nyquist diagram above shows the locus, in blue, of <Font italic="true">K1G8_jw</Font> between <Font family="SymbolPi">w</Font> = 11/sec and <Font encoding="ISO8859-1" family="SymbolPi">w \256 +\245</Font>, which is located at the origin of the diagram.  It also shows the locus, in red, for the portion between <Font family="SymbolPi">w</Font> = -11/sec and <Font encoding="ISO8859-1" family="SymbolPi">w \256</Font> -<Font encoding="ISO8859-1" family="SymbolPi">\245</Font>, also located at the origin.  Reading from the grid, the blue locus shows that the phase margin of <Font italic="true">K1G8_jw</Font> is slightly less than 225 - 180 = 45 degrees.  Using the cursor, it is</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`phase margin` = evalf(arctan(-0.68/(-0.73))*180/Pi*degrees);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">  The remainder of the Nyquist diagram can be inferred from the Bode diagrams.  There, it was seen that the phase margin becomes negative when log10(<Font family="SymbolPi">w</Font><Font italic="true">T</Font>) equals -1.18, 0.38 and 1.81 (i.e. <Font family="SymbolPi">w</Font> = 0.07/sec, 2.4/sec and 64.5/sec).  Hence, the loci of <Font italic="true">K1G8_jw</Font> must cross the negative real axis three times before becoming asymptotic to the imaginary axis.  The highest frequency of crossing has been plotted in the diagram above, while the other two are located far to the left on the negative real axis.  The complete diagram is illustrated in the sketch below, which shows zero encirclements of the Nyquist point.  (The area enclosed always lies to the <Font italic="true">right</Font> of a representative point as it traces the locus from <Font family="SymbolPi">w = </Font>-<Font encoding="ISO8859-1" family="SymbolPi">\245</Font> to <Font encoding="ISO8859-1" family="SymbolPi">+\245</Font>.  Since the Nyquist point is not enclosed within this area, there are zero net encirclements.)  However, if the scale of the axes is either increased, or decreased, while leaving the locus fixed, the locus would encircle the Nyquist point twice in the clockwise sense.  This would violate the stability criterion, and so the system would be unstable.  The change of scale corresponds to a change of the servo's gain, and indicates that the servo is conditionally stable depending upon the gain.  As already noted in the discussion about the Bode diagrams, the case of the gain being too small is unlikely to occur.</Text-field></Input></Group><Group><Input><Text-field alignment="centred"><Image height="298" width="291">TUZOV3RLVWI8b2I8Uj1NRExDZE5GWl1qOmRsOkg6O0JTPUljRz9HQk5TRE9FPHl5eXlZY214YXVgP0ZiQ19iO19kck91Z1FtRz9rVmZiTj5qZVBxdUZ1WGdlV0hibD9bYk9sUFlrOz9sRk9oTk9sZ19cXD5fYlNHeF5OampQX1hPbEV4WklOXztgZF4+Y0hPbE0/c0NvYktGXFx2PnFYUHY/R29vYXZHaXQ+XmVqVmt3UHFEQV9sSHdTd2BTcWhwb2ZEcXlZcG1oeWRPeWdsZnNmUWt1PnBGcVt3aXN1d29YZ3FxUXdPZ3RKR2JEeWVlUF07YXFwQW9rQG94bmpMTl1IV25FQWROaGtrWHdATmJdeV9tPl11cFtpT153aXJ3XlpvZmtUVnZdQVtHYFxcdXZrXFxPXT1HdGdwaGNmeHV4eVFXZ1FIZEleWktHZG1IZGtpeE9Oa0xJaGtQcTtQbmhpb2xWYXBOY2JmW0xHaUFQXXVYbl5nZ0FnY1NGdFVxeWRed0RPcnhoakxoaEpHY1NYbk9vXlVnYnNGZ1xcb2t5bmJFTl9yT2hTX2h2T2tHcHZDTmpNaHZbYGZTTmhEP25iXm1XZ3FmR3hOUGRAXmBHZlpkVmV1V3ddT2RvUVtRR2RASG4+Vlp5UG5LZmI7eXBfXmRRbmtVYV1veVxcRFBhYlZdXVZzYE9aPWh4QV9kXFxGcElBZkFebnlheUBecE1eaDpnbVhXaWNoeHh5blFueV5HcUJAXVxcYHJEcXREcXZHeGBxdndLb3NKZmF0eHdrVnh0V1xcamFpYG5bPVZxYGBxRUdsRT5fTWBsPXZtZFFbTGFeaD94W0B0YV5ea19mQlhzWGlwbXZtd19pdWloWW9zRmBrXFxIaGRnbXdAdlVoaHg+dGJJbWRxb3BAZlBPd3ZhalRIYj9Rd1VveXFRbUxRXlpoZnJJaT5fdU1xXFxgVmFiaVxcRmZ0XUFyXW9tYWBoVVZeZVlhYHZlOm9aTklgRm9ncVF5Uk53Q1dadF93R1ZkV25rU0huYUB2QXF2amZya1Zwd2BqY3hlVWhbRmZeR3BfS0hudD9vREhvUWhfTFBdPHF4UlBnTndwTUlkU3Z3XmBseEZuVz5tT1BkVmZeRXZia15bT3leXVdkV2FqbE5fVllrdE5vYXFaRFBjX0d1TVB1Q1drS2BmQ05gY1ZxYl9jS2l1cE94OlleZHhlW3BiZ1BbU29ob3FtYmBmRnBzSGdtR0ZnQlBacFhcXHBnZmw/b29gakpoXmhgYVRfWlZJclRwXFxPUGtpSW1YUWxDUXhtP2dkYWlqbmVrQHZxSXRHP2xMSWVqUHJGeXFKV3g/Z3FSXnhYXnB1XmVVVmpAPnBnWWdeV2VGYF1kWFxcS3Vod2dYU3U/cWd1eUNRUWRmaUlBQXZra0dnc2JhP1RLd0R0d0hnO2V0a0Q+P2RobUc+d2VpdXd5bUVqYUJfaVY/U1hQV2RPdWllZ3VCc2N5S3VtRVZQc1dibXJqdVRlcXRdb2dwUVVySVdoUVZycVhRdXQ7Z3d0d0luV2VZP2NkWWlGUVlYQUNVZ0J3S1l2UWNmZXdqeVReb1M9P3dZY2JgSVhTXVNSY1NrcWV5cVhsO0RwXVJYc1Vpc3lLO0lWb1lCPXNcXFdHRXdHQF14YlVTb1NXd1lYSENES1ljb2tjWm9jcnVYVGdFaGFmUT1leUFJX2F5SWtyQnlYRHFHb3V0S2NXTllHUEFzQHlkX2tGWGlJYF9kX09ZWHlTYk9JRG1kQ09mRltUeXVTS01obmVkbWlpS2N5X1VmYlNnZl9EQF9UV2NGRUFSS1V0cmd5XWdjREdDUE1kbkdEX01mXmVlYkNTbj9ZUW1zaktoblVDbW1rdkRqRFh5PW1ubDxuXU1rcEhUc1xcUE91VD9QV1R0VXk8bmxJeGplbVpsc209c2RQV1VpbUxlUEVYUz1AUE1UVVl4TUFQSkFRc191cj1Bb2RNSkhtVlthT2d1bV9kS0NgV3hsTT9BV154cl9ITktobG9kVkNxTnVUT0VoWWVAdGJhTnNJcmhZbEdtcElFVWNtTE9kWVtpdm9sWGZNbm1FcUB5bVhxeWtpbW9YdGVobXNNUmhMU2BIdGFwb1hQcXRNWGVlblFcXE5yWU5qbXlBXFx0bVxcWGhpdVVFbEF1d2A9eD5kamtEUnlwdz91UHBFeUZYU0lFcG9AcjthcU90a3NJV2NoeFlgcXBpcnF0Vmh0cWRYeVBZbD08eD5kbmFJdHl0UVJsV3NdeFtYeGJhVl1QbHVpbD9tc2lodXZdU2JJdT1gb3ZYd0FEbUBvdnVQbG9JeDo+d05hYV9IakhnZFhBcEZub0BOeE5WalZpZFJnb2Y/XlhYa05oWnNYY3NXXl9JY2NWXFxsd2xsZ1xcb19uTkhac1hvOkdnSlhxTk9dcWh2d29bP2B1T19uUj52QmZsQGlaP2ZyO3ZiYm5nSVdwUWhiVG5mQlhmWVleaVZvYmBoRHh1OkFnTl5dXUhoVUd1Q3BeO25wSFB3UXdxTk9kaz9qQ3ZaXVZnVUhmYmhob3ddZ0haYEdzXmdbVWheXXB5YD9dZFlaRndpaV5lRmZ4ZVZkTEdwRlhkZGldP19kRk91SEdzeFh2ZkZqSF5hYkF4W0dqSl5uXW9fSHF2WmB5aXd3aE9hVEZlZm5dcl9bcUF0bklxRU9mdG9wcmB2YGFgaFhxSWZsV1lrYFl4RU9hQE95Zk5rVG91V1hpcXFaRFBocHd3Q1lqc3hoSGFwQWdjUlB5UVZ2ZWBkPFl2PEdtUV06RmM/b2M+b288P2Y8Mzw=</Image></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The difficulty with the low frequency portion of the Nyquist diagram can be partly overcome by plotting the inverse diagram, i.e. plotting the reciprocal of <Font italic="true">K1G8_jw</Font>.  Unfortunately, such plots can be difficult to interpret and serve to illustrate the superiority of Bode diagrams in the low frequency region.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Observations</Text-field></Title><Group><Input><Text-field bookmark="Observations" layout="Normal" style="Normal">1.  The choice of compensation and time constants has not been elaborated in this worksheet.  It is, however, worth mentioning that acceleration feedback is sometimes used to increase the load's apparent inertia in a manner analogous to using velocity feedback to increase the load's apparent viscous friction.  In this case, the increased apparent inertia would decrease the throw-off caused by the guns firing.  Having included acceleration feedback, velocity feedback is the best choice for the servo's compensation.  The accelerometer and tacho-generator gains and the filter constants were chosen to cause rapid attenuation of the gain at frequencies above about 70/sec and so provide the high "apparent inertia" to sudden load torques.  They were also chosen to attenuate the large low frequency gain rapidly at frequencies greater than about 0.07/sec, and to produce a long 6 dB/octave attenuation at gain crossover.  These are important for achieving good, constant-velocity performance needed for aiming accurately at moving targets, and for good servo stability with noise rejection.  These choices are usually made with the help of Bode diagrams for the inner feedback loops involved, in this case, the <Font italic="true">G7H1</Font> loop and the <Font italic="true">H3</Font> loop.</Text-field><Text-field layout="Normal" style="Normal"><Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#G7H1" style="Hyperlink">G7H1</Hyperlink> <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#H3" style="Hyperlink">H3</Hyperlink>   </Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">2.  The modern design process includes the classical, open loop, frequency analysis shown in this worksheet together with the closed loop analysis shown in Part II.  However, specialized software for servos is normally used to try out different approaches and to optimize the parameter settings, taking account of non-linearities, backlash, etc.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">3.  Much painstaking work with slide-rule, log-tables, and semi-log charts were needed in the past to perform the calculations and draw the Bode diagrams shown in this worksheet.  Great care was necessary when developing the transfer functions - a minor slip with the algebra could cause immense problems.  Maple does it all quickly, easily and without errors.  Not only that, it can feed optimized code to certain software used for designing servos, thus reducing setting up times and errors.</Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">Part 2</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Synopsis</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">In Part I, transfer functions were developed and open loop responses were analyzed for a multi-input, multi-loop servo representing the training motion of a naval gunmounting.  This worksheet examines the closed loop responses of the servo to standard inputs and the responses to load torques caused by a gun firing.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Introduction</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The purpose of the servo analyzed in this worksheet was described in the introduction to Part I.  Briefly, it was to train a naval gunmounting so that its twin guns were aimed continuously at the target by the ship's gun fire control system, and to resist the load torque created when a gun fired.  Part I, developed the servo's transfer functions and analyzed the open loops.  The closed loop responses to standard signals at the reference input and to approximations for load torques are examined in this worksheet. </Text-field><Text-field layout="Normal" style="Normal">  The servo had three feedback loops - the gunmounting's angular position (i.e. training), angular velocity and angular acceleration.  The block diagram shown below illustrates the inputs, loops and main components of the servo.</Text-field><Text-field layout="Normal" style="Normal">  The ship is assumed to be in flat calm water, so that it does not suffer from sea-induced linear, or angular, motions.  (This avoids the complications that arise when such motions have to be taken into account.)</Text-field><Text-field layout="Normal" style="Normal">  The details of the servo have been adapted from those given in the reference, which are believed to have been used on a gunmounting similar to the one discussed in this worksheet.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field bookmark="Block Diagram" layout="Heading 2" style="Heading 2">Block Diagram</Text-field></Title><Group><Input><Text-field alignment="centred"><Image height="346" width="573">TUZOV3RLVWI8b2I8Uj1NRExDZE5OWndqOmRvOkg6eFZddHFYQjpHQk5TRE9FPHl5eXlZY214YXVgQG5fQ19iO19kck91Z2FtO3dnVWZiTj5qWVF1aFdwR0BcXFhfbkJfblFOcHdScmdkRD9TS2VES2dibFVSZE9GSD1jTUdkcmN4Tm91PFVjakFVREtGUVtSeT1CZGtpTkNYOj1UcmN2SltmZz9kRHNiRD15OklGQT9HV2FUc1t4SmNGVTtma3FmU3V3aldIbW9jZ2l4VHFVSXFTS1dkeEdJZnFXQHF4cXdXWGFpc1VIaWlzaF1IREtTRElWSFtET1tCWztDOz1DPD1EPD9Ebk1uRGRxeUx0XFxYeHZobUtcXFRgTHZmdVZiTU1vZExjQVNlZGtUdW9LTVhkZUtMcGtYQUxrZUw7UEpPTUx3RFRfRXU/dVRfQXdJVHNwWW9pVFFTRFJoSHJIWEx5RHdmYVJCYFlxWWtVPHR5PHlDYXNadVJdYXlSdU5WSU1nXFxmTEFjVUhdTl5ub190QlZvZWZGeVZoPXZYWUlWX3FnRVd3bGpPVU5JYVR2cHl0PXdyaFdlSXhePU9WbVNoPHhlWEtzZXJnbVRtPHdCeUxZUGtwQVFTPFJWYWtaVVN2SVJsVHFLWG9tTWtAQHN5PHFCQFJKTGpYeHJfSVhfUG5iRVlGRVRTQUs8PUw8Pl9PR2hXZmJaZnVTZnBgaFxcZmdiWkhic2ZqblFcXGNQZ1RwalBoa1peWls+O19WUjtmQE1GdFN1W1l0ZUN0PkVIX2FYckV3SD9mYENZXTt2YnFUQTtXeEtkeUtlTVVEbW15Q1lFO0Nybk13YGdJbGtCZUFXc2NDYmlkeEdpTHVzWj1mSW9SSHFjPm9XcUlTYD1UOllkYmVFY21ScT11U1N1Wz1mQUt2VUdVXkt0czt3a01FX2V4eWdjYFVoWVl3PHFmbD9HYl1jQkNScHlUdHN3Q3VVTEdUcGNVQE92THV2anF0cGFWOkVyVVF2XFx3Z1A7aUZZQjpFeXVheVxcP3NJV3U/cVhYeVRRT0Vtb2ZyYVhjZ1RzRWlrU1VHbWhmVWlfZ3dDd2VTSXRJX1U6eUJ5b3R5b2dvX3U7P1d3PVVzVWdQU3dFdVM6X0Naa1ZePVNOV1JvU3ZpbXVRd1JAd3lAc0hNa1hXa1RVXXVxQXd4YVhNO1ZwX3JbPWZjPWZvP3dlZVhoeVdsYWY7U0RbSVltS3lYYVJuX0lgaUdUYWZnU3NpbVdAW2VLO2Zqd3ZuPURTd3ZCUWVGaVk8cXI+d1JXVXJPU3Y9P1lQXVZrbVVaUUVIXVdHYUVEQ1Q6S2VdZXl3bVNyd1NYd0NZUWk6U1RVZURiU2hVc2d2PUVfYXlFd3RzSUdmX3haZVZPcURPd3JfTWNhYXl5aVlbO1lUa3VmT1dUPWh4UUVtdUJgZ1Rnb3VyU1NCRXdVZ1dqXUdVUVl0bUhmZVNHUWJUaVJCU0dpdUlEY1JwR2l4P3VgPVJSeUg/S1VzYWRvY1c9W2hKZXdOSUlQZ2hbdURDcUdlXWdpS2lOZ1VlSXY/UUY6VUhxcUdEVVhtO3ZCXUNLc0RuV2NLPWhLV2N0SUlJUVlOW1hJXXR0TXJ1O2NoY3R1W2lBW0hHdWJ3UXRcXGl5O0lFVUViSkl2clN3U01pdldXdltYXFxvVVhTY3ZNRFdtdW5rSFc/ZVdDd1dtdV1Rd0ZVeV5TSG1pZXBVaE5TaFpBVmB1UmVTdT1TYlhvY15zWXh5Uzp3dV15SVplZHlnZVpndEJddUNdSGZDU0pRVFZFRHBrR3dBRG5NaDtld113aWN5Yl1ndmtNVkh5aWFxUjpdaElVV0FxQlF3VVJLUlJreF9dWTp1RDs9QzxcXEp3eXFIeVFHRXdzQWxeQXBacG1JRFNkcHRLRFQ/TFhcXElVTURTSnV0WGVSWlRLbFxcTENddFdkVEVNajxkVUtQcEVEWEpxTEtdVlNxdHU9UXR5UjtcXFh4WXg9aVJhVE5CeFZsbHd2bXBhbVBXQGtpREt0WHNIdEp4eHBKYUpnQXM7QFlESEpQYXFQQU1ZZWxqQWxBYHNySVk8YU5qbXZpSW1LUW12WVZYdFd0WXI6UHhOaFFdRE10PVJBaWpcXHlPS01zVFFwSWhQb1xcU01ReHZkVD9dTXJJcFNVVE9gVHB0U1pFbXFJbktNV19QUF5AUnZEal9ddTpIcGRVVWRYU0BcXHJ2QU9PaVlNcXB1WHFkPVdFcVF0cGxWSFNYbFg8YGprWHVqZFJIRXVieXheeFRdVE5vcVhAWVBPSXRzTG5iZXJQdXRrTVFUXFxUbmVwRVl2WEhxS0hWO11MdWVvTXVNaW1RV3lYZFBLOk11dVlxSXR2V2BvS0lsUUx1R3FrYXFZdXFPWWV4YEF2WmlSQ0VrUFlyYGhPd3VzcHlMUlFLPjxuW1R1bmxwb0V2dmhWWUVwUnlNVVRqcWVTdV1QZ0hYQFxccGw8dWV1VGB0WUJQbDtYVF9ET0pBVUpFTj5Qal9oWTx5eERZb188TU5pWXd0VXZEVkdNcFRxTndRTktgdEtUbWFpc2JVc0FIc1d4bE9teG49dWNxSkZZS3NsclV4T1NBYz9wZm5hdmlIc1ZAdUpvc18/YEJBcE52dmBxY2VBd2lhX3hualNXWkV5XFxQUG5zZ3dgR11PaG1fX2Q8SW94V2FdYXhheFtLXmpaUmQ/UmF5Z1lrR1RLQ2xVRmp1dmxxZXM/ZUtNRlldZ3NnaVBZVUthUkhbcm1nYmVtV0V5WHBbRjtdSVU9d0g9SVFLZW13WFdNUkFFdkxrZV9VQ0VBdl5BcmJHUnVZZk9Vcl93Q0FBR0lxRExrZEZhaVVhWVxcaUh0XWhkPXdgR3g/eWlTcVRbU2VmV3ZLU2lLW2hwWVZuQ1NKc1lqa1NxZ3NvV1ZnQVNwUUVNPXdLaWhBS1U9cXdwV0ZTX2dZW1lzYVV0R0RcXE9DXkFHPENUT01ES1VWTmFkd2F3cEt3ak1iRkFIaUdWdEFEOz1DPE5iSk5qSl50WV5sTnhaaXFpUGBfPz9kbVhsSndoSmlcXFBvbjpmXnd5d2Zoa3JnYnFfY1VWaXV5WlpOaVRpXk9HZ0lIZ0VpZUNRZ0VxY0JmYkVHcERGXFxNcXhpT2VvRm5oR15jXm5WXlpMWXRmWFxcTWFyblhcXENmaFdXXUtHZ1xcT3U8T19dWF5yPmdLXmpaUkJLWmR5ZHhZWmFodUJxXnNPYlZnOzxKOmBuXFx0TlxcdFRbPFA7Mjs=</Image></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"> <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#H1" style="Hyperlink">Return to H1</Hyperlink></Text-field></Input></Group><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Definitions</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  R</Font> is the reference input (angular position in training required by the gun fire control system).</Text-field><Text-field layout="Normal" style="Normal">  <Font family="SymbolPi">e</Font> is the error between the controlled output and the reference input.</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">L</Font> is an external disturbing torque applied to the motor shaft that drives the load.</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">C</Font> is the controlled output (angular position of the gunmounting in training).</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">C/R</Font> is the closed loop response of the controlled output to the reference input.</Text-field><Text-field layout="Normal" style="Normal">  <Font italic="true">C/L</Font> is the closed loop response of the controlled output to an external disturbing torque applied to the motor shaft that drives the load.</Text-field><Text-field layout="Normal" style="Normal">  </Text-field><Text-field layout="Normal" style="Normal">N.B.  The servo is assumed to be linear throughout the operating region and the principle of super-position therefore applies.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Transfer Functions shown in the Block Diagram</Text-field></Title><Text-field layout="Normal" style="Normal">  The "Bode2" package will be required for plotting the Bode attenuation and phase margin diagrams.  It is therefore necessary to copy bode2.lib, bode2.ind and bode2.hdb to Maple's lib directory, for example, "C:\Program Files\Maple 9.5\LIB".  It is also necessary to call the "plots" package and to change the imaginary unit to <Font bold="true" italic="true">j</Font>, instead of <Font italic="true">I</Font>.  A reduction to 4 decimal places shown makes the display easier to read, so "displayprecision" is set to 4, (computation, however, is at the precision set by "Digits").</Text-field><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">restart; 
interface(imaginaryunit=j,displayprecision=4): 
with(plots):
with(Bode2);
with(inttrans):</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">a.  <Font underline="true">Pre-amplifier</Font> (including gains of synchros and error-phase-discriminating circuits).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K1;</Font><Font italic="true"> </Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  K1</Font> has units of volts/radian.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">b.  <Font underline="true">Power Amplifier and Amplidyne</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K2G2 := K2/((Tf*s + 1)*(Tq*s + 1));</Font><Font italic="true"> </Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  K2</Font> is gain of the power amplifier and amplidyne generator (volt/volt).</Text-field><Text-field layout="Normal" style="Normal"><Font italic="true">  Tf</Font> and <Font italic="true">Tq</Font> are the time constants (sec) of the amplidyne generator field and armature quadrature currents respectively.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">c.  <Font underline="true">Combined Motor and Load</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K3G3 := K3/(s*(Tm*s + 1));</Font>  </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"><Font italic="true">  K3</Font> is the combined motor and load velocity constant (radian/sec/volt).</Text-field><Text-field layout="Normal" style="Normal"><Font italic="true">  Tm</Font> is the motor time constant = the combined motor and load inertias divided by the combined friction and armature current losses (sec).</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">d.  <Font underline="true">Tacho-generator</Font> (gunmounting's angular velocity feedback).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K4G4 := K4*s;</Font> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  K4</Font> is the tacho-generator constant (volt/radian/sec).</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">e.  <Font underline="true">Accelerometer</Font> (gunmounting's angular acceleration feedback).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K5G5 := K5*s^2;</Font> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">  <Font italic="true">K5</Font> is the accelerometer constant (volt/radian/sec squared).</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">f.  <Font underline="true">Filter</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Normal"><Font bold="true" executable="true" family="Monospaced" foreground="[255,0,0]">K6G6 := K6*s^2*(Tc*s + 1)/((Ta*s + 1)*(Tb*s + 1)*(Td*s + 1));</Font> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font italic="true">  K6</Font> is the filter constant (sec squared).</Text-field><Text-field layout="Normal" style="Normal"><Font italic="true">  Ta</Font>, <Font italic="true">Tb</Font>, <Font italic="true">Tc</Font> and <Font italic="true">Td</Font> are time constants (sec) of the filter network.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">g.  <Font underline="true">Applied load torque to motor voltage</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K7;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  <Font italic="true">K7</Font> is the reciprocal of the DC motor torque-per-volt constant, and its units are volt/Nm.  It transforms the external load torque applied at the motor shaft into the equivalent voltage that would be required at the terminals of the motor's armature to produce that torque.  (Note that it is voltage that is summed at the input to the DC motor.)</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">h.  <Font italic="true" underline="true">G7</Font>.</Text-field><Text-field layout="Normal" style="Normal">  Let <Font italic="true">G7</Font> be the transfer function of the forward loop in the block diagram from the input to the power amplifier and amplidyne to the output.</Text-field><Text-field layout="Normal" style="Normal"> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">G7 := K2G2*K3G3;</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">i. <Font italic="true" underline="true">K54G54</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Let <Font italic="true">K54G54</Font> be the ratio of accelerometer transfer function to that of the tacho-generator (sec).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K54G54 := K5G5/K4G4;</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field bookmark="H1" layout="Normal" style="Normal">j.  <Font italic="true" underline="true">H1</Font>.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Let <Font italic="true">H1</Font> be the transfer function of the inner feedback loop (comprising the tacho-generator, accelerometer and filter).  <Hyperlink bold="false" executable="false" family="Times New Roman" italic="false" linktarget="Wks:#Block Diagram" style="Hyperlink">Block Diagram</Hyperlink>  </Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Maple Input">H1 := K4G4*K6G6*(1 + K54G54);</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Component Parameters</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The servo components had the following parameters (expressed in the SI units mentioned above). </Text-field><Text-field layout="Normal" style="Normal"> </Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">Params := {K1 = 14.8, K2 = 2300, K3 = 0.526, K4 = 0.98, K5 = 0.01368, K6 = 0.193, Tf = 1/31.5,Tq = 1/17.6, Tm = 1/9.1, Ta = 1/0.986, Tb = 1/5.26, Tc = 0.016,Td = 0.0016,
K7 = 0.63062};</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Closed-Loop Transfer Functions</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">   In Part I, the loop transfer function from the pre-amplifier input, <Font family="SymbolPi">e</Font>, to the controlled output, <Font italic="true">C</Font>, was found to be,</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K1G8 := K1*K2/(Tf*s+1)/(Tq*s+1)*K3/s/(Tm*s+1)/(1+K2/(Tf*s+1)/(Tq*s+1)*K3*s^2/(Tm*s+1)*K4*K6*(Tc*s+1)/(Ta*s+1)/(Tb*s+1)/(Td*s+1)*(1+K5*s/K4));</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Hence the closed loop transfer function between the reference input, <Font italic="true">R</Font>, and the controlled output, <Font italic="true">C</Font>, is</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/R` := simplify(K1G8/(1+K1G8)):</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  In Part I, the closed loop transfer function between the externally applied load torque, <Font italic="true">L</Font>, and the controlled output, <Font italic="true">C</Font>, was found to be,</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/L` := K7*K3/s/(Tm*s+1)/(1+K3/s/(Tm*s+1)*K2/(Tf*s+1)/(Tq*s+1)*(K1+K4*s^3*K6*(Tc*s+1)/(Ta*s+1)/(Tb*s+1)/(Td*s+1)*(1+K5*s/K4)));</Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Closed-Loop Response to Reference Input</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Frequency Response</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  Inserting the servo's parameters, the closed loop-frequency response is</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/R freq` := eval(`C/R`,Params):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">a := attensemilog(`C/R freq`,.01,100):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">b :=  pointplot([[-3,-3],[2,-3]],connect=true,colour=black,linestyle=2): </Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(a,b,title="C/R Closed Loop gain dB");</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The resonance peak of the closed loop response, <Font italic="true">C/R</Font>, is about 3.5 dB = 1.5 at the damped natural frequency of about <Font family="SymbolPi">w</Font> = 10/sec (since the normalizing time constant <Font italic="true">T</Font> is 1 sec).  The bandwidth of the servo, i.e. the frequency range in which the magnitude of the closed loop response,<Font italic="true"> C/R</Font>, is greater than -3 dB, is approximately 27/sec.  Both can be obtained more accurately by using the "linelog" command to plot <Font italic="true">C/R</Font> instead of "linesemilog" and then using the cursor.  Or, by substituting <Font italic="true">j</Font><Font family="SymbolPi">a</Font> = <Font italic="true">s</Font> (where <Font family="SymbolPi">a</Font> = <Font family="SymbolPi">w</Font><Font italic="true">T</Font>) in <Font italic="true">C/R</Font> <Font italic="true">freq</Font>, a function, <Font italic="true">C/R f</Font>,can be obtained, from which the maximum value of <Font italic="true">C/R</Font> <Font italic="true">freq</Font> can then be found thus, </Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/R f` := map(abs, subs(s=j*alpha,`C/R freq`) ):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/R max` = maximize(`C/R f`,location,alpha=10..20);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The maximize function yields not only the maximum value of <Font italic="true">C/R freq</Font> = 1.5, but also its location, <Font family="SymbolPi">a = w</Font><Font italic="true">T</Font> = 10.72, hence <Font family="SymbolPi">w</Font> = 10.7/sec</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Setting <Font italic="true">C/R f</Font> = -3 dB and solving for <Font family="SymbolPi">a</Font> yields the cutoff frequency, <Equation input-equation="omega[c];" style="2D Comment">NiMmJSZvbWVnYUc2IyUiY0c=</Equation>.  In this case, the bandwidth extends from <Font family="SymbolPi">w</Font> = 0/sec (direct current) to the cutoff frequency, so it has the same value as the cutoff frequency.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">omega[c]*T = fsolve( `C/R f`=10^(-3/20), alpha,0..40);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The bandwidth is an important operational parameter, since fast response requires a wide bandwidth.  However, too wide a bandwidth allows noise to appear in the response.  The bandwidth of 28/sec found here is a satisfactory compromise between these conflicting requirements.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Position Step Input</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  At this point, it is convenient to define some options that will be used in the transient response plots.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">txt := [TIMES,ROMAN,12]:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">opt1 := axesfont=txt,labelfont=txt,titlefont=txt,colour=blue:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">opt2 := style=point,symbol=diamond,axesfont=txt,titlefont=txt,labelfont=[SYMBOL],
        symbolsize=18,colour=red:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  While considering the open loop attenuation and phase margin in Part I, it was noted that a decrease of gain, <Font italic="true">K1</Font>, (the gain of the synchro chain, phase-discrimination circuits and pre-amplifier) could cause a sluggish, oscillatory response of the servo.  It was also suggested that the gain could be increased slightly without affecting stability.  This is examined below, by plotting the closed loop response to a unit step input of position with four different settings of the gain <Font italic="true">K1</Font>, namely 20.92, 14.8, 7.4, and 3.7 volt/radian.  These represent an increase of 3 dB on the proper setting, the proper setting, and reductions of 6 dB and 12 dB from it.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">   First, <Font italic="true">K1</Font> is removed from the set of parameters so that it can be given four different settings, each associated with a colour.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">Params := Params minus {K1=14.8};</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K1s := array([[20.92,14.8,7.4,3.7],[green,blue,magenta,maroon]]);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Next the closed loop transfer functions are obtained for each of the four settings of <Font italic="true">K1</Font>,</Text-field><Text-field layout="Normal" style="Normal"> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/R 4` := eval(`C/R`,Params):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/R 4` := [seq( subs(K1=K1s[1,i],`C/R 4`) ,i=1..4)]:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The Laplace transforms of the four closed loop responses, <Font italic="true">C</Font>, to a unit step of (angular) position input are</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">posn_steps := [seq( `C/R 4`[i]/s ,i=1..4)]:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The inverse Laplace transforms of these four responses are normalized by dividing them by the input (<Font italic="true">R</Font> = 1 radian) and plotted, thus</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">posn_step_resp := seq( invlaplace(posn_steps[i],s,t) ,i=1..4):
</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">a := seq(plot(posn_step_resp[i],t=0..2,
title="Responses to Unit Step of Input Position for\nK1 = 20.92 (green), 14.8 (blue), 7.4 (magenta) &amp; 3.7(maroon)",labels=["t sec","C/R"],opt1,colour=K1s[2,i]),i=1..4):
</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">b:= pointplot([[0,1],[2,1]],connect=true,colour=black,linestyle=2):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(a,b);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The plot shows how the response becomes sluggish and oscillatory as the pre-amplifier gain reduces by 6 dB (halves) from its proper setting and then halves again.  With a setting of <Font italic="true">K1</Font> = 7.4 volts/radian, the settling time more than doubles to <Font italic="true">t</Font> = 1.15 sec, which is probably unacceptable.  On the other hand, a 3 dB increase in <Font italic="true">K1</Font> causes only a slight increase of the overshoot and settling time.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The sensitivity of the closed loop response to gain reductions in the pre-amplifier and power amplifier would need careful evaluation during breadboard and prototype trials.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"> The normalized error (i.e. Error/<Font italic="true">R</Font> = 1 - <Font italic="true">C/R</Font>) is more easily seen in the following plot,</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">plot([seq( (1-posn_step_resp[i]) ,i=1..4)], t=0..3/2,
colour = [seq(K1s[2,i],i=1..4)],title="Error in Response to Unit Step of Input Position for\nK1 = 20.92 (green), 14.8 (blue), 7.4 (magenta) &amp; 3.7(maroon)",labels=["t sec","Error/R"],opt1);
</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"> </Text-field><Text-field layout="Normal" style="Normal">  With the gain setting of <Font italic="true">K1</Font> = 14.8 volts/radian, the error of the response to a step displacement has a single overshoot of 38% and the transient is virtually complete after <Font italic="true">t</Font> = 0.44 sec.  This should be satisfactory for a gun-aiming servo.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Velocity Step Input</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The Laplace of transforms of the responses, <Font italic="true">C</Font>, to a unit step of (angular) velocity input, or ramp input, for each of the 4 settings of <Font italic="true">K1</Font> are</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">vel_steps := [seq( `C/R 4`[i]/s^2 ,i=1..4)]:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The inverse Laplace transforms of these (angular) position responses, <Font italic="true">C</Font>, are then normalized by dividing them by the input's (angular) position, <Font italic="true">R</Font>.  That is, each <Font italic="true">C</Font> is divided by <Font italic="true">R</Font> = <Font italic="true">velocity</Font><Font encoding="ISO8859-1"> \327 </Font><Font italic="true">t</Font><Font encoding="ISO8859-1"> = 1 radian/sec \327 </Font><Font italic="true">t</Font> sec = <Font italic="true">t</Font> radian, thus, </Text-field><Text-field layout="Normal" style="Normal"> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">vel_step_resp := seq( (invlaplace(vel_steps[i],s,t)) / (1*t), i=1..4):
</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">aa := seq(plot(vel_step_resp[i], t=0..3/2,title="Responses to Unit Step of Input Velocity for\nK1 = 20.92 (green), 14.8 (blue), 7.4 (magenta) &amp; 3.7(maroon)",labels=["t sec","C/R"],
         opt1,colour=K1s[2,i]),i=1..4):
</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(aa,b);
</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  With <Font italic="true">K1</Font> = 14.8 volts/radian, <Font italic="true">C/R</Font> reaches unity by <Font italic="true">t</Font> = 0.3 sec and the transient has a very small overshoot of about 2%.  Although the transient has a rather long tail (caused by zeros in the numerator of the closed loop transfer function that are associated with the tacho-generator and accelerometer feedbacks) it is substantially complete by <Font italic="true">t</Font> = 1 sec.  These are important characteristics for gun aiming servos, since they have to respond quickly and accurately to changing input velocities.  (In passing, it should be noted that the accelerations involved are limited by the "g forces" that the target can withstand.)  As <Font italic="true">K1</Font> decreases, the response becomes sluggish and oscillatory, as already seen with a step input of position.</Text-field><Text-field layout="Normal" style="Normal">  A 3 dB increase of <Font italic="true">K1</Font> to 20.92 volts/radian produces a near ideal response.  However, it would increase the servo's damped resonant frequency and its bandwidth, making it more responsive to noise in the input.  For this reason, <Font italic="true">K1</Font> = 14.8 volts/radian, is probably the best compromise.  It would, though, be worth investigating the response to a noisy signal with the higher setting during breadboard and prototype testing.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The response, <Font italic="true">C</Font>, to a unit step of velocity input is simply the normalized response shown above, multiplied by the normalizing factor <Font italic="true">R</Font><Font encoding="ISO8859-1"> = 1 radian/sec \327 t sec, and is shown below with </Font><Font italic="true">K1</Font> = 14.8 volts/radian.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">aaa := plot( vel_step_resp[2]*(1*t),  t=0..1, title="Response (blue) to Unit Step Input of Velocity (red)",labels=["t sec","C and R"],opt1):
</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">bbb := plot(1*t,t=0..1,colour=red):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display([aaa,bbb]);
</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The plot shows that after a transient, the response, <Font italic="true">C</Font>, leads the input, <Font italic="true">R</Font><Font encoding="ISO8859-1"> = 1 radian/sec \327 </Font><Font italic="true">t</Font> sec.  The lead is caused by the velocity feedback and the values chosen for the feedback gain and filter constants.  The very small difference between response and input results from the fairly large value of the velocity error coefficient, <Font italic="true">Kv</Font>, noted in Part I.</Text-field><Text-field layout="Normal" style="Normal">  The difference between the input and the response, <Font italic="true">R - C</Font>, is a position error and is plotted below.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1">plot( 1000*(1*t - vel_step_resp[2]*t), t=0..3/2,title="Error of Response to Unit Step of Velocity", labels=["t sec","1000 \327 Error"],opt1);</Font>
</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The plot shows that while following a constant velocity input of 1 radian/sec, the error is less than 5/1000 radians after <Font italic="true">t</Font>  = 1 sec.  A 500 knot (257 m/sec) target crossing the line of fire perpendicularly at range <Font italic="true">r</Font> would generate an angular rate of <Equation input-equation="257*m/(r*sec);" style="2D Comment">NiMqKCIkZCMiIiIlIm1HRiUqJiUickdGJSUkc2VjR0YlISIi</Equation>, so the gun aiming error would be</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`aiming error` := (257*m/(r*sec))*((5*radian/1000.)/(1*radian/sec))*r;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The gun aiming error of 1.29 m should be acceptable.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The plot shows that the velocity error approaches a very small value asymptotically.  This can be found from the Laplace transform of the velocity by means of the final value theorem.  First, it is necessary to reset <Font italic="true">K1</Font> = 14.8, re-insert it in the servo's parameters, and evaluate <Font italic="true">K1G8</Font> with these parameters.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">Params := Params union {K1=14.8}:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">K1G8 := eval(K1G8,Params):</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Assuming that the input velocity is <Font family="SymbolPi">W</Font>, the transform of the velocity error is given by </Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">epsilon(s) := 'Omega/s^2/(1 + K1G8)';</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The final value theorem states that <Equation input-equation="limit(y(x),x = infinity) = limit(s*y(s),s = 0);" style="2D Comment">NiMvLSUmbGltaXRHNiQtJSJ5RzYjJSJ4Ry9GKiUpaW5maW5pdHlHLUYlNiQqJiUic0ciIiItRig2I0YwRjEvRjAiIiE=</Equation>, hence the steady state error resulting from an input velocity of <Font family="SymbolPi">W</Font> is </Text-field><Text-field layout="Normal" style="Normal"> </Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">ss_error := limit(s*eval(epsilon(s)),s=0);</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"><Font encoding="ISO8859-1">  At full display precision, this is 0.5585019637e-4 \327 </Font><Font family="SymbolPi">W, </Font>which is identical to that found in Part I as <Equation input-equation="Omega/Kv = Omega/17905.04;" style="2D Comment">NiMvKiYlJk9tZWdhRyIiIiUjS3ZHISIiKiZGJUYmLSUmRmxvYXRHNiQiKC8weiIhIiNGKA==</Equation>.  With <Font family="SymbolPi">W</Font> = 1 radian/sec, the steady state velocity error is about 1% of the velocity error (5/1000) when <Font italic="true">t</Font> = 1 sec.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Poles of C/R</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">   The poles of <Font italic="true">C/R</Font> are the roots of the characteristic equation of the <Font italic="true">C/R</Font> transfer function.  This is obtained by removing any common factors from its numerator and denominator and equating the latter to zero.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`char C/R` := denom(normal(eval(`C/R`,Params))) = 0:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/R poles` := [solve(`char C/R`)];</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  These are now plotted, thus</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1">plot(map([Re, Im],`C/R poles`),-16.. 0,-16..16, style=point, symbol=diamond, title="Poles of  C/R", labels=[\302,\301], opt2 );  </Font>
   </Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The plot does not show the distant poles since their effect should be small.  The damping ratio indicated by the complex poles nearest to the imaginary axis is</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">solve( sqrt(1-zeta^2)/zeta = abs( Im(`C/R poles`[6])/Re(`C/R poles`[6]) ),{zeta});</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"> This value of damping ratio suggests a more sluggish response than that obtained for the unit step displacement.  However, the servo is not a simple, </Text-field><Text-field layout="Normal" style="Normal">second-order system, to which the standard response curves apply.</Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Closed Loop Response to Load Torque</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Frequency Response</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The closed loop response to a load torque input is given by</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/L` := eval(`C/L`,Params):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">display(attenlog(`C/L`,.01,100),title="C/L Closed Loop Gain dB");</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The dimensions of <Font italic="true">C/L</Font> = <Equation input-equation="[K7]*[K[3]*G[3]]/(1+[K[3]*G[3]]*[K2*G[2]]*[H[2]]);" style="2D Comment">NiMqKDcjJSNLN0ciIiI3IyomJiUiS0c2IyIiJEYmJiUiR0dGK0YmRiYsJkYmRiYqKEYnRiY3IyomJSNLMkdGJiZGLjYjIiIjRiZGJjcjJiUiSEdGNUYmRiYhIiI=</Equation> are <Equation input-equation="volt/(N*m);" style="2D Comment">NiMqJiUldm9sdEciIiIqJiUiTkdGJSUibUdGJSEiIg==</Equation><Font encoding="ISO8859-1"> \327 </Font><Equation input-equation="radian/volt" style="2D Comment">NiMqJiUncmFkaWFuRyIiIiUldm9sdEchIiI=</Equation>.  This equation is normalized by dividing both sides by 1 <Equation input-equation="radian/(N*m);" style="2D Comment">NiMqJiUncmFkaWFuRyIiIiomJSJOR0YlJSJtR0YlISIi</Equation>, so that the result is non-dimensional and can be expressed in decibels (dB).</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  <Font italic="true">C/L</Font> is the response to torques applied to the <Font underline="true">motor shaft</Font>.  If the torque is applied to the load, i.e. the gunmounting's structure, instead of the motor shaft, it will be reduced by the gearing. </Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The firing load torque can be visualized as an approximately triangular pulse with a duration of about 3/4 second, and can be represented by a Fourier series with frequencies of <Equation input-equation="4*Pi/(3*sec);" style="2D Comment">NiMqKCIiJSIiIiUjUGlHRiUqJiIiJEYlJSRzZWNHRiUhIiI=</Equation>, <Equation input-equation="8*Pi/(3*sec);" style="2D Comment">NiMqKCIiKSIiIiUjUGlHRiUqJiIiJEYlJSRzZWNHRiUhIiI=</Equation>, etc.  In this frequency range, the attenuation of <Font italic="true">C/L</Font> is at least -88 dB.  The gearing ratio of the train drive is about 200:1, so the throw off in the training motion would be less than </Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`throw off` := evalf( (75000*N*m) * (1*radian/(N*m)) * (1/200) * (10^(-88/20)) * (180*degrees/(Pi*radian)) );</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Ramp Input of Load Torque</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The firing load could be considered as a torque that increases steadily for 1/4 second up to a maximum value, and then decreases steadily to zero for 1/2 second.  (The actual shape and duration of the pulse is determined by the design of the gun's recoil and recuperator systems that attach it to the gunmounting.)</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The (angular) position of the output, <Font italic="true">C</Font>, in response to a unit step input of torque rate, i.e. a suddenly applied ramp of 1 Nm/sec of torque input, is given by</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">torque_rate_resp := `C/L`/s^2:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">torque_rate_resp := invlaplace(torque_rate_resp,s,t):</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal"> At this point, it is convenient to obtain a function, <Font italic="true">f</Font>, for the response,<Font italic="true"> C</Font>, to the 1 Nm/sec ramp of torque input thus,</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">f := unapply(torque_rate_resp,t):</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The response is plotted below.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1">plot(10^6*f, 0..0.4, labels=["t sec","C \327 10^6"], title="Response to Unit Ramp of Load Torque", opt1);</Font>
</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  The response, <Font italic="true">C</Font>, to a torque rate of 1 Nm/sec for 1/4 sec, followed by one of  -0.58 Nm/sec, lasting for 1/2 sec is shown in the plot below,</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">plot( piecewise(t&lt;1/4, 10^6*f(t),t&gt;=1/4, 10^6*( f(1/4) - 0.58* f(t- 1/4))),<Font encoding="ISO8859-1">
t=0..3/4,labels=["t sec","C \327 10^6"],</Font>
title="Response to Triangular Pulse of Torque",opt1);
</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Assuming that the firing-load torque increases linearly, it would produce a torque rate of 75 kNm / (0.25 sec) = 300 kNm/sec.  Since the gearing is 200:1, the maximum throw off is given by</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`throw off` := evalf( (300*10^3*N*m/sec) * (1/200) * (simplify(f(1/4)) * radian/((N*m)/sec))
                * (180*degrees/(Pi*radian)) );</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  Thus, the throw-off produced by one gun firing at zero elevation is between 0.6 and 0.9 degrees.  This should be within the amplifier's linear region when controlled by a 36-speed synchro, and so the assumptions of linearity and super-position hold good.  However, this would need to be checked using actual firing pulse data.  The extent of the throw-off and the time taken to recover from it are important operational parameters, and would have to be measured during prototype trials.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Step Input of Load Torque</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">  The (angular) position response, <Font italic="true">C</Font>, of the gunmounting's training servo to a unit step of torque input is shown below for interest.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">torque_step_resp := ( `C/L`/1 )/s:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">torque_step_resp := invlaplace(torque_step_resp,s,t):</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1">plot(10^6*torque_step_resp,t=0..0.5,labels=["t sec","C \327 10^6"],</Font>
       title="Response to Unit Step Input of Load Torque",opt1);
</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">  After <Font italic="true">t</Font> = 0.36 sec, the transient is complete and, in the steady state, the load torque applied at the motor shaft keeps the gunmounting trained away from its correct position by about 18.2 <Font family="SymbolPi">m</Font> radians.  It is analogous to a torque twisting a shaft, or a force compressing a spring.  Thus, the steady state torsional "stiffness" of the training motion from motor shaft to gunmounting position is 18.2 <Font family="SymbolPi">m</Font> radians/Nm (3.75 seconds of arc per Nm).</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Poles of C/L</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">   The poles of <Font italic="true">C/L</Font> are the roots of the characteristic equation of the <Font italic="true">C/L</Font> transfer function.  This is obtained by removing any common factors from its numerator and denominator and equating the latter to zero.</Text-field><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`char C/L` := denom(normal(`C/L`)) = 0:</Text-field><Text-field layout="Normal" prompt="&gt; " style="Maple Input">`C/L poles` := [solve(`char C/L`)];</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1">plot(map([Re, Im],`C/L poles`),-16.. 0,-16..16, title="Poles of  C/L", labels=[\302,\301], opt2 );  </Font>

</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">  The poles are identical to those of <Font italic="true">C/R.</Font></Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Observations</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">1.  The choice of compensation, gains and time constants was discussed in Part I.  The closed loop frequency and transient response analyses in this worksheet demonstrate that the values chosen should have satisfied the operational requirements for a gunmounting's training motion.  The effect of sea-induced ship motion would need to be examined before leaving the theoretical design stages, but that is outside the scope of this analysis.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">2.  The classical analyses shown in this worksheet are part of the modern design process and are often included in specialized software for designing servos.  This is also used to try out different approaches and to optimize the parameter settings, taking account of non-linearities, backlash, etc.  Despite the efficacy of such software, the designer's skill and experience is still the most important factor in achieving a successful design.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">3.  The final hardware will depart from the theoretical design because of the assumptions that have to be made - for example, about the rigidity of the gunmounting's structure and the ship's supporting structure.  Extensive experimental work with breadboard models and a prototype almost invariably follow the theoretical work.  In the case of a naval gunmounting, this work would have included live firings in heavy seas and against drone aircraft.</Text-field><Text-field layout="Normal" style="Normal"/><Text-field layout="Normal" style="Normal">4.  Long, painstaking work with slide-rule, tables of Laplace transforms, log-tables, and semi-log charts (and eraser!) were needed in the past to perform the calculations and draw the plots shown in this worksheet.  The transient responses were particularly arduous.  Maple does it all quickly, easily and without errors.  Not only that, it can feed optimized code to certain software used for designing servos, thus reducing setting up times and errors.</Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">Reference</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">Page 378 excerpted from <Font italic="true">Servomechanisms &amp; Regulating System Design</Font>, <Font italic="true">Vol I</Font>, by Harold Chestnut &amp; Robert W. Mayer, originally published by John Wiley &amp; Sons, 1951 - 1980, copyright by The General Electric Company 1951.  Used with permission of The General Electric Company. </Text-field></Input></Group><Text-field layout="Normal" style="Normal"/></Section><Group><Input><Text-field layout="Normal" style="Text"><Font italic="true">
Legal Notice: The copyright for this application is owned by the author(s). Neither Maplesoft nor the author are responsible for any errors contained within and are not liable for any damages resulting from the use of this material.. This application is intended for non-commercial, non-profit use only. Contact the author for permission if you wish to use this application in for-profit activities.</Font>
</Text-field></Input></Group><Group><Input><Text-field alignment="centred"><Image height="33" width="800">TUZOV3RLVWI8b2I8Uj1NRExDZE5WWlpKOkBMPkg6VEtHeE1rSjo8T2BMb1xcbFF4bFFXZE1XcHNIcVNobVdoWW9lWE9QbVRQbVZgbXZxeXhxPVhqPXhYcXVYYXhuYVhjRVdjPVVSPVV3ZVl3RUxLRExxdFBxPFI6PXJeYXZedVJBdXJaQG5adFZhdVZiPVdiTVl0TXl2YXl2WXl1WVl4bVl4cXl4cVl5dVl5RVlzRVlwbVhweXl5eXlwcXhwPUo6Pjo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6OjpkeTxUeXBDPnFVTENUSmNEWG9YdXNUPGF1cGtjZldNWEBKQ2VVYGROdVRtV3h5eXlwcHVQQ0RTU3VMQ2x1Pjx4VHBRbWxzYl1NaWhVT2BxVGVYU1FPO0BKeFZdd09sOkBzeUZ2PHdcXHRAdHNOblFuXFxWP3c8d1xcP0ZxSmlqWHluWlZ2bnlIRXJtaUJfX3RXaXRbTXl4WVJJSVh2V2d0U1M9O2dRTXdBSUNdSVlyR1hSb2djW0VwcVl0c3huPUJWU1VHdUVBW1d4S3JXYVNIc3NvWUJQa3luS2N0cWdteVVLQVlRWVV3X3JzPXdib1lUV1hJP0lRS3lvW1hAd3lkcXl0WVJHQXlgaXhzW1NseVhhU3lxdXk6bWVsPWRYcXlkSWZ2Z1JJZVNVa1VtVUJHd3VaaXRTO2VRP1M+QWRNYXNua3lTR2JEU3VpbWJTYWJqeXROQXlNdVhsYVRXYUNwO3k/YXQ7X3R4YVR3YXRoP2NqPUdiZ1lWR0NBW2VBa2heaWh5YUlHb1ZkR3h5V2VRYXRhbVZIWXg6U0VJZXd5YWNtY1NCQXZnT3l5c3NFeUJWV0N3UUZ0WVd4WWRNZ2NZX3leVXk/Z2NlW1dYUUNEY3dHdXdITXc/cXd4W2dhY3NjR3J3T3R1S0ZYS3NjW0ZaSUJPcUlySUlda3VJQ2ZSb3NNX3lUU0VXV2NLUXNfcUdIZUlpYVdCc3ZhQVhXb0ZzWVR5dUlZU2RXQ2V0W2ZacE9ZdHZbXFxYU012Tj1YaGx1eGVsXXlsdlVuO1BZc3F2a21tQ3hTRVFQc01PZVVwUUVLTmB5VkFxY3FSUXBZeEhyW3hVXFxBdGdQVmV4bUhIUVlEWHB0TDtleV9cXFhIeHlUcExRPXFKaEprbHFBPXdQeHFPdHBQbXdRPWtXZFNTWWp4aFF0PWxpPFg9UHJcXEhveE1LeHBwZFVQR3hsYDxSYWRXc0VNVWhuTWluYXF2eVxcdF1wSndcXFB0dHQ6bHdfaHk7UHh1RWxXcGZ5cGlReWc8SWJnSHFRP3dSd3ZGZ2NRbm10SV1sWFpvYXV2d1xcXVZpXFw/eXVJakdxeUFfXWpeY2lhXFxedmFZZm1YWXZWX2ZveWRfd1phP3lJUGZOWHBPaW1iSW53aWllUXlaQFtqZltwX2BzP1xcTkBxYXdbPGFfPXFwZEl1XT5nbkhwVWlcXF5hW0FHY1NfeV1wbkhnX29JaT1Ya01gYkteeVVXakZoaENwaWY/bGxoZWxoa0txaz1xZ0NxcUlva0phZFpAXUlPc3BIamdRZ1V2Xk1wXltha1hOb2t4Y0ZheE1YPkVmeD1HSnlZXT11S1dYdWVmY1lDVl9ETztYXW9lRHdJXVVyaElYaEtkdFlndj1zWU14eU1oRUFiZEtkRkVEO01CaW1VWWd2TnNmQnVEZ3F3XnNSWm9pZXlpWUVmRUFzWU9jVTt1Zl9DXjtnPkVJVW1XeV14WltIP1VUaXdoYXliPEVXVUFobWdoVWVlXU9ETHlma1lkT1FETk1zbGVnXW1IR2t5blVyclVoamdidnN0cklDc09pVT91cFVodE1FX2NWVWV5d1dyU2VTdkl3SHFzRVV2d2FTYG12X2tDRWdERUVWT295ZlNGWUdYaFt4ZTt3ZnN5YT9IYmN1X1NpSFVmclN0cXNnSUNVS21SO0lFR0dpRVV4U1Nld2tCUmNpYz9mW0dIc11XQkNlRlNYTWVjQHF3UVlpT0NGaTtiZF9lcGdoQ2NyU0liclVGZktYcE9FPkNkR1VWSF9zcz1HYUVGXFxNaF91REpjWGVXR1NrSUE9VGBbdWhPaUtPeTtJZG9fc0JRZ1BHYmlNeFpJeFs9Uk5RSENVd2xJaFZBcz5NeHY9dDtJZWtlY1tpVG9lQl1ZU1ZzSV1VR2tNZ0M9eE1fY3ZdckNrR2xPeUU9d1ZzeW1vUlBFUkdVV29Lcz4/ZE5HY3FPdkw9RGNnVVVpZD1TZEJZdGFjQmN5VDtzQz8/c1hzQkZFSVBLZHdVaWJVVXVvd3RDeExFUnhHVVBPYz1lZVdXREpfdEJJRmpbUk1XWG9hSW5pRkRZeXZJZkZZSDtFaWZhV0FBZGtRZ1N1SW9ZSFM/c1xcYVlua1ljQ1JYQXk7PXVyU3NVRUdYb3Zta2RVP2JJa3V2SWhmO2hIS1Jtc0lxa0drQ0lFR1NRaVV5P3JbY2h5XURXP1VKd2VvX0hJO0lbaVJQdVlDY2VdeUlRR1NSPVNGY1lASUhOYWJFeWhUO0hcXGdDW2lpRXViWElZWz9GaGtmQWFSeWNjUTtEPE1CTGtzVUd2TV1GT1NXWmFGbm1VVk9CXU1oYGd1XWV3OkNTWFtWVVtkXmlXQ0lUTWtpbmdWbWNZO0V1SWtGWmdldGFTbGtlRF9TbFVkP1NVW1doYF9JSGt1TmFJQkVZQEtoUVtJYlNmbF9DcGdWXUlCZ2NmOkNyT1dXbGlWUFNETXVFa3dCWVFiZ0t4R2lXZmNkZ19jQ29YRHlGb0FGPENZZF9mWlNVS09YbVVFcm12cFdnYVFJZVdHeU1pdU9maGVGWVtVV2dkR3dlWztYQFloPG93c2tUd1VnallkdkVoblRQYExKYXRVbXlvXXhsa1VwZ1BTSG1TT2lTWHRNP0hzSGhXZ2xudT15cE1vc21QV1F0WG1sTERSXmVyYXBwQVBxQFR3dVxcbWY8eXRNb190TlFEbXd1VUJhbFtUS01dVVpcXFZzVVBnXFxPaFhVXWl3PmxUPlR0b2xZVWVNXFxgcTppTkZRa01ldUI8WV55cVtUcXdMeHlZa15tUERoVVRFTFtteGRZVHJVd0hZcHBgUl10c3lobTxcXHJkaE5cXF1WR2VqRXlUQkxsWGhVaWRTa2xWY0lta3VKQVxcT0ZBSnhYVEpcXG9ScFVyXFxxbkVVZjxQT2FvY2lvWHhZVVRSeGhtS0hub1V1QmF2dnh0XUBvcmR5cUlsYHR5Y0V5Zz1TdDxWO0xZYERvREVsQ2hXWWRrcElrU01vcGhuaHFrZU1XPFFYXmRvZ0VtTTxreEFZTT1tcFBLbVRUTW1YZVFMbnVLP0hNZUlVYGBUcU1TZGVOcW14SGVMSz1PVXB4XkBraVlwYHhYVmRvVUBMPVBwckFQSXVSW1FwQFlsdlBXd1FUb01wR2BqT1h5Rmh4QUVUaWVSQURLZ2lvVlBPeVhVbFhUOkl3YzxOZ2VNTnVwXFxYV3JkUUZQUXZsUD1Ub3Nlbz5xWGJpV09cXHlFPVBVaVBBQVNnTHR4WExHPVNUQVNBeGo9QFdpeHdYYFhPQXRIbG9JZW9IaUx2eXVvdU10THRUeUpzQXhCWHJAVHFXWE9zRUtvcHVBRVU8dXlPXFxMVHlQQVhtPXRPVVFuZWFORF1LT1l5THlYYnR4dWhtY1lyWE1raFxceWxMb19lcWB0U2VBT0hdbHFVd2lQbmtQd2xIUGdIcmVoWV5wS2hQd0dQSjs8TzxgcVU9dE14VVVFUFdAUmRJVGZZamphb3dUcU1RalhISlNcXE08RXZhcHBUQG1XTUpAaU9WaHlMUUtxXVQ9RXljPVVocU5hXVBKXFxYXFxMdVtEc1FAT1tYUnc8UmJgUGB0U3VlamNlWVhAVU49ckZleHVIbURta11YUkxhWUVsUm1JUF1QZWNoYHJ4bWE/YXJhYUN4dldRW1xcYVpgeWlGQWo/Z3ZWVmReQG1HeVtoaGp4UXZqSXdNVndQR3lYV19FcGpETm5zeV5FaHZFX2Q6UG5rT2FEQV5DbnhFQW9DaF9ld2M7cGJbSVtad2NVP2twR3d4dmNWVlxcT1dhWUdaV3FiR0dealZrQVFdbVhja2Z3VFZmb3ZaVm5aTHdmb0llUz5lQEh0Y3ZzZ1BuPFlxRE94Y3FiZE5tUHh0cXdoc2ZhZz5teU9lZGhxQ0ZrTldxc3B5XUBfVlFySUl1XW5jTEliPl94ZFFeW3l3XmBeWXFiU3hleWdhPk9rVkBmcFZmZU5obXhlU3duXj9fR09rbGZgUXFnS195Sz95akBweHZ3Ykh0SWB5WWFpP0h2Sl53dlF2WW5nQVZvPVhod2NSZUJJTWZsS1RVX2JgcXJGUUM8VUdSV1k9a1ZXQWl2XVg8Q1N5TXljeXdlb0U+P3R0a3NWZ0JUbXRHSVh2S0RUO0RgYXRwYUdRRVZBPWVmb0hAXVRnc3dzQ2ZXR0ViQ0NMSVl0U3dHO3RSYUM/XWhpW1Rmd1NQVWNTUVlaQ3Vsb0VbS1RuT1NUdURQcWZwUVVfWXhbP1VaPWJgeUN1RVRVZWN0Y3JzYVdJR2hQVVZkQ1hvW0RuO0dUb2Y9QVZCY1lSR2dhYVlic3Z0PVVCdVZJT2VaS2dHbWhIUXJdXXVtc2lmeVRQV3RuZXlaS3lkbUhqb1dSQXNTUUhld0RTPUhqXUM+cWRIW1hISWdrd1RHdXZJX3NnWURnYWJTc2lMWXJiXUljW3VaVXVDZUdOXUlueXlqaVZuTXVKaWJxXUU+PXNIW3RoUURYZ1RcXHFoTndUVm1HZG9TaUtzRF1ERF1VT2tzTz1mWDtYdklkYlV3Umlpc0NFdj90RUFTP2VIW0VIaU95W21jRT9oWTtld0tDclt4O0VDcFVFYUl0Uk1VZU1JQHdGPUd1cUlkcmlYbUFpSG91Ql1VRWt2Ym9EYF1iRGV1XlVIT3N4d0tTb2dWRV9HTlFiQkFkdU1ZUTtZX11YYnFCZVtGRllHRj10WGd4cnlZcEFGRG9pZElSSGdVZj91WEdnXVdndUdpZ11VUlFycDt1PU1IWUlYeGNJYW1zcUVsPHVSPFBNd3R3Tk1xTllNQj9cXGFJaXF2Ym94aGtud0RPdl1ecjphXFxbV2hFeHNuX2NkUW9ATmddb3JMUG5DcHRFP3dKcWk6YWRgP2dqWFxcQm9sOkBkSmlzW3ZlbF5wSz5dVHBjSUhob1Nab1hKT2h3W1dnc2VzdUJmRWddPXV1VVk9cVhaV1ZZTVNaRUNIV0hxZVg8U3VeRXV2WVg7QUZRUUNdXUZsXVNOcUlPPUlMUXdoSXdab2VxRW9PcVZZQFRUcHJXQU5xWXN1eE5BQFdqbHB1YVh5dG1YTVJrZHBJXUtcXExUQD1QZFxcU3hISlNYTmh1bEZZUW10d0poV0k8UXN1UlVwd21cXHJRREx5dU1nTXY+QHBTQHBmdFJpVW5pVFY6dVJSaWw8bFJZPHdsdFNWaUxoSEtEQHZWaVNgRE9mYVR2QXN5TXVLbVFVaHZxbFF1TFdAcWxyYFJkZFJLSW1eUVlBYVh4ZFBcXFR1VmxrdE1ZbXlQQWB4Uml2UlVvTHhLbUFOYWxMYHFWYGVURElPO01ZXFxIb1FpWW5Na0hMTnFoeWxVSlxcdFNedUtKSU1LQVlbcXVmTXJ4QVhmeEp5WHhlYFJQcXhPaW9ybEpXXVhFSFh3XFxsSnFyPVh3TjxUPmBuRlBrbEh2XkxUZF1rdml1Oll3bGhXa1R5RHBMU1VWVXFRQ0F1VFRsaVBvcHVvVEhOU1F5UnRzPklxS1lLaFROUU1zZUFqb2FsclF2YklzbE1wPVxcb2pMVU1EdURReW1hb2lRdWxtUE1FTHdocHVwbG5JdnlwUGBYbENETT5MWUBgcmRxdG95bkBNTEZUVVVQb1xcVVdSXFxXTWV0T0FvRWV3TElVY3RSd0B0XUVSR0BYdHFLdUhRV3FqV0xxWmBMVFVPVHVzbUhQY1lrP0ROPXVUXFxhWFNlTE51S3J0dGZAa0l1blVUWENNdFl5UlVRcGxYd2BYdj1pWHBwdUxtUlVxd1RNbVtdcXhoTEVsdD5sTmlAcVE9UV9sUkw8TmdlcmhoWHdBcnlBTD1pd11JeFlUVXloajtwb3FYUG1VZ0hHXFxnYW5mV2ZGPmhyQXd0d3lbWXM8VnVHWGhTR3hlUGpNXmV4blxcdmFiSE5qVGZmRll3RE5yZUBxb2hlSFdtb1dgXVBcXGdmcV1Ja3h4XFw/dmtubmNcXGdpdXBvdkloTWFaT0lraklkVnF0dj9lZm5oZWBpPU9peFZ1ZVZvcHhqSk91TllgW1dcXGpYXFxTTmtlcXJRX3BVZ2hqTmlOUXRwR1xcQ0llX0lhYllzQHd3QndcXExgeE8/cmBxWmk/Y0BXc1dgXkBmam9nZXBwamtJcG5Ya0tQbmRHYWRHaWRvY0U+bT9GamZfYllmXFxcXD9wXUhpZU5xV2dnZUl1Q0FuaGlad2FlcFlua2dlRnlqdk9odV9bR1FrcGlvU05hP25kaXByVUZqY1ZcXHBRbmd3XVI/XVdGZVd4YD5pX0hAdEF3ZGJueTx4X19PYEZ5Z2dxdWpBdEpoYWlBblNBcz14d3RwXmFZbmxvbG4/ZVlRdEFebUp2d0Q/a1xcUWxdeHFNUGNgX3NqVl1ndnJlT3NJT2twUF5WeV5bVndgT1tnd21McWldTm1aQGhCQXJpUF1PPltASGRtWVp5aXJbTm48WXBlTmZvbnNvXl1kbmZJWXVYd2tFQWNVeW5eQWBdVmV5WXVsUG9nQW47P1xcSz9tdF5ncF5qWEd4Zj55c2Zac2d1PWBzZWJfYUlFU1NKY1dld3RtQ3JFQ2ZnRVJhcUVOQ2hCO2ZeSXZ4WUw9UFNdPXlLWG1HZU1ZTG1yVFNCcExfYFVBbG1YbVhsVVRYRW5eRXNTbW1meVJFWHNERXdlbHZRcWxRYVhAQHRqPHBrVFlrRFNOcXhQUWpsdXNpVEpFTFhRXFxSd2BzUGFTVVlKd1BqZGVzX1FzS2BqQElqX0R1Rm1KbVBMbWxsaDxTU1BLVjxXW2VPYWFUTkB3TGx0dj1xZEBPT0hyYzxLPmh1aFBQPUFwU1VSUF1tYklWU3VybERMcXBLdWFWbGlWPklvT3hKeEx5R1hPaHF0PVFQQlFWSXRSamRWP11QRlBQQ3l2c11ZQl1SWEFzUEx5c1FUXk11TFVPRE11ZURQPVVQcEhzRlV4OlhKYGhObEVZS3lrcVFMUUhTRXVyXmFYX1hKSF1VeXh0Z01SQ1h0anVvP0VRV01MW2FSU2lraWRvZUxzVWR1V0VNdGhZWnlRW3F3eEhUW3RPdTxWR3hxYmBxcDxPUUFXT2VZSUl3XlR2YEhyTnlQO0VLaERMaVRxY1hMcTxOWGVqc0VLc2VUO01ZQTxvc211ZkBVQHR4VU1KWWFNRnV2VmFqVWVsdl14WGBuY3VUaFR4Qlxcd3h0dkNpdUBIc1FVUTptc0p5VVZYTE9lVUFMbWRhWV1UTW91cUVFeFdgeEs9UVFMeUdBeWlIUFxceE9mXXRHPmNKd2BneHdeZl1tSWRKd2dYaXliWF1fXlxcXXhdd1hvb3ZmSmB2Z1FrbFdyaHFgc3hxVGhkX0F1WEhvdGF1eHF2VlBzPmZYUUVHX1lHeXVqR1dxYUNPeUU+V1hbd3VFd3lzTUhzQUNhd1lmc0lpcXZXaVdwV0dvR1ltcXdBZWg7X1hxR1N5W1lRVVc8a0ZhVUdtdWhxZVlFO3hkd2JEVURkV1Y8T1lqbXdjXXJMP1RwdXdGX3NuV3VtaWlhQUlueUJbYVVieXhcXHl5YGNTTG1IeHNJbndZTHdmPW9iX2t0eGdVSldUQl1UdEl2S2tERE1JQ01WWkNIPFdXRjt2WGV1T0dlXlFlTHdpa11Ia0NmclVYdV9EZ29DW09JeXVoX0l5YltlRWhxcnlRP013VGV4SXVOYnVtdjxzT2l3eV11Tz5pZT9vTlhwbkZiXWl5a3l2QHBuTT9eYlFiY09wXUBwTV93T0laXFxpXXRWcEdJdT1QZGJIZk14Y3hYYXQ/YVdQWnN3dz54YUR2djx3cVF2eWtecGlBcl9AZmRZeWZveHNhY3RXX3V2Z0JQbXF2bUtfWk1BclpXWnlBdkNQbXVZZFxcQWJacF1aTmdYd3J5WGF4dmE+d2ZZcGNaZ2VtPnV4aXVbR2lZbnV3UXU8YWlKbnM/XFxVTnBxSGdqZndocVtiYWhiQHhDR2JIVmtrX25UUGVpb2JmeWNVZmBYbmF4aWRsd2lUSGptaGVGP3N3PnFXWHhUV3lnUWJ1cFp0WXBncXBrd3dmV3ZjSFpjQXdbaXVNaXliXm1FZnloX3l5WHNJSW9zWGRKZnh2cV0+eWFSX1pWeHlcXGJTP0ViQXdzXXddd3ZjT0ZvTWh3U1VSYWd5Q1lkaVR3QUJ1QUVHV0Z1U0lHb0VrS1lJR0ZZVVlddXdgdXdYb0d1QUZWV2tHd3F5ZmJAcXJyaWZqP3NZcHU9QF9db249Z1tRQGx0UWJRTlpEZlxcRldlXFx5cXV3WzxwdV4+bHZReFxcWXc8d1xcPFZ4UlBuPXl4aU5bQ05nQl5pck9wd0duRWZ5eVdudHF3Omd3RWZaU3BpX0dcXDw/YFFueFY/d3lnbTxOWl5xeWFHcHh4aU1wa19PaHFZcld4XFx0QHQ/QHZBQVxcZXFfclFxdj51eUB0eWFgV3l5Onh2bXlzWHd5WWZbTVd4b1dtSWd2b0U6O0I6TVRLV0RLV2dKO2VaMTo=</Image></Text-field></Input></Group><Text-field/></Worksheet>