Exemple #1
0
        public void CreateOptiSolidName(ref SAP2000v19.cSapModel Model, ref bool testbool, int arstoix, double xmin, double xmax, double ymin, double ymax, double zmin, double zmax, double diakr, out int OptiSolidNumb, out string[] OptiSolidName)
        {
            int    ret;
            double x_double = 0.0;
            double y_double = 0.0;
            double z_double = 0.0;

            OptiSolidName = new string[arstoix];
            int count   = 0;
            int NumbSel = 0;

            int[] ObjTypeSel;
            ObjTypeSel = new int[1];
            string[] ObjNameSel;
            ObjNameSel = new string[1];
            string voith = "";

            for (x_double = Math.Round(xmax, 3); x_double >= Math.Round(xmin, 3) + Math.Round(diakr / 2, 3); x_double = Math.Round(x_double, 3) - Math.Round(diakr, 3))
            {
                for (z_double = Math.Round(zmax, 3); z_double >= Math.Round(zmin, 3) + Math.Round(diakr / 2, 3); z_double = Math.Round(z_double, 3) - Math.Round(diakr, 3))
                {
                    for (y_double = Math.Round(ymax, 3); y_double >= Math.Round(ymin, 3) + Math.Round(diakr / 2, 3); y_double = Math.Round(y_double, 3) - Math.Round(diakr, 3))
                    {
                        ret = Model.SelectObj.ClearSelection();
                        ret = Model.SelectObj.CoordinateRange(Math.Round(x_double, 3) - Math.Round(diakr, 3), Math.Round(x_double, 3), Math.Round(y_double, 3) - Math.Round(diakr, 3), Math.Round(y_double, 3), Math.Round(z_double, 3) - Math.Round(diakr, 3), Math.Round(z_double, 3), false, "Global", false, false, false, false, true, false);
                        if (ret != 0)
                        {
                            Console.WriteLine("Error: while creating OptiSolidName.");
                        }
                        ret = Model.SelectObj.GetSelected(ref NumbSel, ref ObjTypeSel, ref ObjNameSel);
                        if (ret != 0)
                        {
                            Console.WriteLine("Error: while creating OptiSolidName.");
                        }
                        if (ObjNameSel[0] == voith)
                        {
                            Console.WriteLine("Error: while creating OptiSolidName.");
                        }
                        count = count + 1;
                        OptiSolidName[count - 1] = ObjNameSel[0];
                        voith = ObjNameSel[0];
                    }
                }
            }
            OptiSolidNumb = count;
            if (OptiSolidNumb != arstoix)
            {
                Console.WriteLine("Error: while creating OptiSolidName. OptiSolidNumb is {0} while arstoix is {1}. check if the program didnt read all the elements", OptiSolidNumb, arstoix);
                testbool = false;
            }
        }
