707390346 发表于 2012-3-19 22:22

【求助】关于弯曲薄板简支约束问题

我在计算板振动时(矩形单元),当边界定为一边固支时,可以计算出结果
但定为四边简支时,则最后得到的总刚是奇异的,不知道怎么回事,求教一下大家
一边固支:
disp(1:jdx*jdy,1:3)=1;
constraints=1:jdx:jdx*jdy;
disp(constraints,:)=0;
dof=0;                   %degree of freedom
for ni=1:jdx*jdy
    for nj=1:3
      if disp(ni,nj)~=0
            dof=dof+1;
            disp(ni,nj)=dof;
      end
    end
end
四边简支:
constraints1=1:jdx;    % constraints of x=0
constraints2=1:jdx:jdx*jdy;% constraints of y=0
constraints3=111:jdx*jdy;    % constraints of x=a
constraints4=11:jdx:121;   % constraints of y=b
disp(constraints1,1:2)=0;
disp(constraints2,1:2:3)=0;
disp(constraints3,1:2)=0;
disp(constraints4,1:2:3)=0;
dof=0;                   %degree of freedom
for ni=1:jdx*jdy
    for nj=1:3
      if disp(ni,nj)~=0
            dof=dof+1;
            disp(ni,nj)=dof;
      end
    end
end

ChaChing 发表于 2012-4-18 01:05

回复 1 # 707390346 的帖子

jdx=?, jdy=?
四边简支时,则最后得到的总刚是奇异的? 总刚=disp?
这样的问法可能不够明确, 高手们怎出手帮!?:@)

猜测LZ自由度没对应好, 或许先将原问题简化下, LZ自己就能发现那出错
页: [1]
查看完整版本: 【求助】关于弯曲薄板简支约束问题