如何用约束方程施加预应力

dongd

对于一 有限元梁杆组合结构，下面的例子是一长10米的简支梁 截面0.5*1 在距下缘0.2处设置钢束，只是定性计算 暂时未考虑预应力损失及有效截面积因素。

现在问题是运行下面的命令流后 发现梁体的弯矩和轴力和理论的相差，且出现横向变形，以前没接触到约束方程，希望有了解的dx给提个建议。
问题关键是：考虑预应力对梁体产生的弯矩和轴力
命令流需完善的地方请多多指教，本人还属于摸索阶段，在此 道谢了
/prep7
et,1,beam188
keyopt,1,3,2   
et,2,link10
mp,ex,1,3.45e10
mp,dens,1,1000
mp,prxy,1,0.1667
MP,alpx,1,1.0e-5
mp,ex,2,2e11   
mp,prxy,2,0.3
!mp,dens,2,8500
!MP,alpx,2,1.2e-5
r,2,0.01,0.0001
SECTYPE,   1, BEAM, RECT, s1, 0
SECOFFSET, CENT
SECDATA,0.5,1,5,5,0,0,0,0,0,0   
N,1,0,,,,,,
N,21,10,,,,,,   
n,50,0,3
!*  
FILL,1,21,19, ,1,1,1,1,  
type,1
r,1,
real,1
mat,1
*do,i,1,20,1
e,i,i+1,50
*enddo
n,101,0,-0.3
n,121,10,-0.3
FILL,101,121,19, ,1,1,1,1,
type,2
real,2
mat,2
*do,i,101,120,1
e,i,i+1
*enddo
d,1,ux
d,1,uy
d,1,uz
d,21,uy
d,21,uz
!/eshap,1
*do,i,1,21,1
cp,next,ux,i,i+100  !不考虑滑移
cp,next,uy,i,i+100
cp,next,uz,i,i+100
*enddo
*do,i,2,20,1
ce,next,0,i,rotz,0.3,100+i,ux,-1
*enddo

dongd

个人认为，是约束方程的原因，经计算结构变形明显不对称
ce,next,0,i,rotz,0.3,100+i,ux,-1
！rotz(i)*0.3=ux(j)
其中，i为梁上节点，j为对应钢绞线节点 0.3为偏心距

[ 本帖最后由 dongd 于 2008-2-20 22:07 编辑 ]

dongd

建立杆件与梁体的约束方程时，不仅仅是预应力引起梁体单元的转动，结构自身重力也有影响

dongd

我试试用刚臂的形式来连接预应力杆件和梁体单元看情况如何，我做上面的例子想以此来建立梁杆结构在进行整体分析时简化预应力，但是以前没做过，大家有做过的没呢 ，指点指点～