资源描述:
《abaqus中truss单元预应力的两种施加方法》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、ABAQUS中预应力Truss单元的两种实现方法 例题:100m 长钢缆水平放置从 x=0到 x=100。两端固定。无初始拉力,计算下垂量。截面; A=0.01539m2,Density: r =7800kg/m3,g=9.8m/s2,E=2.1e+11N/m2Analyticalsolutionofmaximumdisplacement(u2)atx=50m: U2_max =-((3*r*g*L^4)/(64*E))^(1/3)=-1.194944005m方法一. 沿 trusselement 加沿长度方向初始拉应力 (see
2、job-1.inp)此文件中使用了initialcondition,type=stress方法加初始拉应力。因工程上无此初应力, 更好的方法是使用降温法。算完后再升温。用降温法。算完后再升温。NOTE: 降温法施加预应力(激活钢绞线)。温度=-力/(膨胀系数*弹模*钢绞线面积)1、第一步,在truss单元中施加一个初始应力,让计算处于初始平衡状态;初始应力设置过小,可能不收敛,应多次试算,找到一个合理的应力值。一般情况下,这个初始值对最终值的影响不大,可以忽略。2、第二步,施加truss单元的自重荷载,打开非线性开关(nlgeom=YES )考虑几何非线性问
3、题;3、本例中初始值采用0.1Mpa。自重作用下缆索的拉应力约为80Mpa。最大位移为 -1.195m,与理论计算值吻合得很好。*HeadingCableapplinggravityloadwithinitailstress ThemaximumAnalyticaldisplacementwithoutinitailstress(atx=50m) U2=-1.194944005meter****Method1.Using*initialcondition,type=stressmethod***Preprint,echo=NO,model=NO,hi
4、story=NO,contact=NO*Node 1, 0., 0. 101, 100., 0.*NGEN,NSET=NALL1,101,1*Element,type=T2D21, 1, 2*ELGEN,elset=ELALL 1,100,1*ELSET,ELSET=EL_OUT1,51,100*SolidSection,elset=ELALL,material=steel 0.01539,***Nset,nset=Left1,*Nset,nset=right101
5、,*Nset,nset=mid51,** **MATERIALS** *Material,name=steel *Density7800.,*Elastic2.1e+11,0.3*initialcondition,type=stress**Note:thesolutionwillnotconvergeastheinitialstress<100,000N/m^2ELALL,100000*BoundaryLeft,1,2Right,1,2*Step,name=Step-0,inc=1000Initialstressequilibrium*Static1,
6、1.,1e-05,1.*Output,field,variable=PRESELECT*Output,history,variable=PRESELECT*Nodeprint,nset=mid,freq=1000U,*ELPRINT,ELSET=EL_OUT,freq=1000S*ENDSTEP**----------------------------------------------------------------** **STEP:Step-1** *Step,name=Step-1,nlgeom=YES,inc=1000Applygravity
7、load*Static0.01,1.,1e-05,1.**Name:GRAVITY-1 Type:Gravity*DloadELALL,GRAV,9.8,0.,-1.****OUTPUTREQUESTS** *Restart,write,numberinterval=1,timemarks=NO *Output,field,variable=PRESELECT*Output,history,variable=PRESELECT*Nodeprint,nset=mid,freq=1000U,*ELPRINT,ELSET=EL_OUT,freq=1000S*En
8、dStep方法二. 使用 STABILIZEpara