[求助]模态分析的振型为什么会这样？

zglecsi

我用ansys对碳纤维加固梁进行模态分析，得出的振型如图所示：
其中齿状的单元是钢筋，怎么会这样呢？请大侠指教！

AaronSpark

把问题描述的具体点

zglecsi

我对一碳纤维加固梁进行模态分析，计算结束后，得到的各阶频率如下表所示，
SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE
1 0.19034E-02 1 1 1
2 0.19273E-02 1 2 2
3 0.19360E-02 1 3 3
4 0.19469E-02 1 4 4
5 0.19532E-02 1 5 5
6 0.19690E-02 1 6 6
7 0.19755E-02 1 7 7
8 0.19955E-02 1 8 8
9 0.20002E-02 1 9 9
10 0.20124E-02 1 10 10
然后我看各阶振型，形状就如图所示，通过单元查看，我发现其中齿状的部分是钢筋单元，这是怎么回事呢？请大侠指教。



附命令流：
finish
/clear,nostart
/config,nres,2000
!碳纤维布加固钢筋混凝土梁有限元分析
/com,structural
/prep7
!定义单元类型
et,1,link8
et,2,solid65
et,3,plane42
et,4,shell41

!定义实常数
r,1,283.8704, , !受拉钢筋1
r,2,199.9996, , !受拉钢筋2
r,3 !混凝土
r,4,0.33 !碳纤维布

!定义材料属性
mp,ex,1,2e5 !受拉钢筋1
mp,prxy,1,0.3
mp,dens,1,7.83e-9
tb,bkin,1,1,2,1
tbdata,,344.75

mp,ex,2,2e5 !受拉钢筋2
mp,prxy,2,0.3
mp,dens,2,7.83e-9
tb,bkin,2,1,2,1
tbdata,,551.6

mp,ex,3,3.6e4 !混凝土
mp,prxy,3,0.2
mp,dens,3,2.5e-9
tb,conc,3,1,9
tbdata,,0.4,1,2.04,-1

mp,ex,4,2.3e5 !碳纤维布
mp,prxy,4,0
mp,dens,4,1.35e-9
tb,bkin,4,1,2,1
tbdata,,3800

/pnum,node,1
/pnum,elem,1

!产生所有的节点
n,1
n,11,380
fill,1,11
ngen,16,11,1,11,1,,32
ngen,123,1000,1,176,1,,,－50
/view,1,1,1,1
nplot

!纵向受拉钢筋1单元建立
type,1
real,1
mat,1
*do,ii,14,121014,1000
e,ii,ii+1000
*enddo
*do,ii,20,121020,1000
e,ii,ii+1000
*enddo
!纵向受拉钢筋2单元建立
type,1
real,2
mat,2
*do,ii,13,121013,1000
e,ii,ii+1000
*enddo
*do,ii,21,121021,1000
e,ii,ii+1000
*enddo
!碳纤维布单元的建立
type,4
real,4
mat,4
*do,ii,3,8,1
*do,jj,ii,121000+ii,1000
e,jj,jj+1,jj+1000
e,jj+1,jj+1000,jj+1000+1
*enddo
*enddo

/view,1,1,1,1
/pnum,elem,0
/pnum,node,0
/eshape,1
eplot

!混凝土单元的建立
type,3
k,1
k,2,380
k,3,380,480
k,4,,480
a,1,2,3,4

lsel,s,loc,y,0
lsel,a,loc,y,480
lesize,all,,,10
lsel,all
lsel,s,loc,x,0
lsel,a,loc,x,380
lesize,all,,,10

amesh,all

type,2
real,3
mat,3
extopt,esize,122
extopt,aclear,1
vext,all,,,,,6100

/pnum,mat,1
/pnum,node,0
/pnum,elem,0
/view,1,1,1,1
eplot
allsel，all
nummrg,all
numcmp,all
eplot

/SOLU
allsel,all
ANTYPE,MODAL
MODOPT,SUBSP,15
nsel,s,loc,y,0
nsel,r,loc,z,150
d,all,uy
d,all,uz
d,all,ux

nsel,s,loc,y,0
nsel,r,loc,z,5950
d,all,uy
d,all,ux

ASEL,S,LOC,X,0
DA,ALL,UX

ASEL,S,LOC,X,380
DA,ALL,UX

allsel,all
MXPAND,15
SOLVE
FINISH

usefully

频率太小了，不正常，你看看单位有没有搞错。

hmilyxfz

频率的大小一看就知道不对劲吧，不管是什么单位。看看是不是中间出错了哦！！

sysh320

　　可能是耦合出了问题

zglecsi

谢谢各位大侠的解答，我又做了一遍，把单位统一成 N，Kg，M。还是不对，结果如下：
SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE
1 0.79338E-04 1 1 1
2 0.79682E-04 1 2 2
3 0.79971E-04 1 3 3
4 0.81474E-04 1 4 4
5 0.82311E-04 1 5 5
6 0.83425E-04 1 6 6
7 0.83660E-04 1 7 7
8 0.84176E-04 1 8 8
9 0.86321E-04 1 9 9
10 0.87037E-04 1 10 10

我分别看了混凝土单元，钢筋单元和碳纤维布单元的一阶振型，如下：
混凝土单元：

钢筋单元：
碳纤维布单元：


到底是怎么回事呢？

zglecsi

哪位大哥大姐能帮小弟解惑啊？万分感谢啊！！
是不是耦合的问题？正如sysh320兄所说。
但是，我的模型好像是共用节点啊？跟耦合有关吗？？

linqus

从图形上看，可能是节点的耦合问题；
建议，
（1）改用弹性材料，做特征值分析看看；
（2）改用弹性材料，并在上表面施加一个均布面力，看看变形情况，可以以此初步判断耦合是否正确；

ps：没有仔细看apdl，探讨。

zglecsi

多谢各位大侠的指点！问题基本解决！
我重新看了一下模型，发现梁横截面混凝土单元的划分和钢筋单元与碳纤维单元所在的位置不一致。重新划分网格后，问题解决。