Exemplo n.º 1
0
        static void Main(string[] args)
        {
            var id = new OpenSees.IDWrapper(10);

            id[1] = 2;
            var mat = new OpenSees.MatrixWrapper(2, 2);

            mat[1, 1] = 3;
            var theDomain = new OpenSees.Components.DomainWrapper();
            var node1     = new OpenSees.Components.NodeWrapper(1, 2, 0, 0);
            var node2     = new OpenSees.Components.NodeWrapper(2, 2, 144.0, 0);
            var node3     = new OpenSees.Components.NodeWrapper(3, 2, 168.0, 0);
            var node4     = new OpenSees.Components.NodeWrapper(4, 2, 72.0, 96.0);

            theDomain.AddNode(new OpenSees.Components.NodeWrapper[] { node1, node2, node3, node4 });

            var theMaterial = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(1, 3000, 0);
            var truss1      = new OpenSees.Elements.TrussWrapper(1, 2, 1, 4, theMaterial, 10.0, 0, 0, 0);
            var truss2      = new OpenSees.Elements.TrussWrapper(2, 2, 2, 4, theMaterial, 5.0, 0, 0, 0);
            var truss3      = new OpenSees.Elements.TrussWrapper(3, 2, 3, 4, theMaterial, 5.0, 0, 0, 0);

            theDomain.AddElement(new OpenSees.Elements.ElementWrapper[] { truss1, truss2, truss3 });

            var sp1 = new OpenSees.Components.Constraints.SP_ConstraintWrapper(1, 0, 0.0, true);
            var sp2 = new OpenSees.Components.Constraints.SP_ConstraintWrapper(1, 1, 0.0, true);
            var sp3 = new OpenSees.Components.Constraints.SP_ConstraintWrapper(2, 0, 0.0, true);
            var sp4 = new OpenSees.Components.Constraints.SP_ConstraintWrapper(2, 1, 0.0, true);
            var sp5 = new OpenSees.Components.Constraints.SP_ConstraintWrapper(3, 0, 0.0, true);
            var sp6 = new OpenSees.Components.Constraints.SP_ConstraintWrapper(3, 1, 0.0, true);

            theDomain.AddSP_Constraint(new OpenSees.Components.Constraints.SP_ConstraintWrapper[] { sp1, sp2, sp3, sp4, sp5, sp6 });

            var theSeries      = new OpenSees.Components.Timeseries.LinearSeriesWrapper();
            var theLoadPattern = new OpenSees.Components.LoadPatterns.LoadPatternWrapper(1);

            theLoadPattern.SetTimeSeries(theSeries);
            theDomain.AddLoadPattern(theLoadPattern);

            var theLoadValues = new OpenSees.VectorWrapper(2);

            theLoadValues[0] = 100;
            theLoadValues[1] = -50;

            var nodalLoad = new OpenSees.Components.Loads.NodalLoadWrapper(1, 4, theLoadValues, false);

            theDomain.AddNodalLoad(nodalLoad, 1);

            var theModel      = new OpenSees.AnalysisModelWrapper();
            var theSolnAlgo   = new OpenSees.Algorithms.LinearWrapper();
            var theIntegrator = new OpenSees.Integrators.Static.LoadControlWrapper(1.0, 1, 1.0, 1.0);
            var theHandler    = new OpenSees.Handlers.PlainHandlerWrapper();
            var theRCM        = new OpenSees.GraphNumberers.RCMWrapper(false);
            var theNumberer   = new OpenSees.Numberers.DOF_NumbererWrapper(theRCM);
            var theSolver     = new OpenSees.Systems.Linears.BandSPDLinLapackSolverWrapper();
            var theSOE        = new OpenSees.Systems.Linears.BandSPDLinSOEWrapper(theSolver);
            var theAnalysis   = new OpenSees.Analysis.StaticAnalysisWrapper(
                theDomain,
                theHandler,
                theNumberer,
                theModel,
                theSolnAlgo,
                theSOE,
                theIntegrator,
                null);

            int numSteps = 1;
            var ret      = theAnalysis.Analyze(numSteps);

            theDomain.CalculateNodalReactions(0);
            var disp = node4.GetCommitDisp();

            Console.WriteLine($"node disps : {{{string.Join(" ,", disp)}}}");
            Console.ReadKey();
        }
