辛苦研究,也来一个吧。COMBIN40单元
双弹簧加阻尼器还有间隙的组合单元,可以模拟隔震支座(我的课题)。该单元单自由度(平动或扭转),模拟时可以在两节点间三个方向分别设置。最后可得到单元双线性滞回曲线。
1.输入:
在两节点间直接设置单元;指定自由度(UX,UY或UZ等);
实常数有:
k1,第一弹簧刚度;
k2,第二弹簧刚度;某个刚度取0.0时则不计该弹簧。
C,单元阻尼系数;
M,单元质量,如TMD质量,也可为0.0;
GAP,间隙,为正值,则有一间隙存在;为负值,则存在初始间隙;为0.0,则无间隙;
FSLIDE,滑移力,相当于最大静摩擦力。为正时,克服滑移力后弹簧1内力保持不变;为负值时一旦克服滑移力,则弹簧1内力降为0;为0.0时则无滑移,即相当于两弹簧并联。
选项:
keyopt(1)=0,标准间隙功能;=1,间隙闭合后则紧锁。
keyopt(3)=0-8,选自由度;
keyopt(4)=0,单元信息输出;=1,不输出。
keyopt(6)=0-3,设置质量分配,是在节点间等分还是在某个节点。
2.输出
计算完成后,除了和其他单元一样得到节点信息外,还可得到该单元计算信息(很重要的)。包括F1(弹簧1内力);F2(弹簧2内力);STAT(间隙闭合信息、节点相对位置信息);OLDST(前一荷载步时STAT值);STR1(弹簧1相对位移);STR2(弹簧2相对位移,也即单元两节点间相对位移);UI,UJ(两节点位移);SLIDE(滑移量)。
3.重要信息提取
(以下假设模型中combin40单元编号为32,两节点号为7和30)
节点位移:NSOL,2,7,U,X,UX_7
节点应力:ESOL,3,32,7,S,X,SX_7
节点力:ESOL,4,32,7,F,X,FX_7
节点弯矩: ESOL,5,32,7,M,Y,MY_7
用sequence numbers方法可得到:
弹簧1内力:ESOL,6,32,,SMIS,1,F1 !必须是1;
弹簧1相对位移:ESOL,7,32,,NMIS,3,STR1 !必须是3;
弹簧2内力:ESOL,8,32,,SMIS,2,F2 !必须是2;
弹簧2相对位移:ESOL,9,32,,NMIS,4,STR2 !必须是4;
简单示例:
(瞬态分析完成后)
finish
/post26
allsel,all
ESOL,2,32,,SMIS,1,F1 !定义变量,弹簧1内力;
ESOL,3,32,,SMIS,2,F2 !定义变量,弹簧2内力;
ESOL,4,32,,NMIS,4,STR2 !定义变量,弹簧2相对位移,即节点相对位移;
*DIM,F1,ARRAY,1024,1,1 !定义数列,1024为地震波点数;
VGET,F1(1),2 !将变量转入数列;
*CFOPEN,F1,DAT, !写入数据文件;
*VWRITE,F1(1)
(e8.4)
*CFCLOS
(另两个数据写法相同)
SAVE
FINISH
上面将要得到的结果写入dat文件里面,通过origin软件可方便处理数据。
F1+F2即为节点内力,STR2即为节点相对位移,则可绘制单元的滞回曲线,分析其与支座性能差异。 |