结构矩阵分析与程序设计钢架vb代码Word文档下载推荐.docx
- 文档编号:15341819
- 上传时间:2022-10-29
- 格式:DOCX
- 页数:16
- 大小:17.45KB
结构矩阵分析与程序设计钢架vb代码Word文档下载推荐.docx
《结构矩阵分析与程序设计钢架vb代码Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《结构矩阵分析与程序设计钢架vb代码Word文档下载推荐.docx(16页珍藏版)》请在冰豆网上搜索。
\juzheng\钢架\tr3.2.11.txt"
ForInputAs#1
\juzheng\钢架\tw3.2.11.txt"
ForOutputAs#2
Callinput1
Callwstiff
Callload
Callbound
Callgauss
Callnqm
Close1
Close2
EndSub
===============================
SUB-1ReadAndPrintIntialData
==============================
Subinput1()
DimintiAsInteger,intjAsInteger,iAsInteger,jAsInteger,kAsInteger
Dimdx,dyAsDouble
Print#2,"
PlaneFramestructuralAnalysis"
*******************************"
inputdata"
====="
Print#2,
structuralcontroldata"
---------------------"
nn"
;
Spc(3);
"
ne"
nf"
nd"
ndf"
Spc
(2);
npj"
npe"
n"
Input#1,nn,ne,nf,nd,ndf,npj,npe
n=3*(nn-nf)
Print#2,nn;
ne;
nf;
nd;
ndf;
npj;
npe;
n
Nodalcoordinates"
Node"
x"
Spc(5);
y"
i=nn
Forinti=1Toi
Input#1,inti,x(inti),y(inti)
Print#2,inti;
x(inti);
y(inti)
Nextinti
Print#2,
Print#2,"
ElementInformation"
Ele.No."
Spc(4);
;
jl"
jr"
Spc(6);
ea"
ei"
al"
i=ne
Input#1,inti,jl(inti),jr(inti),ea(inti),ei(inti)
Forinti=1Toi
Ifjl(inti)>
=jr(inti)ThenStop
Nextinti
j=jl(inti)
k=jr(inti)
dx=x(k)-x(j)
dy=y(k)-y(j)
al(inti)=Sqr(dx*dx+dy*dy)
Print#2,Spc(3);
inti;
jl(inti);
jr(inti);
ea(inti);
ei(inti);
al(inti)
k=npj
Ifk<
>
0Then
NodalLoad"
i"
Spc(13);
mj"
xd"
yd"
md"
Forinti=1Tok
Input#1,inti,mj(inti),qj(inti,1),qj(inti,2),qj(inti,3)
Spc
(1),mj(inti);
Spc
(1);
qj(inti,1);
qj(inti,2);
qj(inti,3)
EndIf
i=npe
Ifi<
Elementloads"
mf"
ind"
aq"
bq"
q1"
q2"
Input#1,inti,mf(inti),ind(inti),aq(inti),bq(inti),q1(inti),q2(inti)
mf(inti);
ind(inti);
aq(inti);
bq(inti);
q1(inti);
q2(inti)
EndIf
j=ndf
Ifj<
Bonundaryconditions"
ibd"
bd"
Forinti=1Toj
Input#1,inti,ibd(inti),bd(inti)
ibd(inti);
bd(inti)
========================================================
sub-2AssemnbleStructuralStiffnessMatrix{R}
Subwstiff()
DimiAsInteger,jAsInteger,ieAsInteger,k1AsInteger,k2AsInteger
Fori=1Ton
Forj=1Ton
r(i,j)=0
Nextj
Nexti
ie=1
DoWhileie<
=ne
Callstiff(ie)
Calllocat(ie)
Fork1=1To6
i=ii(k1)
=nThen
Fork2=k1To6
j=ii(k2)
r(i,j)=r(i,j)+c(k1,k2)
Nextk2
Nextk1
ie=ie+1
Loop
Fori=2Ton
Forj=1To(i-1)
r(i,j)=r(j,i)
sub-3setupStiffnessMatrix[c]
Substiff(ie)
DimiAsInteger,jAsInteger
DimcxAsDouble,cyAsDouble,b1AsDouble,b2AsDouble,b3AsDouble,b4AsDouble
Dims1AsDouble,s2AsDouble,s3AsDouble,s4AsDouble,s5AsDouble,s6AsDouble
i=jl(ie)
j=jr(ie)
cx=(x(j)-x(i))/al(ie)
cy=(y(j)-y(i))/al(ie)
b1=ea(ie)/al(ie)
b2=12#*ei(ie)/al(ie)^3
b3=6#*ei(ie)/al(ie)^2
b4=2#*ei(ie)/al(ie)
s1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 结构 矩阵 分析 程序设计 钢架 vb 代码
![提示](https://static.bdocx.com/images/bang_tan.gif)