A Nonlinear Finite Element Method for Magnetoelectric Composite and the Study on the Influence of Interfacial Bonding

<table>Numerical results on the field distributions of a three-layer laminated composite with an interfacial crack in the middle. (a) Shear stress <svg height="15.1125" id="M156" style="vertical-align:-5.76988pt" version="1.1" viewbox="0 0 22.5 15.1125" width="22.5" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,7.863)"><path d="M471 456q-22 -60 -37 -85q-20 -3 -84 -3q-9 0 -29 2t-29 2q-38 -145 -54 -234q-15 -80 14 -80q23 0 68 38l14 -25q-78 -83 -137 -83q-50 0 -50 64q0 25 7 57l67 266q-57 4 -98 -9t-80 -49l-20 21q14 18 27 31t36.5 30t56.5 26t73 9q26 0 82.5 -2t89.5 -2q21 0 32 5t24 24
z" id="x1D70F"></path></g>
<g transform="matrix(.012,-0,0,-.012,8.275,11.95)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g><g transform="matrix(.012,-0,0,-.012,14.926,11.95)"><path d="M556 393q0 -39 -36 -106q-42 -78 -185 -279q-47 -66 -81 -108t-117 -135l-112 -26l-8 22q150 90 251 219q-6 136 -39 340q-8 53 -21 53q-6 0 -27 -19.5t-38 -42.5l-16 26q80 111 127 111q23 0 35 -28t20 -90q18 -137 27 -263h2q142 200 142 279q0 24 -14 48q-4 7 5 26
q13 28 43 28q18 0 30 -15.5t12 -39.5z" id="x1D466"></path></g>
</svg>, (b) normal stress <svg height="12.325" id="M157" style="vertical-align:-3.3907pt" version="1.1" viewbox="0 0 17 12.325" width="17" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,8.038)"><path d="M548 455q-32 -77 -74 -77q-37 0 -103 10l-3 -2q36 -30 48.5 -62t12.5 -80q0 -103 -75 -179.5t-175 -76.5q-70 0 -113 45t-43 128q0 115 83.5 200.5t209.5 85.5q31 0 77.5 -4t67.5 -4q36 0 62 30zM350 274q0 54 -17 86q-19 35 -57 35q-50 0 -90 -37.5t-59 -91.5t-19 -109
q0 -61 24.5 -96.5t65.5 -35.5q51 0 87.5 46.5t50.5 100.5t14 102z" id="x1D70E"></path></g>
<g transform="matrix(.012,-0,0,-.012,9.663,12.112)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g>
</svg>, (c) electric field <svg height="18.362499" id="M158" style="vertical-align:-5.76988pt" version="1.1" viewbox="0 0 18.375 18.362499" width="18.375" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.113)"><path d="M609 650l-19 -162l-30 -2q2 69 -14 94q-9 18 -28.5 26t-74.5 8h-88q-30 0 -37 -6.5t-12 -36.5l-41 -212h103q48 0 69 5.5t32 19t26 50.5h29l-40 -198h-30q2 58 -11.5 70.5t-85.5 12.5h-99l-30 -167q-17 -87 3 -101q18 -15 111 -15q59 0 89 7t54 27q8 7 15.5 16t16 21.5
l14 20.5t15 23.5t13.5 21.5l29 -10q-52 -129 -71 -163h-500l6 28q66 4 83 17.5t28 73.5l77 409q11 61 -0.5 75.5t-78.5 18.5l10 28h467z" id="x1D438"></path></g>
<g transform="matrix(.012,-0,0,-.012,10.8,15.188)"><path d="M556 393q0 -39 -36 -106q-42 -78 -185 -279q-47 -66 -81 -108t-117 -135l-112 -26l-8 22q150 90 251 219q-6 136 -39 340q-8 53 -21 53q-6 0 -27 -19.5t-38 -42.5l-16 26q80 111 127 111q23 0 35 -28t20 -90q18 -137 27 -263h2q142 200 142 279q0 24 -14 48q-4 7 5 26
q13 28 43 28q18 0 30 -15.5t12 -39.5z" id="x1D466"></path></g>
</svg>, (d) magnetization <svg height="15.4" id="M159" style="vertical-align:-3.3907pt" version="1.1" viewbox="0 0 24.75 15.4" width="24.75" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.112)"><path d="M998 650l-8 -28q-71 -4 -86 -16t-22 -69l-50 -397q-3 -28 -4.5 -44t2 -29t6.5 -18.5t17 -10.5t24.5 -6.5t37.5 -3.5l-8 -28h-271l7 28q63 6 78 22t25 90l60 415h-2l-353 -552h-23l-130 536h-2l-70 -254q-44 -158 -47 -188q-5 -38 9 -51t71 -18l-6 -28h-241l8 28
q45 4 67 18.5t35 45.5q16 38 74 233l52 173q24 79 11.5 98t-89.5 26l6 28h177l136 -508l337 508h172z" id="x1D440"></path></g>
<g transform="matrix(.012,-0,0,-.012,17.412,15.187)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g>
</svg>. Stresses do not concentrate near the tip of the cracks. The applied magnetic field is 630 Oe.</table>

Mathematical Problems in Engineering

fig10

Figure 10

Figure 10: A Nonlinear Finite Element Method for Magnetoelectric Composite and the Study on the Influence of Interfacial Bonding 