中华钢结构论坛                                             

 
*计算钢管混凝土结构的弯矩-曲率曲线的程序或软件
huangqian2005






积分 51
帖子 51
2012-6-30 17:21    顶部


谁有能计算钢管混凝土柱的弯矩-曲率曲线的软件或小程序,麻烦告知一下,急啊!!先谢了,可以发到我qq邮箱:56781421@qq.com
xgzhang1986






积分 8
帖子 8
2012-7-6 20:55    顶部
是啊 我也需要  也给我发一份吧 1139412903@qq.com,谢谢!
游子






积分 1570
帖子 513
2012-10-6 11:19    QQ 顶部
这个程序一共不超过100行。
当时我曾想过放到博士论文里面,可是太简单,怕人笑话,才取消了这个念头。
如果楼上二位是研究生,不妨自己手工编写。你会发现:真的很简单!
KEEROE






积分 3
帖子 5
2012-11-30 11:34    顶部
0 分
最近一直在学钢管混凝土,也 需要这个这个小程序哈,希望柚子大哥发一个,251427352@qq.com,非常感谢哈
zongkun






积分 106
帖子 102
2012-11-30 13:57    顶部
0 分
试试XTRACT
whisperstone






积分 7
帖子 8
2013-4-9 16:40    顶部
0 分
请游子前辈给我发一份吧,万分感谢!wangyawen321@126.com
积润






积分 0
帖子 1
2016-7-2 16:01    顶部
0 分
最近我一直在学钢管混凝土,也 需要这个这个小程序哈,希望柚子大哥发一个,826811076@qq.com,非常感谢哈
游子






积分 1570
帖子 513
2017-3-9 03:10    QQ 顶部
以下是矩形钢管混凝土构件受纯弯曲的程序,Fortran编写。把程序里面的钢材、混凝土的应力-应变关系改动一下,就可以为己所用了;再稍加改动,也能用于圆形钢管混凝土构件和型钢混凝土构件。

程序是十多年前写的,疏漏之处,请多多指正!

致在校的研究生:程序是一个手段,不必过分看重它。

c*******************************************************************
c  This program is for predicting the M-Um response of rectangular
c              Concrete-filled steel tubular columns under pure bending            

c  NOTE: 1. Stress & strain: Compression(+);Tension(-).  
c       2. Units: E--MPa;       fy--MPa;    Strain---1E-6;
c                 H,B,t---mm;   F---kN;  
c                 Moment--kN*m; Curvature---1/mm   
c       3. The cross-section is divided by 100 along the height.
c            
c
c*******************************************************************               
      Program main
        Open(1,fiLe='bending-result.dat',status='unknown')
        Write(*,*)'Please input fy,fcy,Ec,H,B,t'
        Write(1,*)'fy,fu,Es,fcy,H,B,t,L'
        Read(*,*)fy,fcy,Ec,H,B,t,L
c       
        Do 30 i=1,30
        Um=0.5*i
        Curvature=3.14*3.14*Um/(L**2)
c
        Do 20 K=1,10000
              SnCR=-20000+K*2
c**********************************************************
c            Stress of steel in bottom flange
c**********************************************************               
               Sny=fy/Es*1E6  
                         Distflange=0.5*H-0.5*t
                         Snbot=SnCR-Distflange*Curvature*1E6
c
               SNy=fy/Es*1000000.
               SN1=10*SNy
                 SN2=50*SNy
                 Snsn=abs(Snbot)
            If(Snsn.Le.SNy)then
               SSs=Es*Snsn/1000000.
          Elseif(Snsn.GT.SNy.AND.Snsn.LE.SN1)THEN  
               SSs=fy
          Elseif(Snsn.GT.SN1.and.Snsn.LE.SN2)then
               SSs=fu-(fu-fy)*(SN2-Snsn)**2/(SN2-SN1)**2
          Elseif(Snsn.GT.SN2)THEN
               SSs=fu
          Endif
c
          If(Snbot.LT.0)then
                 SSsbot=-SSs
          Endif   
                Fbot=SSsbot*t*B/1000.
            Abot=Fbot*(-Distflange)/1000.
c*****************************************************
c           Stress of steel in top flange
c*****************************************************
              Sntop=SnCR+Distflange*Curvature*1E6
