作业帮matlab 矩阵求方程这是一个矩阵位移法的题,如果6x6的刚度矩阵K我都知道了,但是方程是KU=F,U=[Ux Uy Q Vx Vy P]',F=[Fx Fy M1 Px Py M2]',比如其中Ux Uy Q Px Py M2 都已知,我怎么求其他的数啊,就是对应的位移
来源：学生作业帮助网 编辑：作业帮 时间：2024/05/04 21:05:11
![matlab 矩阵求方程这是一个矩阵位移法的题,如果6x6的刚度矩阵K我都知道了,但是方程是KU=F,U=[Ux Uy Q Vx Vy P]',F=[Fx Fy M1 Px Py M2]',比如其中Ux Uy Q Px Py M2 都已知,我怎么求其他的数啊,就是对应的位移](/uploads/image/z/2542287-39-7.jpg?t=matlab+%E7%9F%A9%E9%98%B5%E6%B1%82%E6%96%B9%E7%A8%8B%E8%BF%99%E6%98%AF%E4%B8%80%E4%B8%AA%E7%9F%A9%E9%98%B5%E4%BD%8D%E7%A7%BB%E6%B3%95%E7%9A%84%E9%A2%98%2C%E5%A6%82%E6%9E%9C6x6%E7%9A%84%E5%88%9A%E5%BA%A6%E7%9F%A9%E9%98%B5K%E6%88%91%E9%83%BD%E7%9F%A5%E9%81%93%E4%BA%86%2C%E4%BD%86%E6%98%AF%E6%96%B9%E7%A8%8B%E6%98%AFKU%3DF%2CU%3D%5BUx+Uy+Q+Vx+Vy+P%5D%27%2CF%3D%5BFx+Fy+M1+Px+Py+M2%5D%27%2C%E6%AF%94%E5%A6%82%E5%85%B6%E4%B8%ADUx+Uy+Q+Px+Py+M2+%E9%83%BD%E5%B7%B2%E7%9F%A5%2C%E6%88%91%E6%80%8E%E4%B9%88%E6%B1%82%E5%85%B6%E4%BB%96%E7%9A%84%E6%95%B0%E5%95%8A%2C%E5%B0%B1%E6%98%AF%E5%AF%B9%E5%BA%94%E7%9A%84%E4%BD%8D%E7%A7%BB)
matlab 矩阵求方程这是一个矩阵位移法的题,如果6x6的刚度矩阵K我都知道了,但是方程是KU=F,U=[Ux Uy Q Vx Vy P]',F=[Fx Fy M1 Px Py M2]',比如其中Ux Uy Q Px Py M2 都已知,我怎么求其他的数啊,就是对应的位移
matlab 矩阵求方程
这是一个矩阵位移法的题,如果6x6的刚度矩阵K我都知道了,但是方程是KU=F,U=[Ux Uy Q Vx Vy P]',F=[Fx Fy M1 Px Py M2]',比如其中Ux Uy Q Px Py M2 都已知,我怎么求其他的数啊,就是对应的位移和等效结点力两者只知道一个,但是未知数与方程数相同,自己手算肯定能算出答案的,matlab里怎么算这种
你这方法不错,但是如果我已知的U不是连续出现的,就不能简单的分块了,要逐个挑出来,形成新的刚度矩阵,我不是要解决具体问题,至少这个方法,我要稍微有点普遍意义啊
matlab 矩阵求方程这是一个矩阵位移法的题,如果6x6的刚度矩阵K我都知道了,但是方程是KU=F,U=[Ux Uy Q Vx Vy P]',F=[Fx Fy M1 Px Py M2]',比如其中Ux Uy Q Px Py M2 都已知,我怎么求其他的数啊,就是对应的位移
你这个有点像有限元里面的东西.
下面的K是我乱写的.
你如果给定的不是前三个位移后三个应力的话,你自己得适当地调整,实际上就是调整下面的K1和K2,你一个命令地看,应当可以OK的.
K=[2,1,0,0,0,0;1,2,1,0,0,0;0,1,2,1,0,0;0,0,1,2,1,0;0,0,0,1,2,1;0,0,0,0,1,2];
u1=[1,2,3]';
F2=[4,5,6]';
K1=K(4:6,1:3);
K2=K(4:6,4:6);
u2=K2\(F2-K1*u1);
u=[u1;u2];
F=K*u;
%
%第一次修改
n=6;
K=[2,1,0,0,0,0;1,2,1,0,0,0;0,1,2,1,0,0;0,0,1,2,1,0;0,0,0,1,2,1;0,0,0,0,1,2];
n_node_ok=3; %三个节点位移己知
n_node=[1,2,4];%第1,2,4个节点的位移己知,这个的长度要与上面匹配
u1=[1,2,3]'; %三个节点(1,2,4节点)的位移值
F2=[4,5,6]'; %其它节点(2,5,6节点)的载荷值
n_f=1:n;
n_f(n_node)=[]; %给定载荷的节点号
K1=K(n_f,n_node);
K2=K(n_f,n_f);
u2=K2\(F2-K1*u1); %给定载荷节点的位移值
%后处理(所有节点的位移值的合成)
u(n_node)=u1;
u(n_f)=u2;
u=u';
F=K*u;
u,F