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

707390346

我在计算板振动时（矩形单元），当边界定为一边固支时，可以计算出结果
但定为四边简支时，则最后得到的总刚是奇异的，不知道怎么回事，求教一下大家
一边固支：
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

回复 1 # 707390346 的帖子

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

猜测LZ自由度没对应好, 或许先将原问题简化下, LZ自己就能发现那出错