Exemple #2
0
        //METHODS
        public void ReadModel(ref SAP2000v19.cSapModel Model, out string[] PointName, out int PointNumb, out string[] SolidName, out int SolidNumb, out string[] FrameName, out int FrameNumb, out string[] AreaName, out int AreaNumb, out string[] LoadName, out int LoadNumb)
        {
            int i;
            int ret;

            //Read points and change their local axes
            PointName = new string[0];
            PointNumb = 0;
            double a = 0.0;
            double b = 0.0;
            double c = 0.0;

            ret = Model.PointObj.GetNameList(ref PointNumb, ref PointName);
            Console.WriteLine("{0} points", PointNumb);
            for (i = 0; i <= PointNumb - 1; i++)
            {
                ret = Model.PointObj.SetLocalAxes(PointName[i], a, b, c, SAP2000v19.eItemType.Objects);
            }
            //Read solid objects and change their local axes
            SolidName = new string[0];
            SolidNumb = 0;
            ret       = Model.SolidObj.GetNameList(ref SolidNumb, ref SolidName);
            Console.WriteLine("{0} solid objects", SolidNumb);
            for (i = 0; i <= SolidNumb - 1; i++)
            {
                ret = Model.SolidObj.SetLocalAxes(SolidName[i], a, b, c, SAP2000v19.eItemType.Objects);
            }
            //Read frame objects and change their local axes
            FrameName = new string[0];
            FrameNumb = 0;
            ret       = Model.FrameObj.GetNameList(ref FrameNumb, ref FrameName);
            Console.WriteLine("{0} Frame objects ", FrameNumb);
            for (i = 0; i <= FrameNumb - 1; i++)
            {
                ret = Model.FrameObj.SetLocalAxes(FrameName[i], 0.0, SAP2000v19.eItemType.Objects);
            }
            //Read areaobject and change their local axes
            AreaName = new string[0];
            AreaNumb = 0;
            ret      = Model.AreaObj.GetNameList(ref AreaNumb, ref AreaName);
            Console.WriteLine("{0} Area objects ", AreaNumb);
            for (i = 0; i <= AreaNumb - 1; i++)
            {
                ret = Model.AreaObj.SetLocalAxes(AreaName[i], 0.0, SAP2000v19.eItemType.Objects);
            }
            //Read loadcases
            LoadNumb = 0;
            LoadName = new string[0];
            ret      = Model.LoadCases.GetNameList(ref LoadNumb, ref LoadName);
        }
Exemple #3
0
        public void UpdateModel(ref SAP2000v19.cSapModel Model, string[] OptiSolidName, int OptiSolidNumb, double[] xkfil)
        {
            int    i;
            int    ret;
            string temp_string;
            double temp_double;

            //update model based on array xkfil
            ret = Model.SetModelIsLocked(false);
            for (i = 0; i <= OptiSolidNumb - 1; i++)
            {
                temp_double = Math.Round(xkfil[i] * xkfil[i] * xkfil[i], 2);
                temp_string = "OptiSolidProp" + string.Format("{0:N2}", temp_double);
                ret         = Model.SolidObj.SetProperty(OptiSolidName[i], temp_string, 0);
            }
        }
Exemple #4
0
        public static void ImportModel(List <int[]> elements, Dictionary <int, double[]> points)

        {
            int ret;

            SAP2000v19.cOAPI   SAPObject = null;
            SAP2000v19.cHelper myHelper  = new SAP2000v19.Helper();

            SAPObject = myHelper.GetObject("CSI.SAP2000.API.SapObject");

            //Get a reference to cSapModel to access all OAPI classes and functions
            SAP2000v19.cSapModel SapModel = SAPObject.SapModel;
            int contadorFrame             = 1;
            int contadorShell             = 1;

            foreach (int[] element in elements)
            {
                double[] firstNode  = points[element[0]];
                double[] secondNode = points[element[1]];

                if (element.Length == 2) // LINE
                {
                    string nameFrame = $"Frame-{contadorFrame++}";
                    ret = SapModel.FrameObj.AddByCoord(firstNode[0], firstNode[1], firstNode[2], secondNode[0], secondNode[1], secondNode[2], ref nameFrame);
                }
                else
                {
                    double[] thirdNode = points[element[2]];
                    double[] xCoords   = new double[element.Length];
                    double[] yCoords   = new double[element.Length];
                    double[] zCoords   = new double[element.Length];

                    for (int i = 0; i < element.Length; i++)
                    {
                        xCoords[i] = points[element[i]][0];
                        yCoords[i] = points[element[i]][1];
                        zCoords[i] = points[element[i]][2];
                    }
                    string nameShell = $"Shell-{contadorShell++}";
                    ret = SapModel.AreaObj.AddByCoord(element.Length, ref xCoords, ref yCoords, ref zCoords, ref nameShell);
                }
            }

            SapModel.View.RefreshView();
        }
