1、juzheng钢架tr3.2.11.txt For Input As #1juzheng钢架tw3.2.11.txt For Output As #2Call input1Call wstiffCall loadCall boundCall gaussCall nqmClose 1Close 2End Sub=SUB-1 Read And Print Intial Data=Sub input1()Dim inti As Integer, intj As Integer, i As Integer, j As Integer, k As IntegerDim dx, dy As Double
2、Print #2, Plane Frame structural Analysis*input data= = = = = Print #2,structural control data-nn; Spc(3); nenfndndf Spc(2);npjnpen Input #1, nn, ne, nf, nd, ndf, npj, npe n = 3 * (nn - nf) Print #2, nn; ne; nf; nd; ndf; npj; npe; nNodal coordinatesNode x Spc(5);yi = nn For inti = 1 To i Input #1, i
3、nti, x(inti), y(inti) Print #2, inti; x(inti); y(inti) Next intiPrint #2,Print #2, Element InformationEle.No. Spc(4); ;jljr Spc(6);eaeiali = ne Input #1, inti, jl(inti), jr(inti), ea(inti), ei(inti)For inti = 1 To iIf jl(inti) = jr(inti) Then StopNext inti j = jl(inti) k = jr(inti) dx = x(k) - x(j)
4、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 If k 0 ThenNodal Loadi Spc(13);mjxdydmd For inti = 1 To k 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,
5、 2); qj(inti, 3) End If i = npe If i Element loadsmfindaqbqq1q2 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)End If j = ndf If j Bonundary conditionsibdbdFor inti = 1 To j Input #1, inti, ibd(inti), bd(inti) ib
6、d(inti); bd(inti)=sub-2 Assemnble Structural Stiffness MatrixRSub wstiff()Dim i As Integer, j As Integer, ie As Integer, k1 As Integer, k2 As IntegerFor i = 1 To n For j = 1 To n r(i, j) = 0 Next jNext iie = 1Do While ie = ne Call stiff(ie) Call locat(ie)For k1 = 1 To 6 i = ii(k1)= n Then For k2 = k
7、1 To 6 j = ii(k2) r(i, j) = r(i, j) + c(k1, k2) Next k2 Next k1 ie = ie + 1LoopFor i = 2 To n For j = 1 To (i - 1) r(i, j) = r(j, i) sub-3 set up Stiffness MatrixcSub stiff(ie)Dim i As Integer, j As IntegerDim cx As Double, cy As Double, b1 As Double, b2 As Double, b3 As Double, b4 As DoubleDim s1 As Double, s2 As Double, s3 As Double, s4 As Double, s5 As Double, s6 As Doublei = 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) 3b3 = 6# * ei(ie) / al(ie) 2b4 = 2# * ei(ie) / al(ie)s1
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1