预应力钢筋如何加在实体ANSYE模型

预应力钢筋加在实体模型的混凝土中
因为钢筋和混凝土是连接在一起的
所以一般是把钢筋自由度与混凝土自由度偶合起来
节点太多了
手工偶合几乎不可能
请问有什么简单的方法实现偶合呢?
我设想编个小命令行程序,
对钢筋的每个节点找到距离最近的混凝土节点
然后偶合起来
大家觉得如何呢?
是不是太笨了?
有更好的方法么?
请指教.
ansys我是才开始用
编此程序还有些问题
请问什么语句可以知道节点坐标呢?
求了节点距离后判断语句是什么格式呢?
请赐教
vmesh之后，混凝土与预应力筋（以及钢筋）是自动耦合在一起的，不需要在做什么耦合了，就这个问题来说，ansys是不麻烦的。

那样的话
你是不是必须把预应力钢筋和钢筋作为混凝土体积的边界线?
如果线数量少还好办
如果数量很多呢?
还会使混凝土的单元划分形状不规则甚至使网格划分失效?
有很多预应力钢筋吗？我的想法是用link8模拟钢筋，solid65模拟混凝土，然后在预应力筋的端部和solid65偶合起来，可以吗？其中link8就是边界or
切出来的线，然后如果是曲线预应力筋的话就在link8上加等效荷载，应为
solid65不能在面上加线载荷，是吗？
大家讨论一下好吗？
是有很多预应力钢筋
而且是曲线的,
你的想法很好
但是仅仅在端部偶合就够了吗?
请教一下:
你的link8上的等效荷载怎么计算出来?
link8是边界或切出来,
钢筋多的时候,怎么样切的又快又好呢?
您是怎么做的?
等效载荷很好算啊，教科书上都有的，其实具体的我也没动手，我现在只见过link8模拟普通钢筋的。

如果是切出来的，都不需要偶合我觉得，应为是一个线，然后同时又是实体的边，
我做的东西钢筋不多，但我想多的话应该也可以8，
你觉得在端部偶合有什么问题吗？
网上大家都说可以用等效载荷or温度来模拟，不知道哪位贴一个等效载荷的例子学习学习，温度的我倒见过，但是嫌它太麻烦。

网上大家都说可以用等效载荷or温度来模拟
呵呵我理解错了，
我以为你说的等效是曲线改直线后再加个等效荷载
原来你是指加预应力的方式啊。

直接加等效应变就可以了。

少数钢筋可以做为混凝土体积边界来切，如果钢筋很多的话当然也可以，
但是复杂性可不是简单的正比例哟，
想象一下很多钢筋和预应力钢筋缠绕在一起，
那么这个切出来的混凝土体将是何等复杂的形状？！
我现在建预应力混凝土梁的模型，准备用SOLID65+LINK8，两者不需要耦合吗？我想的是不是应该这样：
1 VMESH混凝土实体；
2 选择适当实体的节点，建立LINK8的线形，就不用两者耦合了？
这样对吗？急于赶做毕业论文，请各位指点。

谢谢！
只要不是无粘结钢筋，
当然可以选择适当实体的节点，建立LINK8的线形，
问题是连接的线是不是你想要的线？
误差大不大？
To:真诚
是粘结的，在实体里面找预应力筋线形是比较繁，有什么好建议吗？另外，对SOLID65的参数的取值我也不是很清楚，请赐教。

我个人的经验，
还是不要先划网格再连接节点为钢筋单元吧。

传统的方式是把预应力钢筋作为体积的边界，把混凝土体积分割开来，
glue后划分混凝土单元，边界就作为link或beam单元了。

要是嫌那样做复杂，可以考虑一下我前面提议的方法和程序，呵呵。

我的比较发现：只要混凝土单元足够密，那么那样偶合与边界切割法的结果几乎没有差别。

普通钢筋可以用Solid65的分布钢筋模拟，其实常数数据很简单。

混凝土的数据多了点，其说明在solid65理论手册里很全面，我就是看帮助学的。

data table要10个数据，但是简单计算只要输入4个就可以了，
ShrCf-Op: shear transfer coefficients for an open crack 可以取0.3（有人说按照铁路和公路规范是0.5）
ShrCf-Cl: Shear transfer coefficients for a closed crack, 可以取0.5-1.0（有人说按照铁路和公路规范是1）
UnTensSt: 单轴抗拉强度
UnCompSt: 单轴抗压强度
后面6个参数可以不管它，取缺省值，但是如果填写了其中一个，则六个就要全部填写。

帮助里有这些数据的理论意义和缺省值的取值。

solid65的Keyopt(7)=1可以帮助收敛，在cracking is imminent时。

有一篇文章里说觉得solid65单元模拟精度很高，只要各材料参数取的正确就可以十分好的模拟RC构件加载到变形开裂破坏的全过程。

与实验数据十分接近。

但是要注意混凝土的取值，弹性模量抗压强度不能取规范值，要用实验公式。

下面是个混凝土单元数据例子。

ET,2,solid65
Ex,1,3e10
NUXY,1,1/6
Dens,1,2550
TB,CONCR,1,,,
TBMODIF,2,1,0.3
TBMODIF,3,1,0.6
TBMODIF,4,1,2800000
TBMODIF,5,1,28000000
我觉得用初始应变加预应力最好。

下面是论坛里拷贝的一个例子。

