整合整体刚度遇到函数编写的难题

secondye

小弟最近在编写动力有限元的程序，在整合整体刚度矩阵的时候遇到一个难题，由于是初学者，望各位指教：
有x,y分别为2维和6维列向量，

当x=[1 2]'，y=[1 2 3 4 5 6]';

当x=[2 3]'，y=[4 5 6 7 8 9 ]'，

当x=[3 4]'，y=[7 8 9 10 11 12 ]'，……如此类推，求出y的表达式。
谢谢各位了。

ChaChing

个人水平有限, 真看不出y的规律性!:@(

secondye

我的一个哥们帮我搞掂了
y(1:6)=[3*(x(1)-1)+1,3*(x(1)-1)+2,3*(x(1)-1)+3,3*(x(2)-1)+1,3*(x(2)-1)+2,3*(x(2)-1)+3]

ChaChing

第一时间, 我就感觉问题怪怪的, 但一直没能回覆讨论
问题是: 没事搅个x(2)干什麽? 有其他目的吗?

xx=[1,2,3]'; yy=[1,4,7]'; aa=[xx,ones(3,1)]; sol=aa\yy
上面求出sol=[3 -2], 表示yy=3*xx-2, 所以y=3*x(1)-2+(0:5)'即可

还有楼主给的式子或许可以改为
y=[3*(x(1)-1)+(1:3)';3*(x(2)-1)+(1:3)']

secondye

回复 ChaChing 的帖子

哦，谢谢版主的指教。
这个是在用来集成整体刚度矩阵时用的，每个节点有三个自由度。
当x=[1 2]'，y=[1 2 3 4 5 6]'; 表示节点1，2，对应整体坐标中的自由度是1，2，3，4，5，6
当x=[2 3]'，y=[4 5 6 7 8 9 ]'，同理。
   

ChaChing

回复 secondye 的帖子
喔, 终於了解了! 谢谢
那就不能仅使用x(1), 毕竟节点x(2)不一定是x(1)+1
还是得用y=[3*(x(1)-1)+(1:3)';3*(x(2)-1)+(1:3)']

secondye

回复 ChaChing 的帖子

嗯，刚刚把程序调试通了。是关于空间桁架单元的动力有限元程序，把单元类型改一改的话也可以用来梁单元，板单元等

   

zhouyang664

我有点儿不理解,不就是x=[i,i+1],y=3*i-1:3*i+4，不知道你是什么意思？

ChaChing

回复 zhouyang664 的帖子

节点x(2)不一定是x(1)+1 !?:@)

zhouyang664

回复 ChaChing 的帖子

应该是：x=[i,i+1],y=1+3*(i-1):1+3*(i-1)+5=3*i-2:3*i+3   

ChaChing

回复 zhouyang664 的帖子
可能个人未说清楚, LZ的问题系属於FEM
所以节点x(2)不一定是x(1)+1

zhouyang664

回复 ChaChing 的帖子


    专业问题，我理解不上去了！