Exemple #5
0
        public void DivideSolid(ref SAP2000v19.cSapModel Model, ref bool testbool, double xmin, double xmax, double ymin, double ymax, double zmin, double zmax, int nlx, int nly, int nlz, out string SolidNot, out int SolidNumb, out double diakr)
        {
            string[] SolidName2;
            SolidName2 = new string[0];
            SolidNumb  = 0;
            int ret;
            int NumbSel = 0;

            int[] ObjTypeSel;
            ObjTypeSel = new int[1];
            string[] ObjNameSel;
            ObjNameSel = new string[1];
            double diakrx;
            double diakry;
            double diakrz;

            diakrx = Math.Round((xmax - xmin) / nlx, 2);
            Console.WriteLine("diakrx={0}", diakrx);
            diakry = Math.Round((ymax - ymin) / nly, 2);
            Console.WriteLine("diakry={0}", diakry);
            diakrz = Math.Round((zmax - zmin) / nlz, 2);
            Console.WriteLine("diakrz={0}", diakrz);

            if ((diakrx == diakry) && (diakry == diakrz))
            {
                diakr = diakrx;
            }
            else
            {
                Console.WriteLine("Error in divide proccess");
                testbool = false;
                diakr    = Convert.ToDouble((diakrx + diakry + diakrz) / 3);
            }
            ret = Model.SelectObj.CoordinateRange(xmin, xmax, ymin, ymax, zmin, zmax, false, "Global", false, false, false, false, true, false);
            ret = Model.SelectObj.GetSelected(ref NumbSel, ref ObjTypeSel, ref ObjNameSel);
            //Divide
            ret      = Model.EditSolid.Divide(ObjNameSel[0], nlz, nlx, nly, ref SolidNumb, ref SolidName2);
            SolidNot = ObjNameSel[0];
            if (SolidNumb != nlx * nly * nlz)
            {
                Console.WriteLine("Error in divide proccess");
                testbool = false;
            }
            Console.WriteLine("SolidNumb={0}", SolidNumb);
        }
Exemple #6
0
        public void CreateMatSolid(ref SAP2000v19.cSapModel Model, string Region, string Standard, string Grade)
        {
            int    i;
            int    ret;
            string temp_name;

            string[] MatName;
            MatName = new string[101];
            string temp_mat = "";
            double e        = 0.0;
            double amat     = 0.0;
            double u        = 0.0;
            double g        = 0.0;

            //Create 101 materials with names OptiMat0, OptiMat1..
            for (i = 0; i <= 100; i++)
            {
                temp_name  = "OptiMat" + Convert.ToString(i);
                ret        = Model.PropMaterial.AddMaterial(ref temp_mat, SAP2000v19.eMatType.Concrete, Region, Standard, Grade, temp_name);
                MatName[i] = temp_mat;
            }
            //Change of E in each material
            ret = Model.PropMaterial.GetMPIsotropic(MatName[0], ref e, ref u, ref amat, ref g);
            ret = Model.PropMaterial.SetMPIsotropic(MatName[0], e * 0.0000001, u, amat);
            double metr = 0.01;

            for (i = 1; i <= 100; i++)
            {
                ret  = Model.PropMaterial.SetMPIsotropic(MatName[i], e * metr, u, amat);
                metr = metr + 0.01;
                metr = Math.Round(metr, 2);
            }
            //Create 101 solid properties with names OptiSolidProp0.00, OptiSolidProp0.01..
            metr = 0.00;
            for (i = 0; i <= 100; i++)
            {
                temp_name = "OptiSolidProp" + string.Format("{0:N2}", metr);
                ret       = Model.PropSolid.SetProp(temp_name, MatName[i], 0.0, 0.0, 0.0, true, -1, "", "");
                metr      = Math.Round(metr, 2) + 0.01;
                metr      = Math.Round(metr, 2);
            }
        }
