public Tuple <Model, Dictionary <int, INode>, double> Reference2RVEExample10_000withRenumbering_mono_hexa()
{
Model model = new Model();
model.SubdomainsDictionary.Add(1, new Subdomain(1));
Dictionary <int, INode> boundaryNodes = new Dictionary <int, INode>();
// COPY APO: Reference2RVEExample100_000withRenumbering_mono_hexa
double[,] Dq = new double[1, 1];
//Tuple<rveMatrixParameters, grapheneSheetParameters> mpgp;
//rveMatrixParameters mp;
//grapheneSheetParameters gp;
renumbering_vector_path = "..\\..\\..\\RveTemplates\\Input\\RveHomogeneous\\REF_new_total_numbering27.txt";
string Fxk_p_komvoi_rve_path = @"C:\Users\turbo-x\Desktop\notes_elegxoi\REFERENCE_fe2_diafora_check\fe2_tax_me1_arxiko_chol_dixws_me1_OriginalRVEExampleChol_me_a1_REF2_10_000_renu_new_multiple_algorithms_check_stress_27hexa\Fxk_p_komvoi_rve.txt";
int subdiscr1 = 1;
int discr1 = 3;
// int discr2 dn xrhsimopoieitai
int discr3 = 3;
int subdiscr1_shell = 7;
int discr1_shell = 1;
mpgp = FEMMeshBuilder.GetReferenceRveExampleParameters(subdiscr1, discr1, discr3, subdiscr1_shell, discr1_shell);
mp = mpgp.Item1;
gp = mpgp.Item2;
double[][] ekk_xyz = new double[2][] { new double[] { 0, 0, 0 }, new double[] { 0.25 * 105, 0, 0.25 * 40 } };
int graphene_sheets_number = 0; // 0 gra sheets afou exoume mono hexa
o_x_parameters[] model_o_x_parameteroi = new o_x_parameters[graphene_sheets_number];
FEMMeshBuilder.LinearHexaElementsOnlyRVEwithRenumbering_forMS(model, mp, Dq, renumbering_vector_path, boundaryNodes);
double volume = mp.L01 * mp.L02 * mp.L03;
return(new Tuple <Model, Dictionary <int, INode>, double>(model, boundaryNodes, volume));
}
private Tuple <Model, Dictionary <int, Node>, double> Reference2RVEExample10000withRenumberingwithInput_forMS()
{
Model model = new Model();
model.SubdomainsDictionary.Add(1, new Subdomain(1));
Dictionary <int, Node> boundaryNodes = new Dictionary <int, Node>();
//Origin public static void Reference2RVEExample10000withRenumberingwithInput(Model model)
double[,] Dq;
//Tuple<rveMatrixParameters, grapheneSheetParameters> mpgp;
//rveMatrixParameters mp;
//grapheneSheetParameters gp;
var rve_id_data = RVE_id.ToString();
renumbering_vector_path = "..\\..\\..\\RveTemplates\\Input\\RveGrShMultiple\\rve_no_{0}\\REF_new_total_numbering.txt";
renumbering_vector_path = string.Format(renumbering_vector_path, rve_id_data);
string Fxk_p_komvoi_rve_path = "..\\..\\..\\RveTemplates\\Input\\RveGrShMultiple\\rve_no_{0}\\Fxk_p_komvoi_rve.txt";
Fxk_p_komvoi_rve_path = string.Format(Fxk_p_komvoi_rve_path, rve_id_data);
string o_xsunol_input_path_gen = "..\\..\\..\\RveTemplates\\Input\\RveGrShMultiple\\rve_no_{0}\\o_xsunol_gs_";
o_xsunol_input_path_gen = string.Format(o_xsunol_input_path_gen, rve_id_data);
o_xsunol_input_path_gen = o_xsunol_input_path_gen + "{0}.txt";
int subdiscr1 = 2;
int discr1 = 4;
// int discr2 dn xrhsimopoieitai
int discr3 = 8;
int subdiscr1_shell = 6;
int discr1_shell = 1;
mpgp = FEMMeshBuilder.GetReferenceKanonikhGewmetriaRveExampleParametersStiffCase(subdiscr1, discr1, discr3, subdiscr1_shell, discr1_shell);
mp = mpgp.Item1; //mp.hexa1 = 9; mp.hexa2 = 9; mp.hexa3 = 9;
gp = mpgp.Item2;
int graphene_sheets_number = 3;
o_x_parameters[] model_o_x_parameteroi = new o_x_parameters[graphene_sheets_number];
double[][] ekk_xyz = new double[graphene_sheets_number][];
Dq = new double[9, 3 * (((mp.hexa1 + 1) * (mp.hexa2 + 1) * (mp.hexa3 + 1)) - ((mp.hexa1 - 1) * (mp.hexa2 - 1) * (mp.hexa3 - 1)))];
FEMMeshBuilder.HexaElementsOnlyRVEwithRenumbering_forMS(model, mp, Dq, renumbering_vector_path, boundaryNodes);
double volume = mp.L01 * mp.L02 * mp.L03;
int hexaElementsNumber = model.ElementsDictionary.Count();
//IEnumerable<Element> hostGroup = model.ElementsDictionary.Where(x => (x.Key < hexaElementsNumber + 1)).Select(kv => kv.Value);
List <int> EmbeddedElementsIDs = new List <int>();
int element_counter_after_Adding_sheet;
element_counter_after_Adding_sheet = hexaElementsNumber; // initial value before adding first graphene sheet
int shellElementsNumber;
for (int j = 0; j < graphene_sheets_number; j++)
{
string file_no = (j + 1).ToString();
string ox_sunol_input_path = string.Format(o_xsunol_input_path_gen, file_no);
FEMMeshBuilder.AddGrapheneSheet_with_o_x_Input_withRenumberingBondSlip(model, gp, ekk_xyz[j], model_o_x_parameteroi[j], renumbering_vector_path, ox_sunol_input_path);
shellElementsNumber = (model.ElementsDictionary.Count() - element_counter_after_Adding_sheet) / 3; //tha xrhsimefsei
for (int k = shellElementsNumber + element_counter_after_Adding_sheet + 1; k < model.ElementsDictionary.Count() + 1; k++)
{
EmbeddedElementsIDs.Add(model.ElementsDictionary[k].ID);
}
element_counter_after_Adding_sheet = model.ElementsDictionary.Count();
}
int[] EmbElementsIds = EmbeddedElementsIDs.ToArray();
IEnumerable <Element> embdeddedGroup = model.ElementsDictionary.Where(x => (Array.IndexOf(EmbElementsIds, x.Key) > -1)).Select(kv => kv.Value); // dld einai null afth th stigmh
//var embeddedGrouping = new EmbeddedCohesiveGrouping(model, hostGroup, embdeddedGroup);
//var CohesiveGroupings = new EmbeddedCohesiveGrouping[EmbElementsIds.GetLength(0)];
var hostSubGroups = new Dictionary <int, IEnumerable <Element> >();
for (int i1 = 0; i1 < EmbElementsIds.GetLength(0); i1++)
{
hostSubGroups.Add(EmbElementsIds[i1], FEMMeshBuilder.GetHostGroupForCohesiveElement(model.ElementsDictionary[EmbElementsIds[i1]], mp, model, renumbering_vector_path));
//var embeddedGroup_i1 = new List<Element>(1) { model.ElementsDictionary[EmbElementsIds[i1]] };
//CohesiveGroupings[i1] = new EmbeddedCohesiveGrouping(model, hostGroup_i1, embeddedGroup_i1);
}
var CohesiveGroupping = new EmbeddedCohesiveSubGrouping(model, hostSubGroups, embdeddedGroup);
return(new Tuple <Model, Dictionary <int, Node>, double>(model, boundaryNodes, volume));
}
private Tuple <Model, Dictionary <int, Node>, double> Reference2RVEExample10000withRenumberingwithInput_forMS()
{
Model model = new Model();
model.SubdomainsDictionary.Add(1, new Subdomain(1));
Dictionary <int, Node> boundaryNodes = new Dictionary <int, Node>();
//Origin public static void Reference2RVEExample10000withRenumberingwithInput(Model model)
double[,] Dq;
//Tuple<rveMatrixParameters, grapheneSheetParameters> mpgp;
//rveMatrixParameters mp;
//grapheneSheetParameters gp;
var rve_id_data = RVE_id.ToString();
renumbering_vector_path = "..\\..\\..\\RveTemplates\\Input\\RveGrShMultiple\\rve_no_{0}\\REF_new_total_numbering.txt";
renumbering_vector_path = string.Format(renumbering_vector_path, rve_id_data);
string Fxk_p_komvoi_rve_path = "..\\..\\..\\RveTemplates\\Input\\RveGrShMultiple\\rve_no_{0}\\Fxk_p_komvoi_rve.txt";
Fxk_p_komvoi_rve_path = string.Format(Fxk_p_komvoi_rve_path, rve_id_data);
string o_xsunol_input_path_gen = "..\\..\\..\\RveTemplates\\Input\\RveGrShMultiple\\rve_no_{0}\\o_xsunol_gs_";
o_xsunol_input_path_gen = string.Format(o_xsunol_input_path_gen, rve_id_data);
o_xsunol_input_path_gen = o_xsunol_input_path_gen + "{0}.txt";
string subdomainOutputPath_gen = @"C:\Users\turbo-x\Desktop\notes_elegxoi\REFERENCE_kanonikh_gewmetria_fe2_post_dg\REF2_10__000_renu_new_multiple_algorithms_check_develop_gia_fe2_3grsh_4182dofs_multiple2\RVE_database\rve_no_{0}";
string subdomainOutputPath = string.Format(subdomainOutputPath_gen, rve_id_data);
int subdiscr1 = 2;
int discr1 = 4;
// int discr2 dn xrhsimopoieitai
int discr3 = 8;
int subdiscr1_shell = 6;
int discr1_shell = 1;
mpgp = FEMMeshBuilder.GetReferenceKanonikhGewmetriaRveExampleParametersStiffCase(subdiscr1, discr1, discr3, subdiscr1_shell, discr1_shell);
mp = mpgp.Item1; //mp.hexa1 = 9; mp.hexa2 = 9; mp.hexa3 = 9;
gp = mpgp.Item2;
int graphene_sheets_number = 3;
o_x_parameters[] model_o_x_parameteroi = new o_x_parameters[graphene_sheets_number];
double[][] ekk_xyz = new double[graphene_sheets_number][];
Dq = new double[9, 3 * (((mp.hexa1 + 1) * (mp.hexa2 + 1) * (mp.hexa3 + 1)) - ((mp.hexa1 - 1) * (mp.hexa2 - 1) * (mp.hexa3 - 1)))];
FEMMeshBuilder.HexaElementsOnlyRVEwithRenumbering_forMS(model, mp, Dq, renumbering_vector_path, boundaryNodes);
//domain separation ds1
int totalSubdomains = 8;
DdmCalculationsGeneral.BuildModelInterconnectionData(model);
var decomposer = new AutomaticDomainDecomposer2(model, totalSubdomains);
decomposer.UpdateModel();
var subdHexaIds = DdmCalculationsGeneral.DetermineHexaElementsSubdomainsFromModel(model);
double volume = mp.L01 * mp.L02 * mp.L03;
int hexaElementsNumber = model.ElementsDictionary.Count();
//IEnumerable<Element> hostGroup = model.ElementsDictionary.Where(x => (x.Key < hexaElementsNumber + 1)).Select(kv => kv.Value);
List <int> EmbeddedElementsIDs = new List <int>();
int element_counter_after_Adding_sheet;
element_counter_after_Adding_sheet = hexaElementsNumber; // initial value before adding first graphene sheet
int shellElementsNumber;
//ds2
int[] lowerCohesiveBound = new int[graphene_sheets_number]; int[] upperCohesiveBound = new int[graphene_sheets_number]; int[] grShElementssnumber = new int[graphene_sheets_number];
for (int j = 0; j < graphene_sheets_number; j++)
{
string file_no = (j + 1).ToString();
string ox_sunol_input_path = string.Format(o_xsunol_input_path_gen, file_no);
FEMMeshBuilder.AddGrapheneSheet_with_o_x_Input_withRenumberingBondSlip(model, gp, ekk_xyz[j], model_o_x_parameteroi[j], renumbering_vector_path, ox_sunol_input_path);
shellElementsNumber = (model.ElementsDictionary.Count() - element_counter_after_Adding_sheet) / 3; //tha xrhsimefsei
lowerCohesiveBound[j] = shellElementsNumber + element_counter_after_Adding_sheet + 1; //ds3
upperCohesiveBound[j] = 2 * shellElementsNumber + element_counter_after_Adding_sheet;
grShElementssnumber[j] = shellElementsNumber;
for (int k = shellElementsNumber + element_counter_after_Adding_sheet + 1; k < model.ElementsDictionary.Count() + 1; k++)
{
EmbeddedElementsIDs.Add(model.ElementsDictionary[k].ID);
}
element_counter_after_Adding_sheet = model.ElementsDictionary.Count();
}
int[] EmbElementsIds = EmbeddedElementsIDs.ToArray();
IEnumerable <Element> embdeddedGroup = model.ElementsDictionary.Where(x => (Array.IndexOf(EmbElementsIds, x.Key) > -1)).Select(kv => kv.Value); // dld einai null afth th stigmh
//var embeddedGrouping = new EmbeddedCohesiveGrouping(model, hostGroup, embdeddedGroup);
//var CohesiveGroupings = new EmbeddedCohesiveGrouping[EmbElementsIds.GetLength(0)];
var hostSubGroups = new Dictionary <int, IEnumerable <Element> >();
for (int i1 = 0; i1 < EmbElementsIds.GetLength(0); i1++)
{
hostSubGroups.Add(EmbElementsIds[i1], FEMMeshBuilder.GetHostGroupForCohesiveElement(model.ElementsDictionary[EmbElementsIds[i1]], mp, model, renumbering_vector_path));
//var embeddedGroup_i1 = new List<Element>(1) { model.ElementsDictionary[EmbElementsIds[i1]] };
//CohesiveGroupings[i1] = new EmbeddedCohesiveGrouping(model, hostGroup_i1, embeddedGroup_i1);
}
var CohesiveGroupping = new EmbeddedCohesiveSubGrouping(model, hostSubGroups, embdeddedGroup);
//ds4
int[][] subdCohElementIds = DdmCalculationsGeneral.DetermineCoheiveELementsSubdomainsSimple(model, totalSubdomains);
int[][] subdShellElementIds = DdmCalculationsGeneral.DetermineShellELementsSubdomains(model, totalSubdomains, subdCohElementIds,
lowerCohesiveBound, upperCohesiveBound, grShElementssnumber);
int[][] subdElementIds1 = DdmCalculationsGeneral.CombineSubdomainElementsIdsArraysIntoOne(subdHexaIds, subdCohElementIds);
int[][] subdElementIds2 = DdmCalculationsGeneral.CombineSubdomainElementsIdsArraysIntoOne(subdElementIds1, subdShellElementIds);
DdmCalculationsGeneral.UndoModelInterconnectionDataBuild(model);
DdmCalculations.SeparateSubdomains(model, subdElementIds2);
bool print_subdomain_data = false;
if (print_subdomain_data)
{
DdmCalculationsGeneral.PrintSubdomainDataForPostPro(subdHexaIds, subdCohElementIds, subdShellElementIds, subdomainOutputPath);
}
bool get_subdomain_data = true;
if (get_subdomain_data)
{
(hexaPrint, cohePrint, shellPrint) = DdmCalculationsGeneral.GetSubdomainDataForPostPro(subdHexaIds, subdCohElementIds, subdShellElementIds, subdomainOutputPath);
}
return(new Tuple <Model, Dictionary <int, Node>, double>(model, boundaryNodes, volume));
}
public Dictionary <Node, IList <IDofType> > GetModelRigidBodyNodeConstraints(Model model)
{
return(FEMMeshBuilder.GetConstraintsOfDegenerateRVEForNonSingularStiffnessMatrix_withRenumbering(model, mp.hexa1, mp.hexa2, mp.hexa3, renumbering_vector_path));
//TODO: Pithanws na epistrefetai apo GetModelAndBoundaryNodes ... AndConstraints.
}