c
                Snsn=Sntop
          If(Snsn.Le.SNy)then
               SSs=Es*Snsn/1000000.
          Elseif(Snsn.GT.SNy.AND.Snsn.LE.SN1)THEN  
               SSs=fy
          Elseif(Snsn.GT.SN1.and.Snsn.LE.SN2)then
               SSs=fu-(fu-fy)*(SN2-Snsn)**2/(SN2-SN1)**2
          Elseif(Snsn.GT.SN2)THEN
               SSs=fu
          Endif
            If(Sntop.GT.0)then
                 SSstop=SSs
c**************************************************************
c     Attention:The following 2 lines give birth to errors!!
c          Else
c                     SSstop=-SSs
c**************************************************************
                Endif                
c
              Ftop=SSstop*t*B/1000.
              Atop=Ftop*Distflange/1000.
c**************************************************************
c           Stress of steel fibers in webs
c**************************************************************
              Fsweb=0.
              Fcweb=0.
              Asweb=0.
              Acweb=0.
           Do 10 J=1,100
              B0=B-2*t
              H0=H-2*t
                  DistJ=(J-50-0.5)*(H0/100.)
            SnJ=SnCR+DistJ*curvature*1E6
c
              Snsn=abs(SnJ)
          If(Snsn.Le.SNy)then
              SSs=Es*Snsn/1000000.
          Elseif(Snsn.GT.SNy.AND.Snsn.LE.SN1)THEN  
              SSs=fy
          Elseif(Snsn.GT.SN1.and.Snsn.LE.SN2)then
              SSs=fu-(fu-fy)*(SN2-Snsn)**2/(SN2-SN1)**2
          Elseif(Snsn.GT.SN2)THEN
              SSs=fu
          Endif
          If(SnJ.LT.0)then
              SSsweb=-SSs
            Else
                 SSsweb=SSs     
          Endif   
            FsJ=SSsweb*2*t*H0/100/1000
              Fsweb=Fsweb+FsJ
              AsJ=FsJ*DistJ/1000.
              Asweb=Asweb+AsJ
c************************************************************
c          Stress of concrete fibers in the web
c************************************************************
            Ac=(B-2*t)*(H-2*t)
              As=B*H-Ac
              Conf=As*fy/(Ac*0.85*fcy)
           SN0=1300+(1300+18*fcy)*(Conf**0.2)
           fc0=0.85*fcy*(1+(1-0.65*(H/B-1)**3)*(0.37*(conf**0.82)
     $           -0.19*conf)/(1+0.05*(fcy/50)**5.65))
c
            x=SnJ/Sn0
           IF(x.LE.0)then
                ft=0.61*(fcy**0.5)
                Snt=ft/Ec*1000000
              IF(SnJ.GT.(-Snt))then
                    SSc=Ec*SnJ/1000000.
            Else
                SSc=-ft*((-Snt/SnJ)**0.4)           
                  Endif
           Elseif(x.LE.1)THEN
                Y=1.6*X/(0.6+X)
              SSc=fc0*Y
           ELSE
                SSc=fc0
         ENDIF
                    FcJ=SSc*B0*H0/100/1000.
                    Fcweb=Fcweb+FcJ
                AcJ=FcJ*DistJ/1000.
                    Acweb=Acweb+AcJ
10    Continue   
                  F=Fsweb+Fcweb+Fbot+Ftop
              Amoment=Asweb+Acweb+Abot+Atop
              stif=Amoment/(Curvature*1000)
c*****************************************************
                 If(Um.LE.4)then
                    If(abs(F).LE.2)then       
                  Write(1,*)Um,Amoment,stif,SnCR,Sntop
                  Write(*,*)Um,Amoment,stif,SnCR,Sntop
              Endif
           Else
                If(abs(F).LE.0.1)then
                          Write(1,*)Um,Amoment,stif,SnCR,Sntop
                                        Write(*,*)Um,Amoment,Stif,SnCR,Sntop                           
              Endif
                 Endif                
c*****************************************************
20    Continue
30    Continue  
      Write(*,*)'Calculation successful, well done!'
      End
 


上海市通信管理局
沪ICP备020910
联系我们 版权声明 发帖规则 论坛制度 关于我们
版权所有.中华钢结构论坛.永久保留所有权利

[Processing Time] User:0.28, System:0.03, Children of user:0, Children of system:0