| 查看: 5869 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
nayunayunayu铁虫 (小有名气)
|
[求助]
模拟与试验数据差距很大是怎么回事？？ 已有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是我模拟的结果，差的也太多了···· |
回复此楼» 猜你喜欢
一志愿上海交大生物与医药专硕324分，求调剂
已经有4人回复
-
一志愿重庆大学085700资源与环境，总分308求调剂
已经有4人回复
-
0854电子信息求调剂 324
已经有3人回复
-
303求调剂
已经有3人回复
-
336求调剂
已经有3人回复
-
306求调剂
已经有9人回复
-
工科0856求调剂
已经有4人回复
-
工科材料085601 279求调剂
已经有10人回复
-
081700 调剂 267分
已经有4人回复
-
276求调剂。有半年电池和半年高分子实习经历
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
10楼2015-10-23 13:51:21
nayunayunayu
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 672.3
- 红花: 1
- 帖子: 71
- 在线: 28.9小时
- 虫号: 2959745
- 注册: 2014-02-09
- 性别: GG
- 专业: 结构工程
2楼2015-10-22 20:32:23
wzq43798785
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1134.3
- 散金: 80
- 红花: 2
- 帖子: 136
- 在线: 70.9小时
- 虫号: 3736696
- 注册: 2015-03-11
- 性别: GG
- 专业: 防灾工程
3楼2015-10-23 00:33:16
4楼2015-10-23 00:52:24
