[求助] 模拟与试验数据差距很大是怎么回事？？ 已有3人参与

我的命令流如下： FINISH /CLEAR /filname,nayu /title,DSTC-C40-3 /PREP7 ! FRP 厚度, FRP thickness Fthick=4 ! steel 厚度, FRP thickness Sthick=5 ! 混凝土弹模, Concrete Elastic Modulus Ec=3.0e4 ! 混凝土泊松比, Concrete Poisson ratio EMUc=0.2          ! 混凝土抗压抗拉强度 ! Concrete Compressive and Tensile Strength fc=14.3 ft=1.43 ! FRP弹模和泊松比 ! Elastic Modulue and Poisson ratio for FRP Ef=2.2e4 EMUf=0.3 !钢管弹性模量 Es=2.106e5 !钢管泊松比 Esmuf=0.25 !钢管抗压抗拉强度 sc=36.3e4 st=36.3e4 fy=353.2            !常温下屈服应力 exx=2e5             !常温下杨氏模量 ! Height and Radius of Tube ! 圆柱高度和半径 H=750 Rf=121 Rs=72.5 ! 选择混凝土单元, Concrete element ET,1,SOLID65          keyopt,1,5,1                                    keyopt,1,6,3          keyopt,1,7,1 ! 选择FRP单元 ET,2,solid45    ! 输入混凝土材料参数 MPTEMP,,,,,,,,   MPTEMP,1,20 MPDATA,EX,1,,Ec MPDATA,PRXY,1,,EMUc ! 输入应力应变曲线, Concrete stress-stain curve TB,MISO,1,1, TBTEMP,20 tbpt,,164e-6,4.92 tbpt,,328e-6,9.75 tbpt,,493e-6,14.36 tbpt,,657e-6,18.57 tbpt,,821e-6,22.23 tbpt,,985e-6,25.21 tbpt,,1149e-6,27.45 tbpt,,1314e-6,28.96 tbpt,,1478e-6,29.79 tbpt,,1642e-6,30.05 tbpt,,1970e-6,28.75 tbpt,,2463e-6,24.51 ! 输入破坏曲面, Failure envelop surface TB,CONC,1,1,9, TBTEMP,20 TBDATA,,0.4,0.9,ft,-1,, TBDATA,,,,1,,, MPTEMP,,,,,,,, MPTEMP,1,20 ! 输入FRP材料参数 MPDATA,EX,2,,Ef MPDATA,PRXY,2,,EMUf TB,bkin,2,1,2 TBTEMP,20 TBDATA,,430,0 ! 输入钢管材料参数 MPDATA,EX,3,,Es MPDATA,PRXY,3,,Esmuf                                  TB,BISO,3,1,2, 1                                             !双线性各向同性强化模型 TBDATA,,235,2.106E4 ！建立组合柱模型 pcirc,Rs-Sthick,Rs,0,360    !钢管 pcirc,Rf,Rf+Fthick,0,360    !frp管 pcirc,Rs,Rf,0,360          !混凝土 k,100,0,0,0 k,101,0,0,750 l,100,101 vdrag,all,,,,,,25 LDELE,25,,,1 numcmp,all vglue,all numcmp,all !赋予材料属性 VSEL,S,,,1                  !FRP管    2 VATT,2,0,2 VSEL,S,,,3                  !混凝土   1 VATT,1,0,1 VSEL,S,,,2                 !型钢      3 VATT,3,0,2 allsel,all ！划分网格 vsweep,1,1,11,1                 !型钢       vsweep,2,2,20,1                 !FRP管                                 vsweep,3,21,22,1               !混凝土                                        numcmp,all                       !压缩节点编号 nummrg,all ALLSEL,ALL                                                                               !************************* csys,1 et,4,184 keyopt,4,1,1 nsel,s,loc,z,H nsel,r,loc,x,Rs-Sthick,Rf n,350000,0,0,H+10 type,4 *get,nnum,node,0,count *get,ND,node,0,num,min *do,i,2,nnum !!节点个数是nnum，只需要生成nnum个mpc单元 E,350000,ND ND=NDNEXT(ND) *enddo allsel !************************** !施加荷载 !************************** F,350000,FZ,-2702920 !************************** !施加自由度约束 !************************** nsel,s,loc,Z,0 d,all,ux,0,,,,uy,uz allsel,all FINISH !*************************** !进入求解器 !*************************** /SOLU !求解选项设置 ANTYPE,0            !选择分析类型 NLGEOM,1            !计入大变形效应 NSUBST,30,100,20,on   !子步数 AUTOTS,ON            !自动时间步 NEQIT,50            !平衡迭代次数 CNVTOL,U,,0.015    !设置收敛准则 ncnv,2         !如果不收敛时结束而不退出 OUTRES,ALL,ALL      !输出所有子步的结果 LNSRCH,1                  !用Newton-Raphson方法激活线性搜索 cutcontrol,nioterpred,0             !在非线性解析中控制步长的削减 solve 模拟完成！！！！所得结果如下 A1（实验数据）   A2（模拟数据）        荷载        位移         荷载               位移 200        2.9            200           0.096 400        3.91         400                0.193 600        4.76          600           0.293 800        5.56                800          0.398 1000 6.39        1000        0.511 1100 6.81        1100        0.572 1200 7.37        1200        0.639 1360 8.57        1360        0.867 1460 9.75        1460        0.94 1560 11.21      1560        1.091 1660 13.42      1660        1.161 1920  25.74        1920        2.156 试验肯定是没有问题的，但是有限元模拟为什么跟试验差的这么多，问题出在哪里呢？我也检查不出来呀，求教各位大神！！！ A1是师兄试验所得荷载-位移曲线，A2是我模拟的结果，差的也太多了····

回复此楼

» 猜你喜欢

• 有多少人是今天查系统知道结果的？ 已经有17人回复
• 为什么资助数各大高校都创新高，自己申请怎么就这么难 已经有3人回复
• 基金未中，这种答复是模板吗？ 已经有4人回复
• 导师吐槽：我怎么摊上了这么个极品研究生！ 已经有5人回复
• 怎么看青基中了没有啊 已经有5人回复
• 梦想 已经有3人回复
• 面上合作单位盖章 已经有3人回复
• 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… 已经有5人回复
• 基金不中，共勉 已经有7人回复
• 为什么 国际(地区)合作与交流项目 没有放榜？ 已经有11人回复

» 本主题相关商家推荐: (我也要在这里推广)