结构矩阵分析与程序设计钢架vb代码Word文档下载推荐.docx

结构矩阵分析与程序设计钢架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