| 查看: 6121 | 回复: 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是我模拟的结果，差的也太多了···· |
回复此楼» 猜你喜欢
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有4人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
让我中一个面上吧！
已经有17人回复
-
基金不中，共勉
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
wzq43798785
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1134.3
- 散金: 80
- 红花: 2
- 帖子: 136
- 在线: 70.9小时
- 虫号: 3736696
- 注册: 2015-03-11
- 性别: GG
- 专业: 防灾工程
|
一般模型尺寸，加载方式，约束都比较好调，没什么大问题，你这个模型中有两种材料，型钢加混凝土，彼此之间直接共结点，肯定对变形有影响。还有就是混凝土的本构模型是否正确， 发自小木虫Android客户端 |
11楼2015-10-23 19:06:35
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