a.初应力方式
b.温度方式（建议使用）
c.等效荷载法（细部不便模拟）
下面有一小例，使用solid95和link8,或许对你有参考！
!简支梁实体与预应力钢筋分析
/COM, Structural
/PREP7
egjx=2e5 !Ey
agjx=140 !单根钢绞线面积
ehnt=4e4 !Eh
xzxs=1.0e-5 !线胀系数
yjl=200000 !定义预加力
et,1,link8 !定义link8单元
et,2,solid95 !定义solid95单元
r,1,agjx !定义link8单元的面积
r,2 !定义第2种实常数
mp,ex,1,egjx !定义link8单元的弹性模量
mp,prxy,1,0.3 !定义link8单元的泊松系数
mp,alpx,1,1.0e-5 !定义线膨胀系数
mp,ex,2,ehnt !定义solid95单元的弹性模量
mp,prxy,2,0.3 !定义solid95单元的泊松系数
blc4, , ,100,200,3000 !定义梁体
/view,1,1,1,1 !定义ISO查看
/ang,1
vplot !绘制梁体
kwpave,6 !工作平面移动到关键点6
wpoff,-30 !工作平面移动-30mm(X)
wprot,0,0,90 !工作平面旋转
vsbw,1 !分割梁体
wpoff,0,0,-40 !工作平面移动-40mm(Z)
vsbw,2 !分割梁体
wpoff,0,40 !工作平面移动40mm(Y)
wprot,0,90 !工作平面旋转
vsbw,all !分割梁体
wpstyl !关闭工作平面显示
nummrg,all,,,,low !整理
numcmp,all !压缩编号
esize,30 !定义网分时边长控制
lsel,s,,,28,38,10 !定义line28和38为新的选择集latt,1,1,1 !定义选择集的属性
lmesh,all !对线划分单元
allsel,all !新的选择集为所有的实体
gplot !绘制所有的实体
vsel,s,,,all !定义所有体为选择集
vatt,2,2,2 !定义选择集的属性
mshape,0,3d !将体划分单元的形状定位HEX mshkey,1 !采用MAPPED划分器
vmesh,all !对体进行划分单元
finish
/solu
dl,3,,all !对线line7施加约束(UX,UY,UZ)
dl,16,,all !对线line31施加约束(UX,UY,UZ)
dl,23,,all !对线line23施加约束(UX,UY,UZ)
dl,2,,uy !对线line4施加约束(UY)
dl,15,,uy !对线line30施加约束(UY)
dl,22,,uy !对线line23施加约束(UY)
dk,2,,,,,ux,uy !对关键点2约束(UX,UY)
bfl,28,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度bfl,38,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度solve !求解
finish
/post1
plnsol,s,z,0,1 !绘制Z方向的应力
etable,sigi,ls,1 !定义钢筋单元数据表
plls,sigi,sigi,1 !绘制上述应力
!finish
!/exit,nosav
真诚兄：你的例子和见解很好，让我受益匪浅。

我想用ANSYS分析一榀框架的破坏情况，用SOLID65 和LINK8单元，现在有几个问题困扰我：
1.我的框架是4层，高6米，三跨，总长6米。

我担心运算量太大；请问你是否
做过这方面的分析。

有什么经验？
2.我是否需要输入应力－应变曲线？特别是非线形阶段；
3.分析结果都有哪些？如何看到截面的应力和破坏情况？
4.建模时，我想先画出实体，再切割将各个部分的交线定义为钢筋，可以吗？
在进行体网分时，网格密度如何掌握？
问题多多，麻烦你了！
真诚兄：
我还想补充一下，我分析的是普通梁柱模型试验，无预应力钢筋。

SOLID65单元本身就考虑了钢筋，是吗？我在常数定义时，不知道REBAR的几个参数如何填写？我的梁柱截面都是矩形，对称配筋。

如果在65的实常数中定义了钢筋的参数，是不是就可以不用LINK8了？谢谢
To:真诚
你说的：“把预应力钢筋作为体积的边界，把混凝土体积分割开来，”我明白它的意思，但不知“体积的边界”是什么含义，以及如何实现。

是把体的边界面作为工作面，在其上作预应力筋线形吗？
还有，对体单元进行映射划分网格时，我做的箱梁截面很不规则，不能划分，应该再将截面简化吗？
混凝土的开裂及压碎目前还不是我的课题重点，对三个方向的钢筋（主要是各自的角度）该如何取值呢？又如何区查看钢筋的布置呢（图形显示）？POLTCTRL/DEVICE OPTION 命令好像不行。

问题较多，请别介意。

谢谢！
to Cia:
1)运算量不大，我算的箱梁40米。

对称简化，机器是1.5GCPU,260M内存,
运算产生700M左右的硬盘文件。

时间10分钟。

2）钢筋如果考虑材料非线性可以用TB输入非线性模型，混凝土就用isotropic材料吧，它的非线性我没弄过。

3）一般的后处理啊，可以看到应力，plot results可以看到混凝土开列与压碎的区域，
4)可以，有点麻烦罢了。

5）网格是规则的好，密度不必太大，但是钢筋link要是多了，密度只怕小不下来。

据说solid65
单元的最佳尺寸是混凝土粗骨料尺寸的三倍，：）
补充：矩形对称配筋，不必用link了，rebar实常数就是三个正交方向的钢筋的material number,volumn ratio和定义方向，打开GUI一看就明白了。

to 阿蕾:
先输入预应力钢绞线线型，上下平移和左右平移得到两个面，用此二面切割混凝土实体，钢绞线就在混凝土的边界角线上了，
“边界切割法”都是这么做的吧。

：）
箱梁截面很不规则，是不能划分，应该再将截面简化，
能划分规则网的必须是六面体。

比较轻松的办法是先划分面，再extrude出去。

面也要切割成规则一些，还可以用lccat命令concatenates multiple lines into one line for mapped meshing,
前面好象有讨论网格划分的例子。

plot control->style->size and shape->display of moment(on)可以看见solid65 rebar的方向符号示意
rebar的方向是方向余弦值吧，
你输入几个值试一下，plot control->style->size and shape->display of element(on)
一看就明白了。