{"id":3849,"date":"2022-05-06T16:43:27","date_gmt":"2022-05-06T08:43:27","guid":{"rendered":"http:\/\/101.37.83.176\/?p=3849"},"modified":"2022-09-06T08:32:44","modified_gmt":"2022-09-06T00:32:44","slug":"python-concrete-calculation-1","status":"publish","type":"post","link":"http:\/\/101.37.83.176\/?p=3849","title":{"rendered":"Python\u4e4b\u5e38\u7528\u6df7\u51dd\u571f\u6784\u4ef6\u8ba1\u7b9701"},"content":{"rendered":"<p>\u5199\u5728\u524d\u9762\u7684\u8bdd<br \/>\n\u6700\u8fd1\u5728\u7684\u7528Jupyter Notebook\u5b66\u4e60Python\uff0c\u987a\u4fbf\u628a\u4ee5\u524dfx-FD10Pro\u505a\u7684\u6784\u4ef6\u8ba1\u7b97\u6587\u4ef6\u91cd\u5199\u6210Python\uff0c\u5c31\u5f53\u505a\u7ec3\u4e60\u4e86\u3002<\/p>\n<hr \/>\n<p>\u7a0b\u5e8f\u8bf4\u660e\uff1a<br \/>\n\u3010a\u3011\u5df2\u77e5Asp,M\uff0c\u6c42As<br \/>\n<span class=\"katex-eq\" data-katex-display=\"false\">\n\nx=h_0-\\sqrt{h_0^2-\\dfrac{2\\left[\\gamma_{RE}M-f_y^{\\prime}A_s^{\\prime}(h_0-a_s^{\\prime}\\right]}{\\alpha_1f_cb}}\n\n\\\\should,\\ x\\le\\xi_bh_0\n\n\\\\if\\ x\\ge2a^{\\prime},\\ A_s=\\dfrac{\\alpha_1f_cbx+f_y^{\\prime}A_s^{\\prime}}{f_y}\n\n\\\\if\\ x&lt;2a^{\\prime},\\ A_s=\\dfrac{\\gamma_{RE}M}{f_y(h-a_s-a_s^{\\prime})}\n\n<\/span><br \/>\n\u3010b\u3011\u5df2\u77e5As,Asp\uff0c\u6c42Mu<br \/>\n<span class=\"katex-eq\" data-katex-display=\"false\">\n\nx=\\dfrac{f_yA_s-f_y^{\\prime}A_s^{\\prime}}{\\alpha_1f_cb}\n\n\\\\if\\ ,x\\ge2a^{\\prime},\\ M_u=[\\alpha_1f_cbx(h_0-\\dfrac{x}{2})+f_y^{\\prime}A_s^{\\prime}(h_0-a_s^{\\prime})]\/\\gamma_{RE}\n\n\\\\if\\ ,x&lt;2a^{\\prime},\\ M_u=f_yA_s(h-a_s-a_s^{\\prime})\/\\gamma_{RE}\n\n<\/span><\/p>\n<hr \/>\n<pre class=\"brush: python; title: ; notranslate\" title=\"\">\r\n#coding=utf-8\r\n#\u94a2\u7b4b\u6df7\u51dd\u571f\u53c2\u6570\r\ndef c_hrb():\r\n    global fcuk,HRB,Ec,fc,ft,ftk,Es,fy,fyp,fyk\r\n    global a1,epsilon_cu\r\n    fcuk=Ec=fc=ft=ftk=0.0\r\n    HRB=Es=fy=fyp=fyk=0.0\r\n    #\u77e9\u5f62\u5e94\u529b\u56fe\u7cfb\u6570a1\uff0cC50\u4ee5\u4e0b\u4e3a1.0\r\n    a1=1.0\r\n    #\u6b63\u622a\u9762\u6df7\u51dd\u571f\u6781\u9650\u538b\u5e94\u53d8epsilon_cu\uff0cC50\u4ee5\u4e0b\u4e3a0.0033\r\n    epsilon_cu=0.0033    \r\n    fcuk=int(input(&quot;C?30\/35\/40&quot;))\r\n    if (fcuk==30):\r\n        Ec=30000\r\n        fc=14.3\r\n        ft=1.43\r\n        ftk=2.01\r\n    elif (fcuk==35):\r\n        Ec=31500\r\n        fc=16.7\r\n        ft=1.57\r\n        ftk=2.2\r\n    elif (fcuk==40):\r\n        Ec=32500\r\n        fc=19.1\r\n        ft=1.71\r\n        ftk=2.39\r\n    HRB=int(input(&quot;HRB?300\/400\/500&quot;))\r\n    if (HRB==300):\r\n        Es=210000\r\n        fy=270\r\n        fyp=270\r\n        fyk=300\r\n    elif (HRB==400):\r\n        Es=200000\r\n        fy=360\r\n        fyp=360\r\n        fyk=400\r\n    elif (HRB==500):\r\n        Es=200000\r\n        fy=435\r\n        fyp=410\r\n        fyk=500\r\n    print('fcuk=',fcuk)\r\n    print('Ec=',Ec)\r\n    print('fc=',fc)\r\n    print('ft=',ft)\r\n    print('ftk=',ftk)\r\n    print('HRB=',HRB)\r\n    print('Es=',Es)\r\n    print('fy=',fy)\r\n    print('fyp=',fyp)\r\n    print('fyk=',fyk)\r\n#\u51fd\u6570\u7ed3\u675f\r\n\r\n#1.\u53d7\u5f2f\u6784\u4ef6\u6b63\u622a\u9762\r\n#\u77e9\u5f62\u53cc\u7b4b\uff08\u975e\u6297\u9707\uff0cgamma_RE=1.0;\u6297\u9707\uff0c\u300a\u6df7\u89c4\u300b11.1.6\uff0cgamma_RE=0.75\uff09\r\n#\u3010a\u3011\u5df2\u77e5Asp,M\uff0c\u6c42As\r\ndef given_asp_m_solve_as():\r\n    import math\r\n    c_hrb()\r\n    gamma_RE=float(input(&quot;\u6297\u97070.75,gamma_RE=&quot;))\r\n    nsp=float(input(&quot;\u53d7\u538b\u7b4b\u4e2a\u6570nsp=&quot;))\r\n    dsp=float(input(&quot;\u53d7\u538b\u7b4b\u76f4\u5f84dsp=&quot;))\r\n    Asp=nsp*0.25*3.14*(dsp**2)\r\n    M=float(input(&quot;kN*m,M=&quot;))\r\n    h=float(input(&quot;mm,h=&quot;))\r\n    b=float(input(&quot;mm,b=&quot;))\r\n    ast=float(input(&quot;mm,\u53d7\u62c9\u94a2\u7b4b\u5408\u529b\u70b9\u81f3\u53d7\u62c9\u8fb9\u7f18\u7684\u8ddd\u79bbast=&quot;))\r\n    asp=float(input(&quot;mm,\u53d7\u538b\u94a2\u7b4b\u5408\u529b\u70b9\u81f3\u53d7\u538b\u8fb9\u7f18\u7684\u8ddd\u79bbasp=&quot;))\r\n    h0=h-ast\r\n    #\u754c\u9650\u53d7\u538b\u533a\u9ad8\u5ea6ksi_b\r\n    ksi_b=round(0.8\/(1 + fy \/ (Es*epsilon_cu)),3)\r\n    #\u53d7\u538b\u533a\u9ad8\u5ea6x\r\n    x=round(h0-math.sqrt(h0**2 - 2*(gamma_RE*M*(10**6) - fyp*Asp*(h0-asp)) \/ (a1*fc*b)),1)\r\n    print('h0=',h0)\r\n    print('ksi_b=',ksi_b)\r\n    print('x=',x)\r\n    if (x&lt;=ksi_b*h0):\r\n        #\u5f2f\u66f2\u7834\u574f\u65f6\uff0c\u622a\u9762\u7684\u5ef6\u6027\u53d6\u51b3\u4e8e\u53d7\u538b\u533a\u9ad8\u5ea6\u7684\u5927\u5c0f\uff0c\u53d7\u538b\u533a\u9ad8\u5ea6\u8d8a\u5c0f\uff0c\u622a\u9762\u8f6c\u52a8\u5c31\u8d8a\u5927\uff0c\u5ef6\u6027\u5c31\u8d8a\u597d\r\n        print(x,&quot;\u975e\u6297\u9707\u622a\u9762\uff0c&lt;=ksi_b*h0&quot;,ksi_b*h0)\r\n        print(x,&quot;\u6297\u9707\u7aef\u622a\u9762\uff0c\u4e00\u7ea7x&lt;=0.25h0&quot;,0.25*h0)\r\n        print(x,&quot;\u6297\u9707\u7aef\u622a\u9762\uff0c\u4e8c\u3001\u4e09\u7ea7x&lt;=0.35h0&quot;,0.35*h0)        \r\n    else:\r\n        #\u8d85\u7b4b\u6881\uff0c\u8106\u6027\u7834\u574f\r\n        print(x,&quot;&gt;ksi_b*h0,\u622a\u9762\u8fc7\u5c0f\uff01&quot;,ksi_b*h0)\r\n    if (x&gt;=2*asp):\r\n        print(x,&quot;&gt;=2*asp,\u6df7\u89c46.2.10-4\u5f0f&quot;,2*asp)\r\n        As=(a1*fc*b*x+fyp*Asp)\/fy\r\n    else:\r\n        print(x,&quot;&lt;2*asp,\u6df7\u89c46.2.14\u6761&quot;,2*asp)\r\n        As=gamma_RE*M*(10**6) \/ (fy*(h - ast - asp))\r\n    As=round(As,1)\r\n    print('\u53d7\u62c9\u7b4b\u9762\u79efAs=',As)\r\n    rho=round(As\/(b*h)*100,3)\r\n    print('\u53e6\u9700\u6ee1\u8db3\u6700\u5c0f\u914d\u7b4b\u7387\uff1b\u914d\u7b4b\u7387rho=',rho,'%')\r\n#\u51fd\u6570\u7ed3\u675f\r\n\r\n#\u3010b\u3011\u5df2\u77e5As,Asp\uff0c\u6c42Mu\r\ndef given_as_asp_solve_mu():\r\n    import math\r\n    c_hrb()\r\n    gamma_RE=float(input(&quot;\u6297\u97070.75,gamma_RE=&quot;))\r\n    nst=float(input(&quot;\u53d7\u62c9\u7b4b\u4e2a\u6570nst=&quot;))\r\n    dst=float(input(&quot;\u53d7\u62c9\u7b4b\u76f4\u5f84dst=&quot;))\r\n    As=nst*0.25*3.14*(dst**2)\r\n    nsp=float(input(&quot;\u53d7\u538b\u7b4b\u4e2a\u6570nsp=&quot;))\r\n    dsp=float(input(&quot;\u53d7\u538b\u7b4b\u76f4\u5f84dsp=&quot;))\r\n    Asp=nsp*0.25*3.14*(dsp**2)\r\n    h=float(input(&quot;mm,h=&quot;))\r\n    b=float(input(&quot;mm,b=&quot;))\r\n    ast=float(input(&quot;mm,\u53d7\u62c9\u94a2\u7b4b\u5408\u529b\u70b9\u81f3\u53d7\u62c9\u8fb9\u7f18\u7684\u8ddd\u79bbast=&quot;))\r\n    asp=float(input(&quot;mm,\u53d7\u538b\u94a2\u7b4b\u5408\u529b\u70b9\u81f3\u53d7\u538b\u8fb9\u7f18\u7684\u8ddd\u79bbasp=&quot;))\r\n    h0=h-ast\r\n    #\u754c\u9650\u53d7\u538b\u533a\u9ad8\u5ea6ksi_b\r\n    ksi_b=round(0.8\/(1 + fy \/ (Es*epsilon_cu)),3)\r\n    #\u53d7\u538b\u533a\u9ad8\u5ea6x\r\n    x=(fy*As-fyp*Asp)\/(a1*fc*b)\r\n    print('h0=',h0)\r\n    print('ksi_b=',ksi_b)\r\n    print('x=',x)\r\n    if (x&lt;=ksi_b*h0):\r\n        #\u5f2f\u66f2\u7834\u574f\u65f6\uff0c\u622a\u9762\u7684\u5ef6\u6027\u53d6\u51b3\u4e8e\u53d7\u538b\u533a\u9ad8\u5ea6\u7684\u5927\u5c0f\uff0c\u53d7\u538b\u533a\u9ad8\u5ea6\u8d8a\u5c0f\uff0c\u622a\u9762\u8f6c\u52a8\u5c31\u8d8a\u5927\uff0c\u5ef6\u6027\u5c31\u8d8a\u597d\r\n        print(x,&quot;\u975e\u6297\u9707\u622a\u9762\uff0c&lt;=ksi_b*h0&quot;,ksi_b*h0)\r\n        print(x,&quot;\u6297\u9707\u7aef\u622a\u9762\uff0c\u4e00\u7ea7x&lt;=0.25h0&quot;,0.25*h0)\r\n        print(x,&quot;\u6297\u9707\u7aef\u622a\u9762\uff0c\u4e8c\u3001\u4e09\u7ea7x&lt;=0.35h0&quot;,0.35*h0)        \r\n    else:\r\n        x=ksi_b*h0\r\n        print(x,&quot;\u53d6x=ksi_b*h0&quot;,ksi_b*h0)\r\n    if (x&gt;=2*asp):\r\n        print(x,&quot;&gt;=2*asp,\u6df7\u89c46.2.10-4\u5f0f&quot;,2*asp)\r\n        Mu=(a1*fc*b*x*(h0-x\/2) + fyp*Asp*(h0-asp)) \/ gamma_RE\/(10**6)\r\n    else:\r\n        print(x,&quot;&lt;2*asp,\u6df7\u89c46.2.14\u6761&quot;,2*asp)\r\n        Mu=fy*As*(h - ast - asp) \/ gamma_RE\/(10**6)\r\n    Mu=round(Mu,3)\r\n    print('\u6781\u9650\u6297\u5f2f\u627f\u8f7d\u529bMu=',Mu)\r\n#\u51fd\u6570\u7ed3\u675f\r\n<\/pre>\n","protected":false},"excerpt":{"rendered":"<p>\u7528Jupyter Notebook\u5b66\u4e60Python\uff0c\u987a\u4fbf\u628a\u4ee5\u524dfx-FD10Pro\u505a\u7684\u6784\u4ef6\u8ba1\u7b97\u6587\u4ef6\u91cd\u5199\u6210Python\u3002<\/p>\n","protected":false},"author":1,"featured_media":4222,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"ngg_post_thumbnail":0,"footnotes":""},"categories":[203],"tags":[222,204,223],"class_list":["post-3849","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-python","tag-jupyter-notebook","tag-python","tag-223"],"_links":{"self":[{"href":"http:\/\/101.37.83.176\/index.php?rest_route=\/wp\/v2\/posts\/3849","targetHints":{"allow":["GET"]}}],"collection":[{"href":"http:\/\/101.37.83.176\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/101.37.83.176\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/101.37.83.176\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/101.37.83.176\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=3849"}],"version-history":[{"count":0,"href":"http:\/\/101.37.83.176\/index.php?rest_route=\/wp\/v2\/posts\/3849\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"http:\/\/101.37.83.176\/index.php?rest_route=\/wp\/v2\/media\/4222"}],"wp:attachment":[{"href":"http:\/\/101.37.83.176\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=3849"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/101.37.83.176\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=3849"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/101.37.83.176\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=3849"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}