mdlcommand.bas

用visual basic编写的有限元程序！

Attribute VB_Name = "mdlCommand"
Public Function Do_APDL(ARG$()) As Long
    Rem APDL命令处理,By Sdh,2005-11-9
    Dim nNODE&, sLab$, uValue#
    Dim num&, Point#(2), Value1#
    Dim P1%, P2%, Pc%, R#
    Dim Lines%(10), LCount%
    Dim I&, J&, K&
    Select Case ARG(0)
        Case "N"
            num = Val(ARG(1)): Point(0) = Val(ARG(2)): Point(1) = Val(ARG(3))
            If (num > 0) Then
                Call Add_Node(Point, num)
            End If
        Case "E"
            I = Val(ARG(1)): J = Val(ARG(2)): K = Val(ARG(3))
            If (I > 0) And (J > 0) And (K > 0) Then
                Call Add_Elem(I, J, K)
            End If
        Case "D"
            nNODE = Val(ARG(1)): sLab = ARG(2): uValue = Val(ARG(3))
            Call Add_Dof(nNODE, sLab, uValue)
        Case "F"
            nNODE = Val(ARG(1)): sLab = ARG(2): uValue = Val(ARG(3))
            Call Add_Force(nNODE, sLab, uValue)
        Case "MPDATA"
            sLab = ARG(1): num = Val(ARG(2)): uValue = Val(ARG(4))
            Call Add_MpData(num, sLab, uValue)
        Case "TYPE":   DEF_ELEM.TYP = Val(ARG(1))
        Case "MAT":    DEF_ELEM.MAT = Val(ARG(1))
        Case "REAL":   DEF_ELEM.REL = Val(ARG(1))
        Case "ESYS":   DEF_ELEM.ESY = Val(ARG(1))
        Case "SECNUM": DEF_ELEM.SEC = Val(ARG(1))
        Case "R"
            num = Val(ARG(1)):  uValue = Val(ARG(2))
            Call Add_Real(num, uValue)
        Case "NMODIF"
            num = Val(ARG(1)):  Point(0) = dtr(Val(ARG(5)))
            Call Node_Theta(Point, num)
        Case Else
    End Select
End Function


Function COM_NMODIF(NODE&, X, Y, Z, THXY, THYZ, THZX)
    If Not IsEmpty(X) Then
        
    End If
End Function