Exemple #7
0
        public void DeleteSolid(ref SAP2000v19.cSapModel Model, string[] OptiSolidName, int OptiSolidNumb, double[] xkfil)
        {
            int    i;
            int    ret;
            double temp_double;
            string temp_string;

            ret = Model.SetModelIsLocked(false);
            for (i = 0; i <= OptiSolidNumb - 1; i++)
            {
                if (xkfil[i] > 0.6)
                {
                    temp_double = Math.Round(xkfil[i] * xkfil[i] * xkfil[i], 2);//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!PROSOXH DEN EXEI ANANEW8EI O PINAKAS X WS TRITH RIZA KTL
                    temp_string = "OptiSolidProp" + string.Format("{0:N2}", temp_double);
                    ret         = Model.SolidObj.SetProperty(OptiSolidName[i], temp_string, 0);
                }
                else
                {
                    ret = Model.SolidObj.Delete(OptiSolidName[i], 0);
                }
            }
        }
Exemple #8
0
        public void CalcDerivative(ref SAP2000v19.cSapModel Model, ref bool testbool, string[] OptiSolidName, int OptiSolidNumb, double[] xkfil, out double[] der, out double comp2)
        {
            //declaration for solid joint force
            int nresultS = 0;

            string[] ObjeS;
            ObjeS = new string[8];
            string[] elmS;
            elmS = new string[8];
            string[] LoadCaseS;
            LoadCaseS = new string[8];
            string[] StepTypeS;
            StepTypeS = new string[8];
            double[] StepNumS;
            StepNumS = new double[8];
            double[] F1S;
            F1S = new double[8];
            double[] F2S;
            F2S = new double[8];
            double[] F3S;
            F3S = new double[8];
            double[] M1S;
            M1S = new double[8];
            double[] M2S;
            M2S = new double[8];
            double[] M3S;
            M3S = new double[8];
            string[] pointelmS;
            pointelmS = new string[8];

            //declaration for Joint Displacements
            int nresult = 0;

            string[] LoadCase;
            LoadCase = new string[1];
            string[] StepType;
            StepType = new string[1];
            double[] StepNum;
            StepNum = new double[1];
            double[] U1;
            U1 = new double[1];
            double[] U2;
            U2 = new double[1];
            double[] U3;
            U3 = new double[1];
            double[] R1;
            R1 = new double[1];
            double[] R2;
            R2 = new double[1];
            double[] R3;
            R3 = new double[1];
            string[] Obje;
            Obje = new string[1];
            string[] elm;
            elm = new string[1];

            //rest of declarations
            int i;
            int j;
            int ret;

            der = new double[OptiSolidNumb];
            double compder;

            comp2 = 0.0;
            for (i = 0; i <= OptiSolidNumb - 1; i++)
            {
                compder = 0.0;
                ret     = Model.Results.SolidJointForce(OptiSolidName[i], SAP2000v19.eItemTypeElm.ObjectElm, ref nresultS, ref ObjeS, ref elmS, ref pointelmS, ref LoadCaseS, ref StepTypeS, ref StepNumS, ref F1S, ref F2S, ref F3S, ref M1S, ref M2S, ref M3S);
                if (nresultS == 8)
                {
                    for (j = 0; j <= nresultS - 1; j++)
                    {
                        ret = Model.Results.JointDispl(pointelmS[j], SAP2000v19.eItemTypeElm.ObjectElm, ref nresult, ref Obje, ref elm, ref LoadCase, ref StepType, ref StepNum, ref U1, ref U2, ref U3, ref R1, ref R2, ref R3);

                        compder = compder + (U1[0] * F1S[j]) + (U2[0] * F2S[j]) + (U3[0] * F3S[j]) + (R1[0] * M1S[j]) + (R2[0] * M2S[j]) + (R3[0] * M3S[j]);
                    }
                }
                else
                {
                    Console.WriteLine("Error: while calculating the compliance of a solid element (Method CalcDerivatives). check nresultS");
                    Console.WriteLine("number of results: {0}", nresultS);
                    testbool = false;
                }
                der[i] = -3 * compder / Math.Max(xkfil[i], Math.Pow(10, -3));
                if (der[i] > 0)
                {
                    Console.WriteLine("Error: positive derivative = {0} in solid element {1}", der[i], OptiSolidName[i]);
                    testbool = false;
                }
                comp2 = comp2 + compder;
            }
        }
