FOPI_1_part.ipynb

{
 "metadata": {
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.2"
  },
  "orig_nbformat": 2,
  "kernelspec": {
   "name": "python392jvsc74a57bd0d20511710d0ef1b8f61e9561ea7a1dacb95f86736785c369355f9020fcd31a57",
   "display_name": "Python 3.9.2 64-bit"
  },
  "metadata": {
   "interpreter": {
    "hash": "d20511710d0ef1b8f61e9561ea7a1dacb95f86736785c369355f9020fcd31a57"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2,
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      " Sm= 4.155e-04\n Sδ= 4.712e-04\n Sm/Sδ= 8.817e-01\n Bm= 5.898e-01\n\n Bm_Hm= 1.072e-05\n\n Lm= 8.681e-03\n Dm= 2.300e-02\n\n"
     ]
    }
   ],
   "source": [
    "import math as m\n",
    "###3.1.1\n",
    "#///условия////\n",
    "a=2*1e-3\n",
    "δ=1.5*1e-3\n",
    "Dδ=25*1e-3\n",
    "h=6*1e-3\n",
    "Bδ=0.4\n",
    "Dm=Dδ-a\n",
    "k=0.3\n",
    "μ0=m.pi*4*1e-7\n",
    "#///условия////\n",
    "\n",
    "###3.2.1\n",
    "\n",
    "Sm=m.pi*Dm**2/4\n",
    "Sδ=h*m.pi*Dδ\n",
    "Sm_Sδ=Sm/Sδ\n",
    "Bm=Bδ*1/Sm_Sδ*(1+k)\n",
    "print(' Sm= %.3e\\n' %Sm,'Sδ= %.3e\\n' %Sδ,'Sm/Sδ= %.3e\\n' %Sm_Sδ,'Bm= %.3e\\n' %Bm)\n",
    "\n",
    "Hm=55*1e3\n",
    "Bm_Hm=Bm/Hm\n",
    "print(' Bm_Hm= %.3e\\n' %Bm_Hm)\n",
    "Lm=Bm_Hm*(δ/(μ0*(1+k)))*Sm_Sδ\n",
    "\n",
    "print(' Lm= %.3e\\n' %Lm,'Dm= %.3e\\n' %Dm)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      " Lm= 1.472e-02\n Dm= 1.800e-02\n Sm_Sδ= 5.400e-01\n\n"
     ]
    }
   ],
   "source": [
    "#Если не сошлось\n",
    "Bm=0.95 #вводим желаемое значение исходя из рабочей точки графика\n",
    "Hm=32*1e3 #эквивалентное ей значение напряженности\n",
    "Dm=m.sqrt((Bδ*4*Dδ*h*(1+k))/(Bm))\n",
    "Dm=18*1e-3 # округляем\n",
    "Sm=m.pi*Dm**2/4\n",
    "Sm_Sδ=Sm/Sδ\n",
    "Bm_Hm=Bm/Hm\n",
    "Lm=Bm_Hm*(δ/(μ0*(1+k)))*Sm_Sδ\n",
    "print(' Lm= %.3e\\n' %Lm,'Dm= %.3e\\n' %Dm,'Sm_Sδ= %.3e\\n' %Sm_Sδ)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      " Φm= 2.417e-04\n Φδ= 1.885e-04\n Bδ= 4.000e-01\n\n"
     ]
    }
   ],
   "source": [
    "###3.2.2 проверка \n",
    "Bm_Hm=μ0*(1+k)*Sδ/Sm*Lm/δ\n",
    "Hm2=20*1e3 #указываем значение точки по y через которую будем проводить прямую\n",
    "Bm2=Bm_Hm*Hm2\n",
    "#проводим прямую и смотрим значения пересечения\n",
    "Bm_izm=0.97\n",
    "Hm_izm=31*1e3\n",
    "Φm=Bm*Sm # ищем потоки и сравниваем с условием индукцию зазора\n",
    "Φδ=Bδ*Sδ\n",
    "Bδ=Φδ/Sδ\n",
    "print(' Φm= %.3e\\n' %Φm,'Φδ= %.3e\\n' %Φδ,'Bδ= %.3e\\n' %Bδ)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "import math as m\n",
    "###4.1.1\n",
    "#///условия////\n",
    "a=10 \n",
    "b=10\n",
    "c=20\n",
    "dδ=0.2\n",
    "h=50\n",
    "S0=1500\n",
    "Fizm=20\n",
    "f=a\n",
    "#///условия////\n",
    "d0=0.2\n",
    "d1=0.22\n",
    "lamb=1.5\n",
    "eps=3\n",
    "μ0=m.pi*4*1e-7"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "w=  8508\n"
     ]
    }
   ],
   "source": [
    "f0_list=[0.6,0.55,0.5]\n",
    "f0=f0_list[0] # Выбираем сечение провода(у всех дефолтное)\n",
    "w=4*f0*(c-eps-2*lamb)*(h-lamb-f)/(m.pi*d1**2) #считаем количесвто витков, округляем до целого\n",
    "print('w= ','%.0f'%w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "output_type": "execute_result",
     "data": {
      "text/plain": [
       "42.72344373373544"
      ]
     },
     "metadata": {},
     "execution_count": 20
    }
   ],
   "source": [
    "qq=(365.502+(133591+(181184-0.04*(1.636802e+03)**2))**0.5)/19.22\n",
    "qq\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "output_type": "execute_result",
     "data": {
      "text/plain": [
       "310488.4067404736"
      ]
     },
     "metadata": {},
     "execution_count": 21
    }
   ],
   "source": [
    "133591+(181184-(1.636802e+03*0.04)**2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "δ0min1= 1.0 \nδ0max1= 1.333\n\nδ0min2= 1.0 \nδ0max2= 2.0\n"
     ]
    }
   ],
   "source": [
    "δ0max1=dδ/0.15 #определяем границы начального зазора 1 соотношение\n",
    "δ0min1=dδ/0.2\n",
    "print('δ0min1=',δ0min1,'\\nδ0max1=','%.3f'%δ0max1)\n",
    "\n",
    "δ0min2=a*0.1 #определяем границы начального зазора 2  соотношение\n",
    "δ0max2=a*0.2\n",
    "print('\\nδ0min2=',δ0min2,'\\nδ0max2=',δ0max2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "δ0= 1.0\n"
     ]
    }
   ],
   "source": [
    "#Выбрать наименьшее значение, которые находится сразу в обоих множествах\n",
    "min = 1e10\n",
    "if δ0min1<δ0min2: \n",
    "    δ0=round(δ0min2,2)\n",
    "else:\n",
    "    δ0=round(δ0min1,2)\n",
    "print('δ0=',δ0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "δ0max= 1.2 \nδ0min= 0.8\n"
     ]
    }
   ],
   "source": [
    "δ0max=δ0+dδ # добавляем ход\n",
    "δ0min=δ0-dδ\n",
    "print('δ0max=',δ0max,'\\nδ0min=',δ0min)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "def p4_1_3_2():\n",
    "    global m1,m2\n",
    "    #4.1.3.2.\n",
    "    #Для расчета имеем п-образный сердечник. Для него следует выбирать      индукцию B=0.2\n",
    "    B=0.3 #0,3 для Ш-образного\n",
    "    #Построили схему замещения\n",
    "    #Определение магнитного сопротивления примаксимальном воздушном         зазоре\n",
    "    X_=c/2\n",
    "    print('X_=',X_)\n",
    "    #по рис. 7 по известной X'/δmax определяем Zb/δmax\n",
    "    print('X_/δ0max=',X_/δ0max) \n",
    "    #По графику ищем соответствие (для значений  <35 коэффициент 1.3)\n",
    "    Zb_δmax=X_/δ0max /1.3\n",
    "    print('Zb_δmax=',Zb_δmax) \n",
    "    Zb=Zb_δmax*δ0max\n",
    "    print('Zb=',Zb)\n",
    "    Z_a=a\n",
    "    print('Z_a=',Z_a)\n",
    "    Z__b=Z_a\n",
    "    print('Z__b=',Z__b)\n",
    "    x__δmax= Z__b/δ0max*1.3\n",
    "    x__=x__δmax*δ0max\n",
    "    print('x__=','%.3f'%x__)\n",
    "    # НОВЫЙ ШАГ\n",
    "    m1=(Z__b+(x__-δ0max/2))/2\n",
    "    m2=(Zb+(X_-δ0max))/2\n",
    "    print('m1=%.2f'%m1,'\\nm2=%.2f'%m2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "X_= 10.0\nX_/δ0max= 8.333333333333334\nZb_δmax= 6.410256410256411\nZb= 7.6923076923076925\nZ_a= 10\nZ__b= 10\nx__= 13.000\nm1=11.20 \nm2=8.25\n"
     ]
    }
   ],
   "source": [
    "p4_1_3_2()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "def x_δ (δ0max,m1=m1,m2=m2,w=w,k1=1):\n",
    "    global Ze_max,Ze_min\n",
    "    B=0.3#0,3 для Ш-образного\n",
    "    #Построили схему замещения\n",
    "    G1=μ0*((2*b*1e-3)/(m.pi*((δ0max/m2)+0.5)))\n",
    "    print('m2=',m2,'\\n\\nG1=','%.3f e-9'%(G1*1e9))\n",
    "\n",
    "    G2=μ0*0.52*b*1e-3\n",
    "    G3=μ0*(a*b*1e-3/δ0max)\n",
    "    G4=μ0*0.26*b*1e-3\n",
    "    print('G2=','%.3f e-9'%(G2*1e9))\n",
    "    print('G3=','%.3f e-9'%(G3*1e9))\n",
    "    print('G4=','%.3f e-9'%(G4*1e9))\n",
    "\n",
    "    \n",
    "    G5=μ0*((2*b*1e-3)/(m.pi*((δ0max/m1)+1)))\n",
    "    print('m1=',m1,'\\n\\nG5=','%.3fe-9'%(G5*1e9))\n",
    "    #до зазора δ2\n",
    "\n",
    "    G6=μ0*0.077*δ0max*1e-3\n",
    "    G7=μ0*m1*1e-3/4\n",
    "    G8=2*G6+2*G7\n",
    "    print('G6=','%.3f e-9'%(G6*1e9))\n",
    "    print('G7=','%.3f e-9'%(G7*1e9))\n",
    "    print('G8=','%.3f e-9'%(G8*1e9))\n",
    "\n",
    "    Gb4=G1+G2+G3+3*(G4+G5)+2*(G6+G7+G8)\n",
    "    print('\\nGb4=','%.3f e-9'%(Gb4*1e9))\n",
    "\n",
    "    Rb4=1/Gb4\n",
    "    print('\\nRb4=','%.3f e6'%(Rb4*1e-6))\n",
    "\n",
    "    Gb2=G1+G2+G3+2*(G4+G5+G8)\n",
    "    print('\\nGb2=','%.3f e-9'%(Gb2*1e9))\n",
    "    Rb2=1/Gb2\n",
    "    print('\\nRb2=','%.3f e6'%(Rb2*1e-6))\n",
    "\n",
    "    R_ms=(Rb2+Rb4)/2\n",
    "\n",
    "    print('\\nR_ms','%.3f e6'%(R_ms*1e-6))\n",
    "    #4.1.3.3.\n",
    "    ro=0.0174\n",
    "    \n",
    "    # lsr=2*(2*(c-eps-lamb)+3*a+2*eps+4*lamb) для Ш-образного\n",
    "    lsr=4*(h-a-lamb)/k1+8*lamb+2*eps+4*a  # Для П-образного\n",
    "\n",
    "    q=m.pi*d0**2/4\n",
    "    ow=2*m.pi*50\n",
    "    Ra_r=ro*lsr*1e-3*w/q\n",
    "    Ra_im=ow*w**2/R_ms\n",
    "    Ze=(Ra_r**2+Ra_im**2)**0.5\n",
    "    if δ0max==1.2:\n",
    "        Ze_min=(Ra_r**2+Ra_im**2)**0.5\n",
    "        print('\\nZe_min= %.2f e3 '%(Ze_min*1e-3))\n",
    "    else:\n",
    "        Ze_max=(Ra_r**2+Ra_im**2)**0.5\n",
    "        print('\\nZe_max= %.2f e3 '%(Ze_max*1e-3))\n",
    "    print('lsr=%.2f'%lsr)\n",
    "\n",
    "    print('\\nRa_r=%.2f'%Ra_r)\n",
    "    print('Ra_im=%.2f'%Ra_im)\n",
    "\n",
    "    \n",
    "    #4.1.3.4.\n",
    "    I=(B*2*a*b*R_ms*10**(-6))/(w*(2**0.5))\n",
    "    print('I= %.4f'%I)\n",
    "    j=I/q\n",
    "    print('j<=2..2.5 j= %.2f'%j)\n",
    "    U=I*Ze\n",
    "    print('U= %.2f'%U)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "m2= 8.246153846153847 \n\nG1= 12.393 e-9\nG2= 6.535 e-9\nG3= 104.720 e-9\nG4= 3.267 e-9\nm1= 11.2 \n\nG5= 7.226e-9\nG6= 0.116 e-9\nG7= 3.519 e-9\nG8= 7.269 e-9\n\nGb4= 176.935 e-9\n\nRb4= 5.652 e6\n\nGb2= 159.172 e-9\n\nRb2= 6.283 e6\n\nR_ms 5.967 e6\n\nZe_min= 3.94 e3 \nlsr=212.00\n\nRa_r=998.94\nRa_im=3810.59\nI= 0.0298\nj<=2..2.5 j= 0.95\nU= 117.23\n\n\n\n\n\n\nm2= 8.246153846153847 \n\nG1= 13.400 e-9\nG2= 6.535 e-9\nG3= 157.080 e-9\nG4= 3.267 e-9\nm1= 11.2 \n\nG5= 7.467e-9\nG6= 0.077 e-9\nG7= 3.519 e-9\nG8= 7.192 e-9\n\nGb4= 230.792 e-9\n\nRb4= 4.333 e6\n\nGb2= 212.866 e-9\n\nRb2= 4.698 e6\n\nR_ms 4.515 e6\n\nZe_max= 5.13 e3 \nlsr=212.00\n\nRa_r=998.94\nRa_im=5035.79\nI= 0.0225\nj<=2..2.5 j= 0.72\nU= 115.60\ndZe=1194.57 \n"
     ]
    }
   ],
   "source": [
    "x_δ (δ0max)\n",
    "print('\\n\\n\\n\\n\\n')\n",
    "x_δ (δ0min)\n",
    "print('dZe=%.2f '%((Ze_max-Ze_min)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "dZe=1.194567e+03\nS=2.986417e+03\nS0=1.500000e+03\nk1=1.411008\nw_new=6029.415082\n"
     ]
    }
   ],
   "source": [
    "dZe=Ze_max-Ze_min\n",
    "print('dZe=%e'%dZe)\n",
    "S=dZe/(2*dδ)\n",
    "print('S=%e'%S)\n",
    "print('S0=%e'%S0)\n",
    "k1=(S/S0)**0.5\n",
    "print('k1=%f'%k1)\n",
    "w_new=w/k1\n",
    "print('w_new=%f'%w_new)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "output_type": "execute_result",
     "data": {
      "text/plain": [
       "8507.555139821316"
      ]
     },
     "metadata": {},
     "execution_count": 37
    }
   ],
   "source": [
    "w"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "m2= 8.246153846153847 \n\nG1= 13.400 e-9\nG2= 6.535 e-9\nG3= 157.080 e-9\nG4= 3.267 e-9\nm1= 11.2 \n\nG5= 7.467e-9\nG6= 0.077 e-9\nG7= 3.519 e-9\nG8= 7.192 e-9\n\nGb4= 230.792 e-9\n\nRb4= 4.333 e6\n\nGb2= 212.866 e-9\n\nRb2= 4.698 e6\n\nR_ms 4.515 e6\n\nZe_max= 2.59 e3 \nlsr=167.14\n\nRa_r=558.12\nRa_im=2529.00\nI= 0.0318\nj<=2..2.5 j= 1.01\nU= 82.29\nm2= 8.246153846153847 \n\nG1= 12.393 e-9\nG2= 6.535 e-9\nG3= 104.720 e-9\nG4= 3.267 e-9\nm1= 11.2 \n\nG5= 7.226e-9\nG6= 0.116 e-9\nG7= 3.519 e-9\nG8= 7.269 e-9\n\nGb4= 176.935 e-9\n\nRb4= 5.652 e6\n\nGb2= 159.172 e-9\n\nRb2= 6.283 e6\n\nR_ms 5.967 e6\n\nZe_min= 1.99 e3 \nlsr=167.14\n\nRa_r=558.12\nRa_im=1913.70\nI= 0.0420\nj<=2..2.5 j= 1.34\nU= 83.71\n\n\n\n\n\n\ndZe=5.964306e+02\nS=1.491077e+03\nk1=0.997021\nw_new=6047.013424\n\n\n\n\n\n\nm2= 8.246153846153847 \n\nG1= 13.400 e-9\nG2= 6.535 e-9\nG3= 157.080 e-9\nG4= 3.267 e-9\nm1= 11.2 \n\nG5= 7.467e-9\nG6= 0.077 e-9\nG7= 3.519 e-9\nG8= 7.192 e-9\n\nGb4= 230.792 e-9\n\nRb4= 4.333 e6\n\nGb2= 212.866 e-9\n\nRb2= 4.698 e6\n\nR_ms 4.515 e6\n\nZe_max= 2.64 e3 \nlsr=212.46\n\nRa_r=711.57\nRa_im=2544.12\nI= 0.0317\nj<=2..2.5 j= 1.01\nU= 83.69\nm2= 8.246153846153847 \n\nG1= 12.393 e-9\nG2= 6.535 e-9\nG3= 104.720 e-9\nG4= 3.267 e-9\nm1= 11.2 \n\nG5= 7.226e-9\nG6= 0.116 e-9\nG7= 3.519 e-9\nG8= 7.269 e-9\n\nGb4= 176.935 e-9\n\nRb4= 5.652 e6\n\nGb2= 159.172 e-9\n\nRb2= 6.283 e6\n\nR_ms 5.967 e6\n\nZe_min= 2.05 e3 \nlsr=212.46\n\nRa_r=711.57\nRa_im=1925.14\nI= 0.0419\nj<=2..2.5 j= 1.33\nU= 85.93\n\n\n\n\n\n\ndZe=5.893232e+02\nS=1.473308e+03\nk1=0.991063\nw_new=6101.531106\n\n\n\n\n\n\nm2= 8.246153846153847 \n\nG1= 13.400 e-9\nG2= 6.535 e-9\nG3= 157.080 e-9\nG4= 3.267 e-9\nm1= 11.2 \n\nG5= 7.467e-9\nG6= 0.077 e-9\nG7= 3.519 e-9\nG8= 7.192 e-9\n\nGb4= 230.792 e-9\n\nRb4= 4.333 e6\n\nGb2= 212.866 e-9\n\nRb2= 4.698 e6\n\nR_ms 4.515 e6\n\nZe_max= 2.69 e3 \nlsr=213.39\n\nRa_r=721.18\nRa_im=2590.61\nI= 0.0314\nj<=2..2.5 j= 1.00\nU= 84.42\nm2= 8.246153846153847 \n\nG1= 12.393 e-9\nG2= 6.535 e-9\nG3= 104.720 e-9\nG4= 3.267 e-9\nm1= 11.2 \n\nG5= 7.226e-9\nG6= 0.116 e-9\nG7= 3.519 e-9\nG8= 7.269 e-9\n\nGb4= 176.935 e-9\n\nRb4= 5.652 e6\n\nGb2= 159.172 e-9\n\nRb2= 6.283 e6\n\nR_ms 5.967 e6\n\nZe_min= 2.09 e3 \nlsr=213.39\n\nRa_r=721.18\nRa_im=1960.32\nI= 0.0415\nj<=2..2.5 j= 1.32\nU= 86.66\n\n\n\n\n\n\ndZe=6.003542e+02\nS=1.500885e+03\nk1=1.000295\nw_new=6100.199920\n\n\n\n\n\n\nz=3\nw_new=6100.000000\n"
     ]
    }
   ],
   "source": [
    "z=0\n",
    "while m.fabs(S-S0)>1:\n",
    "    z=z+1\n",
    "    w_new=round(w_new)\n",
    "    x_δ (δ0max=δ0min,w=w_new,k1=k1)\n",
    "    x_δ (δ0max=δ0max,w=w_new,k1=k1)\n",
    "    print('\\n\\n\\n\\n\\n')\n",
    "    dZe=m.fabs(Ze_max-Ze_min)\n",
    "    print('dZe=%e'%dZe)\n",
    "    S=dZe/(2*dδ)\n",
    "    print('S=%e'%S)\n",
    "    k1=(S/S0)**0.5\n",
    "    print('k1=%f'%k1)\n",
    "    w_new=w_new/k1\n",
    "    print('w_new=%f'%w_new)\n",
    "    print('\\n\\n\\n\\n\\n')\n",
    "print('z=%d'%z)   \n",
    "m.fabs(S-S0)\n",
    "print('w_new=%f'%round(w_new))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "output_type": "stream",
     "name": "stdout",
     "text": [
      "Fsum=7.12\nFizm=20.00\nFizm>Fsum*5..10\n20.00>35.59..71.18\n"
     ]
    }
   ],
   "source": [
    "w_new=round(w_new)\n",
    "#4.1.3.7. Определение измерительного усилия\n",
    "Iδmin=0.0415\n",
    "Gb4_max= 230.792e-9\n",
    "Gb4_min= 176.935e-9\n",
    "Iδmax= 0.0314\n",
    "Fsum=w_new**2*(Iδmin**2*Gb4_min/(δ0min*1e-3)-Iδmax**2*Gb4_max/(δ0max*1e-3))\n",
    "print('Fsum=%.2f'%Fsum)\n",
    "print('Fizm=%.2f'%Fizm)\n",
    "print('Fizm>Fsum*5..10')\n",
    "print('%.2f>%.2f..%.2f'%(Fizm,Fsum*5,Fsum*10))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ]
}