A Numerical Method of the Euler-Bernoulli Beam with Optimal Local Kelvin-Voigt Damping

<table>The optimal controls for <svg height="16.275" id="M212" style="vertical-align:-3.13504pt" version="1.1" viewbox="0 0 52.099998 16.275" width="52.099998" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.137)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,5.944,12.137)"><path d="M381 665q131 0 226.5 -93t95.5 -239q0 -158 -96 -253t-238 -95q-137 0 -231 95t-94 238q0 142 92.5 244.5t244.5 102.5zM359 629q-89 0 -151 -74.5t-62 -208.5q0 -141 69 -233t175 -92q90 0 150.5 74t60.5 211q0 152 -69 237.5t-173 85.5z" id="x4F"></path></g><g transform="matrix(.017,-0,0,-.017,18.931,12.137)"><path d="M614 175l29 -10q-33 -109 -57 -154q-121 -26 -184 -26q-90 0 -160.5 29t-112.5 77t-63.5 105.5t-21.5 119.5q0 157 108 253t277 96q36 0 71.5 -5t69 -13.5t36.5 -8.5q15 -102 20 -150l-29 -8q-20 79 -66.5 114t-128.5 35q-119 0 -187.5 -86t-68.5 -207
q0 -140 73.5 -227.5t188.5 -87.5q73 0 119.5 37.5t86.5 116.5z" id="x43"></path></g><g transform="matrix(.017,-0,0,-.017,30.235,12.137)"><path d="M46 650h247q112 0 169 -41q64 -49 64 -141q0 -91 -56.5 -143t-138.5 -62q-23 -3 -45 -1l-74 19v-157q0 -63 15.5 -77t83.5 -19v-28h-271v28q60 5 74 20t14 77v400q0 65 -12.5 79t-69.5 18v28zM212 559v-245q31 -12 73 -12q59 0 102 38t43 124q0 77 -42.5 115.5
t-108.5 38.5q-44 0 -56 -12q-11 -8 -11 -47z" id="x50"></path></g><g transform="matrix(.017,-0,0,-.017,39.806,12.137)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g>
<g transform="matrix(.012,-0,0,-.012,45.7,16.212)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g>
</svg> and <svg height="16.275" id="M213" style="vertical-align:-3.13504pt" version="1.1" viewbox="0 0 52.099998 16.275" width="52.099998" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.137)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,5.944,12.137)"><path d="M381 665q131 0 226.5 -93t95.5 -239q0 -158 -96 -253t-238 -95q-137 0 -231 95t-94 238q0 142 92.5 244.5t244.5 102.5zM359 629q-89 0 -151 -74.5t-62 -208.5q0 -141 69 -233t175 -92q90 0 150.5 74t60.5 211q0 152 -69 237.5t-173 85.5z" id="x4F"></path></g><g transform="matrix(.017,-0,0,-.017,18.931,12.137)"><path d="M614 175l29 -10q-33 -109 -57 -154q-121 -26 -184 -26q-90 0 -160.5 29t-112.5 77t-63.5 105.5t-21.5 119.5q0 157 108 253t277 96q36 0 71.5 -5t69 -13.5t36.5 -8.5q15 -102 20 -150l-29 -8q-20 79 -66.5 114t-128.5 35q-119 0 -187.5 -86t-68.5 -207
q0 -140 73.5 -227.5t188.5 -87.5q73 0 119.5 37.5t86.5 116.5z" id="x43"></path></g><g transform="matrix(.017,-0,0,-.017,30.235,12.137)"><path d="M46 650h247q112 0 169 -41q64 -49 64 -141q0 -91 -56.5 -143t-138.5 -62q-23 -3 -45 -1l-74 19v-157q0 -63 15.5 -77t83.5 -19v-28h-271v28q60 5 74 20t14 77v400q0 65 -12.5 79t-69.5 18v28zM212 559v-245q31 -12 73 -12q59 0 102 38t43 124q0 77 -42.5 115.5
t-108.5 38.5q-44 0 -56 -12q-11 -8 -11 -47z" id="x50"></path></g><g transform="matrix(.017,-0,0,-.017,39.806,12.137)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g>
<g transform="matrix(.012,-0,0,-.012,45.7,16.212)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g>
</svg>.</table>

Journal of Applied Mathematics

fig2

Figure 2

Figure 2: A Numerical Method of the Euler-Bernoulli Beam with Optimal Local Kelvin-Voigt Damping 