Exemplo n.º 2
0
        static void Main(string[] args)
        {
            var theDomain = new OpenSees.Components.DomainWrapper();

            #region paramaters

            var length = 1400.0;
            var width_Timoshenko_long      = 1100.0;
            var width_Timoshenko_trans     = 1400.0;
            var width_inside_long          = 385.0;
            var width_inside_trans         = 465.0;
            var width_Perimeter_long_right = 220.0;
            var width_Perimeter_long_left  = 110.0;
            var width_Perimeter_trans_up   = 250.0;
            var width_Perimeter_trans_down = 140.0;
            var tickness = 40.0;

            var Tension       = 666.0;
            var IP            = 8590.0;
            var OOP           = 7177.0;
            var about_Tension = 1000000000000000.0;
            var about_IP      = 1000000000000000.0;;
            var about_OOP     = 93000.0;

            var E_Timber_Long  = 13200.0;
            var E_Timber_Trans = 2200.0;
            var poisson_ratio  = 0.4;
            var G_Timber       = 820.0;

            var Area_Timoshenko_long        = width_Timoshenko_long * tickness;
            var Iz_Timoshenko_long          = (tickness * width_Timoshenko_long * width_Timoshenko_long * width_Timoshenko_long) / 12.0;
            var Iy_Timoshenko_long          = (width_Timoshenko_long * tickness * tickness * tickness) / (12.0);
            var Av_Timoshenko_long          = (5 * width_Timoshenko_long * tickness) / (6);
            var Torsional_J_Timoshenko_long = (width_Timoshenko_long * tickness * (width_Timoshenko_long * width_Timoshenko_long + tickness * tickness)) / 12.0;

            var Area_Timoshenko_trans        = width_Timoshenko_trans * tickness;
            var Iz_Timoshenko_trans          = (tickness * width_Timoshenko_trans * width_Timoshenko_trans * width_Timoshenko_trans) / (12.0);
            var Iy_Timoshenko_trans          = (width_Timoshenko_trans * tickness * tickness * tickness) / (12.0);
            var Av_Timoshenko_trans          = (5 * width_Timoshenko_trans * tickness) / (6);
            var Torsional_J_Timoshenko_trans = (width_Timoshenko_trans * tickness * (width_Timoshenko_trans * width_Timoshenko_trans + tickness * tickness)) / 12.0;

            var Area_Perimeter_long_right        = width_Perimeter_long_right * tickness;
            var Iz_Perimeter_long_right          = 1.0 * (tickness * width_Perimeter_long_right * width_Perimeter_long_right * width_Perimeter_long_right) / (12.0);
            var Iy_Perimeter_long_right          = 1.0 * (width_Perimeter_long_right * tickness * tickness * tickness) / (12.0);
            var Av_Perimeter_long_right          = (5 * width_Perimeter_long_right * tickness) / (6);
            var Torsional_J_Perimeter_long_right = (width_Perimeter_long_right * tickness * (width_Perimeter_long_right * width_Perimeter_long_right + tickness * tickness)) / 12.0;

            var Area_Perimeter_long_left        = width_Perimeter_long_left * tickness;
            var Iz_Perimeter_long_left          = 1.0 * (tickness * width_Perimeter_long_left * width_Perimeter_long_left * width_Perimeter_long_left) / (12.0);
            var Iy_Perimeter_long_left          = 1.0 * (width_Perimeter_long_left * tickness * tickness * tickness) / (12.0);
            var Av_Perimeter_long_left          = (5 * width_Perimeter_long_left * tickness) / (6);
            var Torsional_J_Perimeter_long_left = (width_Perimeter_long_left * tickness * (width_Perimeter_long_left * width_Perimeter_long_left + tickness * tickness)) / 12.0;

            var Area_Perimeter_trans_up        = width_Perimeter_trans_up * tickness;
            var Iz_Perimeter_trans_up          = 1.0 * (tickness * width_Perimeter_trans_up * width_Perimeter_trans_up * width_Perimeter_trans_up) / (12.0);
            var Iy_Perimeter_trans_up          = 1.0 * (width_Perimeter_trans_up * tickness * tickness * tickness) / (12.0);
            var Av_Perimeter_trans_up          = (5 * width_Perimeter_trans_up * tickness) / (6);
            var Torsional_J_Perimeter_trans_up = (width_Perimeter_trans_up * tickness * (width_Perimeter_trans_up * width_Perimeter_trans_up + tickness * tickness)) / 12.0;

            var Area_Perimeter_trans_down        = width_Perimeter_trans_down * tickness;
            var Iz_Perimeter_trans_down          = 1.0 * (tickness * width_Perimeter_trans_down * width_Perimeter_trans_down * width_Perimeter_trans_down) / (12.0);
            var Iy_Perimeter_trans_down          = 1.0 * (width_Perimeter_trans_down * tickness * tickness * tickness) / (12.0);
            var Av_Perimeter_trans_down          = (5 * width_Perimeter_trans_down * tickness) / (6);
            var Torsional_J_Perimeter_trans_down = (width_Perimeter_trans_down * tickness * (width_Perimeter_trans_down * width_Perimeter_trans_down + tickness * tickness)) / 12.0;

            var Area_inside_long        = width_inside_long * tickness;
            var Iz_inside_long          = (tickness * width_inside_long * width_inside_long * width_inside_long) / (12.0);
            var Iy_inside_long          = (width_inside_long * tickness * tickness * tickness) / (12.0);
            var Av_inside_long          = (5 * width_inside_long * tickness) / (6);
            var Torsional_J_inside_long = (width_inside_long * tickness * (width_inside_long * width_inside_long + tickness * tickness)) / 12.0;

            var Area_inside_trans        = width_inside_trans * tickness;
            var Iz_inside_trans          = (tickness * width_inside_trans * width_inside_trans * width_inside_trans) / (12.0);
            var Iy_inside_trans          = (width_inside_trans * tickness * tickness * tickness) / (12.0);
            var Av_inside_trans          = (5 * width_inside_trans * tickness) / (6);
            var Torsional_J_inside_trans = (width_inside_trans * tickness * (width_inside_trans * width_inside_trans + tickness * tickness)) / 12.0;


            var rigidmattag      = 1;
            var semi_rigidmattag = 101;

            var Freemattag = 2;

            var timbermattag_long  = 30;
            var timbermattag_trans = 31;

            var Axial_perimeter_long_right = 40;
            var Axial_perimeter_long_left  = 41;
            var Axial_perimeter_trans_up   = 42;
            var Axial_perimeter_trans_down = 43;
            var Axial_inside_long          = 44;
            var Axial_inside_trans         = 45;

            var Shear_1_perimeter_long_right = 50;
            var Shear_1_perimeter_long_left  = 51;
            var Shear_1_perimeter_trans_up   = 52;
            var Shear_1_perimeter_trans_down = 53;
            var Shear_1_inside_long          = 54;
            var Shear_1_inside_trans         = 55;

            var Shear_2_perimeter_long  = 60;
            var Shear_2_perimeter_trans = 61;
            var Shear_2_inside_long     = 62;
            var Shear_2_inside_trans    = 63;

            var Trans_Shear_Timoshenko_Eq = 70;
            var Long_Shear_Timoshenko_Eq  = 71;

            var Stiffness_Tension       = 10;
            var Stiffness_IP            = 11;
            var Stiffness_OOP           = 12;
            var Stiffness_about_Tension = 13;
            var Stiffness_about_IP      = 14;
            var Stiffness_about_OOP     = 15;
            #endregion

            #region materials

            var rigidmattagMat      = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(rigidmattag, 1.0e+15, 0);                                                                              //uniaxialMaterial Elastic $$rigidmattag[expr 1.0e+15]
            var semi_rigidmattagMat = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(semi_rigidmattag, 1.0e+4, 0);                                                                          //uniaxialMaterial Elastic $semi_rigidmattag[expr 1.0e+4]
            var FreemattagMat       = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Freemattag, 1.0e-1, 0);                                                                                // uniaxialMaterial Elastic $Freemattag[expr 1.0e-1]

            var timbermattag_longMat          = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(timbermattag_long, E_Timber_Long, 0);                                                        //  uniaxialMaterial Elastic $timbermattag_long          $E_Timber_Long
            var timbermattag_transMat         = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(timbermattag_trans, E_Timber_Trans, 0);                                                      //  uniaxialMaterial    Elastic $timbermattag_trans         $E_Timber_Trans
            var Axial_perimeter_long_rightMat = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Axial_perimeter_long_right, 100.0 * E_Timber_Long * Area_Perimeter_long_right / length, 0);  // uniaxialMaterial    Elastic $Axial_perimeter_long_right[expr 100.0 *$E_Timber_Long *$Area_Perimeter_long_right /$length];         # for LT2
            var Axial_perimeter_long_leftMat  = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Axial_perimeter_long_left, 100.0 * E_Timber_Long * Area_Perimeter_long_left / length, 0);    // uniaxialMaterial Elastic $Axial_perimeter_long_left[expr 100.0 *$E_Timber_Long *$Area_Perimeter_long_left /$length];      # for LT2
            var Axial_perimeter_trans_upMat   = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Axial_perimeter_trans_up, 100.0 * E_Timber_Trans * Area_Perimeter_trans_up / length, 0);     // uniaxialMaterial Elastic $Axial_perimeter_trans_up[expr 100.0 *$E_Timber_Trans *$Area_Perimeter_trans_up /$length];      # for LT2
            var Axial_perimeter_trans_downMat = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Axial_perimeter_trans_down, 100.0 * E_Timber_Trans * Area_Perimeter_trans_down / length, 0); //uniaxialMaterial Elastic $Axial_perimeter_trans_down[expr 100.0 *$E_Timber_Trans *$Area_Perimeter_trans_down /$length];        # for LT2
            var Axial_inside_longMat          = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Axial_inside_long, 1.0 * E_Timber_Long * Area_inside_long / length, 0);                      // uniaxialMaterial Elastic $Axial_inside_long[expr 1.*$E_Timber_Long *$Area_inside_long /$length];             # for LT2
            var Axial_inside_transMat         = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Axial_inside_trans, 1.0 * E_Timber_Trans * Area_inside_trans / length, 0);                   //  uniaxialMaterial Elastic $Axial_inside_trans[expr 1.*$E_Timber_Trans *$Area_inside_trans /$length];           # for LT2


            var Shear_1_perimeter_long_rightMat = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_1_perimeter_long_right, 100.0 * (5 * width_Perimeter_long_right * tickness) / (6), 0); // uniaxialMaterial Elastic $Shear_1_perimeter_long_right[expr 100.0 * (5 *$width_Perimeter_long_right *$tickness)/ (6)];
            var Shear_1_perimeter_long_leftMat  = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_1_perimeter_long_left, 100.0 * (5 * width_Perimeter_long_left * tickness) / (6), 0);   // uniaxialMaterial Elastic $Shear_1_perimeter_long_left[expr 100.0 * (5 *$width_Perimeter_long_left *$tickness)/ (6)];
            var Shear_1_perimeter_trans_upMat   = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_1_perimeter_trans_up, 100.0 * (5 * width_Perimeter_trans_up * tickness) / (6), 0);     // uniaxialMaterial Elastic $Shear_1_perimeter_trans_up[expr 100.0 * (5 *$width_Perimeter_trans_up *$tickness)/ (6)];

            var Shear_1_perimeter_trans_downMat = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_1_perimeter_trans_down, 100.0 * (5 * width_Perimeter_trans_down * tickness) / (6), 0); // uniaxialMaterial Elastic $Shear_1_perimeter_trans_down[expr 100.0 * (5 *$width_Perimeter_trans_down *$tickness)/ (6)];
            var Shear_1_inside_longMat          = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_1_inside_long, 1.0 * (5 * width_inside_long * tickness) / (6), 0);                     // uniaxialMaterial Elastic $Shear_1_inside_long[expr 1.0 * (5 *$width_inside_long *$tickness)/ (6)];
            var Shear_1_inside_transMat         = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_1_inside_trans, 1.0 * (5 * width_inside_trans * tickness) / (6), 0);                   // uniaxialMaterial Elastic $Shear_1_inside_trans[expr 1.0 * (5 *$width_inside_trans *$tickness)/ (6)];

            var Shear_2_perimeter_longMat  = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_2_perimeter_long, 1.0e+15, 0);                                                              //uniaxialMaterial Elastic $Shear_2_perimeter_long     1.0e+15;
            var Shear_2_perimeter_transMat = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_2_perimeter_trans, 1.0e+15, 0);                                                             // uniaxialMaterial Elastic $Shear_2_perimeter_trans    1.0e+15;
            var Shear_2_inside_longMat     = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_2_inside_long, 1.0e+15, 0);                                                                 // uniaxialMaterial Elastic $Shear_2_inside_long        1.0e+15;
            var Shear_2_inside_transMat    = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Shear_2_inside_trans, 1.0e+15, 0);                                                                // uniaxialMaterial Elastic $Shear_2_inside_trans       1.0e+15;


            var Trans_Shear_Timoshenko_EqMat = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Trans_Shear_Timoshenko_Eq, 34787.878, 0); // uniaxialMaterial Elastic $Trans_Shear_Timoshenko_Eq[expr 34787.878]
            var Long_Shear_Timoshenko_EqMat  = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Long_Shear_Timoshenko_Eq, 16081.971, 0);  // uniaxialMaterial Elastic $Long_Shear_Timoshenko_Eq[expr 16081.971]


            var Stiffness_TensionMat       = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Stiffness_Tension, Tension, 0);             //  uniaxialMaterial Elastic $Stiffness_Tension[expr   $Tension]
            var Stiffness_IPMat            = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Stiffness_IP, IP, 0);                       // uniaxialMaterial Elastic $Stiffness_IP[expr   $IP]
            var Stiffness_OOPMat           = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Stiffness_OOP, OOP, 0);                     // uniaxialMaterial Elastic $Stiffness_OOP[expr   $OOP]
            var Stiffness_about_TensionMat = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Stiffness_about_Tension, about_Tension, 0); // uniaxialMaterial Elastic $Stiffness_about_Tension[expr   $about_Tension]
            var Stiffness_about_IPMat      = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Stiffness_about_IP, about_IP, 0);           // uniaxialMaterial Elastic $Stiffness_about_IP[expr   $about_IP]
            var Stiffness_about_OOPMat     = new OpenSees.Materials.Uniaxials.ElasticMaterialWrapper(Stiffness_about_OOP, about_OOP, 0);         //uniaxialMaterial Elastic $Stiffness_about_OOP[expr   $about_OOP]

            #endregion


            #region nodes
            var node1    = new OpenSees.Components.NodeWrapper(1, 6, 1100, 0, 700);
            var node2    = new OpenSees.Components.NodeWrapper(2, 6, 1100, 0, 466.7);
            var node3    = new OpenSees.Components.NodeWrapper(3, 6, 1100, 0, 765);
            var node4    = new OpenSees.Components.NodeWrapper(4, 6, 0, 0, 700);
            var node5    = new OpenSees.Components.NodeWrapper(5, 6, 0, 0, 933.4);
            var node6    = new OpenSees.Components.NodeWrapper(6, 6, 0, 0, 655);
            var node7    = new OpenSees.Components.NodeWrapper(7, 6, 1100, 0, 1155);
            var node8    = new OpenSees.Components.NodeWrapper(8, 6, 1100, 0, 1400);
            var node9    = new OpenSees.Components.NodeWrapper(9, 6, 1100, 0, 933.4);
            var node10   = new OpenSees.Components.NodeWrapper(10, 6, 1100, 0, 355);
            var node11   = new OpenSees.Components.NodeWrapper(11, 6, 1100, 0, 0);
            var node12   = new OpenSees.Components.NodeWrapper(12, 6, 0, 0, 1045);
            var node13   = new OpenSees.Components.NodeWrapper(13, 6, 0, 0, 1400);
            var node14   = new OpenSees.Components.NodeWrapper(14, 6, 0, 0, 466.7);
            var node15   = new OpenSees.Components.NodeWrapper(15, 6, 0, 0, 245);
            var node16   = new OpenSees.Components.NodeWrapper(16, 6, 0, 0, 0);
            var node17   = new OpenSees.Components.NodeWrapper(17, 6, 1140, 0, 355);
            var node18   = new OpenSees.Components.NodeWrapper(18, 6, 1140, 0, 765);
            var node19   = new OpenSees.Components.NodeWrapper(19, 6, 1140, 0, 1155);
            var node20   = new OpenSees.Components.NodeWrapper(20, 6, 255.207, 0, 1440);
            var node21   = new OpenSees.Components.NodeWrapper(21, 6, 255.207, 0, 1400);
            var node22   = new OpenSees.Components.NodeWrapper(22, 6, 672.394, 0, 1440);
            var node23   = new OpenSees.Components.NodeWrapper(23, 6, 672.394, 0, 1400);
            var node24   = new OpenSees.Components.NodeWrapper(24, 6, -40, 0, 1045);
            var node25   = new OpenSees.Components.NodeWrapper(25, 6, -40, 0, 655);
            var node26   = new OpenSees.Components.NodeWrapper(26, 6, -40, 0, 245);
            var node27   = new OpenSees.Components.NodeWrapper(27, 6, 426.25, 0, -40);
            var node28   = new OpenSees.Components.NodeWrapper(28, 6, 426.25, 0, 0);
            var node29   = new OpenSees.Components.NodeWrapper(29, 6, 838.75, 0, -40);
            var node30   = new OpenSees.Components.NodeWrapper(30, 6, 838.75, 0, 0);
            var node31   = new OpenSees.Components.NodeWrapper(31, 6, 733.34, 0, 1400);
            var node32   = new OpenSees.Components.NodeWrapper(32, 6, 366.67, 0, 1400);
            var node33   = new OpenSees.Components.NodeWrapper(33, 6, 550, 0, 1400);
            var node34   = new OpenSees.Components.NodeWrapper(34, 6, 733.34, 0, 0);
            var node35   = new OpenSees.Components.NodeWrapper(35, 6, 550, 0, 0);
            var node36   = new OpenSees.Components.NodeWrapper(36, 6, 366.67, 0, 0);
            var node1601 = new OpenSees.Components.NodeWrapper(1601, 6, 0, 0, 0);
            var node1602 = new OpenSees.Components.NodeWrapper(1602, 6, 0, 0, 0);
            var node1301 = new OpenSees.Components.NodeWrapper(1301, 6, 0, 0, 1400);
            var node1302 = new OpenSees.Components.NodeWrapper(1302, 6, 0, 0, 1400);
            var node1101 = new OpenSees.Components.NodeWrapper(1101, 6, 1100, 0, 0);
            var node1102 = new OpenSees.Components.NodeWrapper(1102, 6, 1100, 0, 0);
            var node802  = new OpenSees.Components.NodeWrapper(802, 6, 1100, 0, 1400);
            var node801  = new OpenSees.Components.NodeWrapper(801, 6, 1100, 0, 1400);
            var node3201 = new OpenSees.Components.NodeWrapper(3201, 6, 366.67, 0, 1400);
            var node3101 = new OpenSees.Components.NodeWrapper(3101, 6, 733.34, 0, 1400);
            var node3601 = new OpenSees.Components.NodeWrapper(3601, 6, 366.67, 0, 0);
            var node3401 = new OpenSees.Components.NodeWrapper(3401, 6, 733.34, 0, 0);
            var node1402 = new OpenSees.Components.NodeWrapper(1402, 6, 0, 0, 466.7);
            var node502  = new OpenSees.Components.NodeWrapper(502, 6, 0, 0, 933.4);
            var node902  = new OpenSees.Components.NodeWrapper(902, 6, 1100, 0, 933.4);
            var node202  = new OpenSees.Components.NodeWrapper(202, 6, 1100, 0, 466.7);

            var nodes = new OpenSees.Components.NodeWrapper[] { node1, node2, node3, node4, node5, node6, node7, node8, node9, node10,
                                                                node11, node12, node13, node14, node15, node16, node17, node18, node19, node20,
                                                                node21, node22, node23, node24, node25, node26, node27, node28, node29, node30,
                                                                node31, node32, node33, node34, node35, node36,
                                                                node1601, node1602, node1301, node1302, node1101, node1102, node802, node801, node3201, node3101,
                                                                node3601, node3401, node1402, node502, node902, node202 };
            theDomain.AddNode(nodes);

            foreach (var tag in new[] { 27, 29 /*, 24,25,26*/ })
            {
                for (var dof = 0; dof < 6; dof++)
                {
                    theDomain.AddSP_Constraint(new OpenSees.Components.Constraints.SP_ConstraintWrapper(tag, dof, 0.0, true));
                }
            }


            #endregion

            #region elements
            var xzplaneVec = new OpenSees.VectorWrapper(new double[] { 0, 1, 0 });
            var lineartag_TimoshenkoCrdT = new OpenSees.Elements.CrdTransfs.LinearCrdTransf3dWrapper(xzplaneVec);
            var lineartag_VerticalCrdT   = new OpenSees.Elements.CrdTransfs.LinearCrdTransf3dWrapper(xzplaneVec);
            var lineartag_HorizontalCrdT = new OpenSees.Elements.CrdTransfs.LinearCrdTransf3dWrapper(xzplaneVec);

            {
                var dirIds1        = new OpenSees.IDWrapper(new int[] { 1 }); // 2-- for wrapper
                var ele_link_tnl1  = new OpenSees.Elements.TwoNodeLinkWrapper(1, 3, 35, 33, dirIds1, new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Long_Shear_Timoshenko_EqMat }, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl58 = new OpenSees.Elements.TwoNodeLinkWrapper(58, 3, 1, 4, dirIds1, new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Trans_Shear_Timoshenko_EqMat }, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl1);
                theDomain.AddElement(ele_link_tnl58);
            }


            //puts "The Timoshenko beams are defined"
            //# Defining Link Elements Between Zero-Length Nodes
            var dirIds2 = new OpenSees.IDWrapper(new int[] { 0, 1, 2, 3, 4, 5 });


            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_perimeter_long_leftMat, Shear_1_perimeter_long_leftMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl2       = new OpenSees.Elements.TwoNodeLinkWrapper(2, 3, 16, 1601, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl2);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_perimeter_trans_downMat, Shear_1_perimeter_trans_downMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl3       = new OpenSees.Elements.TwoNodeLinkWrapper(3, 3, 1602, 16, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(new double[] { -1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl3);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_perimeter_long_leftMat, Shear_1_perimeter_long_leftMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl4       = new OpenSees.Elements.TwoNodeLinkWrapper(4, 3, 1301, 13, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl4);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_perimeter_trans_upMat, Shear_1_perimeter_trans_upMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl5       = new OpenSees.Elements.TwoNodeLinkWrapper(5, 3, 1302, 13, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(new double[] { -1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl5);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_perimeter_long_rightMat, Shear_1_perimeter_long_rightMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl6       = new OpenSees.Elements.TwoNodeLinkWrapper(6, 3, 11, 1101, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl6);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_perimeter_trans_downMat, Shear_1_perimeter_trans_downMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl7       = new OpenSees.Elements.TwoNodeLinkWrapper(7, 3, 11, 1102, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(new double[] { -1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl7);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_perimeter_long_rightMat, Shear_1_perimeter_long_rightMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl8       = new OpenSees.Elements.TwoNodeLinkWrapper(8, 3, 801, 8, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl8);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_perimeter_trans_upMat, Shear_1_perimeter_trans_upMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl9       = new OpenSees.Elements.TwoNodeLinkWrapper(9, 3, 8, 802, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(new double[] { -1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl9);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_inside_transMat, Shear_1_inside_transMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl10      = new OpenSees.Elements.TwoNodeLinkWrapper(10, 3, 2, 202, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(new double[] { -1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl11      = new OpenSees.Elements.TwoNodeLinkWrapper(11, 3, 9, 902, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(new double[] { -1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl12      = new OpenSees.Elements.TwoNodeLinkWrapper(12, 3, 502, 5, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(new double[] { -1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl13      = new OpenSees.Elements.TwoNodeLinkWrapper(13, 3, 1402, 14, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(new double[] { -1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl10);
                theDomain.AddElement(ele_link_tnl11);
                theDomain.AddElement(ele_link_tnl12);
                theDomain.AddElement(ele_link_tnl13);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Axial_inside_longMat, Shear_1_inside_longMat, rigidmattagMat, FreemattagMat, FreemattagMat, FreemattagMat };
                var ele_link_tnl14      = new OpenSees.Elements.TwoNodeLinkWrapper(14, 3, 36, 3601, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl15      = new OpenSees.Elements.TwoNodeLinkWrapper(15, 3, 34, 3401, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl16      = new OpenSees.Elements.TwoNodeLinkWrapper(16, 3, 3201, 32, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl17      = new OpenSees.Elements.TwoNodeLinkWrapper(17, 3, 3101, 31, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl14);
                theDomain.AddElement(ele_link_tnl15);
                theDomain.AddElement(ele_link_tnl16);
                theDomain.AddElement(ele_link_tnl17);
            }

            //# Defining Vertical Perimeter Elements LEFT
            {
                var ele_timbeam18 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(18, 1601, 15, E_Timber_Long, G_Timber, Area_Perimeter_long_left, Torsional_J_Perimeter_long_left, Iy_Perimeter_long_left, Iz_Perimeter_long_left, Av_Perimeter_long_left, Av_Perimeter_long_left, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam19 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(19, 15, 14, E_Timber_Long, G_Timber, Area_Perimeter_long_left, Torsional_J_Perimeter_long_left, Iy_Perimeter_long_left, Iz_Perimeter_long_left, Av_Perimeter_long_left, Av_Perimeter_long_left, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam20 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(20, 14, 6, E_Timber_Long, G_Timber, Area_Perimeter_long_left, Torsional_J_Perimeter_long_left, Iy_Perimeter_long_left, Iz_Perimeter_long_left, Av_Perimeter_long_left, Av_Perimeter_long_left, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam21 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(21, 6, 4, E_Timber_Long, G_Timber, Area_Perimeter_long_left, Torsional_J_Perimeter_long_left, Iy_Perimeter_long_left, Iz_Perimeter_long_left, Av_Perimeter_long_left, Av_Perimeter_long_left, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam22 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(22, 4, 5, E_Timber_Long, G_Timber, Area_Perimeter_long_left, Torsional_J_Perimeter_long_left, Iy_Perimeter_long_left, Iz_Perimeter_long_left, Av_Perimeter_long_left, Av_Perimeter_long_left, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam23 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(23, 5, 12, E_Timber_Long, G_Timber, Area_Perimeter_long_left, Torsional_J_Perimeter_long_left, Iy_Perimeter_long_left, Iz_Perimeter_long_left, Av_Perimeter_long_left, Av_Perimeter_long_left, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam24 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(24, 12, 1301, E_Timber_Long, G_Timber, Area_Perimeter_long_left, Torsional_J_Perimeter_long_left, Iy_Perimeter_long_left, Iz_Perimeter_long_left, Av_Perimeter_long_left, Av_Perimeter_long_left, lineartag_VerticalCrdT, 0, 0);
                theDomain.AddElement(ele_timbeam18);
                theDomain.AddElement(ele_timbeam19);
                theDomain.AddElement(ele_timbeam20);
                theDomain.AddElement(ele_timbeam21);
                theDomain.AddElement(ele_timbeam22);
                theDomain.AddElement(ele_timbeam23);
                theDomain.AddElement(ele_timbeam24);
            }

            //# Defining Vertical Perimeter Elements RIGHT
            {
                var ele_timbeam25 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(25, 1101, 10, E_Timber_Long, G_Timber, Area_Perimeter_long_right, Torsional_J_Perimeter_long_right, Iy_Perimeter_long_right, Iz_Perimeter_long_right, Av_Perimeter_long_right, Av_Perimeter_long_right, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam26 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(26, 10, 2, E_Timber_Long, G_Timber, Area_Perimeter_long_right, Torsional_J_Perimeter_long_right, Iy_Perimeter_long_right, Iz_Perimeter_long_right, Av_Perimeter_long_right, Av_Perimeter_long_right, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam27 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(27, 2, 1, E_Timber_Long, G_Timber, Area_Perimeter_long_right, Torsional_J_Perimeter_long_right, Iy_Perimeter_long_right, Iz_Perimeter_long_right, Av_Perimeter_long_right, Av_Perimeter_long_right, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam28 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(28, 1, 3, E_Timber_Long, G_Timber, Area_Perimeter_long_right, Torsional_J_Perimeter_long_right, Iy_Perimeter_long_right, Iz_Perimeter_long_right, Av_Perimeter_long_right, Av_Perimeter_long_right, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam29 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(29, 3, 9, E_Timber_Long, G_Timber, Area_Perimeter_long_right, Torsional_J_Perimeter_long_right, Iy_Perimeter_long_right, Iz_Perimeter_long_right, Av_Perimeter_long_right, Av_Perimeter_long_right, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam30 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(30, 9, 7, E_Timber_Long, G_Timber, Area_Perimeter_long_right, Torsional_J_Perimeter_long_right, Iy_Perimeter_long_right, Iz_Perimeter_long_right, Av_Perimeter_long_right, Av_Perimeter_long_right, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam31 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(31, 7, 801, E_Timber_Long, G_Timber, Area_Perimeter_long_right, Torsional_J_Perimeter_long_right, Iy_Perimeter_long_right, Iz_Perimeter_long_right, Av_Perimeter_long_right, Av_Perimeter_long_right, lineartag_VerticalCrdT, 0, 0);
                theDomain.AddElement(ele_timbeam25);
                theDomain.AddElement(ele_timbeam26);
                theDomain.AddElement(ele_timbeam27);
                theDomain.AddElement(ele_timbeam28);
                theDomain.AddElement(ele_timbeam29);
                theDomain.AddElement(ele_timbeam30);
                theDomain.AddElement(ele_timbeam31);
            }

            //# Defining Vertical Perimeter Elements DOWN
            {
                var ele_timbeam32 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(32, 1102, 30, E_Timber_Trans, G_Timber, Area_Perimeter_trans_down, Torsional_J_Perimeter_trans_down, Iy_Perimeter_trans_down, Iz_Perimeter_trans_down, Av_Perimeter_trans_down, Av_Perimeter_trans_down, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam33 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(33, 30, 34, E_Timber_Trans, G_Timber, Area_Perimeter_trans_down, Torsional_J_Perimeter_trans_down, Iy_Perimeter_trans_down, Iz_Perimeter_trans_down, Av_Perimeter_trans_down, Av_Perimeter_trans_down, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam34 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(34, 34, 35, E_Timber_Trans, G_Timber, Area_Perimeter_trans_down, Torsional_J_Perimeter_trans_down, Iy_Perimeter_trans_down, Iz_Perimeter_trans_down, Av_Perimeter_trans_down, Av_Perimeter_trans_down, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam35 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(35, 35, 28, E_Timber_Trans, G_Timber, Area_Perimeter_trans_down, Torsional_J_Perimeter_trans_down, Iy_Perimeter_trans_down, Iz_Perimeter_trans_down, Av_Perimeter_trans_down, Av_Perimeter_trans_down, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam36 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(36, 28, 36, E_Timber_Trans, G_Timber, Area_Perimeter_trans_down, Torsional_J_Perimeter_trans_down, Iy_Perimeter_trans_down, Iz_Perimeter_trans_down, Av_Perimeter_trans_down, Av_Perimeter_trans_down, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam37 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(37, 36, 1602, E_Timber_Trans, G_Timber, Area_Perimeter_trans_down, Torsional_J_Perimeter_trans_down, Iy_Perimeter_trans_down, Iz_Perimeter_trans_down, Av_Perimeter_trans_down, Av_Perimeter_trans_down, lineartag_HorizontalCrdT, 0, 0);

                theDomain.AddElement(ele_timbeam32);
                theDomain.AddElement(ele_timbeam33);
                theDomain.AddElement(ele_timbeam34);
                theDomain.AddElement(ele_timbeam35);
                theDomain.AddElement(ele_timbeam36);
                theDomain.AddElement(ele_timbeam37);
            }

            //# Defining Vertical Perimeter Elements UP
            {
                var ele_timbeam38 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(38, 802, 31, E_Timber_Trans, G_Timber, Area_Perimeter_trans_up, Torsional_J_Perimeter_trans_up, Iy_Perimeter_trans_up, Iz_Perimeter_trans_up, Av_Perimeter_trans_up, Av_Perimeter_trans_up, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam39 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(39, 31, 23, E_Timber_Trans, G_Timber, Area_Perimeter_trans_up, Torsional_J_Perimeter_trans_up, Iy_Perimeter_trans_up, Iz_Perimeter_trans_up, Av_Perimeter_trans_up, Av_Perimeter_trans_up, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam40 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(40, 23, 33, E_Timber_Trans, G_Timber, Area_Perimeter_trans_up, Torsional_J_Perimeter_trans_up, Iy_Perimeter_trans_up, Iz_Perimeter_trans_up, Av_Perimeter_trans_up, Av_Perimeter_trans_up, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam41 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(41, 33, 32, E_Timber_Trans, G_Timber, Area_Perimeter_trans_up, Torsional_J_Perimeter_trans_up, Iy_Perimeter_trans_up, Iz_Perimeter_trans_up, Av_Perimeter_trans_up, Av_Perimeter_trans_up, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam42 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(42, 32, 21, E_Timber_Trans, G_Timber, Area_Perimeter_trans_up, Torsional_J_Perimeter_trans_up, Iy_Perimeter_trans_up, Iz_Perimeter_trans_up, Av_Perimeter_trans_up, Av_Perimeter_trans_up, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam43 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(43, 21, 1302, E_Timber_Trans, G_Timber, Area_Perimeter_trans_up, Torsional_J_Perimeter_trans_up, Iy_Perimeter_trans_up, Iz_Perimeter_trans_up, Av_Perimeter_trans_up, Av_Perimeter_trans_up, lineartag_HorizontalCrdT, 0, 0);

                theDomain.AddElement(ele_timbeam38);
                theDomain.AddElement(ele_timbeam39);
                theDomain.AddElement(ele_timbeam40);
                theDomain.AddElement(ele_timbeam41);
                theDomain.AddElement(ele_timbeam42);
                theDomain.AddElement(ele_timbeam43);
            }

            //# Defining Vertical Inside Elements
            {
                var ele_timbeam44 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(44, 3601, 3201, E_Timber_Long, G_Timber, Area_inside_long, Torsional_J_inside_long, Iy_inside_long, Iz_inside_long, Av_inside_long, Av_inside_long, lineartag_VerticalCrdT, 0, 0);
                var ele_timbeam45 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(45, 3401, 3101, E_Timber_Long, G_Timber, Area_inside_long, Torsional_J_inside_long, Iy_inside_long, Iz_inside_long, Av_inside_long, Av_inside_long, lineartag_VerticalCrdT, 0, 0);

                theDomain.AddElement(ele_timbeam44);
                theDomain.AddElement(ele_timbeam45);
            }

            //# Defining Vertical Inside Elements
            {
                var ele_timbeam46 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(46, 202, 1402, E_Timber_Trans, G_Timber, Area_inside_trans, Torsional_J_inside_trans, Iy_inside_trans, Iz_inside_trans, Av_inside_trans, Av_inside_trans, lineartag_HorizontalCrdT, 0, 0);
                var ele_timbeam47 = new OpenSees.Elements.ElasticTimoshenkoBeam3dWrapper(47, 902, 502, E_Timber_Trans, G_Timber, Area_inside_trans, Torsional_J_inside_trans, Iy_inside_trans, Iz_inside_trans, Av_inside_trans, Av_inside_trans, lineartag_HorizontalCrdT, 0, 0);

                theDomain.AddElement(ele_timbeam46);
                theDomain.AddElement(ele_timbeam47);
            }

            //# Defining Joint Elements
            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Stiffness_TensionMat, Stiffness_IPMat, Stiffness_OOPMat, Stiffness_about_TensionMat, Stiffness_about_OOPMat, Stiffness_about_IPMat };
                var ele_link_tnl48      = new OpenSees.Elements.TwoNodeLinkWrapper(48, 3, 27, 28, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl49      = new OpenSees.Elements.TwoNodeLinkWrapper(49, 3, 29, 30, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl50      = new OpenSees.Elements.TwoNodeLinkWrapper(50, 3, 23, 22, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl51      = new OpenSees.Elements.TwoNodeLinkWrapper(51, 3, 21, 20, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 1, 0, 0 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                theDomain.AddElement(ele_link_tnl48);
                theDomain.AddElement(ele_link_tnl49);
                theDomain.AddElement(ele_link_tnl50);
                theDomain.AddElement(ele_link_tnl51);
            }

            {
                var jointElementMatsSet = new OpenSees.Materials.Uniaxials.UniaxialMaterialWrapper[] { Stiffness_TensionMat, Stiffness_IPMat, Stiffness_OOPMat, Stiffness_about_TensionMat, Stiffness_about_OOPMat, Stiffness_about_IPMat };
                var ele_link_tnl52      = new OpenSees.Elements.TwoNodeLinkWrapper(52, 3, 19, 7, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl53      = new OpenSees.Elements.TwoNodeLinkWrapper(53, 3, 18, 3, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl54      = new OpenSees.Elements.TwoNodeLinkWrapper(54, 3, 17, 10, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);

                var ele_link_tnl55 = new OpenSees.Elements.TwoNodeLinkWrapper(55, 3, 12, 24, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl56 = new OpenSees.Elements.TwoNodeLinkWrapper(56, 3, 6, 25, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);
                var ele_link_tnl57 = new OpenSees.Elements.TwoNodeLinkWrapper(57, 3, 15, 26, dirIds2, jointElementMatsSet, new OpenSees.VectorWrapper(new double[] { 0, 0, 1 }), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), new OpenSees.VectorWrapper(0), 0, 0);

                theDomain.AddElement(ele_link_tnl52);
                theDomain.AddElement(ele_link_tnl53);
                theDomain.AddElement(ele_link_tnl54);

                theDomain.AddElement(ele_link_tnl55);
                theDomain.AddElement(ele_link_tnl56);
                theDomain.AddElement(ele_link_tnl57);
            }

            #endregion

            #region recorders
            var savepath = System.Environment.CurrentDirectory + @"\opsnet_results\";
            if (!System.IO.Directory.Exists(savepath))
            {
                System.IO.Directory.CreateDirectory(savepath);
            }


            {
                var nodeTags = new int[] { 21, 23, 33, 13, 8, 7, 3, 10, 19, 18, 17, 20, 22, 28 };
                foreach (var tag in nodeTags)
                {
                    var opsstream = new OpenSees.Handlers.DataFileStreamWrapper(savepath + $@"\OPSNET_Disp_{tag}.txt");
                    var recorder  = new OpenSees.Recorders.NodeRecorderWrapper(dirIds2, new OpenSees.IDWrapper(new int[] { tag }), 0, "disp", theDomain, opsstream);
                    theDomain.AddRecorder(recorder);
                }
            }

            {
                var nodeTags = new int[] { 27, 29 };
                foreach (var tag in nodeTags)
                {
                    var opsstream = new OpenSees.Handlers.DataFileStreamWrapper(savepath + $@"\OPSNET_Nodal_Reactions_{tag}.txt");
                    var recorder  = new OpenSees.Recorders.NodeRecorderWrapper(dirIds2, new OpenSees.IDWrapper(new int[] { tag }), 0, "reaction", theDomain, opsstream);
                    theDomain.AddRecorder(recorder);
                }
            }

            {
                var eleTags = new int[] { 1, 2 };

                foreach (var tag in eleTags)
                {
                    var opsstream = new OpenSees.Handlers.DataFileStreamWrapper(savepath + $@"\OPSNET_Element_Forces_{tag}.txt");
                    var recorder  = new OpenSees.Recorders.ElementRecorderWrapper(new OpenSees.IDWrapper(new int[] { tag }), new string[] { "force" }, true, theDomain, opsstream, 0, new OpenSees.IDWrapper(0));
                    theDomain.AddRecorder(recorder);
                }
            }
            #endregion

            #region loading
            var theSeries      = new OpenSees.Components.Timeseries.LinearSeriesWrapper();
            var theLoadPattern = new OpenSees.Components.LoadPatterns.LoadPatternWrapper(1);
            theLoadPattern.SetTimeSeries(theSeries);
            theDomain.AddLoadPattern(theLoadPattern);

            var theLoadValues = new OpenSees.VectorWrapper(6);
            theLoadValues.Zero();
            theLoadValues[0] = 9871.8;

            theDomain.AddNodalLoad(new OpenSees.Components.Loads.NodalLoadWrapper(1, 20, theLoadValues, false), 1);
            theDomain.AddNodalLoad(new OpenSees.Components.Loads.NodalLoadWrapper(2, 22, theLoadValues, false), 1);

            var theModel      = new OpenSees.AnalysisModelWrapper();
            var theSolnAlgo   = new OpenSees.Algorithms.NewtonRaphsonWrapper();
            var theIntegrator = new OpenSees.Integrators.Static.LoadControlWrapper(0.0005, 1, 0.0005, 0.0005);
            var theHandler    = new OpenSees.Handlers.PlainHandlerWrapper();
            var theRCM        = new OpenSees.GraphNumberers.RCMWrapper(false);
            var theNumberer   = new OpenSees.Numberers.DOF_NumbererWrapper(theRCM);
            var theSolver     = new OpenSees.Systems.Linears.BandGenLinLapackSolverWrapper();
            var theSOE        = new OpenSees.Systems.Linears.BandGenLinSOEWrapper(theSolver);
            var test          = new OpenSees.ConvergenceTests.CTestNormDispIncrWrapper(1e-12, 20, 2, 2, 1.0e10);
            var theAnalysis   = new OpenSees.Analysis.StaticAnalysisWrapper(
                theDomain,
                theHandler,
                theNumberer,
                theModel,
                theSolnAlgo,
                theSOE,
                theIntegrator,
                test);

            theAnalysis.Analyze(2000);


            Console.ReadKey();
            #endregion
        }