function y = QuadraticQuadAssemble(K,k,i,j,m,p,q,r,s,t)
%QuadraticQuadAssemble This function assembles the element
% stiffness matrix k of the quadratic
% quadrilateral element with nodes i, j,
% m, p, q, r, s, and t into the global
% stiffness matrix K.
% This function returns the global
% stiffness matrix K after the element
% stiffness matrix k is assembled.
K(2*i-1,2*i-1) = K(2*i-1,2*i-1) + k(1,1);
K(2*i-1,2*i) = K(2*i-1,2*i) + k(1,2);
K(2*i-1,2*j-1) = K(2*i-1,2*j-1) + k(1,3);
K(2*i-1,2*j) = K(2*i-1,2*j) + k(1,4);
K(2*i-1,2*m-1) = K(2*i-1,2*m-1) + k(1,5);
K(2*i-1,2*m) = K(2*i-1,2*m) + k(1,6);
K(2*i-1,2*p-1) = K(2*i-1,2*p-1) + k(1,7);
K(2*i-1,2*p) = K(2*i-1,2*p) + k(1,8);
K(2*i-1,2*q-1) = K(2*i-1,2*q-1) + k(1,9);
K(2*i-1,2*q) = K(2*i-1,2*q) + k(1,10);
K(2*i-1,2*r-1) = K(2*i-1,2*r-1) + k(1,11);
K(2*i-1,2*r) = K(2*i-1,2*r) + k(1,12);
K(2*i-1,2*s-1) = K(2*i-1,2*s-1) + k(1,13);
K(2*i-1,2*s) = K(2*i-1,2*s) + k(1,14);
K(2*i-1,2*t-1) = K(2*i-1,2*t-1) + k(1,15);
K(2*i-1,2*t) = K(2*i-1,2*t) + k(1,16);
K(2*i,2*i-1) = K(2*i,2*i-1) + k(2,1);
K(2*i,2*i) = K(2*i,2*i) + k(2,2);
K(2*i,2*j-1) = K(2*i,2*j-1) + k(2,3);
K(2*i,2*j) = K(2*i,2*j) + k(2,4);
K(2*i,2*m-1) = K(2*i,2*m-1) + k(2,5);
K(2*i,2*m) = K(2*i,2*m) + k(2,6);
K(2*i,2*p-1) = K(2*i,2*p-1) + k(2,7);
K(2*i,2*p) = K(2*i,2*p) + k(2,8);
K(2*i,2*q-1) = K(2*i,2*q-1) + k(2,9);
K(2*i,2*q) = K(2*i,2*q) + k(2,10);
K(2*i,2*r-1) = K(2*i,2*r-1) + k(2,11);
K(2*i,2*r) = K(2*i,2*r) + k(2,12);
K(2*i,2*s-1) = K(2*i,2*s-1) + k(2,13);
K(2*i,2*s) = K(2*i,2*s) + k(2,14);
K(2*i,2*t-1) = K(2*i,2*t-1) + k(2,15);
K(2*i,2*t) = K(2*i,2*t) + k(2,16);
K(2*j-1,2*i-1) = K(2*j-1,2*i-1) + k(3,1);
K(2*j-1,2*i) = K(2*j-1,2*i) + k(3,2);
K(2*j-1,2*j-1) = K(2*j-1,2*j-1) + k(3,3);
K(2*j-1,2*j) = K(2*j-1,2*j) + k(3,4);
K(2*j-1,2*m-1) = K(2*j-1,2*m-1) + k(3,5);
K(2*j-1,2*m) = K(2*j-1,2*m) + k(3,6);
K(2*j-1,2*p-1) = K(2*j-1,2*p-1) + k(3,7);
K(2*j-1,2*p) = K(2*j-1,2*p) + k(3,8);
K(2*j-1,2*q-1) = K(2*j-1,2*q-1) + k(3,9);
K(2*j-1,2*q) = K(2*j-1,2*q) + k(3,10);
K(2*j-1,2*r-1) = K(2*j-1,2*r-1) + k(3,11);
K(2*j-1,2*r) = K(2*j-1,2*r) + k(3,12);
K(2*j-1,2*s-1) = K(2*j-1,2*s-1) + k(3,13);
K(2*j-1,2*s) = K(2*j-1,2*s) + k(3,14);
K(2*j-1,2*t-1) = K(2*j-1,2*t-1) + k(3,15);
K(2*j-1,2*t) = K(2*j-1,2*t) + k(3,16);
K(2*j,2*i-1) = K(2*j,2*i-1) + k(4,1);
K(2*j,2*i) = K(2*j,2*i) + k(4,2);
K(2*j,2*j-1) = K(2*j,2*j-1) + k(4,3);
K(2*j,2*j) = K(2*j,2*j) + k(4,4);
K(2*j,2*m-1) = K(2*j,2*m-1) + k(4,5);
K(2*j,2*m) = K(2*j,2*m) + k(4,6);
K(2*j,2*p-1) = K(2*j,2*p-1) + k(4,7);
K(2*j,2*p) = K(2*j,2*p) + k(4,8);
K(2*j,2*q-1) = K(2*j,2*q-1) + k(4,9);
K(2*j,2*q) = K(2*j,2*q) + k(4,10);
K(2*j,2*r-1) = K(2*j,2*r-1) + k(4,11);
K(2*j,2*r) = K(2*j,2*r) + k(4,12);
K(2*j,2*s-1) = K(2*j,2*s-1) + k(4,13);
K(2*j,2*s) = K(2*j,2*s) + k(4,14);
K(2*j,2*t-1) = K(2*j,2*t-1) + k(4,15);
K(2*j,2*t) = K(2*j,2*t) + k(4,16);
K(2*m-1,2*i-1) = K(2*m-1,2*i-1) + k(5,1);
K(2*m-1,2*i) = K(2*m-1,2*i) + k(5,2);
K(2*m-1,2*j-1) = K(2*m-1,2*j-1) + k(5,3);
K(2*m-1,2*j) = K(2*m-1,2*j) + k(5,4);
K(2*m-1,2*m-1) = K(2*m-1,2*m-1) + k(5,5);
K(2*m-1,2*m) = K(2*m-1,2*m) + k(5,6);
K(2*m-1,2*p-1) = K(2*m-1,2*p-1) + k(5,7);
K(2*m-1,2*p) = K(2*m-1,2*p) + k(5,8);
K(2*m-1,2*q-1) = K(2*m-1,2*q-1) + k(5,9);
K(2*m-1,2*q) = K(2*m-1,2*q) + k(5,10);
K(2*m-1,2*r-1) = K(2*m-1,2*r-1) + k(5,11);
K(2*m-1,2*r) = K(2*m-1,2*r) + k(5,12);
K(2*m-1,2*s-1) = K(2*m-1,2*s-1) + k(5,13);
K(2*m-1,2*s) = K(2*m-1,2*s) + k(5,14);
K(2*m-1,2*t-1) = K(2*m-1,2*t-1) + k(5,15);
K(2*m-1,2*t) = K(2*m-1,2*t) + k(5,16);
K(2*m,2*i-1) = K(2*m,2*i-1) + k(6,1);
K(2*m,2*i) = K(2*m,2*i) + k(6,2);
K(2*m,2*j-1) = K(2*m,2*j-1) + k(6,3);
K(2*m,2*j) = K(2*m,2*j) + k(6,4);
K(2*m,2*m-1) = K(2*m,2*m-1) + k(6,5);
K(2*m,2*m) = K(2*m,2*m) + k(6,6);
K(2*m,2*p-1) = K(2*m,2*p-1) + k(6,7);
K(2*m,2*p) = K(2*m,2*p) + k(6,8);
K(2*m,2*q-1) = K(2*m,2*q-1) + k(6,9);
K(2*m,2*q) = K(2*m,2*q) + k(6,10);
K(2*m,2*r-1) = K(2*m,2*r-1) + k(6,11);
K(2*m,2*r) = K(2*m,2*r) + k(6,12);
K(2*m,2*s-1) = K(2*m,2*s-1) + k(6,13);
K(2*m,2*s) = K(2*m,2*s) + k(6,14);
K(2*m,2*t-1) = K(2*m,2*t-1) + k(6,15);
K(2*m,2*t) = K(2*m,2*t) + k(6,16);
K(2*p-1,2*i-1) = K(2*p-1,2*i-1) + k(7,1);
K(2*p-1,2*i) = K(2*p-1,2*i) + k(7,2);
K(2*p-1,2*j-1) = K(2*p-1,2*j-1) + k(7,3);
K(2*p-1,2*j) = K(2*p-1,2*j) + k(7,4);
K(2*p-1,2*m-1) = K(2*p-1,2*m-1) + k(7,5);
K(2*p-1,2*m) = K(2*p-1,2*m) + k(7,6);
K(2*p-1,2*p-1) = K(2*p-1,2*p-1) + k(7,7);
K(2*p-1,2*p) = K(2*p-1,2*p) + k(7,8);
K(2*p-1,2*q-1) = K(2*p-1,2*q-1) + k(7,9);
K(2*p-1,2*q) = K(2*p-1,2*q) + k(7,10);
K(2*p-1,2*r-1) = K(2*p-1,2*r-1) + k(7,11);
K(2*p-1,2*r) = K(2*p-1,2*r) + k(7,12);
K(2*p-1,2*s-1) = K(2*p-1,2*s-1) + k(7,13);
K(2*p-1,2*s) = K(2*p-1,2*s) + k(7,14);
K(2*p-1,2*t-1) = K(2*p-1,2*t-1) + k(7,15);
K(2*p-1,2*t) = K(2*p-1,2*t) + k(7,16);
K(2*p,2*i-1) = K(2*p,2*i-1) + k(8,1);
K(2*p,2*i) = K(2*p,2*i) + k(8,2);
K(2*p,2*j-1) = K(2*p,2*j-1) + k(8,3);
K(2*p,2*j) = K(2*p,2*j) + k(8,4);
K(2*p,2*m-1) = K(2*p,2*m-1) + k(8,5);
K(2*p,2*m) = K(2*p,2*m) + k(8,6);
K(2*p,2*p-1) = K(2*p,2*p-1) + k(8,7);
K(2*p,2*p) = K(2*p,2*p) + k(8,8);
K(2*p,2*q-1) = K(2*p,2*q-1) + k(8,9);
K(2*p,2*q) = K(2*p,2*q) + k(8,10);
K(2*p,2*r-1) = K(2*p,2*r-1) + k(8,11);
K(2*p,2*r) = K(2*p,2*r) + k(8,12);
K(2*p,2*s-1) = K(2*p,2*s-1) + k(8,13);
K(2*p,2*s) = K(2*p,2*s) + k(8,14);
K(2*p,2*t-1) = K(2*p,2*t-1) + k(8,15);
K(2*p,2*t) = K(2*p,2*t) + k(8,16);
K(2*q-1,2*i-1) = K(2*q-1,2*i-1) + k(9,1);
K(2*q-1,2*i) = K(2*q-1,2*i) + k(9,2);
K(2*q-1,2*j-1) = K(2*q-1,2*j-1) + k(9,3);
K(2*q-1,2*j) = K(2*q-1,2*j) + k(9,4);
K(2*q-1,2*m-1) = K(2*q-1,2*m-1) + k(9,5);
K(2*q-1,2*m) = K(2*q-1,2*m) + k(9,6);
K(2*q-1,2*p-1) = K(2*q-1,2*p-1) + k(9,7);
K(2*q-1,2*p) = K(2*q-1,2*p) + k(9,8);
K(2*q-1,2*q-1) = K(2*q-1,2*q-1) + k(9,9);
K(2*q-1,2*q) = K(2*q-1,2*q) + k(9,10);
K(2*q-1,2*r-1) = K(2*q-1,2*r-1) + k(9,11);
K(2*q-1,2*r) = K(2*q-1,2*r) + k(9,12);
K(2*q-1,2*s-1) = K(2*q-1,2*s-1) + k(9,13);
K(2*q-1,2*s) = K(2*q-1,2*s) + k(9,14);
K(2*q-1,2*t-1) = K(2*q-1,2*t-1) + k(9,15);
K(2*q-1,2*t) = K(2*q-1,2*t) + k(9,16);
K(2*q,2*i-1) = K(2*q,2*i-1) + k(10,1);
K(2*q,2*i) = K(2*q,2*i) + k(10,2);
K(2*q,2*j-1) = K(2*q,2*j-1) + k(10,3);
K(2*q,2*j) = K(2*q,2*j) + k(10,4);
K(2*q,2*m-1) = K(2*q,2*m-1) + k(10,5);
K(2*q,2*m) = K(2*q,2*m) + k(10,6);
K(2*q,2*p-1) = K(2*q,2*p-1) + k(10,7);
K(2*q,2*p) = K(2*q,2*p) + k(10,8);
K(2*q,2*q-1) = K(2*q,2*q-1) + k(10,9);
K(2*q,2*q) = K(2*q,2*q) + k(10,10);
K(2*q,2*r-1) = K(2*q,2*r-1) + k(10,11);
K(2*q,2*r) = K(2*q,2*r) + k(10,12);
K(2*q,2*s-1) = K(2*q,2*s-1) + k(10,13);
K(2*q,2*s) = K(2*q,2*s) + k(10,14);
K(2*q,2*t-1) = K(2*q,2*t-1) + k(10,15);
K(2*q,2*t) = K(2*q,2*t) + k(10,16);
K(2*r-1,2*i-1) = K(2*r-1,2*i-1) + k(11,1);
K(2*r-1,2*i) = K(2*r-1,2*i) + k(11,2);
K(2*r-1,2*j-1) = K(2*r-1,2*j-1) + k(11,3);
K(2*r-1,2*j) = K(2*r-1,2*j) + k(11,4);
K(2*r-1,2*m-1) = K(2*r-1,2*m-1) + k(11,5);
K(2*r-1,2*m) = K(2*r-1,2*m) + k(11,6);
K(2*r-1,2*p-1) = K(2*r-1,2*p-1) + k(11,7);
K(2*r-1,2*p) = K(2*r-1,2*p) + k(11,8);
K(2*r-1,2*q-1) = K(2*r-1,2*q-1) + k(11,9);
K(2*r-1,2*q) = K(2*r-1,2*q) + k(11,10);
K(2*r-1,2*r-1) = K(2*r-1,2*r-1) + k(11,11);
K(2*r-1,2*r) = K(2*r-1,2*r) + k(11,12);
K(2*r-1,2*s-1) = K(2*r-1,2*s-1) + k(11,13);
K(2*r-1,2*s) = K(2*r-1,2*s) + k(11,14);
K(2*r-1,2*t-1) = K(2*r-1,2*t-1) + k(11,15);
K(2*r-1,2*t) = K(2*r-1,2*t) + k(11,16);
K(2*r,2*i-1) = K(2*r,2*i-1) + k(12,1);
K(2*r,2*i) = K(2*r,2*i) + k(12,2);
K(2*r,2*j-1) = K(2*r,2*j-1) + k(12,3);
K(2*r,2*j) = K(2*r,2*j) + k(12,4);
K(2*r,2*m-1) = K(2*r,2*m-1) + k(12,5);
K(2*r,2*m) = K(2*r,2*m) + k(12,6);
K(2*r,2*p-1) = K(2*r,2*p-1) + k(12,7);
K(2*r,2*p) = K(2*r,2*p) + k(12,8);
K(2*r,2*q-1) = K(2*r,2*q-1) + k(12,9);
K(2*r,2*q) = K(2*r,2*q) + k(12,10);
K(2*r,2*r-1) = K(2*r,2*r-1) + k(12,11);
K(2*r,2*r) = K(2*r,2*r) + k(12,1
K.zip （75个子文件） 
K 
SpaceTrussElementLength.m 398B SpringElementForces.m 261B BeamAssemble.m 1013B QuadraticQuadElementArea.m 665B PlaneTrussElementLength.m 366B PlaneFrameAssemble.m 2KB BilinearQuadElementStresses.m 2KB SpaceFrameElementMomentZDiagram.m 464B PlaneFrameElementAxialDiagram.m 432B QuadTriangleElementArea.m 441B QuadraticBarAssemble.m 695B SpaceFrameElementShearZDiagram.m 453B LinearBarElementForces.m 303B LinearTriangleElementArea.m 450B SpaceFrameElementShearYDiagram.m 453B SpaceFrameElementAxialDiagram.m 432B SpringElementStiffness.m 305B GridElementStiffness.m 912B SpaceTrussElementForce.m 595B SpaceFrameAssemble.m 6KB SpaceFrameElementMomentYDiagram.m 464B PlaneTrussElementStiffness.m 650B GridAssemble.m 2KB SpaceFrameElementStiffness.m 2KB QuadraticBarElementStiffness.m 498B QuadraticQuadElementPStresses.m 510B QuadTriangleAssemble.m 6KB PlaneFrameElementMomentDiagram.m 443B QuadraticBarElementForces.m 319B BeamElementShearDiagram.m 367B PlaneFrameInclinedSupport.m 534B PlaneFrameElementStiffness.m 885B PlaneFrameElementForces.m 831B GridElementLength.m 329B BilinearQuadAssemble.m 3KB SpringAssemble.m 481B SpaceTrussAssemble.m 2KB TetrahedronElementStresses.m 2KB GridElementForces.m 823B BilinearQuadElementArea.m 655B QuadTriangleElementStiffness.m 3KB QuadraticQuadElementStresses.m 3KB SpaceFrameElementLength.m 398B LinearTriangleElementStresses.m 1KB SpaceTrussElementStress.m 549B BilinearQuadElementStiffness2.m 2KB PlaneTrussElementForce.m 460B SpaceFrameElementForces.m 2KB QuadraticBarElementStresses.m 401B PlaneFrameElementLength.m 366B BeamElementMomentDiagram.m 377B LinearTriangleElementPStresses.m 515B TetrahedronElementStiffness.m 2KB LinearBarElementStiffness.m 430B LinearBarAssemble.m 529B LinearBarElementStresses.m 351B SpaceFrameElementTorsionDiagram.m 437B LinearTriangleElementStiffness.m 1KB BeamElementForces.m 253B QuadraticQuadElementStiffness.m 3KB PlaneTrussAssemble.m 1KB QuadraticQuadAssemble.m 11KB BilinearQuadElementPStresses.m 505B PlaneTrussInclinedSupport.m 434B PlaneTrussElementStress.m 439B QuadTriangleElementStresses.m 3KB BilinearQuadElementStiffness.m 2KB TetrahedronElementPStresses.m 615B TetrahedronElementVolume.m 631B LinearTriangleAssemble.m 2KB TetrahedronAssemble.m 6KB SpaceTrussElementStiffness.m 736B QuadTriangleElementPStresses.m 505B BeamElementStiffness.m 499B PlaneFrameElementShearDiagram.m 433B
鹰忍
- 粉丝: 99
最新资源
- PLC舞台灯光设计方案.doc
- 学生信息管理系统-C语言课程方案设计书.doc
- 实验六教学板自检程序设计方案.doc
- 基于单片机大屏幕显示研究设计.doc
- web协同商务系统研究与原型开发.doc
- 钢结构CAD软件STS的功能及应用.docx
- 嵌入式单片机PPP协议的应用研究.doc
- 公路造价师考试辅导：流动资金扩大指标估算法试题.docx
- 用于预测性维护与健康管理的大型语言模型（故障诊断大模型；剩余使用寿命预测大模型）
- 2017年软件实施工程师笔试面试题及答案.docx
- 住宅小区海康网络监控系统方案.doc
- 结合电气工程及其自动化剖析机器人设计.docx
- 《信息系统分析与设计》第3章：通信与计算机网络.ppt
- Python编程作图物理仿真项目进阶设计.docx
- 基于区块链技术的电子轮机日志系统.docx
- 基于51单片机用LCD1602显示的DS18B20课程设计-键控上下限报警功能.doc
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
评论4