马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
您需要 登录 才可以下载或查看，没有账号？我要加入
x
本帖最后由 fflyer 于 2011-9-14 22:28 编辑
首先，我对悬索进行找形分析，之后进行模态分析，但是得到的频率值确非常小，不知道错在哪儿了，还请大侠们指点！
找形命令流：
/CLEAR /FILNAME,EX820 /PREP7 ! 1.定义儿何参数和荷载参数等，定义单元类型和材料性质 L0=120 XH=20 AREA=7.016E-4 EM=7.89E10 !定义几何参数、面积、弹性模量 Q0=65 QF=10000 !定义索单位重量（N/m）和集中荷载（N） H0=9000 !定义自重作用下的水平张力(已知) ERR0=1/1000 ENUM=60 ISTRAN=1.0E-6 !定义迭代条件，单元数目、初应变 ET,1,LINK10 R,1,AREA,ISTRAN !定义单元和实常数 MP,EX,1,EM MP,PRXY，1，0.3 !定义材料性质 MP,DENS,1,Q0/AREA !采用换算密度，且为N/m3单位 ! 2.在弦线位置创建模型，施加约束和白重荷载 K,1 K,2,L0,-XH L,1,2 !创建几何模型 LESIZE,ALL,,,ENUM LMESH,ALL !生成有限元模型 D,NODE(0,0,0),ALL D,NODE(L0,-XH,0),ALL !施加约束 NODE1=NELEM(ENUM/2,1) !获得中间单元的两个节点号 NODE2=NELEM(ENUM/2,2) !以备后面使用 ACEL,,1.0 FINISH !施加值为1的加速度(自重) ! 求解、进人后处理获得索内.力、更新有限元模型、判别收敛条件是否满足 PASS1=1 *DOWHILE,PASS1
/SOLU ANTYPE,0 NLGEOM,ON SSTIF,ON !定义静态求解、打开大变形与应力刚度选顶 NSUBST,20 OUTRES,ALL,ALL SOLVE FINISH !定义子步数、输出结果、求解 /POST1
SET,LAST,LAST !进入后处理，选择最后荷载步的最后子步 *GET,NFOR,ELEM,ENUM/2,SMISC,1 !获得跨中单元的索力并计算其余弦 COSREF=(NX(NODE2)-NX(NODE1))/DISTND(NODE1,NODE2) NFOR=NFOR*ABS(COSREF) !计算跨中单元的水平张力（可用其他单元) ERR1=ABS(NFOR-H0)/H0 FINISH !计算当前索水平张力误差 /PREP7 !进人前处理，更新有限元模型 *IF,ERR1,LT,0.05,THEN !如果误差小于5%时 UPGEOM,0.1,LAST,LAST,EX820,RST *ELSE !模型更新系数采用较小数值 UPGEOM,1,LAST,LAST,EX820,RST *ENDIF !否则模型更新系数采用较大数值 *IF,ERR1,LT,ERR0,EXIT *ENDDO ！如满足迭代条件，退出循环
找形后进行模态分析：
/PREP7
ANTYPE,MODAL
MODOPT,SUBSP,10
MXPAND,10
/SOLU
SOLVE
FINISH
得到的一阶频率值为0.8E-06，后面的也是负的5次数量级的，不明白是什么原因，高手帮忙看看啊！！！
|
[复制链接]
发表于 2011-9-14 22:24
提升卡
变色卡
显身卡
楼主
发表于 2011-9-29 23:59