Exemple #9
0
        public double CalcCompliance(ref SAP2000v19.cSapModel Model, ref bool testbool, string[] SolidName, string SolidNot, int SolidNumb, string[] FrameName, int FrameNumb, string[] AreaName, int AreaNumb)
        {
            //Declaration for Joint Displacements
            int nresult = 0;

            string[] LoadCase;
            LoadCase = new string[1];
            string[] StepType;
            StepType = new string[1];
            double[] StepNum;
            StepNum = new double[1];
            double[] U1;
            U1 = new double[1];
            double[] U2;
            U2 = new double[1];
            double[] U3;
            U3 = new double[1];
            double[] R1;
            R1 = new double[1];
            double[] R2;
            R2 = new double[1];
            double[] R3;
            R3 = new double[1];
            string[] Obje;
            Obje = new string[1];
            string[] elm;
            elm = new string[1];

            //Declaration for solid joint force
            int nresultS = 0;

            string[] ObjeS;
            ObjeS = new string[8];
            string[] elmS;
            elmS = new string[8];
            string[] LoadCaseS;
            LoadCaseS = new string[8];
            string[] StepTypeS;
            StepTypeS = new string[8];
            double[] StepNumS;
            StepNumS = new double[8];
            double[] F1S;
            F1S = new double[8];
            double[] F2S;
            F2S = new double[8];
            double[] F3S;
            F3S = new double[8];
            double[] M1S;
            M1S = new double[8];
            double[] M2S;
            M2S = new double[8];
            double[] M3S;
            M3S = new double[8];
            string[] pointelmS;
            pointelmS = new string[8];

            //Declaration for frame joint forces
            int nresultF = 0;

            string[] LoadCaseF;
            LoadCaseF = new string[2];
            string[] StepTypeF;
            StepTypeF = new string[2];
            double[] StepNumF;
            StepNumF = new double[2];
            string[] pointelmF;
            pointelmF = new string[2];
            double[] F1F;
            F1F = new double[2];
            double[] F2F;
            F2F = new double[2];
            double[] F3F;
            F3F = new double[2];
            double[] M1F;
            M1F = new double[2];
            double[] M2F;
            M2F = new double[2];
            double[] M3F;
            M3F = new double[2];
            string[] ObjeF;
            ObjeF = new string[2];
            string[] elmF;
            elmF = new string[2];

            //Declaration for AreaJointForceShell and AreaJointForcePlane
            int nresultA = 0;

            string[] LoadCaseA;
            LoadCaseA = new string[4];
            string[] StepTypeA;
            StepTypeA = new string[4];
            double[] StepNumA;
            StepNumA = new double[4];
            string[] pointelmA;
            pointelmA = new string[4];
            double[] F1A;
            F1A = new double[4];
            double[] F2A;
            F2A = new double[4];
            double[] F3A;
            F3A = new double[4];
            double[] M1A;
            M1A = new double[4];
            double[] M2A;
            M2A = new double[4];
            double[] M3A;
            M3A = new double[4];
            string[] ObjeA;
            ObjeA = new string[4];
            string[] elmA;
            elmA = new string[4];
            int    proptype = 0;
            string propname = "";

            //COMPLIANCE
            int    i;
            int    j;
            int    ret;
            double comp = 0.0;

            if (SolidNumb != 0)
            {
                for (i = 0; i <= SolidNumb - 1; i++)
                {
                    if (SolidName[i] != SolidNot)
                    {
                        ret = Model.Results.SolidJointForce(SolidName[i], SAP2000v19.eItemTypeElm.ObjectElm, ref nresultS, ref ObjeS, ref elmS, ref pointelmS, ref LoadCaseS, ref StepTypeS, ref StepNumS, ref F1S, ref F2S, ref F3S, ref M1S, ref M2S, ref M3S);
                        if (nresultS == 8)
                        {
                            for (j = 0; j <= nresultS - 1; j++)
                            {
                                ret  = Model.Results.JointDispl(pointelmS[j], SAP2000v19.eItemTypeElm.ObjectElm, ref nresult, ref Obje, ref elm, ref LoadCase, ref StepType, ref StepNum, ref U1, ref U2, ref U3, ref R1, ref R2, ref R3);
                                comp = comp + U1[0] * F1S[j] + U2[0] * F2S[j] + U3[0] * F3S[j] + R1[0] * M1S[j] + R2[0] * M2S[j] + R3[0] * M3S[j];
                            }
                        }
                        else
                        {
                            Console.WriteLine("Error: while calculating the compliance of a solid element (Method CalcCompliance). check nresultS");
                            Console.WriteLine("number of results: {0}", nresultS);
                            testbool = false;
                        }
                    }
                }
            }
            if (FrameNumb != 0)
            {
                for (i = 0; i <= FrameNumb - 1; i++)
                {
                    ret = Model.Results.FrameJointForce(FrameName[i], SAP2000v19.eItemTypeElm.ObjectElm, ref nresultF, ref ObjeF, ref elmF, ref pointelmF, ref LoadCaseF, ref StepTypeF, ref StepNumF, ref F1F, ref F2F, ref F3F, ref M1F, ref M2F, ref M3F);
                    if (nresultF == 2)
                    {
                        for (j = 0; j <= nresultF - 1; j++)
                        {
                            ret  = Model.Results.JointDispl(pointelmF[j], SAP2000v19.eItemTypeElm.ObjectElm, ref nresult, ref Obje, ref elm, ref LoadCase, ref StepType, ref StepNum, ref U1, ref U2, ref U3, ref R1, ref R2, ref R3);
                            comp = comp + U1[0] * F1F[j] + U2[0] * F2F[j] + U3[0] * F3F[j] + R1[0] * M1F[j] + R2[0] * M2F[j] + R3[0] * M3F[j];
                        }
                    }
                    else
                    {
                        Console.WriteLine("Error: while calculating the compliance of a Frame element (Method CalcCompliance). check nresultF");
                        Console.WriteLine("number of results: {0}", nresultF);
                        testbool = false;
                    }
                }
            }
            if (AreaNumb != 0)
            {
                for (i = 0; i <= AreaNumb - 1; i++)
                {
                    ret = Model.AreaObj.GetProperty(AreaName[i], ref propname);
                    ret = Model.PropArea.GetTypeOAPI(propname, ref proptype); //1 = Shell 2 = Plane 3 = Asolid
                    if (proptype == 1)
                    {
                        ret = Model.Results.AreaJointForceShell(AreaName[i], SAP2000v19.eItemTypeElm.ObjectElm, ref nresultA, ref ObjeA, ref elmA, ref pointelmA, ref LoadCaseA, ref StepTypeA, ref StepNumA, ref F1A, ref F2A, ref F3A, ref M1A, ref M2A, ref M3A);
                        if (nresultA == 4)
                        {
                            for (j = 0; j <= nresultA - 1; j++)
                            {
                                ret = Model.Results.JointDispl(pointelmA[j], SAP2000v19.eItemTypeElm.ObjectElm, ref nresult, ref Obje, ref elm, ref LoadCase, ref StepType, ref StepNum, ref U1, ref U2, ref U3, ref R1, ref R2, ref R3);

                                comp = comp + U1[0] * F1A[j] + U2[0] * F2A[j] + U3[0] * F3A[j] + R1[0] * M1A[j] + R2[0] * M2A[j] + R3[0] * M3A[j];
                            }
                        }
                        else
                        {
                            Console.WriteLine("Error: while calculating the compliance of an area-shell element (Method CalcCompliance). check nresultA");
                            Console.WriteLine("number of results: {0}", nresultA);
                            testbool = false;
                        }
                    }
                    else if (proptype == 2)
                    {
                        ret = Model.Results.AreaJointForcePlane(AreaName[i], SAP2000v19.eItemTypeElm.ObjectElm, ref nresultA, ref ObjeA, ref elmA, ref pointelmA, ref LoadCaseA, ref StepTypeA, ref StepNumA, ref F1A, ref F2A, ref F3A, ref M1A, ref M2A, ref M3A);
                        if (nresultA == 4)
                        {
                            for (j = 0; j <= nresultA - 1; j++)
                            {
                                ret  = Model.Results.JointDispl(pointelmA[j], SAP2000v19.eItemTypeElm.ObjectElm, ref nresult, ref Obje, ref elm, ref LoadCase, ref StepType, ref StepNum, ref U1, ref U2, ref U3, ref R1, ref R2, ref R3);
                                comp = comp + U1[0] * F1A[j] + U2[0] * F2A[j] + U3[0] * F3A[j] + R1[0] * M1A[j] + R2[0] * M2A[j] + R3[0] * M3A[j];
                            }
                        }
                        else
                        {
                            Console.WriteLine("Error: while calculating the compliance of an area-plane element(Method CalcCompliance). check nresultA");
                            Console.WriteLine("number of results: {0}", nresultA);
                            testbool = false;
                        }
                    }
                    else
                    {
                        Console.WriteLine("Error: while calculating the compliance of an area element(Method CalcCompliance). check proptype ");
                        Console.WriteLine("proptype = {0}", proptype);
                        testbool = false;
                    }
                }
            }
            return(comp);
        }
