<?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" linebreak="newline" name="Text Output"/><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]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" name="Text" opaque="false" size="12" 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" underline="false"/><Font background="[0,0,0]" executable="false" family="Times New Roman" foreground="[0,0,0]" name="2D Math" opaque="false" size="12"/><Font background="[0,0,0]" bold="true" family="Arial" foreground="[0,0,0]" italic="false" name="Heading 2" opaque="false" size="16" underline="false"/><Font background="[0,0,0]" bold="true" family="Arial" foreground="[0,0,0]" italic="false" name="Heading 1" opaque="false" size="18" underline="false"/><Font background="[0,0,0]" family="Times New Roman" name="Author" opaque="false" size="12"/><Font background="[0,0,0]" family="Monospaced" foreground="[0,0,255]" name="Text Output" opaque="false" readonly="true" size="12"/><Font background="[0,0,0]" bold="true" family="Arial" foreground="[0,0,0]" italic="false" name="Title" opaque="false" size="36" 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 Protection Current Transformer Under Steady-State Conditions Of Operation</Text-field></Input></Group><Group><Input><Text-field layout="Author" style="Author"><Font bold="false" foreground="[0,0,0]" underline="false">Dr. Wlodzislaw Kostecki<Font executable="false" italic="false">
Senior Lecturer</Font>
PNG University of Technology
Department of Electrical and Communication Engineering
Lae, Morobe Province, Papua New Guinea</Font></Text-field><Text-field layout="Author" style="Author"><Font encoding="ISO8859-1">Copyright \251 2001 by Wlodzislaw Kostecki All rights reserved</Font></Text-field></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 describe a current transformer from an over-current protection system.</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 style="Text">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">Analyzing 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 waveforms generated by the current transformer</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">1. Problem Statement</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">1.1 System Description and Data</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">A through-type protection current transformer <Font executable="false">(CT)</Font> with a rated transformation ratio of <Font executable="false">800</Font>/<Font executable="false">5</Font> has an annular ferromagnetic core with a mean diameter of <Font executable="false">60 </Font><Equation input-equation="mm" style="2D Math">NiMlI21tRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and an effective cross-sectional area of 100 </Font><Equation input-equation="mm^2" style="2D Math">NiMqJCUjbW1HIiIj</Equation><Font bold="false" italic="false" style="2D Math" underline="false">. The material of the core has a relative permeability of 5000 for small values of magnetic flux density. The primary consists of a single bus bar passing through the window of the core. A high-voltage ac transmission line current, which flows through the bar primary has a rated frequency of 50 </Font><Equation input-equation="Hz" style="2D Math">NiMlI0h6Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">. A load (burden) connected across the secondary terminals of the CT is the operating coil of an overcurrent relay, which controls the trip coil of a circuit breaker. The resistance of the operating coil is 0.05 </Font><Equation input-equation="Omega" style="2D Math">NiMlJk9tZWdhRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and inductance 0.15 </Font><Equation input-equation="mH" style="2D Math">NiMlI21IRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">. The combined resistance in the secondary is 0.07 </Font><Equation input-equation="Omega" style="2D Math">NiMlJk9tZWdhRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and inductance 0.2 </Font><Equation input-equation="mH" style="2D Math">NiMlI21IRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">.</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">1.2 Objectives</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">A. General objectives</Text-field></Title><Group><Input><Text-field><Font background="[0,0,0]" bold="false" encoding="UTF-8" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" size="12" underline="false">\342\200\242</Font><Font background="[0,0,0]" family="Times New Roman">  Provide a general description of the <Font background="[0,0,0]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" size="12" underline="false">CT</Font><Font background="[0,0,0]" family="Times New Roman"> as a component of the over-current protection system.</Font></Font></Text-field></Input></Group><Group><Input><Text-field><Font background="[0,0,0]" bold="false" encoding="UTF-8" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" size="12" underline="false">\342\200\242</Font><Font background="[0,0,0]" family="Times New Roman">  Describe the <Font background="[0,0,0]" bold="false" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" size="12" underline="false">CT</Font><Font background="[0,0,0]" family="Times New Roman"> in physical and mathematical terms using fundamental laws in their symbolic form.</Font></Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Develop a mathematically tractable model of the electrical-circuit configuration of the <Font executable="false">CT</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Produce an equivalent circuit diagram of the <Font executable="false">CT</Font> model including designations of its parameters and terminals.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Obtain symbolic equations describing operation of the <Font executable="false">CT</Font> in terms of the magnetic flux and electrical quantities.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Solve the set of equations to obtain symbolic forms of the functions describing the basic time-varying parameters.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Derive the symbolic functions describing the supplementary time-varying parameters.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">B. Specific objectives</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">Use the numerical values given in Section 1.1 to</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  obtain plots of waveforms of the <Font executable="false">CT</Font> variables,</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  obtain numerical forms of the functions describing the time-varying parameters,</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  compute numerical values of time-varying and time-invariant parameters.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  verify Kirchhoff<Font encoding="UTF-8" executable="false">\342\200\231</Font>s current law <Font executable="false">(KCL)</Font> at a node of the <Font executable="false">CT</Font>-model equivalent-circuit diagram for an arbitrary value of the independent variable.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">C. Additional objective</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Restore the symbolic forms of the functions describing the time-varying parameters.</Text-field></Input></Group></Section></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">2. Description of Current Transformer: Functional, Physical, and Mathematical.</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">2.1 Fundamentals of application and operation</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">A current transformer, known alternatively as a <Font executable="false" family="Arial Narrow">series transformer</Font>, is an instrument transformer for the transformation of current. Its primary winding is in series with a conductor of the power system whose current is to be measured or controlled. In through-type <Font executable="false">CT</Font>s, the primary winding is a single-turn winding that is the line conductor and is not an integral part of the <Font executable="false">CT</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The protection <Font executable="false">CT</Font> is a device that acts as a buffer or interface between the power system and the protection gear in the event of a major fault in the system. The essence of a <Font executable="false">CT</Font> is that it should reproduce in its secondary the exact replica of the current waveform present in its primary. This must be ensured over a wide range of magnitudes of the primary current.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">To specify the highest magnitude, relevant standards use the concept of the <Font executable="false" family="Arial Narrow">rated accuracy limit factor</Font>, which is expressed as a multiple of the rated effective value of the primary current. This multiple ranges from <Font executable="false">20</Font> to <Font executable="false">30</Font> and the <Font executable="false">CT</Font> is required to replicate the primary-current waveform over such a range as faithfully as possible.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">This implies that the magnetization characteristic, or the <Font executable="false">B-H</Font> curve of the <Font executable="false">CT</Font>-core material must be <Font executable="false" family="Arial Narrow">linear</Font> in the entire range. In other words, the core must on no account saturate until the primary current is above the accuracy limit factor. For this reason, the operating point of a <Font executable="false">CT</Font> at the <Font executable="false" family="Arial Narrow">rated</Font> primary current is chosen well below the saturation region, typically at <Font executable="false">B = 0.1 T</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">If the <Font executable="false">CT</Font> iron-core became prematurely saturated, the <Font executable="false">CT</Font> would no longer reproduce in its secondary the exact replica of the current waveform in its primary. Instead, it would produce a distorted <Font executable="false">(</Font>flattened<Font executable="false">)</Font> version of the fault-current waveform. This translates into a false <Font executable="false">(</Font>lower<Font executable="false">)</Font><Font encoding="UTF-8"> magnitude of the secondary current, which \342\200\223 in turn \342\200\223 delays the protection gear to respond to the fault, or it may not operate at all to clear the fault. The likely result is a damage to the system and its components.</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">2.2 Assumptions</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">To develop a mathematically tractable model of the <Font executable="false">CT</Font>, the following assumptions are made.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(1)</Font> The secondary winding is uniformly distributed on the <Font executable="false">CT</Font> core and there is no leakage of the magnetic flux from either winding. All the magnetic flux is confined to the core and links all the turns of either winding. In consequence, all the magnetic flux is a <Font executable="false" family="Arial Narrow">mutual</Font> flux.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(2)</Font> The effects of hysteresis and eddy-currents in the <Font executable="false">CT</Font> core are <Font executable="false" family="Arial Narrow">negligible</Font>. In consequence, the active power loss in the core is negligible and so is the core-loss current.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(3)</Font> The magnetization characteristic, <Font executable="false">B = f(H)</Font>, or the <Font executable="false">B-H</Font> curve of the <Font executable="false">CT</Font>-core ferromagnetic material is linear for the primary current ranging from its rated effective value up to the value resulting from the given rated accuracy limit factor. In consequence, the relative permeability of the material is <Font executable="false" family="Arial Narrow">constant</Font> in this range.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(4)</Font> All the electrical parameters of either winding and of the burden may be represented as <Font executable="false" family="Arial Narrow">lumped</Font> parameters. All these parameters are temperature-independent and time-invariant.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(5)</Font> The primary-current waveform is an ideal sinusoid.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">2.3 Equivalent schematic diagram</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">An equivalent schematic diagram of the <Font executable="false">CT</Font> is instrumental in carrying out its physical and mathematical analyses. It is convenient to produce such a diagram as viewed from the <Font executable="false" family="Arial Narrow">secondary</Font> side of the <Font executable="false">CT</Font>. The pertinent diagram is shown in Fig. 1.</Text-field></Input></Group><Group><Input><Text-field alignment="centred"><Image height="286" width="573">TUZOV3RLVWI8b2I8Uj1NRExDZE5OWndqOnRLPkg6VG9JTU9Vcjo8T2BMb1xcanl5eXlLPlxcSj48S0I8S0ZcXEtGPExKPExOXFxMTjxNUjxNVlxcTVY8Tlo8Tl5cXE5ePE9iPE9mXFxPZjxQajxQblxcUG48UXI8UXZcXFF2PFI6PVI+XVI+PVNCPVNGXVNGPVRKPVROXVROPVVSPVVWXVVWPVZaPVZeXVZePVdiPVdmXVdmPVhqPVhuXVhuPVlyPXl5eUk7Sjo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6OnJHU3ZxUUw6XFxrVGBOXFxATmRcXFFncXhoWGp1ZFRgVWtITEo8TE5EaHFlWXl5eXFkWXRWZVBAUWtCQFN5V2NOQXVKZmI7WXBRaWVLVmg9QHNGcXZuSWE9SGZCPlpyWFxcaz5dPEJmSUpwQUo8ZUo8bGpaaU1KTVhYXFx2eD13eHBQSGxUbVlNWmhOZ3RucmxwPmRzb01MP21wQ1xcTHNAanBhaztFb0s8V2hYdEppVmVMbG1wU1NtTEBIVUNJS1Jlc0lhcj5heURgdmJJT2NkUGFgdmBJVVp1c21EVEdxcFxcWFE+aHleQXE7bFFRWUxDWVBMTUxaXXU6YHdzeEteZE1XSUpuQFRzRWx5PW1YbW9BWXZPPVZJSWxFRWpBdWx3VHZyQFRTdFFkWXJkPU87QW1xWVlwQFlKYXM9QFdbeFNeSXNbYW5beVRwYXBeRE1DRVk+cEp0TFN1TG1YZVZVaVA/RFlvRXZ5PXFiTVVHSVRAdU5hZWt3VXk+bExXeGtTbWxPcFVFPU5TUUxkUWxvVHh0SU5PTE5iRWVed2tpX3llX1t5X3lkcHg9eGhUQGNiYXhcXHZlZFFpdmFcXE92d1JxWmBvcFFedWd4WkxIb2JAWz5gb11eckBBaW12X3N2cWI+XWJwYElUVktTW1tWRFVlZWFCSFlzYl9mbWlSQFVyRFVDWXVFQ1VkRWtEWV1UXFxJU2NNdm9LclQ/clBHSV5DY2JLVmZJVVJrYmhhSU09Y0tjY1Q/YkpVVkx3RmJDRjxbd2tbV2JteGVvV0JBZ29ZZFBPdl1fRT89eVRvSVVVVmdPdHRzaEZnc1I/ZVxcaXdOcWY8T2RkPWhdT0ViZ3U9XVdmQXJSU2lOWUZiWWlVW0M7aWdOdVlJV2Y/dXJKdWNwS0heY1dqcUVJZXRQV0NfbURmaUdBS1dxQVdTSUNsO0V1SWRvU2RdU0hsXXdba0ZTVVRdeUZ4cWVFc3RQSXddc0ZtbXZzU1lAV2ZePXJgQVJaZ0RUW0ZdXVZ1ZVRQZ1JNRUlXeWNUQ1k+T2JAYVlVVVNweXZPdVU6Z1JrS1ZLYUNOZVhbRVJsd2JMPURaSVViZUNQPWNcXGFzOlV4d0lSVGlJc1NmPlNkdmFIbldkQkNYUD1oP1FGO01zWE9oP3d4Wm1kXV9CRktHTUVCW0FoVFdGP29GcXd2WmFUP1l0QXlIdz1ETV9VQTtFVUNjR0lCZEdoY2d4QV9UU2FocnFyXT1EZ0V1TW1US0tmWV9mZ1NIUFN1d2loYkFlV3dDU11UY0diYVNXSmFIS2tTUnFlVG10bHdDQlV5VlF5TGFIX2FZWUlTckVlUGtEXXVHYEFWRWtCPEVjeT1GdW1lOm1la1FSUGdobFlUWG91UGFjbGtYRT9GRm1TWElpTW9WRF10dj9kSnNjSV9ESEF2Xl9oTm9FY0NTV3NHWmNjXnVic2NkYmNmXktGc1lkUXNkb1dDWD1ZW01GQ1FoWGVEXFxzc1hrdW5LYnd3aGRBdEdRRkJpeGFdY2dvZUlXU3NpQ2lTVHZdVUdVU2VZdmZTeFdvVl5ZdUNZUkhjV2dlRk9jVm1jZGxnUlNRaHJPRl5fV3ltV1lbR11ZVHh3R1pHWT9vRVl3Vl1FdEdrU2RNY0tDWXBBREVdcmBZWDtrSV9TeG9dSVVhVXhDRmhnc2BteUNDQ25vVkZVVj9xWHJTZmJxdVNVVGBncm5vVHZneGhTZWVLdV89ZjppR01xZ3RJRE5VZmJRV2RDUk09c1tXWVprckhRWF1PdVl1dVRvZXJlU2B5dl5ndU1lWURdVGlFaGxdQko9YmNHZGdlWE9DcnlreW1nY051cktbRWg9VlN3WEZjZXZnSFhVZVtTRD5PU1Rzdk5hVW5jZmV3RV0/VF15VXRrZ0k7ZTx1Vnc9aEVnZ2x5VkVVSEdHd2BBYnR3ZmpbY2tpV2hjZVZfc19HRXltQkt1eHhtQjtTd2JRc0ljUkxJZVxcSUdfcUllTWhtS1Vdc3hyR0Zsb3hJSVI8Z1hfb0JzZXhdPVdkVVk+aVdjS2RUQWRBY2ZxR0lJX2huVXdaUVk8X3d5PVVIVWM8a1ZqYWJTYVNUeUNSc1Nwd1Vac1VZb2hEVUhnR3VJc2lIX1VPVWRNY1dXW1ZkcWdZa2RvVVJ4V0dpS1V0WVl1PWhVS3lXTXJmYXlKZ1RxT2hBc2ldcURtdVNsV1liTVJ0TXJKdUZmU1RNU2JubUNkVUdRUUdCd0lJWUZLT2RGX1VZd2ddcUlcXG1HYD9Gb2VXaXNzZUtnWFtSUHVTZkNHaGlHbFNndW12PXl5PEl5OlN4VlVIbztWU1d0RVlJZGNiaGtYWHdYQ0t0bmd3WGl3PT9FZXdkUmdncldGZWdEZHRwcllsb11sQU1uclFuSkBzQ1VvU3huQDxPTDxPZVh4TmR5TmF1VE1rV2B3Pl1TSFRMcXVyX0V1QnVYPUBxP1FueHVXbkBQYXlVXXhSdVBPeFlNS3VTSE1RYXBWU2lVTkhRa3RTR3RwQUlYXzxuP2xLSGhYblFwaFxcS0txbHdESk5gTG9QdWF4a1l5TUk8bHc8cFhYc2FMdWFtaz1UUmVweHRYbUFseHZNUVRgWT9cXHdYPHJXQFBheHhBZUs7bXh2TVF4aVBNPXJYPHJ3cFhNcGZwYWU+Z15EaXFFYWk8P3dFP3JYWGtSZ2tHPm5DSV1ySXdnYHhTbmRtdnJ4WVpUYG07VlpYVnlPPmVaTmlYUGpJQGZJZ2lWWW89YWQ9VmloWVpmUGE8bnhWeWtpTltgTnlmbmpCX2g9aHNmZmw7Vm9TPmhVTmdQYFxcZXByZndmPEdybnduRD9qR09cXEU/Zkk+dnBeeXJmcHRXaXdfcGNpXl9pcWVwaF5oYFFnXk9QdE1HX2o+bUlHcmd3c3R4bF1AcnJnYzxueGZZYk1eWldxeUFQX1xcTmhRbnFKSWVOUGNQSHdvd2NmUV1Hd2lsdm9SaXBHWXA9YHldX246aGF1QV1zVnZLaHlAZ3lUeXlEcF9zX292b2JQWGNtdmZtQHBQQWFJR2VZXnNYcHBnT19peXBQeGxESWFpSVxcT0FnVWlhTkZuTXBjX195dGBnRWFcXEp2dGBgWltnXVdnclFeXFx4SXJmSWJyZmlSYG5yeG5uaF9UcV11QWs6X3hpR2xnbl5fSGBeX1xcZk9da2FhaU9gQ0ZieWFmcGFoTWBtaWF4R19yXWlhU0FpdHd3R09iXFw+XFxgb15SX21cXEBfZk9hXUdjPj54P0BrdV93d09naUdlRndxQ2hwS19nVmF3cmh0WmljPmhwW295O2hhcVZuXVdha1luZHdmcXFaOk5tRVltTXZzZEFgY0FedW93b1huW09zO29YO1ZJO0JSWUpmSTpxRTo7aEE8dT1KOjxyWDpueDo+SjxKUlk6Vmk6Wjo+dT1aWHFbSlliXVlqSlFkQ3hfTmF0Sm50PT9vZEdfSHFxVGZoPUZ2eUFcXGRAXk52dXlXdm5hZmNQYUVmcHVBclg+clRubkhPXj5ob09pY2ZncUJ4dXh3amdoZEpXaXlpc3B5X2hgd1R2XFxnaWFsP3lfcGNeUF54YGtcXGhuWWFxST5lPWd1WXh0VW94QklncG9hUHF3Y2duYVBrW2BbbT9bandfdHZtZWdgSXFmaFlzYnZbZj5hQ2hkTV5paFhbZklhdGdaQVZ5Wldubk9tUG5vcElxZmBqQXl2ZE5bc1loZF5hd3dzYHhxVVlobk55SFhsTVh2YklwYUFwZVliUWZ0bmBzblF4U05lcj91XWdtaWl2X3hcXHZ4eFRIZ3d5cHM+XmlBdVdPcGZ2aTpgX2RRbnZvc1FZYmxxY1RYZ2JwYW53XFxbcWB1SGFCP2pIUXV2dmNHWHZZQWxgaWZVeWpGdlt1T19jcXI9d2FpdmxWV2d3YFtqcWZHVmR3bnlHd3VyYG9VUHlbSGZJZ2duVmlmcWZedmVTd2J0YW5yd2FHb21UXmNLPmVRblp1eG5VaGluSG50aXRsaG1vcWtiWXNQVnRHQXdGb3hZQXRqb3N0YF9Xb29mSVxceHFjUFd4XXlgXUZnWVdbUHBlP3lxYmZodkl5bW5cXE95bWVGbk5HX3ZuYVFZYFp5cE0/ZW1vc25ZXnReYE9GXWdWaHhRZEw/bF9mcVxceVpoVl5VdmM/aWxfZndtPmx0cWRrcV1eQXVSZl9FVmJQd25zSGdiT3lkcXNgSHdfX2pyUXROd3FxaWxsV3RMcWxTdmhuR109YWs/cWRYaWtDeGhsSXFeYWo6QGR3T3htTmVwV3g9R3JvcXVlWG5pR3hveFtiX2xeaHNLP2FwWHg/X29QbndLYG1PQXlNSHR2SWZyX2lodnhhQHNAaXZbT2JNSF1xcWdAYWV2WWxxVmxyb2B4aWRfcXJReGtteXVKV29xb2p5WGV5YWZKQGY8YV55TltBSGVMd1p3VnJRWGVrcWpmSW9haGdDSHU8P11FR2Q/R2A7R2dpQHlIZl9FZnBUV2dVX3NLaGlnXnZ4SGdbaHhlWGhSUVtbVmlUcW1fYWNrXmZ4dmxaR11MTmJzSFp1WWZ5SHlbcWlSSG9Nd15qZmRlPmk6WGVkeF54WWhkUWFWVmxrT3JCaXRYZ25VSF1Pb2RrYGZ5RmJITm93dnJ2YVxcQU5dRFBabHlkXFxoamA/dV53c1ZJa0pnclBXd3JYcHJQcENZdGJpeF5JYmk+Y2VJYVtBXUVhbW1fcFRnb3Q/b0lJX1NIXmNAd3NnYWtweVhPU3lmPD1zT290a3FYcGVTXWtZXktXZElJaUNpQHloRm9WTl90PT9JPl1zU3l2YGtzQDt5P2tSSEFoZml3aFdiTEdITElkP190YElTR1doXFxHRDxteFJ3UjxTSG1TZkljUklnd21xR28/aG13RlhdSXJNaFVxc3BdVz51R0J5V0o9UmJTWTxFZ2ljV3U7ZkZFVEZvV2tVR05rWE47RztXSV5zU2BjY3hpRURTZj53RXBhSFhVVlVrWW9xQ0VjWUBDVT5dQlVlZ09zY01DeFlBR19lTmpka2ZEUDxFS0dQT3lVT0Rga1BUckw9S3NRaldcXHNwQVRrdFlyaHhbVFBfYEpEUG5QcFNzdGtIVVh5aFB5dFRsYGpySW9iTU5HXVdsWFZQRE54aHJndW5PbFJ3eHdGbXd4WE1KTGtraHE+eHk8QVBgcWp2UFByZFRqXFxqZ3lPWGlqTFxccXNVeGdMS0VYbllddG5gbXVNcGhdb1t1cVtobFxcTGxpZUxmSHZCUXZGPE94aU9GTFlfUGpcXGRqaUxWdnlWWT13TnFZYD1PPmVUa3lXSWxPbFVyT0BsaVRVeXlYUWVycWRXWXVZYHh0S11LcER0VVFObmlwVHRZRGhxakx2UlFYS1xcU2RlWG1wVHh0c1hEc1BtWWRITD1FU1plU049dHJZVklBU3llVVNFTGN0UGVwbm15eWdBVUdgTVFZcVBtVUA8S21kbT5ZT2hEeXFdTGZwVjpZWGpcXFB2aVJLUVY6aW9JSXJRSVJXTFc7UG5seHRWRXVQeE1ZWFNKPFZCWW5lPFlKTG1WRVVneVhHYVc7RHdsSXNiPHZPeU5kbU86XW5KYVFgUXNRRG46VFk6dHA7UXlNTE1eZXhFUXJBaHhidXhOaXdlPGtSZUpXQE9AUW9wQVY8TE5AcFhkUVhTQVlVWUx1TVhiRVZUUUpDQWxKXU9BcHJqTHFLUVJmeWtmQG9DbXFQZGtWXFxZWG1ZP0hWbnFuckxuYFBVQ21WQlROWERPQ1VYdllrQEFsd2BVaVlVc2hxXXFXS0FTeXF0T0F1QWVwT0V4QWxqdkxzajxYcFFwPEh1SVFsbHFXQ2hsP3B4VEB0eFFVR2lURHVNWWxMX2VYQVFNT1VqW2xrSmxLPHhZc0xqVFhTWExUWlFZWEVUTlBqWEBNV0B4c1BzTFBQc1lYPGFUU11sW11zXVRqXWBPY3RvZmxvV0hwWWRTXmhXaGBUQFxcdWZhakFxVVBNdGF0U0lVT1txbkldeGV5bFR5c1ZdT3d4V3Bha3hATExgdD1gdUBxbHdMT3NNVUBoWEJkUE5ZVnVJSj9JUztsdVxcRVdTPVdZQFZJRUxXWHBJeGpiXFxZZ3h0VEhrck10YXVOZV1TSHl0XXRybXhPSWl1dG1LbWBuP0lUYEFwUWx0VGhRR0VYbG14PWxOXFxVUF5Vak5BTVVodl5Bc1l1SkdxcWRlczphdnVcXHNMcGxASW1vPUthQFVzPHBFXFxSaFFyPHhMaGlNUGlUckV5SExzXFxET2JRbFlpdD1YVj5tcXU9clxcdG5dYFJZeHBSeHdabU9xZE5WRHNNVVB4aHA/SWpmSVVYTU1geG1xPGdDeG5Tdl91aGM+SV1oUGJ1cXFiVnhoaV15WHZtZmw6X25MPmxhP2xPaW1gdmpHSG50WWVmYHdseG5zSHB4T21mZmZdeWt4eGt3cGVvZ2hBZlo+Z2Rob3lYPnlxeWBfeWF5WXVFeFxcWFltd2ZpXFxuXFw9X15zV2h1P2tgeGZXX3FhPnJFTmZmX2dtWGY6R2hCRmVTT290WHRWV252WHVUWG5oXmZWeF9cXF5dYVdibXl3SnFgblhqZFhyVVZteFh2a0lmO3ZtdXlsdz5bdXd4SFlkX3BfR3BycXhcXFM/ZVh4d2FocWlxYWdYZUJQb3BIXFxbcF5PSHhLX2FeSHBSd153TmBtWXNTaWtsUGpCSGZLV148aGM8aGBQVmNwcXJYZmt3UXhyP11JaXhnWWRSZnRUd2JcXHl3dFFxcEl5WndgWWBab15hakF1Vz53VVBcXGNxY3ZReDxfZHZHYWE+XFxrQGk+cHBnWHBobmdiPmhlR1tedmhtV2pJeXd2R19pV2VMQG1rbmJSb3Rfdms7X192eGl1SG54cXdfRl9eT2VeZ2FmUHB3WXV1XnNhbnhSWXhVSWN2UXA6WHdMPmVWRmxyWV5QcHJ3WXNAcVxcdm9scklxUE5lY2duQF55SXhqQHlkd3lteFh5Ul5eQF9sVV53Tnl2a3dvdV5acFZtR194cFldVl9wY2ZicF9gbldvc1dkdVBrQWd4SFZ5RE9kTm93QldoZ3Fdbl5lZGZnPFFjb195TkdpSmBqSWZnTUd3W2ZvdnBxTD9rYld1aHd2X3lxZHZzU0l5V0dcXE9XcF9YbnRwa2pfX3NueFJxeWpQdVFHXFxwcXk7UHdbb21weGZicGBOeW1eWGZJYWtUeXNAaF9hYFh5dUNTUjtXWT1HRGRLd2xpaXdRdU9JSUBBV2VbUlRLST9HWW11aXVlTkR4alpgcUNpUGVxbm5NUnVoVm9ZU0VBWWp1eHhYbFBZS3VAVWl1d0ptUU1MU0Jhc0xIeEFsdmBRa01gcz90b2hAVkpUbHV1TEBEcHJMWF9AV0dlcnZwSnhhV3VIS0l5al51a0RgbHhUcUNpTW9kUXFUVV54d1p1U3FIdkxNcWRkbldNTFN1cmVkdlNIblBNWVJ1SnFxWVZVT1dNdU9pa1NgdVdhdkdRUlVkS3Z1dWhIWHNAcG55SlVtbmd4Tl9RcGlAc2R0cHVZbGBgTUc9dHZQcVxcPEs8QW1peG1zZVZdXFx4cUVPRUhZT1BNW1BvVlByXWlOYWh1PFVyPFF2cXFwd2luTnBSdm1qWU10ZExqVGVVa0V1dU1wX2RwU1BVcFlsX0VscmBScmx4UzxuP1xcd2BEd0ZZVXVxbWFNSm1pWXFVd1hJVj5JeF1pT3JkbV5wcHZtdXh5d2lNcmVFVkl1cFJOcz95dmRYd3FodktgaXFWaXJxZ1tvbHNAaXRAZl5oXmhmdndYXT1Ad3BWcmxJaj5OW0FWaU5WdnBAa1tJZmlBcEE+cGFnd2tYaFFgeVNhdXNwdmxgd0BJWkZfcndpZ29QYEJJW1FfeWNJalJWdHJQZ0pZbEx3dGN2eG9mc2tfZVRGWkhJcGFRbj5wZVhYeEBmanROYlNPa1xcaVxcbHFlQ3ZlYWBebmlzXUhieE9jT0hkR1BuZF9oYFZjZ1hiUWZaZEFyXFxXXWxfdTxwW3JfZ0ZXZmd3X2tgdlZgZnlWbFhmblVfYFFAZFpOY2FfYGFXW11xcHB2dktHb0BPalBmbD9gZ1JWbz5hY3dYZGtnc0dWW21eZ1VhaFtHXUVwdEpHdlxcYXM+aG1pYF0/R2VQWHhlWFpJRmJKVmtuaV1wSGZlcXd4bltHbmd4SFxcblFmZEdsdHFrPHZcXFRec3hQeUlvc21JbWVnZTxBbV9pX3h2ZHlWaWhveVFIcG1JcFBGZ2ludmpBW0ZxZm9od0xOYFZRZV5uW2ZmbWBgYXFud0JxYFZZcVs+Yl1nYD5Zckt5eGRobllJcUlGeUpxaEFQW054b0I+dHBWZDtXb2FYc1lObVRJYXdHZFlPY0d4WkleYmtIak5Ra2dZX19YdHZuckxYdD1IYFRYXUNpYkpYbFxcX15CYF1EQGFgX3FrQHFmP3RLd3A/WWl0Z3BNaXBjWGxnUVo7cWxIYHNUTnJ5V2pxQWl0SXBnXlxcVGlaRndmRl9lYE90clZqamZpRm9mcVFadl95SnhbcFZmY1dhPmZsQ190UlF4b18/b1NjXVlTZ0lXWXVSWWJ0Z1NHR1NZY2I6U1VtZ1RwdVJRU3RzeWhRQ1Z4PXlAd2M9cVZhTVQ9WXRQUXlEZ3N3Z2g8Y0NDX1NBc2l5X1Z5X0dfdUNoT3dlVUZud2g6SVRcXHFna1lmTmVZX1F4QmdpSW12P1lyPG11XFw/c0JTaUlzV3FbWG53R3RtZW11V2VTVlVPY3NpaUxlUllfRkBzZ1JTVzxlSW5fWVxca2hmU2lob3NOY3ZEXWdfUXM9U2hpX1lAW2ZAWXhVUVhzUUNJWVdecWlhR0Q9UUJicXhdTXR5TWh0Z05NTVlHTXFseXVbbWtyQG1CVFFZbUtwPVVkXU9uWHdPPE9saW9seHJbeGtneXVFWW1rTGxjXFxYckFycVhXUFl4RW1MYEBXdHlOYEBWYnlvdD1RP2lYVmBQc1RSPGhRZFRsTGFMX2F2ZlFYPHlub0hVdEVqaXhNVW1VX1hSXXVNZ3lWQ0VxU01vdWR4Z2RtSXR1SDxURXFZcnl2XlRqXFx1WUZka2hZV3R4UW9VbXRQc2B0VUtNV11Ybm9gdUhheVRNdGZRVkxATmN4bVxcRFFkPE12WFJZbG1ebUxbXFxsPEVxP3RSdjxWQnV5XkF1V0x0akF5UWlQZ3RMZkhvUVBNUERrOklvVUVtYlRRUGF3eFhLVz1sc0hKT3FzXk1ZT3RUQVF2RWBURF1LRXRWcFVRZmRUUUhNbHFOcVxcSnFBc1hYU3A9ckVxV3U9a2JlUEVhWD1cXHFcXFVSPE1YbFFSaElYU2F0SkVTPV1OUHhMPFlwRT12SE11blFzbFBUcmhyUmBxXFxpbz5obmVpdm5kamlMalZJUD1UcmtVdj5kT2xRdU5kT2RtUlBEdk1FT0NYdVJta2ZAUGhBVHRpTz5pc155d3BJU2dgbkV0c0VpcUhZc1hsWV11V2FwamlNeFZhUEZcXHBCTHhrZG9ZQG9mXFxTSkFLSVh2X11XO2Vya3RUW2hQR3lMRGBWeTxNbkBNS1BPaV10UVVxaGxzRGFubkx4ZEV1WHhWSWFXYV5xQXhwY1BbZD9wb3BpeVFvXFxAaUBveHFmaD5hcElmbDs/X0phd1ZOXXJYak5WcVhBYm15ZVBGbHlBXFxRaWFoUF53UGRVWF1UcGlHUW5OYFt0QGRubmpZWHd1eGZab2o/UHJBcHJrSF5xUGZfdmVgdmlCT2llZ3ZOT2hMPnlIP255UXhhYXNgV3V4d21dQGNad3lbPl9JWGBYVnN2UXJbVmlsQWVkP0pZcjtFU1ZxR1hveT1Fcz9zZ09bZFs7ZHdjWHdNV3dfZmtBWTtFZUFVY2ZzYktkU11NcVZZcT1Ram9ZamdIdU1VeUd1VkhQckhxbWRAa05kdl90VFhgVl1gT3FNbVxccGxVeW9WXFxydGh2b1lQVGFNWGFVPFxceXhhd2tAcFFAS3ZIUTxBb2J0eHR0c0JEa2NJZEJ5aHNfbFNQdWhfZmtxdjxgc3F5cHVPeWZZX3BWZlFpY3BGXFxzTnhudmhrSXI8ZmNgdnhkbmtVYHBCcGpjXnhXX1xcYndwaUlaeWhfYT5fP3B1aEl3XFxQXFxPXmlCUVxcaVluQklyVHFseWBkXmF3O19lalh2aFhmU1huWFdzR1ZaQ0hhaFdwTGhzQkBqQXFhOldzdHlhPj5pVmBhaU5bc09meUBpZWhpalFzR05tcWhqbUBmcXdwSU5iZndwaz9ySWBmY15xeWl2T0ByWHhbPFdpST9qSUhhRE9hTV9fQ1dlQ0h4YHBuQ0hrOnB2PF5hO2ZdTlFgR05hb09uTFhvTVFgcWhkXnhlOnFdOnFzcVBhQk93VHFiSG9mdk9fSU5pbVF5Tz50VnddPnhhU09kPlhkUVhrQFlnOnFdOklhalhbbFFoZ1dtPE5tRT5oVGZvS29eYHlveEZbUll4Oz9hakhqVU92OnloWVllSmFlWU5rPEBjZnBvTXhqPkZ1PT55PVB0TGFeaUdvUlB3PVZpOlZlQWhlVUFwW1BebXdbU29dQU5tZVlqalhzTHB4Q0h4Z1dvaU9tWV5fdU9lYE5gVD53RT5nbXlaXV5wRFlscV5kS15feT5lcEFtdHBnY1FpTT5zbGhtZmlhXT5pYGdvZD5uc3BfdD9vSGl1W2ldWmhzQFZ5U09kbm5rYz9yWFhvcXdqO1ZndkFncnhbYHdzUW5qXFw/ZnBhYXFJd0ZAYEJeaj1hXXlRaUtmYFp3a1FeaGo+cHRveHdAW2l4Zj54cWRudVFOaTpJbXc+ZWd2bm53bm9wb1pJcGp2dldxdkdReV9XeGB5eFdgeFNhdztedztWcFBoc0hOb2dmdVN3Z0Rwbk1mc0tfW05Xd1xccHJMaXhoXndxbnNsWVp1eV9WZmp5YHFkUWhad254eW1hcHVnZ3RyVm1DPmlHZ2JTb1xcPFNcXHF3V29CbD1laWlnSVNmZ1NZbFFIPEVTYnllamFSYV9GRGF2a21IX19VT3dpW01lO1lWcXFHVHdiSFFUQmVpWU1iY21CcEFZQHFUZ3djW01leVtVbWVTaUFXUk9pZ01nYG11a21nSnNCaF1SQGdJU1V0dUlJUHFkd0NodW1IT2NyPztGQV1UUllYUU1oRE91ZndGQ3dYW21lO3F3a1FjXFxRWU9rSFU9YmNtQnB2eWxvW3dPbHNwYEZXeEFYZUd5aFtYZXF2aj9eXURfYDpoYXFJY1leYjtPbU5eYFE+WlJZWjpCOjtCandKUllKOmx4Ojw6O2A6PnU9WjpEOjt1PTtCOlNZOjpbPz5qdzprQDw6OmhBPjpKOjs7OnZcXEtIdjtoYTpYaTtxXztoYTpYcUs+aFFIZ1JJZmNOY1JOX1pWaTpGZjs+c1hAaF1IZmJvZG9pZnU+Xjxgd0NGYWZWc3N5cEBPX0Vedzs+aTp5bG14cF1GcV9PXFxAXnc7PmluYlJ5b15Banc8Skk8d1JZWHJIdURFbj5tc11Bclg8alFcXHJdVXVnYExGaUpReUpwQUp2QG53WG08SW5UPFBzPWp3RExqXFxtWXVYYz1WSTxyUUhzWERWbkBUeWhMa2l1c1BVSFB4PGx4dkRKR0huPnB5RVV4V0BLWVFWd3hvWlxcbFdQamBRWXZEWWBMbHdZVFJocFxcWFlFaHBLWGxRZG1eRWp3dEtBQVlbcXFCVHRicVFZVFZTaVlXaFZZSWt3bFReXFxPPERXbT1wRmBPdV1Lc0F4T1lMWDxQQ0F1Y2B0Qm1VRmByUnFvcHhKbl1UU1lWZ1V3cz1YbHFrbUh5O1lUUllMT0hwV2FxTFBZTnV3alBxV0lXeEROeWB2dElUcF5nRV9yTUlsd1hiPFh1QnBtTllnaHhrQT9lcGFsUVFxZXZhSGB3V05bX1dib3d3bj94VFh3QUdlbUd1THlycD5tTVduYnZ1W05jVF53dlZ5R0Bqdz5qdT9kWVBjaVZxcF52c0Rqd0pSOz5aWnc7a0BKOlpzYWh2UHBMeFtsYWVbV2Jzb2hfeVtqR1tSWVo6ZnM/PlpwQVo6YkFCdmdDWld3Xz1TS0dIeXF5aEN0QWt4OmtXO1NWd0dyXk1iQXFDOlNpQTtCZjxaOz5qSWlaQkZaVFddVUlxPV53Oz5pQFBedmBkPUBqSWlaOlZxQT5aZkU6PlI8Wj0+aVhueUlecHJQcD55d1lfZGM+OjxqXlBOYUxOUUVOakQ1Qg==</Image></Text-field></Input></Group><Group><Input><Text-field alignment="centred" layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">Fig. 1</Font>  An approximate equivalent circuit diagram of the <Font executable="false">CT</Font> as viewed from its secondary side.</Text-field></Input></Group><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">Designations used in the diagram:</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="i[p_s]" style="2D Math">NiMmJSJpRzYjJSRwX3NH</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  primary current referred to the secondary,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="v[p_s]" style="2D Math">NiMmJSJ2RzYjJSRwX3NH</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  primary-winding voltage drop referred to the secondary </Font><Font bold="false" italic="false" style="2D Math" underline="false">(represented as an  ac  voltage source),</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="i[ep_s]" style="2D Math">NiMmJSJpRzYjJSVlcF9zRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  primary-winding excitation current referred to the secondary,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="i[cp_s]" style="2D Math">NiMmJSJpRzYjJSVjcF9zRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  primary-winding core-loss current referred to the secondary,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="i[mp_s]" style="2D Math">NiMmJSJpRzYjJSVtcF9zRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  primary-winding magnetizing current referred to the secondary,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="e[p_s]" style="2D Math">NiMmJSJlRzYjJSRwX3NH</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  primary-winding induced electromotive force referred to the secondary </Font><Font bold="false" italic="false" style="2D Math" underline="false">(equal to the electromotive force induced in the secondary winding),</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="i[s]" style="2D Math">NiMmJSJpRzYjJSJzRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  secondary current,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="v[s]" style="2D Math">NiMmJSJ2RzYjJSJzRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  secondary voltage defined as the voltage across the burden,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="R[cp_s]" style="2D Math">NiMmJSJSRzYjJSVjcF9zRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  equivalent core-loss resistance referred to the secondary,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="L[mp_s]" style="2D Math">NiMmJSJMRzYjJSVtcF9zRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  magnetizing inductance referred to the secondary,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="R[es]" style="2D Math">NiMmJSJSRzYjJSNlc0c=</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  equivalent resistance of both windings as viewed from the secondary terminals,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="L[les]" style="2D Math">NiMmJSJMRzYjJSRsZXNH</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  equivalent leakage inductance of both windings as viewed from the secondary terminals,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="R[b]" style="2D Math">NiMmJSJSRzYjJSJiRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  resistance of the burden,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="L[b]" style="2D Math">NiMmJSJMRzYjJSJiRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  inductance of the burden,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="Z[b]" style="2D Math">NiMmJSJaRzYjJSJiRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  impedance of the burden,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">k</Font>  and  <Font executable="false">l</Font><Font encoding="UTF-8">  \342\200\223  secondary terminals.</Font></Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">2.4 Physical and mathematical analysis</Text-field></Title><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  The following analysis uses designations found in the equivalent circuit diagram of the <Font executable="false">CT</Font> of Fig. 1 and other designations that are introduced in the course of the analysis.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">In analyzing the operation of a <Font executable="false">CT</Font>, the well-established theory of <Font executable="false" family="Arial Narrow">on-load</Font> transformer operation is applicable with some specific modifications inherent to the <Font executable="false" family="Arial Narrow">protection</Font> <Font executable="false">CT</Font>.</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 italic="false" underline="false">restart :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The <Font executable="false" family="Arial Narrow">rated transformation ratio</Font> of the <Font executable="false">CT</Font>, <Equation input-equation="a[i][n]" style="2D Math">NiMmJiUiYUc2IyUiaUc2IyUibkc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, is the ratio of the rated primary current, </Font><Equation input-equation="I_p[n]" style="2D Math">NiMmJSRJX3BHNiMlIm5H</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, to the rated secondary current, </Font><Equation input-equation="I_s[n]" style="2D Math">NiMmJSRJX3NHNiMlIm5H</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, or</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">a[i][n] := I_p[n]/I_s[n]  :  'a[i][n]' = a[i][n] ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The <Font executable="false" family="Arial Narrow">turns ratio</Font> of the <Font executable="false">CT</Font>, <Equation input-equation="a[t]" style="2D Math">NiMmJSJhRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, is the ratio of the number of primary turns, </Font><Equation input-equation="N[p]" style="2D Math">NiMmJSJORzYjJSJwRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, to the number of secondary turns, </Font><Equation input-equation="N[s]" style="2D Math">NiMmJSJORzYjJSJzRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, or</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">a[t] := N[p]/N[s]  :  'a[t]' = a[t] ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Throughout this worksheet, the <Font executable="false" family="Arial Narrow">turns ratio</Font> of the <Font executable="false">CT</Font> is used for convenience.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  The turns ratio is equal to the reciprocal of the <Font executable="false" family="Arial Narrow">rated</Font> transformation ratio of the <Font executable="false">CT</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">A current in the line conductor, <Equation input-equation="i[L](t)" style="2D Math">NiMtJiUiaUc2IyUiTEc2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, is the primary current, </Font><Equation input-equation="i[p](t)" style="2D Math">NiMtJiUiaUc2IyUicEc2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, of the CT. The impedance of the primary winding of a through-type CT is negligibly small and does not affect the line current.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">In a bar-primary, ring-type core <Font executable="false">CT</Font>, the leakage flux of both the primary and secondary is negligible and so are both leakage inductances. The primary resistance is also negligibly small. In consequence, with negligible impedance of the primary, the voltage drop across the primary, <Equation input-equation="v[p](t)" style="2D Math">NiMtJiUidkc2IyUicEc2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, is very small. Since the primary voltage is small, the core-loss current, </Font><Equation input-equation="i[cp](t)" style="2D Math">NiMtJiUiaUc2IyUjY3BHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, flowing through the primary is also small and may be omitted.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The excitation current, <Equation input-equation="i[ep](t)" style="2D Math">NiMtJiUiaUc2IyUjZXBHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, drawn by the CT from the ac voltage source is the sum of the core-loss current, </Font><Equation input-equation="i[cp](t)" style="2D Math">NiMtJiUiaUc2IyUjY3BHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, and the magnetizing current, </Font><Equation input-equation="i[mp](t)" style="2D Math">NiMtJiUiaUc2IyUjbXBHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, or</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[ep](t) := i[cp](t) + i[mp](t)  :  'i[ep](t)' = i[ep](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">With the negligible core-loss current, or</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[cp](t) := 0  :  'i[cp](t)' = 0 ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">the excitation current becomes practically the magnetizing current, or</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[ep](t) := eval(i[ep](t))  :  'i[ep](t)' = i[ep](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The primary-winding magnetizing current produces a magnetomotive force <Font executable="false">(mmf)</Font>, <Equation input-equation="f[mp](t)" style="2D Math">NiMtJiUiZkc2IyUjbXBHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, which equals the corresponding current linkage, or</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">f[mp](t) := N[p]*i[mp](t)  :  'f[mp](t)' = f[mp](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">This <Font executable="false">mmf</Font> produces a co-phasal, mutual magnetic flux, <Equation input-equation="phi[m](t)" style="2D Math">NiMtJiUkcGhpRzYjJSJtRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, which is confined to the CT core and links both windings. The flux is given by magnetic Ohm<Font encoding="UTF-8">\342\200\231</Font>s law, viz.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">phi[m](t) := f[mp](t)/R[m]  :  'phi[m](t)' = phi[m](t)  ;  phi[m](t) := 'phi[m](t)' :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">where <Equation input-equation="R[m]" style="2D Math">NiMmJSJSRzYjJSJtRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> is the reluctance of the magnetic core.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  A formula for the reluctance is given and briefly discussed later.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The formula for <Equation input-equation="phi[m](t)" style="2D Math">NiMtJiUkcGhpRzYjJSJtRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> yields equation <Font family="Arial Narrow">Eq[1]</Font> in two unknowns, </Font><Equation input-equation="phi[m](t)" style="2D Math">NiMtJiUkcGhpRzYjJSJtRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and </Font><Equation input-equation="i[mp](t)" style="2D Math">NiMtJiUiaUc2IyUjbXBHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Eq[1] := phi[m](t) - f[mp](t)/R[m] = 0  :  Eq[1]  ;  phi[m](t) := 'phi[m](t)' :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The mutual magnetic flux, <Equation input-equation="phi[m](t)" style="2D Math">NiMtJiUkcGhpRzYjJSJtRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, induces in the secondary an electromotive force (emf), </Font><Equation input-equation="e[s](t)" style="2D Math">NiMtJiUiZUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, given by Faraday<Font encoding="UTF-8">\342\200\231</Font>s law of electromagnetic induction, viz.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s](t) := N[s]*diff(phi[m](t), t)  :  'e[s](t)' = e[s](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">This <Font executable="false">emf</Font> circulates a current, <Equation input-equation="i[s](t)" style="2D Math">NiMtJiUiaUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, in the secondary circuit, which is obtained from Kirchhoff<Font encoding="UTF-8">\342\200\231</Font>s current law (KCL), viz.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[s](t) := i[p_s](t) - i[mp_s](t)  :  'i[s](t)' = i[s](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">where <Equation input-equation="i[p_s](t)" style="2D Math">NiMtJiUiaUc2IyUkcF9zRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and </Font><Equation input-equation="i[mp_s](t)" style="2D Math">NiMtJiUiaUc2IyUlbXBfc0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> are the primary and magnetizing current, respectively, referred (reflected) to the secondary and given by the formulae</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p_s](t) := a[t]*i[p](t)  :  i[mp_s](t) := a[t]*i[mp](t)  :  'i[p_s](t)' = i[p_s](t)  ;  'i[mp_s](t)' = i[mp_s](t)  ;  i[p](t) := 'i[p](t)'  :  i[mp](t) := 'i[mp](t)'  :  i[s](t) := 'i[s](t)' :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The formula obtained from <Font executable="false">KCL</Font> yields equation <Font executable="false" family="Arial Narrow">Eq[2]</Font> in two unknowns, <Equation input-equation="i[mp](t)" style="2D Math">NiMtJiUiaUc2IyUjbXBHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and </Font><Equation input-equation="i[s](t)" style="2D Math">NiMtJiUiaUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Eq[2] := a[t]*i[p](t) - a[t]*i[mp](t) - i[s](t) = 0  :  Eq[2] ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Application of Kirchhoff<Font encoding="UTF-8" executable="false">\342\200\231</Font>s voltage law <Font executable="false">(KVL)</Font> to the secondary circuit gives the formula</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s](t) := R[Ts] * i[s](t) + L[Ts] * diff(i[s](t), t)  :  'e[s](t)' = e[s](t)  ;  e[s](t) := 'e[s](t)' :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">where <Equation input-equation="R[Ts]" style="2D Math">NiMmJSJSRzYjJSNUc0c=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and </Font><Equation input-equation="L[Ts]" style="2D Math">NiMmJSJMRzYjJSNUc0c=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> are the total resistance and total inductance of the secondary circuit, respectively.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The formula for <Equation input-equation="e[s](t)" style="2D Math">NiMtJiUiZUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> yields equation <Font family="Arial Narrow">Eq[3]</Font> in two unknowns, </Font><Equation input-equation="i[mp](t)" style="2D Math">NiMtJiUiaUc2IyUjbXBHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and </Font><Equation input-equation="i[s](t)" style="2D Math">NiMtJiUiaUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Eq[3] := e[s](t) - R[Ts] * i[s](t) - L[Ts] * diff(i[s](t), t) = 0  :  Eq[3] ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Substitution of <Equation input-equation=" e[s](t)" style="2D Math">NiMtJiUiZUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> given by Faraday<Font encoding="UTF-8">\342\200\231</Font>s law of electromagnetic induction yields the desired form of equation <Font family="Arial Narrow">Eq[3]</Font> in two unknowns, </Font><Equation input-equation="phi[m](t)" style="2D Math">NiMtJiUkcGhpRzYjJSJtRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and </Font><Equation input-equation="i[s](t)" style="2D Math">NiMtJiUiaUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, viz.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Eq[3] := subs(e[s](t) = N[s]*diff(phi[m](t), t), lhs(Eq[3])) = 0  :  Eq[3] ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">3. Solution of Simultaneous Equations for Time-varying Parameters</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">The system of the three equations, <Font executable="false" family="Arial Narrow">Eq[1]</Font>, <Font executable="false" family="Arial Narrow">Eq[2]</Font>, and <Font executable="false" family="Arial Narrow">Eq[3]</Font>, will be solved symbolically for the three unknowns, <Equation input-equation="i[mp](t)" style="2D Math">NiMtJiUiaUc2IyUjbXBHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, </Font><Equation input-equation="phi[m](t)" style="2D Math">NiMtJiUkcGhpRzYjJSJtRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, and </Font><Equation input-equation="i[s](t)" style="2D Math">NiMtJiUiaUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> with the primary current, </Font><Equation input-equation="i[p](t)" style="2D Math">NiMtJiUiaUc2IyUicEc2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, given as</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p](t) := I_p_mx * sin(omega*t + theta[p])  :  'i[p](t)' = i[p](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">where:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="I_p_mx" style="2D Math">NiMlJ0lfcF9teEc=</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  </Font><Font bold="false" family="Arial Narrow" italic="false" style="2D Math" underline="false">amplitude</Font><Font bold="false" italic="false" style="2D Math" underline="false"> of the primary current,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="omega" style="2D Math">NiMlJm9tZWdhRw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  </Font><Font bold="false" family="Arial Narrow" italic="false" style="2D Math" underline="false">angular frequency</Font><Font bold="false" italic="false" style="2D Math" underline="false"> of the system voltage,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="theta[p]" style="2D Math">NiMmJSZ0aGV0YUc2IyUicEc=</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  initial phase angle, or phase shift of the primary current with respect to the system voltage.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  It is convenient to adopt the primary current as a <Font executable="false" family="Arial Narrow">reference</Font> waveform with <Font executable="false">zero</Font> initial phase angle. However, all the computations found in this worksheet are valid for a general case when the primary current has an arbitrary <Font executable="false">non-zero</Font> initial phase angle.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Execution of the following command line yields the solution of the three simultaneous equations in the form of an unordered solution set.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">sln := dsolve({Eq[1], Eq[2], Eq[3]}, {i[mp](t), phi[m](t), i[s](t)}) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The formulae for the magnetizing current <Equation input-equation="i[mp](t)" style="2D Math">NiMtJiUiaUc2IyUjbXBHNiMlInRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, mutual magnetic flux </Font><Equation input-equation="phi[m](t)" style="2D Math">NiMtJiUkcGhpRzYjJSJtRzYjJSJ0Rw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, and secondary current </Font><Equation input-equation="i[s](t)" style="2D Math">NiMtJiUiaUc2IyUic0c2IyUidEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, extracted from the set of solutions are as follows.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">assign(sln) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  The magnetizing current flowing through the primary:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[mp](t) := simplify(subs(_C1=0, i[mp](t)))  :  'i[mp](t)' = i[mp](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">A more transparent form of the above formula is</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[mp](t) := collect(collect(i[mp](t), sin(omega*t + theta[p])), cos(omega*t + theta[p])) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'i[mp](t)' = i[mp](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  The mutual magnetic flux in the core:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">phi[m](t) := simplify(subs(_C1=0, phi[m](t)))  :  'phi[m](t)' = phi[m](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">A more transparent form of the above formula is</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">phi[m](t) := collect(collect(phi[m](t), sin(omega*t + theta[p])), cos(omega*t + theta[p])) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'phi[m](t)' = phi[m](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  The secondary current:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[s](t) := simplify(subs(_C1=0, i[s](t)))  :  'i[s](t)' = i[s](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">A more transparent form of the above formula is</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[s](t) := collect(collect(i[s](t), sin(omega*t + theta[p])), cos(omega*t + theta[p])) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'i[s](t)' = i[s](t) ;</Font></Text-field></Input></Group><Text-field layout="Text Output" style="Text Output"/></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">4. Derivation of Formulae for Supplementary Time-varying Parameters</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">Using the above solutions, formulae for three supplementary electrical quantities are obtained as follows.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  The secondary <Font executable="false">emf</Font><Font encoding="UTF-8">  \342\200\223  obtainable from the formula</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred" bullet="none" firstindent="0.0" leftmargin="0.0" linebreak="space" linespacing="0.0" rightmargin="0.0" spaceabove="0.0" spacebelow="0.0"><Equation input-equation="e[s](t) = N[s]*diff(phi[m](t),t)" style="2D Math">NiMvLSYlImVHNiMlInNHNiMlInRHKiYmJSJOR0YnIiIiLSUlZGlmZkc2JC0mJSRwaGlHNiMlIm1HRilGKkYu</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s](t) := N[s]*diff(phi[m](t), t)  :  'e[s](t)' = e[s](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">A more transparent form of the above formula is</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s](t) := collect(collect(simplify(e[s](t)), cos(omega*t + theta[p])), sin(omega*t + theta[p])) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'e[s](t)' = e[s](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  The primary <Font executable="false">emf</Font><Font encoding="UTF-8">  \342\200\223  obtainable from the formula</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="e[p](t) = a[t]*e[s](t)" style="2D Math">NiMvLSYlImVHNiMlInBHNiMlInRHKiYmJSJhR0YpIiIiLSZGJjYjJSJzR0YpRi4=</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[p](t) := a[t]*e[s](t)  :  'e[p](t)' = e[p](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">A more transparent form of the above formula is</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[p](t) := collect(collect(simplify(e[p](t)), cos(omega*t + theta[p])), sin(omega*t + theta[p])) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'e[p](t)' = e[p](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font><Font encoding="UTF-8">  The secondary voltage defined as the voltage across the burden  \342\200\223  obtainable from </Font><Font executable="false">KVL</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="v[s](t) = R[b]*i[s](t)+L[b]*diff(i[s](t),t)" style="2D Math">NiMvLSYlInZHNiMlInNHNiMlInRHLCYqJiYlIlJHNiMlImJHIiIiLSYlImlHRidGKUYxRjEqJiYlIkxHRi9GMS0lJWRpZmZHNiRGMkYqRjFGMQ==</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">v[s](t) := R[b] * i[s](t) + L[b] * diff(i[s](t), t)  :  'v[s](t)' = v[s](t) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">A more transparent form of the above formula is</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">v[s](t) := collect(collect(simplify(v[s](t)), sin(omega*t + theta[p])), cos(omega*t + theta[p])) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'v[s](t)' = v[s](t) ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">5. Transformation of Formulae for Time-varying Parameters from Time Domain to Angular Domain</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">It is convenient to replace the variable <Equation input-equation="omega*t" style="2D Math">NiMqJiUmb21lZ2FHIiIiJSJ0R0Yl</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> with </Font><Equation input-equation="x" style="2D Math">NiMlInhH</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> by putting</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">omega*t = x ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">whence</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">t = x/omega ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Therefore, the relevant formulae obtained earlier in the <Font executable="false" family="Arial Narrow">time domain</Font> assume the following forms in the new domain, which may be called the <Font executable="false" family="Arial Narrow">angular domain</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p](x) := subs(t=x/omega, i[p](t))  :  i_p_f := i[p](x) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[mp](x) := subs(t=x/omega, i[mp](t))  :  i_mp_f := i[mp](x) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">phi[m](x) := subs(t=x/omega, phi[m](t))  :  phi_m_f := phi[m](x) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s](x) := subs(t=x/omega, e[s](t))  :  e_s_f := e[s](x) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[s](x) := subs(t=x/omega, i[s](t))  :  i_s_f := i[s](x) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[p](x) := subs(t=x/omega, e[p](t))  :  e_p_f := e[p](x) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">v[s](x) := subs(t=x/omega, v[s](t))  :  v_s_f := v[s](x) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">In the above formulae, there are several quantities that must be defined at this point of computations. The respective defining formulae follow.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Amplitude of the primary current</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_p_mx := sqrt(2)*I_p  :  I_p[mx] = I_p_mx ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">where <Font executable="false">I_p</Font> is the effective value of the current defined as the root-mean-square, or <Font executable="false">rms</Font> value.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Reluctance of the <Font executable="false">CT</Font> magnetic circuit</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">R[m] := l[av]/(mu[0]*mu[r]*A)  :  'R[m]' = R[m] ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">where:</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="l[av]" style="2D Math">NiMmJSJsRzYjJSNhdkc=</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  mean length of the annular magnetic circuit with mean diameter </Font><Equation input-equation="d" style="2D Math">NiMlImRH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, given by</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">l[av] := Pi*d  :  'l[av]' = l[av] ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">A</Font><Font encoding="UTF-8">  \342\200\223  active cross-sectional area of the magnetic circuit,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="mu[0]" style="2D Math">NiMmJSNtdUc2IyIiIQ==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  absolute permeability given as </Font><Equation input-equation="mu[0] = 4*Pi*10^`-7`" style="2D Math">NiMvJiUjbXVHNiMiIiEqKCIiJSIiIiUjUGlHRiopIiM1JSMtN0dGKg==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> H/m,</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="mu[r]" style="2D Math">NiMmJSNtdUc2IyUickc=</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  relative permeability of the magnetic-circuit material.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Heading 2" style="Heading 2">Notes:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">1. An <Font executable="false" family="Arial Narrow">active</Font> cross-sectional area of the magnetic core is the <Font executable="false" family="Arial Narrow">net</Font> cross-sectional area of the magnetic material. The gross, or geometric cross-sectional area, <Equation input-equation="A[g]" style="2D Math">NiMmJSJBRzYjJSJnRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, of the core is reduced by the lamination stacking factor, </Font><Equation input-equation="k[s]" style="2D Math">NiMmJSJrRzYjJSJzRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, which depends on the <Font family="Arial Narrow">thickness</Font> of the inter-laminar <Font family="Arial Narrow">insulation</Font> and the <Font family="Arial Narrow">tightness</Font> of the core clamping. Thus, the active cross-sectional area may be expressed by the formula</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="A = k[s]*A[g]" style="2D Math">NiMvJSJBRyomJiUia0c2IyUic0ciIiImRiQ2IyUiZ0dGKg==</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">with <Equation input-equation="k[s]" style="2D Math">NiMmJSJrRzYjJSJzRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> being typically between 0.9 and 0.95.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">2. In general, the relative permeability is <Font executable="false" family="Arial Narrow">not</Font> constant and depends on the magnetic flux density, <Font executable="false">B</Font>, in the core. If the <Font executable="false" family="Arial Narrow">hysteresis</Font> phenomenon in the magnetic material is neglected, the relative permeability may be expressed as</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="mu[r](B) = [1/mu[0]]*diff(B(H),H)" style="2D Math">NiMvLSYlI211RzYjJSJyRzYjJSJCRyomNyMqJiIiIkYuJkYmNiMiIiEhIiJGLi0lJWRpZmZHNiQtRio2IyUiSEdGOEYu</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The adoption of a constant value of <Equation input-equation="mu[r]" style="2D Math">NiMmJSNtdUc2IyUickc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> is equivalent to the adoption of the average value of the slope of the B-H curve in the range </Font><Equation input-equation="[-H[max], H[max]]" style="2D Math">NiM3JCwkJiUiSEc2IyUkbWF4RyEiIkYl</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, which is usually expressed as follows</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="mu[r][av] = B[max]/(mu[0]*H[max])" style="2D Math">NiMvJiYlI211RzYjJSJyRzYjJSNhdkcqJiYlIkJHNiMlJG1heEciIiIqJiZGJjYjIiIhRjAmJSJIR0YuRjAhIiI=</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  Angular frequency of the system voltage</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">omega := 2*Pi*f[n]  :  'omega' = omega ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">where <Equation input-equation="f[n]" style="2D Math">NiMmJSJmRzYjJSJuRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> is the rated frequency of the system voltage.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">6. Specification of Input Data</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2"><Font encoding="UTF-8">1. Entering the data names with assigned numerical values that are expressed in units belonging to the Syst\303\250me International d\342\200\231Unit\303\251s (SI) for use in computations:</Font></Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">N[p] := 1  :  N[s]:=160  :  f[n] := 50  :  I_p := 800  :  A := 100*10^(-6)  :  d :=60*10^(-3)  :  R[b] := 0.05  :  R[Ts] := 0.07  :  L[b] := 0.15*10^(-3)  :  L[Ts] := 0.2*10^(-3)  :  mu[r] := 5000  :  mu[0] := 4*Pi*10^(-7) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal"><Font background="[0,0,0]" bold="true" family="Times New Roman" italic="false" size="12" underline="false">IMPORTANT:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The next parameter to be entered is the <Font executable="false" family="Arial Narrow">initial phase angle</Font> of the primary current. This parameter is <Font executable="false" family="Arial Narrow">critical</Font> for the numerical computations to follow. Read first the <Font executable="false" family="Arial Narrow">NOTE</Font> before entering its value.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal"><Font background="[0,0,0]" bold="true" family="Arial Narrow" foreground="[0,0,0]" italic="false" size="12" underline="false">NOTE:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">If the primary current is to be referenced to the system voltage, the initial phase angle <Equation input-equation="theta[p]" style="2D Math">NiMmJSZ0aGV0YUc2IyUicEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> will be different than zero. Its value must be given and entered in <Font family="Arial Narrow">electrical degrees</Font> as</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  a <Font executable="false" family="Arial Narrow">negative</Font><Font encoding="UTF-8"> number  \342\200\223  for current </Font><Font executable="false" family="Arial Narrow">lagging</Font> the voltage,</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  a <Font executable="false" family="Arial Narrow">positive</Font><Font encoding="UTF-8"> number  \342\200\223  for current </Font><Font executable="false" family="Arial Narrow">leading</Font> the voltage.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font style="Text">The initial phase angle of the primary current:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[p] := 0 :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[p[deg]] := theta[p] :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[p[rad]] := Pi*theta[p]/180  :  theta[p] := theta[p[rad]] :</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">2. Displaying numerical values of the data in the 'engineering form' together with practical units of measure:</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'N[p]' = N[p]  ;  'N[s]' = N[s]  ;  'f[n]' = f[n]*Hz  ;  'I_p' = I_p*'A'  ;  'A' = A*10^6*mm^2  ;  'd' = d*10^3*mm  ;  'R[b]' = R[b]*Omega  ;  'R[Ts]' = R[Ts]*Omega  ;  'L[b]' = evalf(L[b]*10^3, 2)*mH  ;  'L[Ts]' = evalf(L[Ts]*10^3, 2)*mH  ;  'mu[r]' = mu[r]  ;  'mu[0]' = mu[0]*`H/m` ;<Font encoding="ISO8859-1">
if abs(frac(theta[p[deg]])) = 0 then print('theta[p]' = theta[p[deg]] * `\260`) else print('theta[p]' = evalf(theta[p[deg]], 3) * `\260`) end if :</Font></Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">7. Plots of Waveforms of Time-varying Parameters</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">According to the common engineering practice, the horizontal-axis scale will be expressed in <Font executable="false" family="Arial Narrow">electrical degrees</Font>. The preceding functions represent time-varying quantities in the <Font executable="false" family="Arial Narrow">angular domain</Font>, which uses the <Font executable="false" family="Arial Narrow">radian scale</Font>. To obtain plots of these functions against the variable <Equation input-equation="x" style="2D Math">NiMlInhH</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> expressed in <Font family="Arial Narrow">degrees</Font>, this variable must be divided by 180/</Font><Equation input-equation="Pi" style="2D Math">NiMlI1BpRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">. Thus, a new independent variable is introduced for plotting purposes <Font family="Arial Narrow">only</Font>, which is labeled </Font><Equation input-equation="xn" style="2D Math">NiMlI3huRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and defined as follows:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">xn := x/(180/Pi)  :  'xn' = xn ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The adoption of <Font executable="false" family="Arial Narrow">electrical degrees</Font> for the <Equation input-equation="x" style="2D Math">NiMlInhH</Equation><Font bold="false" italic="false" style="2D Math" underline="false">-axis scale requires re-defining the functions for plotting purposes as follows:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p](x) := subs(x=xn, i[p](x)) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[mp](x) := subs(x=xn, i[mp](x)) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">phi[m](x) := subs(x=xn, phi[m](x)) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s](x) := subs(x=xn, e[s](x)) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[s](x) := subs(x=xn, i[s](x)) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[p](x) := subs(x=xn, e[p](x)) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">A procedure for plotting waveforms of the time-varying parameters along with the necessary designations of the axes is as follows.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">with(plots, display, textplot) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">plot_i[p] := plot(i[p](x), x = 0..660, y = -1300..2200, thickness = 2, color = red, labels = [``,``]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">plot_i[mp] := plot(i[mp](x)*10^2, x = 0..660, color = blue) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">plot_phi[m] := plot(phi[m](x)*10^7, x = 0..660, color = green) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">plot_e[s] := plot(e[s](x)*10^3, x = 0..660, color = cyan) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">plot_i[s] := plot(i[s](x)*10^2, x = 0..660, color = brown) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">plot_e[p] := plot(e[p](x)*10^5, x = 0..660, color = magenta) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">plot_hax := plot(0, x = 0..720, color = black) :</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">ha_s := textplot([648, 200, `x (\260 el.)`], align = RIGHT, font = [HELVETICA, 11]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">s_i := textplot([15, 2100, `i_p (A), i_s (A/100), i_mp (A/100)`], align = RIGHT, font = [HELVETICA, 10]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">s_e := textplot([15, 1800, `e_s (mV), e_p (mV/100)`], align = RIGHT, font = [HELVETICA, 10]):</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">s_phi_s := textplot([15,1500, `f`], align = RIGHT, font = [SYMBOL,12]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">s_phi_u1 := textplot([30, 1490, `_m (`], align = RIGHT, font = [HELVETICA, 10]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">s_phi_u2 := textplot([71, 1490, `m`], align = RIGHT, font = [SYMBOL,11]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">s_phi_u3 := textplot([85, 1490, `Wb/10)`], align = RIGHT, font = [HELVETICA, 10]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">display([ plot_i[p], plot_i[mp], plot_phi[m], plot_e[s], plot_i[s], plot_e[p], </Font>plot_hax,<Font italic="false" underline="false"> ha_s, s_i, s_e, s_phi_s, s_phi_u1, s_phi_u2, s_phi_u3 ], title=`WAVEFORMS OF MAGNETIC FLUX AND ELECTRICAL QUANTITIES`, titlefont=[HELVETICA, BOLD, 12]) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">Fig. 2</Font>  Waveforms of magnetic flux and electrical quantities characteristic of the <Font executable="false">CT</Font>, plotted as functions of the variable <Equation input-equation="x=omega*t" style="2D Math">NiMvJSJ4RyomJSZvbWVnYUciIiIlInRHRic=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> expressed in <Font family="Arial Narrow">electrical degrees</Font>:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  primary current, <Font executable="false">i_p</Font>, <Font executable="false">(</Font>thick <Font executable="false" family="Arial Narrow">red</Font> curve<Font executable="false">)</Font>;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  magnetizing current flowing through the primary, <Font executable="false">i_mp</Font>, multiplied by <Equation input-equation="10^2" style="2D Math">NiMqJCIjNSIiIw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> (<Font family="Arial Narrow">blue</Font> curve);</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  mutual magnetic flux in the core, <Equation input-equation="phi" style="2D Math">NiMlJHBoaUc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">_m, multiplied by </Font><Equation input-equation="10^7" style="2D Math">NiMqJCIjNSIiKA==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> (<Font family="Arial Narrow">green</Font> curve);</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  secondary <Font executable="false">emf</Font>, <Font executable="false">e_s</Font>, multiplied by <Equation input-equation="10^3" style="2D Math">NiMqJCIjNSIiJA==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> (<Font family="Arial Narrow">cyan</Font> curve);</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  secondary current, <Font executable="false">i_s</Font>, multiplied by <Equation input-equation="10^2" style="2D Math">NiMqJCIjNSIiIw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> (<Font family="Arial Narrow">brown</Font> curve);</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  primary <Font executable="false">emf</Font>, <Font executable="false">e_p</Font>, multiplied by <Equation input-equation="10^5" style="2D Math">NiMqJCIjNSIiJg==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> (<Font family="Arial Narrow">magenta</Font> curve).</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal"><Font background="[0,0,0]" executable="false" family="Arial Narrow" foreground="[0,0,0]" italic="false" size="12">Observations from Fig. 2<Font background="[0,0,0]" executable="false" family="Arial Narrow" foreground="[0,0,0]" italic="false" size="12" underline="false">:</Font></Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(1)</Font> The primary and secondary currents seem to be <Font executable="false" family="Arial Narrow">co-phasal</Font>. In fact, pertinent computation provided later reveals that the secondary current <Font executable="false" family="Arial Narrow">leads</Font> the primary current <Font executable="false">(</Font>referred to the secondary<Font executable="false">)</Font> by a very small angle called the <Font executable="false" family="Arial Narrow">phase error</Font> of the <Font executable="false">CT</Font>. <Font executable="false">(</Font>Refer also to Fig. 3.<Font executable="false">)</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(2)</Font> The magnetizing current and magnetic flux are <Font executable="false" family="Arial Narrow">co-phasal</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(3)</Font> The primary and secondary <Font executable="false">emf</Font>s are <Font executable="false" family="Arial Narrow">co-phasal</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(4)</Font> The magnetizing current is <Font executable="false" family="Arial Narrow">in quadrature</Font> with the primary-induced <Font executable="false">emf</Font>, i.e., lags the <Font executable="false">emf</Font> by <Font encoding="ISO8859-1" executable="false">90\272</Font>. <Font executable="false">(</Font>Refer also to pertinent computation of the <Font executable="false" family="Arial Narrow">phase difference</Font> between the primary <Equation input-equation="emf" style="2D Math">NiMlJGVtZkc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and magnetizing current.)</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(5)</Font> The magnetizing current <Font executable="false" family="Arial Narrow">lags</Font> the primary current. <Font executable="false">(</Font>Refer also to pertinent computation of the <Font executable="false" family="Arial Narrow">phase difference</Font> between the primary current and magnetizing current.<Font executable="false">)</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false">(6)</Font> The secondary current <Font executable="false" family="Arial Narrow">lags</Font> the secondary-induced <Font executable="false">emf</Font>. <Font executable="false">(</Font>Refer also to pertinent computation of the <Font executable="false" family="Arial Narrow">phase difference</Font> between the secondary current and secondary <Font executable="false">emf</Font>.<Font executable="false">)</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">If parts of waveforms of the primary and secondary current are plotted within a very narrow region where they intersect the horizontal axis, the phase lead of the secondary current can be clearly seen. To this end, a region of length of <Font encoding="ISO8859-1" executable="false">2\272</Font> is chosen.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">To ensure selection of a proper region of the same length irrespective of the initial phase angle of the primary current, a reference abscissa, <Equation input-equation="x[r]" style="2D Math">NiMmJSJ4RzYjJSJyRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, is defined. It corresponds to the zero-crossing of the primary-current waveform, or</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">x[r] := 180 - signum(theta[p[deg]]) * abs(theta[p[deg]]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">if abs(frac(theta[p[deg]])) = 0 then print('x[r]'=x[r] * `\260`) else print('x[r]'=evalf(x[r], 4) * `\260`) end if :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Three more ancillary abscissae are defined for the same purpose, viz.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">x[b] := x[r] - 1  :  x[e] := x[r] + 0.7  :  x[hax] := x[r] + 0.99 :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The plotting procedure proper follows.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">plot_i[p] := plot(i[p](x), x = x[b]..x[e], y = -14..24, color = red, labels = [``,``]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">plot_i[s] := plot(i[s](x)*10^2, x = x[b]..x[e], color = blue) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input">plot_hax := plot(0, x = x[b]..x[hax], color = black) :</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">ha_s := textplot([x[hax] - 0.6, 2, `x (\260 el.)`], align = RIGHT, font = [HELVETICA, 11]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">s_i := textplot([x[b]+0.04, 23, `i_p (A), i_s (A/100)`], align = RIGHT, font = [HELVETICA, 10]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">display([ plot_i[p], plot_i[s], ha_s, s_i ], title=`PARTS OF WAVEFORMS OF PRIMARY AND SECONDARY CURRENT`, titlefont=[HELVETICA, BOLD, 12]);</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">Fig. 3</Font>  Parts of waveforms of primary and secondary current to show the phase lead of the secondary current:</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  primary current, <Font executable="false">i_p</Font>, <Font executable="false">(<Font family="Arial Narrow">red</Font></Font> curve<Font executable="false">)</Font>;</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  secondary current, <Font executable="false">i_s</Font>, multiplied by <Equation input-equation="10^2" style="2D Math">NiMqJCIjNSIiIw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> (<Font family="Arial Narrow">blue</Font> curve).</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Before proceeding to the next section, all the re-defined functions must be 'forgotten' and their original symbolic forms in the <Font executable="false" family="Arial Narrow">radian</Font> <Font executable="false" family="Arial Narrow">angular domain</Font> restored.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">To unassign the new definitions from the names of the functions re-defined above only for plotting purposes, each function name is enclosed in single forward quotes:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p](x) := 'i[p](x)'  :  i[mp](x) := 'i[mp](x)'  :  phi[m](x) := 'phi[m](x)'  :  e[s](x) := 'e[s](x)'  :  i[s](x) := 'i[s](x)'  :  e[p](x) := 'e[p](x)' :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">8. Numerical Evaluation of Formulae for Time-varying Parameters</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">It should be noticed that the time-varying quantities obtained earlier have the following forms in the <Font executable="false" family="Arial Narrow">angular domain</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  magnetizing current <Equation input-equation="i[mp]" style="2D Math">NiMmJSJpRzYjJSNtcEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, magnetic flux </Font><Equation input-equation="phi[m]" style="2D Math">NiMmJSRwaGlHNiMlIm1H</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, and secondary current </Font><Equation input-equation="i[s]" style="2D Math">NiMmJSJpRzYjJSJzRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="quantity[1](x) = a[1]*sin(x+theta[p])+a[2]*cos(x+theta[p])" style="2D Math">NiMvLSYlKXF1YW50aXR5RzYjIiIiNiMlInhHLCYqJiYlImFHRidGKC0lJHNpbkc2IywmRipGKCYlJnRoZXRhRzYjJSJwR0YoRihGKComJkYuNiMiIiNGKC0lJGNvc0dGMUYoRig=</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  secondary <Font executable="false">emf</Font> <Equation input-equation="e[s]" style="2D Math">NiMmJSJlRzYjJSJzRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false"> and secondary voltage </Font><Equation input-equation="v[s]" style="2D Math">NiMmJSJ2RzYjJSJzRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">:</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="quantity[2](x) = a[1]*cos(x+theta[p])+a[2]*sin(x+theta[p])" style="2D Math">NiMvLSYlKXF1YW50aXR5RzYjIiIjNiMlInhHLCYqJiYlImFHNiMiIiJGMC0lJGNvc0c2IywmRipGMCYlJnRoZXRhRzYjJSJwR0YwRjBGMComJkYuRidGMC0lJHNpbkdGM0YwRjA=</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Either form may be transformed to an equivalent single-trigonometric-function form, which is more suitable for the analysis of such quantities. The equivalent forms are, respectively,</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="quantity[1](x) = A*sin(x+theta[p]+alpha[1])" style="2D Math">NiMvLSYlKXF1YW50aXR5RzYjIiIiNiMlInhHKiYlIkFHRigtJSRzaW5HNiMsKEYqRigmJSZ0aGV0YUc2IyUicEdGKCYlJmFscGhhR0YnRihGKA==</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="quantity[2](x) = A*sin(x+theta[p]+alpha[2])" style="2D Math">NiMvLSYlKXF1YW50aXR5RzYjIiIjNiMlInhHKiYlIkFHIiIiLSUkc2luRzYjLChGKkYtJiUmdGhldGFHNiMlInBHRi0mJSZhbHBoYUdGJ0YtRi0=</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">where the function parameters denote:</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="A" style="2D Math">NiMlIkFH</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  </Font><Font bold="false" family="Arial Narrow" italic="false" style="2D Math" underline="false">amplitude</Font><Font bold="false" italic="false" style="2D Math" underline="false">, which is determined from the formula</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="A = sqrt(a[1]^2+a[2]^2)" style="2D Math">NiMvJSJBRy0lJXNxcnRHNiMsJiokJiUiYUc2IyIiIiIiI0YtKiQmRis2I0YuRi5GLQ==</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Equation input-equation="alpha[1]" style="2D Math">NiMmJSZhbHBoYUc2IyIiIg==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, </Font><Equation input-equation="alpha[2]" style="2D Math">NiMmJSZhbHBoYUc2IyIiIw==</Equation><Font bold="false" encoding="UTF-8" italic="false" style="2D Math" underline="false">  \342\200\223  initial </Font><Font bold="false" family="Arial Narrow" italic="false" style="2D Math" underline="false">phase angles</Font><Font bold="false" italic="false" style="2D Math" underline="false"> in <Font family="Arial Narrow">radians</Font>, which are determined from the formulae</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="alpha[1] = sign(a[2])*arctan(abs(a[2]/a[1]))" style="2D Math">NiMvJiUmYWxwaGFHNiMiIiIqJi0lJXNpZ25HNiMmJSJhRzYjIiIjRictJSdhcmN0YW5HNiMtJSRhYnNHNiMqJkYsRicmRi1GJiEiIkYn</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field alignment="centred"><Equation input-equation="alpha[2] = sign(a[2])*arctan(abs(a[1]/a[2]))" style="2D Math">NiMvJiUmYWxwaGFHNiMiIiMqJi0lJXNpZ25HNiMmJSJhR0YmIiIiLSUnYXJjdGFuRzYjLSUkYWJzRzYjKiYmRi02I0YuRi5GLCEiIkYu</Equation></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">All the time-varying parameters of the <Font executable="false">CT</Font> will be computed by numerical evaluation of their respective symbolic formulae. Where appropriate, they will be transformed to the form containing the <Font executable="false" family="Arial Narrow">sin</Font> function. The transformed form will be given with the phase angle in electrical <Font executable="false" family="Arial Narrow">radians</Font> and <Font executable="false" family="Arial Narrow">degrees</Font>. The latter form will be given with rounded-off amplitude and phase angle, and with a practical unit of measure, according to the commonly adopted engineering practice.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  In all the formulae that follow, the variable <Font executable="false">x</Font> is in <Font executable="false" family="Arial Narrow">radians</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">1. The primary current</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  With <Font executable="false" family="Arial Narrow">negative</Font> values of the initial phase angle, <Equation input-equation="theta[p]" style="2D Math">NiMmJSZ0aGV0YUc2IyUicEc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, of the primary current, the latter is returned in the form containing the <Font family="Arial Narrow">cos</Font> function. This form poses a slight inconvenience in the further analysis. It is, therefore, transformed into an equivalent form containing the <Font family="Arial Narrow">sin</Font> function. A programming construct <Font family="Arial Narrow">if ... then ... end if</Font> is used to control the final form of the equation for the primary current.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p_v](x) := evalf(eval(i_p_f)) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">if signum(theta[p[deg]]) = -1 then i[p_v](x) := -op(1, i[p_v](x)) * sin(x + evalf(theta[p])) end if :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">In effect, irrespective of the value <Font executable="false">(</Font>negative, <Font executable="false">zero</Font>, or positive<Font executable="false">)</Font> of the initial phase angle, the primary current is input and displayed with the sinusoidal term, viz.</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p](x) = i[p_v](x) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_p_mx := op(1, i[p_v](x)) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[i_p] := evalf(theta[p]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][i_p] := evalf(180*theta[i_p]/Pi)  :  theta[deg][i_p_r] := evalf(theta[deg][i_p], 3) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">i[p](x) = evalf(I_p_mx, 5) * sin(180*x/Pi + theta[deg][i_p_r] * `\260`) * 'A' ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">2. The primary current referred (reflected) to the secondary</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p_s](x) = 'a[t]' * i[p](x)  ;  i[p_s_v](x) := a[t]*i[p_v](x)  :  i[p_s](x) = i[p_s_v](x) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_p_s_mx := evalf(a[t]*I_p_mx) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][i_p_s] := theta[deg][i_p]  :  theta[deg][i_p_s_r] := evalf(theta[deg][i_p_s], 3) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">i[p_s](x) = evalf(I_p_s_mx, 3) * sin(180*x/Pi + theta[deg][i_p_s_r] * `\260`) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the primary current referred to the secondary:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_p_s := evalf(I_p_s_mx/sqrt(2))  :  'I_p_s' = evalf(I_p_s, 2) * 'A' ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">3. The magnetizing current flowing through the primary</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  original form:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[mp_v](x) := evalf(eval(i_mp_f))  :  i[mp](x) = i[mp_v](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  transformed forms:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">a[1] := op(1, op(1, i[mp_v](x)))  :  a[2] := op(1, op(2, i[mp_v](x))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_mp_mx := sqrt(a[1]^2 + a[2]^2) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">alpha[1] := sign(a[2]) * evalf(arctan(abs(a[2]/a[1]))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[i_mp] := evalf(theta[p] + alpha[1]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][i_mp] := evalf(180*theta[i_mp]/Pi) : theta[deg][i_mp_r] := evalf(theta[deg][i_mp], 3) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[mp_v](x) := I_mp_mx*sin(x + theta[i_mp])  :  i[mp](x) = i[mp_v](x) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">i[mp](x) = evalf(I_mp_mx, 3) * sin(180*x/Pi + theta[deg][i_mp_r] * `\260`) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The angular displacement of the magnetizing current with respect to the primary current:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">delta[1] := theta[deg][i_mp] - theta[deg][i_p] :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">'theta[deg][i_mp] - theta[deg][i_p]' = evalf(delta[1], 3) * `\260` ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  This result implies that the magnetizing current <Font executable="false" family="Arial Narrow">lags</Font> the primary current.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The actual effective value of the magnetizing current flowing through the primary:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_mp := evalf(I_mp_mx/sqrt(2))  :  'I_mp' = evalf(I_mp, 2) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The percentage <Font executable="false" family="Arial Narrow">current error</Font> of the <Font executable="false">CT</Font> at the adopted magnitude of the primary current, defined as the effective value of the magnetizing current expressed as a percentage of the primary current:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">epsilon[i] = '100*I_mp/I_p'  ;  epsilon[i] := 100*I_mp/I_p  :  epsilon[i] := evalf(epsilon[i], 2)  :  'epsilon[i]' = epsilon[i] * convert(convert(convert([37], 'bytes'), string), symbol) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  As expected and desired, the percentage <Font executable="false" family="Arial Narrow">current error</Font> of the <Font executable="false">CT</Font> at the adopted magnitude of the primary current is <Font executable="false" family="Arial Narrow">very</Font> low.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the magnetic field strength in the core:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">H = 'I_mp*N[p]/l[av]'  ;  H := I_mp*N[p]/l[av]  :  'H' = evalf(H, 3) * `A/m`  ;  H_v := H :</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">4. The primary-winding magnetizing current referred (reflected) to the secondary</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[mp_s](x) = 'a[t]' * i[mp](x)  ;  i[mp_s_v](x) := a[t]*i[mp_v](x)  :  i[mp_s](x) = i[mp_s_v](x) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_mp_s_mx := evalf(a[t]*I_mp_mx) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][i_mp_s] := theta[deg][i_mp]  :  theta[deg][i_mp_s_r] := evalf(theta[deg][i_mp_s], 3) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">i[mp_s](x) = evalf(I_mp_s_mx*10^3, 3) * sin(180*x/Pi + theta[deg][i_mp_s_r] * `\260`) * mA ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the magnetizing current referred to the secondary:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_mp_s := evalf(I_mp_s_mx/sqrt(2))  :  'I_mp_s' = evalf(I_mp_s*10^3, 3) * mA ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">5. The mutual magnetic flux in the core</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  original form:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">phi[m_v](x) := evalf(eval(phi_m_f))  :  phi[m](x) = phi[m_v](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  transformed forms:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">a[1] := op(1, op(1, phi[m_v](x)))  :  a[2] := op(1, op(2, phi[m_v](x))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Phi_m_mx := sqrt(a[1]^2 + a[2]^2) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">alpha[1] := sign(a[2]) * evalf(arctan(abs(a[2]/a[1]))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[phi_m] := evalf(theta[p] + alpha[1])  :  theta[deg][phi_m] := evalf(180*theta[phi_m]/Pi) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][phi_m_r] := evalf(theta[deg][phi_m], 3) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">phi[m_v](x) := Phi_m_mx*sin(x + theta[phi_m])  :  phi[m](x) = phi[m_v](x) * Wb ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">phi[m](x) = evalf(Phi_m_mx*10^6, 3) * sin(180*x/Pi + theta[deg][phi_m_r] * `\260`) * mu*Wb;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The maximum, or peak value of the magnetic flux in the core:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Phi_m_mx := Phi_m_mx  :  Phi[m][mx] = evalf(Phi_m_mx * 10^6, 3) * mu*Wb ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the magnetic flux in the core:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Phi[m] := evalf(Phi_m_mx/sqrt(2))  :  'Phi[m]' = evalf(Phi[m] * 10^6, 2) * mu*Wb ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The peak value of the magnetic flux density in the core:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">B[mx] = 'Phi[m][mx]/A'  ;  B[mx] := Phi_m_mx/A  :  'B[mx]' = evalf(B[mx], 2) * T ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  The peak value of the magnetic flux density is required if the <Font executable="false" family="Arial Narrow">hysteresis</Font> loss is to be calculated.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the magnetic flux density in the core:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">B := evalf(B[mx]/sqrt(2))  :  'B' = evalf(B, 2) * T ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  The effective value of the magnetic flux density is required if the <Font executable="false" family="Arial Narrow">eddy-current</Font> loss is to be calculated.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  Recall that the values of the magnetic field strength and flux density correspond to the operation of the <Font executable="false">CT</Font> at the <Font executable="false" family="Arial Narrow">rated</Font> primary current. If a prospective <Font executable="false" family="Arial Narrow">short-circuit</Font> current is <Font executable="false">k</Font> times greater than the rated primary current, both the corresponding magnetic field strength and flux density will also be <Font executable="false">k</Font> times greater at a constant relative permeability, <Equation input-equation="mu[r]" style="2D Math">NiMmJSNtdUc2IyUickc=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, of the CT-core material. Exemplarily, if the multiple k is</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">k := 25  :  'k' = k ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">then the magnetic field strength and the flux density at <Font executable="false" family="Arial Narrow">short-circuit</Font> would be, respectively</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'H[sc]' = 'k * H'  ;  H[sc] := k*H  :  H[sc_r] := evalf(H[sc], 4)  :  'H[sc]' = H[sc_r] * `A/m` ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'B[sc]' = 'k * B'  ;  B[sc] := k*B  :  B[sc_r] := evalf(B[sc], 3)  :  'B[sc]' = B[sc_r] * T ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">Comment:</Font>  With the assumed value of the multiple <Font executable="false">k</Font>, magnitudes of the primary current vary in the range</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">range_I_p = '[ - k * I_p, k * I_p ]' ;
if frac(k*I_p/10^3) = 0 then print(range_I_p = [round(-k*I_p/10^3), round(k*I_p/10^3)] * `kA`) else print(range_I_p = evalf([-k*I_p/10^3, k*I_p/10^3], 3) * `kA`) fi :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">This implies that the values of the magnetic field strength and magnetic flux density will vary in the following ranges, respectively</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">range_H = '[ -H[sc], H[sc] ]'  ;  range_H = [ -H[sc_r], H[sc_r] ] * `A/m` ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">range_B = '[ -B[sc], B[sc] ]'  ;  range_B = [ -B[sc_r], B[sc_r] ] * T ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">This calls for <Font executable="false">B</Font> to be a linear function of <Font executable="false">H</Font> for a given material at least within the above-determined range of magnetic field strength.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The angular displacement of the magnetic flux with respect to the magnetizing current:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">delta[2] := theta[deg][phi_m] - theta[deg][i_mp] :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">theta[deg][phi*_m] - 'theta[deg][i_mp]' = round(delta[2]) * `\260`;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  The above <Font executable="false">zero</Font> angular displacement verifies that the magnetic flux is <Font executable="false" family="Arial Narrow">co-phasal</Font> with the magnetizing current.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">6. The  emf  induced in the secondary</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  original form:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s_v](x) := evalf(eval(e_s_f))  :  e[s](x) = e[s_v](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  transformed forms:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">a[1] := op(1, op(2, e[s_v](x)))  :  a[2] := op(1, op(1, e[s_v](x))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">E_s_mx := sqrt(a[1]^2 + a[2]^2) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">alpha[2] := sign(a[2]) * evalf(arctan(abs(a[2]/a[1]))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[e_s] := evalf(theta[p] + alpha[2]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][e_s] := evalf(180*theta[e_s]/Pi)  :  theta[deg][e_s_r] := evalf(theta[deg][e_s], 3) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s_v](x) := E_s_mx*sin(x + theta[e_s])  :  e[s](x) = e[s_v](x) * V ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">e[s](x) = evalf(E_s_mx, 2) * sin(180*x/Pi + theta[deg][e_s_r] * `\260`) * V ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the secondary <Font executable="false">emf</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(a) computed from the maximum value of the secondary <Font executable="false">emf</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">E[s] := evalf(E_s_mx/sqrt(2))  :  'E[s]' = evalf(E[s], 2) * V ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(b) computed from the transformer effective <Font executable="false">emf</Font> equation for the secondary:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'E[s]' = 'omega*N[s]*Phi[m][mx]/sqrt(2)'  ;  E[s] := evalf(omega*N[s]*Phi_m_mx/sqrt(2))  :  'E[s]' = evalf(E[s], 2) * V ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">7. The secondary current</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  original form:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[s_v](x) := evalf(eval(i_s_f))  :  i[s](x) = i[s_v](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  transformed forms:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">a[1] := op(1, op(1, i[s_v](x)))  :  a[2] := op(1, op(2, i[s_v](x))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_s_mx := sqrt(a[1]^2 + a[2]^2) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">alpha[1] := sign(a[2]) * evalf(arctan(abs(a[2]/a[1]))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[i_s] := evalf(theta[p] + alpha[1])  :  theta[deg][i_s] := evalf(180*theta[i_s]/Pi) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">if theta[p] = 0 then rd := 2 else rd := 4 end if :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][i_s_r] := evalf(theta[deg][i_s], rd) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[s_v](x) := I_s_mx*sin(x + theta[i_s])  :  i[s](x) = i[s_v](x) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">i[s](x) = evalf(I_s_mx, 3) * sin(180*x/Pi + theta[deg][i_s_r] * `\260`) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The <Font executable="false" family="Arial Narrow">phase error</Font> of the <Font executable="false">CT</Font> at the adopted magnitude of the primary current, defined as the angular displacement between the secondary current and the primary current referred to the secondary:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">epsilon[theta] = abs(abs('theta[deg][i_s]') - abs('theta[deg][i_p_s]'))  ;  epsilon[theta] := abs(abs(theta[deg][i_s]) - abs(theta[deg][i_p_s]))  :  epsilon[theta] := evalf(epsilon[theta], 2)  :  'epsilon[theta]' = epsilon[theta] * `\260`;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  As expected and desired, the <Font executable="false" family="Arial Narrow">phase error</Font> of the <Font executable="false">CT</Font> is a small angle, by which the secondary current <Font executable="false" family="Arial Narrow">leads</Font> the primary current referred to the secondary.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The angular displacement of the secondary current with respect to the secondary <Font executable="false">emf</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">delta[3] := theta[deg][i_s] - theta[deg][e_s] :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">'theta[deg][i_s] - theta[deg][e_s]' = evalf(delta[3], 3) * `\260` ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  This result verifies that the secondary current <Font executable="false" family="Arial Narrow">lags</Font> the secondary <Font executable="false">emf</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the secondary current:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">I_s := evalf(I_s_mx/sqrt(2))  :  'I_s' = evalf(I_s, 3) * 'A' ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  Disregarding the magnetizing current flowing through the primary, the product of the turns ratio with the effective value of the primary current should be <Font executable="false" family="Arial Narrow">approximately equal</Font> to the effective value of the secondary current computed above. Thus,</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'a[t] * I_p' = 'I_s'  ;  'a[t] * I_p' = round(a[t]*I_p) * 'A' ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">8. The  emf  induced in the primary</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[p](x) = 'a[t]*e[s](x)'  ;  e[p_v](x) := a[t]*e[s_v](x)  :  e[p](x) = e[p_v](x) * V ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">E_p_mx := evalf(a[t]*E_s_mx) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][e_p] := theta[deg][e_s]  :  theta[deg][e_p_r] := evalf(theta[deg][e_p], 3) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">e[p](x) = evalf(E_p_mx*10^3, 3) * sin(180*x/Pi + theta[deg][e_p_r] * `\260`) * mV ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The angular displacement of the magnetizing current with respect to the primary <Font executable="false">emf</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">delta[4] := theta[deg][i_mp] - theta[deg][e_p] :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">'theta[deg][i_mp]' - 'theta[deg][e_p]' = evalf(delta[4], 3) * `\260`;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  The magnetizing current is, as expected, <Font executable="false" family="Arial Narrow">in quadrature</Font> with the primary <Font executable="false">emf</Font>, i.e., <Font executable="false" family="Arial Narrow">lags</Font> the <Font executable="false">emf</Font> by <Font encoding="ISO8859-1" executable="false">90\260</Font>.</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the primary <Font executable="false">emf</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(a) computed from the maximum value of the primary <Font executable="false">emf</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">E[p] := evalf(E_p_mx/sqrt(2))  :  'E[p]' = evalf(E[p]*10^3, 3) * mV ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(b) computed from the transformer effective <Font executable="false">emf</Font> equation for the primary:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'E[p]' = 'omega*N[p]*Phi[m][mx]/sqrt(2)'  ;  E[p] := evalf(omega*N[p]* Phi_m_mx/sqrt(2))  :  'E[p]' = evalf(E[p]*10^3, 3) * mV ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">9. The secondary voltage</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  original form:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">v[s_v](x) := evalf(eval(v_s_f))  :  'v[s](x)' = v[s_v](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font encoding="UTF-8" executable="false">\342\200\242</Font>  transformed forms:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">a[1] := op(1, op(2, v[s_v](x)))  :  a[2] := op(1, op(1, v[s_v](x))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">V_s_mx := sqrt(a[1]^2 + a[2]^2) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">alpha[2] := sign(a[2]) * evalf(arctan(abs(a[2]/a[1]))) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[v_s] := evalf(theta[p] + alpha[2]) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[deg][v_s] := evalf(180*theta[v_s]/Pi)  :  theta[deg][v_s_r] := evalf(theta[deg][v_s], 3) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">v[s_v](x) := V_s_mx*sin(x + theta[v_s])  :  'v[s](x)' = v[s_v](x) * V ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">'v[s](x)' = evalf(V_s_mx, 2) * sin(180*x/Pi + theta[deg][v_s_r] * `\260`) * V ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The angular displacement of the secondary current with respect to the secondary voltage:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">delta[5] := theta[deg][i_s] - theta[deg][v_s] :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">'theta[deg][i_s] - theta[deg][v_s]' = evalf(delta[5], 3) * `\260` ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  This result verifies that the secondary current <Font executable="false" family="Arial Narrow">lags</Font> the secondary voltage. The obtained value of the phase difference is close to a typical value of <Font encoding="UTF-8" executable="false">\342\200\22345\302\272.</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">The effective value of the secondary voltage:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">V[s] := evalf(V_s_mx/sqrt(2))  :  'V[s]' = evalf(V[s], 2) * V ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal"><Font background="[0,0,0]" bold="false" family="Times New Roman" foreground="[0,0,0]" italic="false" size="12" underline="false">Verification of KCL at a node of the equivalent circuit diagram of the CT model</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Kirchhoff<Font encoding="UTF-8" executable="false">\342\200\231</Font>s current law applied to a selected node has the following form in the <Font executable="false" family="Arial Narrow">angular domain</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">KCL := i[p_s](x) - i[s](x) - i[mp_s](x) = 0  :  KCL ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Substitution of transformed forms of the currents results in the following form of <Font executable="false">KCL</Font>:</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">KCL:=subs(i[p_s](x) = i[p_s_v](x), i[s](x) = i[s_v](x), i[mp_s](x) = i[mp_s_v](x), KCL)  :  KCL ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Choosing <Font executable="false">x<Font encoding="ISO8859-1"> = 60\272</Font></Font> and substituting its radian equivalent in the left-hand side of the <Font executable="false">KCL</Font> equation give</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">lhs_KCL := subs(x = 60*Pi/180, lhs(KCL))  :  'lhs_KCL' = lhs_KCL ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Floating-point evaluation of the above result yields</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">lhs_KCL := evalf(lhs_KCL)  :  'lhs_KCL' = lhs_KCL ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Text Output" style="Text Output"/></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">Since the deviation from the expected <Font executable="false">zero</Font> is practically <Font executable="false" family="Arial Narrow">negligible</Font>, verification of <Font executable="false">KCL</Font> is considered satisfactory.</Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">9. Numerical Evaluation of Time-invariant Parameters</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">1. The actual value of the magnetic-core reluctance (as viewed from the primary-winding terminals)</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'R[m]' = 'l[av]/(mu[0]*mu[r]*A)'  ;  R[m] := R[m]  :  'R[m]' = round(R[m]) * `1/H` ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">2. The actual value of the magnetizing inductance (as viewed from the primary-winding terminals)</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">L[mp] = 'N[p]^2/R[m]'  ;  L[mp]:=N[p]^2/R[m]  :  'L[mp]' = evalf(L[mp]*10^6, 3) * mu*'H' ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">3. The actual value of the magnetizing reactance (as viewed from the primary-winding terminals)</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">X[mp] = 'omega*L[mp]'  ;  X[mp] := omega*L[mp]  :  'X[mp]' = evalf(X[mp]*10^3, 3) * m*Omega ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">4. The magnetizing inductance referred (reflected) to the secondary</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">L[mp_s] = 'L[mp]/a[t]^2'  ;  L[mp_s_v] := L[mp]/a[t]^2  :  L[mp_s] = evalf(L[mp_s_v]*10^3, 3) * mH ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">5. The magnetizing reactance referred (reflected) to the secondary</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">X[mp_s] = 'X[mp]/a[t]^2'  ;  X[mp_s_v] := X[mp]/a[t]^2  :  X[mp_s] = evalf(X[mp_s_v], 3)*Omega ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">6. The effective value of the  mmf  that establishes the magnetic flux in the core</Text-field></Title><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3">(a) computed from the primary-current linkage:</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'F[m]' = 'N[p]*I_mp'  ;  F[m] := N[p]*I_mp  :  'F[m]' = evalf(F[m], 3) * 'A' ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3"><Font encoding="UTF-8">(b) computed from Amp\303\250re\342\200\231s circuital law:</Font></Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal"/></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'F[m]' = H * 'l[av]'  ;  H := H_v  :  F[m] := evalf(H*l[av])  :  'F[m]' = evalf(F[m], 3) * 'A' ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 3" style="Heading 3"><Font encoding="UTF-8">(c) computed from magnetic Ohm\342\200\231s law:</Font></Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'F[m]' = 'Phi[m]*R[m]'  ;  F[m] := evalf(Phi[m]*R[m])  :  'F[m]' = evalf(F[m], 3) * 'A' ;</Font></Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">7. The modulus of the total impedance of the secondary circuit</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Z[Ts] := evalf(sqrt(R[Ts]^2 + (omega*L[Ts])^2))  :  'Z[Ts]' = evalf(Z[Ts]*10^3, 3) * m*Omega ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">8. The ratio of the magnetizing reactance referred (reflected) to the secondary to the modulus of the total impedance of the secondary circuit</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">ratio[Xmp/Z] := X[mp_s_v]/Z[Ts]  :  'X[mp_s]/Z[Ts]' = evalf(ratio[Xmp/Z], 4) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  This ratio is, as expected, equal to the reciprocal of the ratio of effective values of the primary-winding magnetizing current referred to the secondary and the secondary current, i.e.,</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">reciprocal_ratio[Imp_s/Is] := 1/(I_mp_s/I_s) :</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'1/(I_mp_s/I_s)' = evalf(reciprocal_ratio[Imp_s/Is], 4) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  Since the total impedance of the secondary circuit, <Equation input-equation="Z[Ts]" style="2D Math">NiMmJSJaRzYjJSNUc0c=</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, is <Font family="Arial Narrow">very small</Font> compared with the magnetizing reactance referred to the secondary, </Font><Equation input-equation="X[mp_s]" style="2D Math">NiMmJSJYRzYjJSVtcF9zRw==</Equation><Font bold="false" italic="false" style="2D Math" underline="false">, or</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">'Z[Ts]' = evalf('X[mp_s]'/ratio[Xmp/Z], 2) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">almost all the primary current referred to the secondary passes through the <Font executable="false" family="Arial Narrow">burden</Font> and the <Font executable="false">CT</Font> operates similarly as a voltage transformer on <Font executable="false" family="Arial Narrow">short-circuit</Font>.</Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">9. The impedance angle of the secondary circuit</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">theta[s] = 'arccos(R[Ts]/Z[Ts])'  ;  theta[s] := evalf(arccos(R[Ts]/Z[Ts])*180/Pi)  :  'theta[s]' = evalf(theta[s], 3) * `\260` ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  The impedance angle of the secondary circuit is, as expected, equal to the absolute value of the angular displacement of the secondary current with respect to the secondary <Font executable="false">emf</Font>, which was computed earlier, i.e.,</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">abs('theta[deg][i_s] - theta[deg][e_s]') = abs(evalf(delta[3], 3)) * `\260` ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">10. The modulus of impedance of the burden (overcurrent-relay coil)</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">Z[b] := evalf(sqrt(R[b]^2 + (omega*L[b])^2))  :  'Z[b]' = evalf(Z[b]*10^3, 3)*m*Omega ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">11. The impedance angle of the burden (overcurrent-relay coil)</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">theta[b] = 'arccos(R[b]/Z[b])'  ;  theta[b] := evalf(arccos(R[b]/Z[b])*180/Pi)  :  'theta[b]' = evalf(theta[b], 3) * `\260` ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal"><Font executable="false" family="Arial Narrow">N.B.</Font>  The impedance angle of the burden is, as expected, equal to the absolute value of the angular displacement of the secondary current with respect to the secondary voltage, which was computed earlier, i.e.,</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">abs('theta[deg][i_s] - theta[deg][v_s]') = abs(evalf(delta[5], 3)) * `\260` ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">12. The power factor of the secondary circuit</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">PF[s] = 'cos(R[Ts]/Z[Ts])'  ;  PF[s] := cos(R[Ts]/Z[Ts])  :  'PF[s]' = evalf(PF[s], 2) ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">13. The power factor of the burden (overcurrent-relay coil)</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">PF[b] = 'cos(R[b]/Z[b])'  ;  PF[b] := cos(R[b]/Z[b])  :  'PF[b]' = evalf(PF[b], 2) ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">14. The active power dissipated in the total resistance of the secondary circuit</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">P[s] = 'I_s^2 * R[Ts]'  ;  P[s] := I_s^2 * R[Ts]  :  'P[s]' = evalf(P[s], 3) * W ;</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">15. The apparent power absorbed by the secondary circuit</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font encoding="ISO8859-1" italic="false" underline="false">S[s] = 'I_s^2 * Z[Ts]'  ;  S[s] := I_s^2 * Z[Ts]  :  'S[s]' = evalf(S[s], 3) * `V\267A` ;</Font></Text-field></Input></Group></Section></Section><Section collapsed="true"><Title><Text-field layout="Heading 1" style="Heading 1">10. Restoring Symbolic Forms of Formulae for Time-varying Parameters</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">The <Font executable="false" family="Arial Narrow">symbolic</Font> forms of the formulae for the functions obtained earlier in the <Font executable="false" family="Arial Narrow">angular domain</Font> may be required for some further analyses. These forms are made available at this point of the computational algorithm through a procedure, which is provided hereunder in steps.</Text-field></Input></Group><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Step 1. Save the numerical values of all the parameters found in these functions, using some other names (e.g., by adding the suffix _v to the names 'bearing' numerical values):</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[p_v] := theta[p]  :  I_p_mx_v := I_p_mx  :  L[b_v] := L[b]  :  L[Ts_v] := L[Ts]  :  N[p_v] := N[p]  :  N[s_v] := N[s]  :  omega_v := omega  :  R[m_v] := R[m]  :  R[b_v] := R[b]  :  R[Ts_v] := R[Ts] :</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Step 2. Unassign the numerical values from the names of all the parameters (or, 'unevaluate' the names 'bearing' numerical values) to restore the symbolic names of the parameters (by enclosing each name in single forward quotes):</Text-field></Title><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">theta[p] := 'theta[p]'  :  I_p_mx := 'I_p_mx'  :  L[b] := 'L[b]'  :  L[Ts] := 'L[Ts]'  :  N[p] := 'N[p]'  :  N[s] := 'N[s]'  :  omega := 'omega'  :  R[m] := 'R[m]'  :  R[b] := 'R[b]'  :  R[Ts] := 'R[Ts]' :</Font></Text-field></Input></Group></Section><Section collapsed="true"><Title><Text-field layout="Heading 2" style="Heading 2">Step 3. Re-enter and display the required functions in the angular domain:</Text-field></Title><Group><Input><Text-field layout="Normal" style="Normal">(a) primary current</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[p](x) := i_p_f  :  'i[p](x)' = i[p](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(b) magnetizing current flowing in the primary</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[mp](x) := i_mp_f  :  'i[mp](x)' = i[mp](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(c) magnetic flux</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">phi[m](x) := phi_m_f  :  'phi[m](x)' = phi[m](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(d) secondary <Font executable="false">emf</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[s](x) := e_s_f  :  'e[s](x)' = e[s](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(e) secondary current</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">i[s](x) := i_s_f  :  'i[s](x)' = i[s](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(f) primary <Font executable="false">emf</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">e[p](x) := e_p_f  :  'e[p](x)' = e[p](x) ;</Font></Text-field></Input></Group><Group><Input><Text-field layout="Normal" style="Normal">(g) secondary voltage</Text-field></Input></Group><Group><Input><Text-field layout="Normal" prompt="&gt; " style="Maple Input"><Font italic="false" underline="false">v[s](x) := v_s_f  :  'v[s](x)' = v[s](x) ;</Font></Text-field></Input></Group></Section><Group><Input><Text-field layout="Normal" style="Text">All the formulae describing the time-varying parameters are again available in their symbolic forms and may be used for a further analysis, if required.
</Text-field></Input></Group></Section><Group><Input><Text-field layout="Normal" style="Text"/><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>