原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
桁架结构计算 第四章 P56
*******************************************************************************
function y=plane_truss_element_stiffness(E,A,L,theta) % 平面桁架单元刚度
x=theta*pi/180;
C=cos(x);
S=sin(x);
y=E*A/L*[
C*C C*S -C*C -C*S;
C*S S*S -C*S -S*S;
-C*C -C*S C*C C*S;
-C*S -S*S C*S S*S]; %平面桁架刚度矩阵
*******************************************************************************
function y=plane_truss_assemble(K,k,i,j) %平面桁架组装
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,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*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,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);
y=K;
*******************************************************************************
function y=plane_truss_element_force(E,A,L,theta,u) %力的表达式
x=theta*pi/180;
C=cos(x);
S=sin(x);
y=E*A/L*[-C -S C S]*u;
*******************************************************************************
function y=plane_truss_element_stress(E,L,theta,u) %应力表达式
x=theta*pi/180;
C=cos(x);
S=sin(x);
y=E/L*[-C -S C S]*u;
*****************************************************************************************************
*****************************************************************************************************
clear;
clc;
E= ;
A= ;
L1= ;%杆长度
L2= ;
L3= ;
L4= ;%杆长度
L5= ;
L6= ;
L7= ;
L8= ;
t1= ;%杆角度
t2= ;
t3= ;
t4= ;%杆角度
t5= ;
t6= ;
t7= ;
t8= ;
k1=plan
文档评论（0）