Exemple #10
0
        public void CreateMapDist(ref SAP2000v19.cSapModel Model, int OptiSolidNumb, int nlx, int nly, int nlz, double diakr, double rmin, out int[,] Map, out double[,] Dist)
        {
            int    i;
            int    j;
            int    x_int;
            int    y_int;
            int    z_int;
            double x_double_max;
            double y_double_max;
            double z_double_max;
            double x_double_min;
            double y_double_min;
            double z_double_min;
            double x_double;
            double y_double;
            double z_double;

            Map  = new int[OptiSolidNumb, 27];
            Dist = new double[OptiSolidNumb, 27];
            int    e_int;
            double e_double;
            double H;

            for (i = 0; i <= OptiSolidNumb - 1; i++)
            {
                for (j = 0; j <= 26; j++)
                {
                    Map[i, j]  = 0;
                    Dist[i, j] = 0.0;
                }
            }
            for (x_int = 1; x_int <= nlx; x_int++)
            {
                for (z_int = 1; z_int <= nlz; z_int++)
                {
                    for (y_int = 1; y_int <= nly; y_int++)
                    {
                        e_int        = (x_int - 1) * nlz * nly + (z_int - 1) * nly + y_int;
                        x_double_min = Math.Truncate(Math.Max(x_int - rmin, 0.0)) + 1;
                        z_double_min = Math.Truncate(Math.Max(z_int - rmin, 0.0)) + 1;
                        y_double_min = Math.Truncate(Math.Max(y_int - rmin, 0.0)) + 1;
                        x_double_max = Math.Min(Math.Truncate(x_int + rmin), nlx);
                        z_double_max = Math.Min(Math.Truncate(z_int + rmin), nlz);
                        y_double_max = Math.Min(Math.Truncate(y_int + rmin), nly);
                        j            = 0;
                        for (x_double = x_double_min; x_double <= x_double_max; x_double++)
                        {
                            for (z_double = z_double_min; z_double <= z_double_max; z_double++)
                            {
                                for (y_double = y_double_min; y_double <= y_double_max; y_double++)
                                {
                                    e_double           = (x_double - 1) * nlz * nly + (z_double - 1) * nly + y_double;
                                    H                  = (rmin - Math.Sqrt(Math.Pow(x_int - Convert.ToInt32(x_double), 2) + Math.Pow(y_int - Convert.ToInt32(y_double), 2) + Math.Pow(z_int - Convert.ToInt32(z_double), 2))) * diakr;
                                    Map[e_int - 1, j]  = Convert.ToInt32(e_double);
                                    Dist[e_int - 1, j] = Math.Max(H, 0);
                                    j                  = j + 1;
                                }
                            }
                        }
                    }
                }
            }
        }