ABAQUS中Truss单元预应力的两种施加方法.docVIP

ABAQUS中Truss单元预应力的两种施加方法 ABAQUS中Truss单元预应力的两种施加方法 PAGE / NUMPAGES ABAQUS中Truss单元预应力的两种施加方法 v1.0 可编写可改正 ABAQUS中预应力 Truss 单元的两种实现方法 例题： 100 m 长钢缆水平搁置从 x = 0 到 x = 100 。两头固定。无初始拉力，计算 下垂量。 截面； A = m 2 ，Density: r =7800 kg/m3, g = m/s 2,E=+11 N/m2 Analytical solution of maximum displacement (u 2 ) at x = 50m : U2_max = -((3*r*g*L^4)/(64*E))^(1/3) = m 方法一 . 沿 truss element 加沿长度方向初始拉应力 (see 此文件中使用了 initial condition, type = stress 方法加初始拉应力。因工 程上无此初应力 , 更好的方法是使用降温法。算完后再升温。用降温法。算完后再升温。 NOTE: 降温法 施加预应力（激活钢绞线）。温度 =- 力 / （膨胀系数 * 弹模 * 钢绞线 面积） 1、第一步，在 truss 单元中施加一个初始应力，让计算处于初始均衡状态；初 始应力设置过小，可能不收敛，应多次试算，找到一个合理的应力值。一般状况 下，这个初始值对最后值的影响不大，能够忽视。 2、第二步，施加 truss 单元的自重荷载，翻开非线性开关（ nlgeom=YES ）考 虑几何非线性问题； 3、本例中初始值采纳。自重作用下缆索的拉应力约为 80Mpa。最大位移为 m， 与理论计算值符合得很好。 *Heading Cable appling gravity load with initail stress The maximum Analytical displacement without initail stress (at x = 50 m) U2 = meter ** Method 1. Using * initial condition,type = stress method *Preprint, echo=NO, model=NO, history=NO, contact=NO *Node 1, 0., 0. 1 v1.0 可编写可改正 101, 100., 0. *NGEN, NSET = NALL 1, 101, 1 *Element, type=T2D2 1, 1, 2 *ELGEN, elset = ELALL 1, 100, 1 *ELSET,ELSET=EL_OUT 1, 51, 100 *Solid Section, elset=ELALL, material=steel , ** *Nset, nset=Left 1, *Nset, nset=right 101, *Nset, nset=mid 51, ** MATERIALS *Material, name=steel *Density 7800., *Elastic +11, *initial condition, type = stress Note: the solution will not converge as the initial stress 100,000 N/m^2 ELALL, 100000 *Boundary Left, 1, 2 Right, 1,2 *Step, name=Step-0, inc=1000 2 v1.0 可编写可改正 Initial stress equilibrium *Static 1, 1., 1e-05, 1. *Output, field, variable=PRESELECT *Output, history, variable=PRESELECT *Node print, nset = mid, freq = 1000 U, *EL PRINT, ELSET=EL_OUT, freq = 1000 S *END STEP STEP: Step-1 ** *Step, name=Step-1, nlgeom=YES, inc=1000 Apply gravity load *Static , 1., 1e-05, 1. ** Name: GRAVITY-1 Type: Gravity *Dload ELALL, GRAV, , 0., -1. ** OUTPUT REQUESTS *Restart, write, number interval=1, time marks=NO *Output, fiel