private void button_create_Geom_Click(object sender, RoutedEventArgs e)
{
#region a test example, not important
m = 9;
width = 8 / 1000.0;
thickness = 1 / 1000.0;
l1 = 30 / 1000.0;
l2 = 16 / 1000.0;
l3 = 21 / 1000.0;
F1 = 1;
F2 = 37;
F3 = 49;
l4 = l2;
l1_ = l1 - 2 * width;
l2_ = l2;
l3_ = l3 + 2 * width;
l4_ = l4;
cell_length = l1 + l3;
cell_height = l2 + width;
#endregion
Count = int.Parse(textBox_count.Text);
Random randomSeed = new Random();
Random random_m = new Random(randomSeed.Next());
Random random_width = new Random(randomSeed.Next());
Random random_thickness = new Random(randomSeed.Next());
Random random_l1 = new Random(randomSeed.Next());
Random random_l2 = new Random(randomSeed.Next());
Random random_l3 = new Random(randomSeed.Next());
Random random_F1 = new Random(randomSeed.Next());
Random random_F2 = new Random(randomSeed.Next());
Random random_F3 = new Random(randomSeed.Next());
for (int j = 0; j < Count; j++)
{
Console.WriteLine(string.Format("---------------------------------now No.{0}", j));
#region generate random Data
try
{
/* change to double
* ran_m = random_m.Next(int.Parse(textBox_hm_cells_min.Text), int.Parse(textBox_hm_cellls_max.Text));
* ran_width = random_width.Next(int.Parse(textBox_width_min.Text), int.Parse(textBox_width_max.Text));
* ran_thickness = random_thickness.Next(int.Parse(textBox_thickness_min.Text), int.Parse(textBox_thickness_max.Text));
* ran_l1 = random_l1.Next(int.Parse(textBox_l1_min.Text), int.Parse(textBox_l1_max.Text));
* ran_l2 = random_l2.Next(int.Parse(textBox_l2_min.Text), int.Parse(textBox_l2_max.Text));
* ran_l3 = random_l3.Next(int.Parse(textBox_l3_min.Text), int.Parse(textBox_l3_max.Text));
* ran_F1 = random_F1.Next(int.Parse(textBox_F1_min.Text), int.Parse(textBox_F1_max.Text));
* ran_F2 = random_F2.Next(int.Parse(textBox_F2_min.Text), int.Parse(textBox_F2_max.Text));
* ran_F3 = random_F3.Next(int.Parse(textBox_F3_min.Text), int.Parse(textBox_F3_max.Text));
*/
ran_m = random_m.Next(int.Parse(textBox_hm_cells_min.Text), int.Parse(textBox_hm_cellls_max.Text));
ran_width = random_width.NextDouble() * (int.Parse(textBox_width_max.Text) - int.Parse(textBox_width_min.Text)) + double.Parse(textBox_width_min.Text);
ran_thickness = random_thickness.NextDouble() * (int.Parse(textBox_thickness_max.Text) - int.Parse(textBox_thickness_min.Text)) + double.Parse(textBox_thickness_min.Text);
ran_l1 = random_l1.NextDouble() * (int.Parse(textBox_l1_max.Text) - int.Parse(textBox_l1_min.Text)) + double.Parse(textBox_l1_min.Text);
ran_l2 = random_l2.NextDouble() * (int.Parse(textBox_l2_max.Text) - int.Parse(textBox_l2_min.Text)) + double.Parse(textBox_l2_min.Text);
ran_l3 = random_l3.NextDouble() * (int.Parse(textBox_l3_max.Text) - int.Parse(textBox_l3_min.Text)) + double.Parse(textBox_l3_min.Text);
ran_F1 = random_F1.NextDouble() * (int.Parse(textBox_F1_max.Text) - int.Parse(textBox_F1_min.Text)) + double.Parse(textBox_F1_min.Text);
ran_F2 = random_F2.NextDouble() * (int.Parse(textBox_F2_max.Text) - int.Parse(textBox_F2_min.Text)) + double.Parse(textBox_F2_min.Text);
ran_F3 = random_F3.NextDouble() * (int.Parse(textBox_F3_max.Text) - int.Parse(textBox_F3_min.Text)) + double.Parse(textBox_F3_min.Text);
}
catch (Exception)
{
Console.WriteLine("Input Format Error");
return;
}
m = ran_m;
width = (double)(ran_width / 1000.0);
thickness = (double)(ran_thickness / 1000.0);
l1 = (double)(ran_l1 / 1000.0);
l2 = (double)(ran_l2 / 1000.0);
l3 = (double)(ran_l3 / 1000.0);
F1 = (double)ran_F1;
F2 = (double)ran_F2;
F3 = (double)ran_F3;
F3 = 0;
l4 = l2;
l1_ = l1 - 2 * width;
l2_ = l2;
l3_ = l3 + 2 * width;
l4_ = l4;
cell_length = l1 + l3;
cell_height = l2 + width;
#endregion
Console.WriteLine(string.Format("{0},{1},{2},{3},{4},{5}", m, width, thickness, l1, l2, l3));
#region geometrie
Console.WriteLine("Start creating new Germetrie");
try
{
swModel = swApp.NewPart();
}
catch (Exception)
{
Console.WriteLine("starting SolidWorks");
try
{
swApp = (SldWorks)Marshal.GetActiveObject("SldWorks.Application");
}
catch (Exception)
{
swApp = new SldWorks();
swApp.Visible = false;
}
//get MaterialLib
strMaterialLib = swApp.GetExecutablePath() + "\\lang\\english\\sldmaterials\\solidworks materials.sldmat";
Console.WriteLine("SolidWorks successfully started");
j--;
continue;
}
swModel.Extension.SelectByID("前视基准面", "PLANE", 0, 0, 0, false, 1, null);
swModel.InsertSketch2(true);
#region sketch
for (int i = 0; i < m; i++)
{
if (i == 0)
{
swModel.SketchManager.CreateLine(0, 0, 0, l1 - width, 0, 0);
swModel.SketchManager.CreateLine(l1 - width, 0, 0, l1 - width, -l2, 0);
swModel.SketchManager.CreateLine(l1 - width, -l2, 0, l1 + l3 - width, -l2, 0);
swModel.SketchManager.CreateLine(l1 + l3 - width, -l2, 0, l1 + l3 - width, 0, 0);
swModel.SketchManager.CreateLine(0, 0, 0, 0, -width, 0);
swModel.SketchManager.CreateLine(0, -width, 0, l1_, -width, 0);
swModel.SketchManager.CreateLine(l1_, -width, 0, l1_, -(width + l2_), 0);
swModel.SketchManager.CreateLine(l1_, -(width + l2_), 0, l1_ + l3_, -(width + l2_), 0);
swModel.SketchManager.CreateLine(l1_ + l3_, -(width + l2_), 0, l1_ + l3_, -width, 0);
}
else
{
x_o = i * cell_length - width;
y_o = 0;
x_o_ = x_o + width;
y_o_ = y_o - width;
swModel.SketchManager.CreateLine(x_o, y_o, 0, x_o + l1, 0, 0);
swModel.SketchManager.CreateLine(x_o + l1, 0, 0, x_o + l1, -l2, 0);
swModel.SketchManager.CreateLine(x_o + l1, -l2, 0, x_o + l1 + l3, -l2, 0);
swModel.SketchManager.CreateLine(x_o + l1 + l3, -l2, 0, x_o + l1 + l3, 0, 0);
swModel.SketchManager.CreateLine(x_o_, y_o_, 0, x_o_ + l1_, y_o_, 0);
swModel.SketchManager.CreateLine(x_o_ + l1_, y_o_, 0, x_o_ + l1_, y_o_ - l2_, 0);
swModel.SketchManager.CreateLine(x_o_ + l1_, y_o_ - l2_, 0, x_o_ + l1_ + l3_, y_o_ - l2_, 0);
swModel.SketchManager.CreateLine(x_o_ + l1_ + l3_, y_o_ - l2_, 0, x_o_ + l1_ + l3_, y_o_, 0);
}
}
swModel.SketchManager.CreateLine(x_o_ + l1_ + l3_, y_o_, 0, x_o_ + l1_ + l3_, 0, 0);
swModel.SketchManager.CreateLine(x_o_ + l1_ + l3_, 0, 0, x_o + l1 + l3, 0, 0);
#endregion
swModel.FeatureManager.FeatureExtrusion2(
true, false, false, 0, 0, thickness, 0, false, false, false, false, 0, 0, false, false, false, false, true, true, true, 0, 0, true
);
swModel.SaveAsSilent("D:\\TUD\\7.Semeter\\SA\\SA_code\\c#\\W_Form_analyse_get_TrainingDaten\\Geometrie.sldprt", true);
swApp.CloseAllDocuments(true);
swApp.OpenDoc("D:\\TUD\\7.Semeter\\SA\\SA_code\\c#\\W_Form_analyse_get_TrainingDaten\\Geometrie.sldprt", (int)swOpenDocOptions_e.swOpenDocOptions_Silent);
Console.WriteLine("Geometrie success");
#endregion
#region simulaiton
string path_to_cosworks_dll = @"C:\Program Files\SOLIDWORKS Corp\SOLIDWORKS\Simulation\cosworks.dll";
errors = swApp.LoadAddIn(path_to_cosworks_dll);
COSMOSObject = (CwAddincallback)swApp.GetAddInObject("SldWorks.Simulation");
try
{
COSMOSWORKS = (CosmosWorks)COSMOSObject.CosmosWorks;
}
catch (Exception)
{
Console.WriteLine("something wrong in Simulaiton Add In, start a new one");
continue;
}
COSMOSWORKS = COSMOSObject.CosmosWorks;
//Get active document
ActDoc = (CWModelDoc)COSMOSWORKS.ActiveDoc;
//Create new static study
StudyMngr = (CWStudyManager)ActDoc.StudyManager;
Study = (CWStudy)StudyMngr.CreateNewStudy("static study", (int)swsAnalysisStudyType_e.swsAnalysisStudyTypeStatic, 0, out errCode);
//Add materials
SolidMgr = Study.SolidManager;
SolidComp = SolidMgr.GetComponentAt(0, out errCode);
SolidBody = SolidComp.GetSolidBodyAt(0, out errCode);
intStatus = SolidBody.SetLibraryMaterial(strMaterialLib, "AISI 1020");
//fixed restraints
LBCMgr = Study.LoadsAndRestraintsManager;
swModel = (ModelDoc2)swApp.ActiveDoc;
swModel.ShowNamedView2("", (int)swStandardViews_e.swIsometricView);
selectionMgr = (SelectionMgr)swModel.SelectionManager;
isSelected = swModel.Extension.SelectByID2("", "FACE", 0, -width / 2.0, thickness / 2.0, false, 0, null, 0);
if (isSelected)
{
object selectedFace = (object)selectionMgr.GetSelectedObject6(1, -1);
object[] fixedFaces = { selectedFace };
CWRestraint restraint = (CWRestraint)LBCMgr.AddRestraint((int)swsRestraintType_e.swsRestraintTypeFixed, fixedFaces, null, out errCode);
}
swModel.ClearSelection2(true);
//add force
selectionMgr = (SelectionMgr)swModel.SelectionManager;
isSelected = swModel.Extension.SelectByID2("", "FACE", x_o + l1 + l3 + (width / 2.0), 0, thickness / 2.0, false, 0, null, 0);
if (isSelected)
{
object selectedFace = (object)selectionMgr.GetSelectedObject6(1, -1);
object[] forceAdd = { selectedFace };
selectionMgr = (SelectionMgr)swModel.SelectionManager;
swModel.Extension.SelectByID2("", "FACE", x_o + l1 + l3 + width, -(l4 + width) / 2.0, thickness / 2.0, false, 0, null, 0);
object selectedFaceToForceDir = (object)selectionMgr.GetSelectedObject6(1, -1);
double[] distValue = null;
double[] forceValue = null;
double[] Force = { F2, F3, F1 };
cwForce = (CWForce)LBCMgr.AddForce3((int)swsForceType_e.swsForceTypeForceOrMoment, (int)swsSelectionType_e.swsSelectionFaceEdgeVertexPoint,
2, 0, 0, 0,
(distValue),
(forceValue),
false, true,
(int)swsBeamNonUniformLoadDef_e.swsTotalLoad,
0, 7, 0.0,
Force,
false, false,
(forceAdd),
(selectedFaceToForceDir),
false, out errCode); //i have tried to figure out these arguments for one day, keep them and dont't change them.
//the way to check cwForce : cwForce.GetForceComponentValues()
//ForceComponet: [int b1, // 1 if x-direction hat Force-komponent, else 0
// int b2, // 1 if y-direction hat Force-komponent, else 0
// int b3, // 1 if z-direction hat Force-komponent, else 0
// double d1, // Force-komponent in x
// double d2, // Force-komponent in y
// double d3 // Force-komponent in z
// ]
//PS: the definition of xyz seems like not the same as the global XYZ system in SW.
swModel.ClearSelection2(true);
//meshing
CWMesh CWMeshObj = default(CWMesh);
CWMeshObj = Study.Mesh;
CWMeshObj.MesherType = (int)swsMesherType_e.swsMesherTypeStandard;
CWMeshObj.Quality = (int)swsMeshQuality_e.swsMeshQualityDraft;
errCode = Study.CreateMesh(0, MeshEleSize, MeshTol);
CWMeshObj = null;
//run analysis
errCode = Study.RunAnalysis();
if (errCode != 0)
{
Console.WriteLine(string.Format("RunAnalysis errCode = {0}", errCode));
Console.WriteLine(string.Format("RunAnalysis failed"));
swApp.CloseAllDocuments(true);
errors = swApp.UnloadAddIn(path_to_cosworks_dll);
Console.WriteLine(string.Format("ready to start a new one"));
continue;
}
Console.WriteLine("RunAnalysis successed, ready to get results");
//get results
CWFeatobj = Study.Results;
//get max von Mieses stress
Stress = (object[])CWFeatobj.GetMinMaxStress((int)swsStressComponent_e.swsStressComponentVON,
0, 0, null,
(int)swsStrengthUnit_e.swsStrengthUnitNewtonPerSquareMillimeter,
out errCode);
maxStress = (float)Stress[3]; //Stress: {node_with_minimum_stress, minimum_stress, node_with_maximum_stress, maximum_stress}
Console.WriteLine(maxStress);
/*
* if (maxStress >= 351.6)
* {
* Console.WriteLine("out of yield stress, start a new example");
* errors = swApp.UnloadAddIn(path_to_cosworks_dll);
* swApp.CloseAllDocuments(true);
* j--;
* continue;
* }
*/
//get max URES displacement
Disp = (object[])CWFeatobj.GetMinMaxDisplacement((int)swsDisplacementComponent_e.swsDisplacementComponentURES,
0, null,
(int)swsLinearUnit_e.swsLinearUnitMillimeters,
out errCode);
maxDisp = (float)Disp[3]; //Disp: {node_with_minimum_displacement, minimum_displacement, node_with_maximum_displacement, maximum_displacement}
CWFeatobj = null;
Console.WriteLine(string.Format("max Displacement: {0:f4} mm", maxDisp));
//output to Excel
exlSheet.Cells[j + 2, 1] = m;
exlSheet.Cells[j + 2, 2] = width * 1000;
exlSheet.Cells[j + 2, 3] = thickness * 1000;
exlSheet.Cells[j + 2, 4] = l1 * 1000;
exlSheet.Cells[j + 2, 5] = l2 * 1000;
exlSheet.Cells[j + 2, 6] = l3 * 1000;
exlSheet.Cells[j + 2, 7] = F1;
exlSheet.Cells[j + 2, 8] = F2;
exlSheet.Cells[j + 2, 9] = F3;
exlSheet.Cells[j + 2, 10] = maxStress;
exlSheet.Cells[j + 2, 11] = maxDisp;
if (j % 5 == 0)
{
exlBook.Save();//C:\Users\Zhao\Documents\工作簿1.xlsx
}
errors = swApp.UnloadAddIn(path_to_cosworks_dll);
swApp.CloseAllDocuments(true);
}
#endregion
}
exlBook.SaveCopyAs("D:\\TUD\\7.Semeter\\SA\\SA_code\\c#\\W_Form_analyse_get_TrainingDaten\\W_Form_simulationDaten_" + GetTimeStamp() + "_F3=0" + ".xlsx");
exlBook.Save();//C:\Users\zhaojie\Documents\工作簿1.xlsx
exlApp.Quit();
}
// another thread work
void backgroundWorker_DoWork(object sender, DoWorkEventArgs e)
{
BackgroundWorker worker = sender as BackgroundWorker;
// banning empty fields
if (string.IsNullOrWhiteSpace(strBox.Text))
{
MessageBox.Show("Enter force value", "Error",
MessageBoxButtons.OK, MessageBoxIcon.Warning);
return;
}
processNameLabel.Invoke((MethodInvoker) delegate {
processNameLabel.Text = "Launching SolidWorks...";
});
// check for manual cancelling of analysis
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate {
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
// opening SW 2017
Guid Guid1 = new Guid("6B36082E-677B-49E8-BCB2-76698EBD2906");
object processSW = System.Activator.CreateInstance(System.Type.GetTypeFromCLSID(Guid1));
swApp = (SldWorks)processSW;
swApp.Visible = true;
processNameLabel.Invoke((MethodInvoker) delegate {
processNameLabel.Text = "Preparing Simulation...";
});
CosmosWorks COSMOSWORKS = default(CosmosWorks);
CwAddincallback CWAddinCallBack = default(CwAddincallback);
CWModelDoc ActDoc = default(CWModelDoc);
CWStudyManager StudyMngr = default(CWStudyManager);
CWStudy Study = default(CWStudy);
CWStaticStudyOptions CWOptions = default(CWStaticStudyOptions);
CWLoadsAndRestraintsManager LBCMgr = default(CWLoadsAndRestraintsManager);
CWRestraint CWBind = default(CWRestraint);
CWForce CWStr = default(CWForce);
CWMesh CWMesh = default(CWMesh);
CWResults CWResult = default(CWResults);
CWPlot plot = default(CWPlot);
ModelDoc2 Part = default(ModelDoc2);
ModelDoc2 swDoc = null;
int longstatus = 0;
int longwarnings = 0;
int errCode = 0;
object[] results = null;
// analysis name
string study_name = studyNameBox.Text;
if (string.IsNullOrWhiteSpace(study_name))
{
study_name = "Static_beam";
}
bwStudy.ReportProgress(5);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate {
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
processNameLabel.Invoke((MethodInvoker) delegate {
processNameLabel.Text = "Creating beam model...";
});
// opening SLDPRT file containing beam model
Part = swApp.OpenDoc6("C:\\Beam.SLDPRT", 1, 0, "", ref longstatus, ref longwarnings);
if (Part == null)
{
MessageBox.Show("Failed to open model");
Application.Exit();
}
string path_to_cosworks_dll = swApp.GetExecutablePath() + @"\Simulation\cosworks.dll";
longstatus = swApp.LoadAddIn(path_to_cosworks_dll);
// isometric view
swDoc = swApp.IActiveDoc2;
swDoc.ShowNamedView2("*Isometric", 7);
swDoc.ViewZoomtofit2();
bwStudy.ReportProgress(10);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Starting Simulation...";
});
// ready-check Simulation add-on
CWAddinCallBack = (CwAddincallback)swApp.GetAddInObject("CosmosWorks.CosmosWorks");
if (CWAddinCallBack == null)
{
MessageBox.Show("Failed to start Simulation");
Application.Exit();
}
bwStudy.ReportProgress(12);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
// ready-check add-on object
COSMOSWORKS = CWAddinCallBack.CosmosWorks;
if (COSMOSWORKS == null)
{
MessageBox.Show("Failed to create CosmosWorks object");
Application.Exit();
}
bwStudy.ReportProgress(14);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
// connecting add-on to document
ActDoc = COSMOSWORKS.ActiveDoc;
if (ActDoc == null)
{
MessageBox.Show("Failed to connect to document");
Application.Exit();
}
bwStudy.ReportProgress(16);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Creating analysis...";
});
// creating static analysis
StudyMngr = ActDoc.StudyManager;
if (StudyMngr == null)
{
MessageBox.Show("Failed to start analysis manager");
Application.Exit();
}
bwStudy.ReportProgress(18);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
Study = StudyMngr.CreateNewStudy3(study_name, (int)swsAnalysisStudyType_e.swsAnalysisStudyTypeStatic, (int)swsMeshType_e.swsMeshTypeSolid, out errCode);
if (Study == null)
{
MessageBox.Show("Failed to create static analysis");
Application.Exit();
}
bwStudy.ReportProgress(20);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Adding attachment object...";
});
// loading indicator manager
LBCMgr = Study.LoadsAndRestraintsManager;
if (LBCMgr == null)
{
MessageBox.Show("Failed to load indicator manager");
Application.Exit();
}
bwStudy.ReportProgress(25);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
// selecting face for attachment and adding it to indicator manager
ISelectionMgr selectionMgr = (ISelectionMgr)swDoc.SelectionManager;
bool isSelected = swDoc.Extension.SelectByRay(-0.023870052576342005, 0.064003232713616853, 0, -0.57738154519998108, -0.57728771208554841, -0.57738154519997897, 0.0037366684901689287, 2, true, 0, 0);
if (isSelected)
{
object face_bind = selectionMgr.GetSelectedObject6(1, -1);
object[] fixedpart = { face_bind };
CWBind = (CWRestraint)LBCMgr.AddRestraint(0, (fixedpart), null, out errCode);
if (errCode != 0)
{
MessageBox.Show("Failed to add attachment object");
Application.Exit();
}
}
bwStudy.ReportProgress(35);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Adding force indicator...";
});
swDoc.ClearSelection2(true);
isSelected = swDoc.Extension.SelectByRay(0.00094849719312151137, 0.029999999999915872, 0.97601905600942018, -0.57738154519998108, -0.57728771208554841, -0.57738154519997897, 0.0037366684901689287, 2, true, 0, 0);
if (isSelected)
{
object face_force = selectionMgr.GetSelectedObject6(1, -1);
object[] fixedpart = { face_force };
object[] DistanceValues = { };
object[] ForceValues = { };
object[] ComponentValues = { };
CWStr = (CWForce)LBCMgr.AddForce3(1, 0, -1, 0, 0, 0, DistanceValues, ForceValues, false, false, 0, 0, 0, Convert.ToInt32(strBox.Text), ComponentValues, false, false, (fixedpart), null, false, out errCode);
if (errCode != 0)
{
MessageBox.Show("Failed to add force indicator");
Application.Exit();
}
CWStr.ForceBeginEdit();
CWStr.Unit = 2;
CWStr.ForceEndEdit();
}
bwStudy.ReportProgress(50);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Creating mesh object...";
});
swDoc.ClearSelection2(true);
// mesh object parameters
CWMesh = ((CWMesh)(Study.Mesh));
CWMesh.Quality = 1;
CWMesh.UseJacobianCheckForSolids = 1;
CWMesh.MesherType = 0;
CWMesh.AutomaticTransition = 0;
CWMesh.AutomaticLooping = 0;
CWMesh.NumberOfLoops = 3;
CWMesh.SaveSettingsWithoutMeshing = 0;
CWMesh.Unit = 0;
bwStudy.ReportProgress(70);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
// creating mesh object
errCode = Study.CreateMesh(0, 13.3821, 0.669103);
if (errCode != 0)
{
MessageBox.Show("Failed to create mesh object");
Application.Exit();
}
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Starting analysis...";
});
// starting analysis
errCode = Study.RunAnalysis();
if (errCode != 0)
{
MessageBox.Show("Failed to start analysis");
Application.Exit();
}
bwStudy.ReportProgress(90);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Receiving analysis results...";
});
// obtaining results
CWResult = (CWResults)Study.Results;
if (CWResult == null)
{
MessageBox.Show("Failed to receive results");
return;
}
// stress graph results
plot = CWResult.GetPlot("Stress1", out errCode);
if (plot == null)
{
MessageBox.Show("Failed to receive stress results");
return;
}
plot.SetComponentUnitAndValueByElem(9, 3, false);
results = (object[])plot.GetMinMaxResultValues(out errCode);
max_tens = Convert.ToDouble(results[3].ToString());
min_tens = Convert.ToDouble(results[1].ToString());
bwStudy.ReportProgress(93);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
// displacement graph results
plot = CWResult.GetPlot("Displacement1", out errCode);
if (plot == null)
{
MessageBox.Show("Failed to receive displacement results");
return;
}
plot.SetComponentUnitAndValueByElem(3, 0, false);
results = (object[])plot.GetMinMaxResultValues(out errCode);
max_move = Convert.ToDouble(results[3].ToString());
min_move = Convert.ToDouble(results[1].ToString());
bwStudy.ReportProgress(96);
if (bwStudy.CancellationPending == true)
{
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Analysis has been cancelled";
});
e.Cancel = true;
return;
}
// strain graph results
plot = CWResult.GetPlot("Strain1", out errCode);
if (plot == null)
{
MessageBox.Show("Failed to receive strain results");
return;
}
plot.SetComponentUnitAndValueByElem(6, 3, false);
results = (object[])plot.GetMinMaxResultValues(out errCode);
bwStudy.ReportProgress(99);
max_deform = Convert.ToDouble(results[3].ToString());
min_deform = Convert.ToDouble(results[1].ToString());
maxTensBox.Invoke((MethodInvoker) delegate
{
maxTensBox.Text = max_tens.ToString();
});
minTensBox.Invoke((MethodInvoker) delegate
{
minTensBox.Text = min_tens.ToString();
});
maxMoveBox.Invoke((MethodInvoker) delegate
{
maxMoveBox.Text = max_move.ToString();
});
minMoveBox.Invoke((MethodInvoker) delegate
{
minMoveBox.Text = min_move.ToString();
});
maxDeformBox.Invoke((MethodInvoker) delegate
{
maxDeformBox.Text = max_deform.ToString();
});
minDeformBox.Invoke((MethodInvoker) delegate
{
minDeformBox.Text = min_deform.ToString();
});
bwStudy.ReportProgress(100);
processNameLabel.Invoke((MethodInvoker) delegate
{
processNameLabel.Text = "Done";
});
}
private void Button_generateFEM_Click(object sender, RoutedEventArgs e)
{
#region Variable definition
int m; //how many cells
int Count; // how many data
double width, thickness, l1, l2, l3, F1, F2, F3;
double l4, l1_, l2_, l3_, l4_, cell_length, cell_height;
double x_o = 0.0;
double y_o = 0.0;
double x_o_ = 0.0;
double y_o_ = 0.0;
const double MeshEleSize = 2.0;
const double MeshTol = 0.1;
string strMaterialLib = null;
object[] Disp = null;
object[] Stress = null;
ModelDoc2 swModel = null;
SelectionMgr selectionMgr = null;
CosmosWorks COSMOSWORKS = null;
CwAddincallback COSMOSObject = default(CwAddincallback);
CWModelDoc ActDoc = default(CWModelDoc);
CWStudyManager StudyMngr = default(CWStudyManager);
CWStudy Study = default(CWStudy);
CWSolidManager SolidMgr = default(CWSolidManager);
CWSolidBody SolidBody = default(CWSolidBody);
CWSolidComponent SolidComp = default(CWSolidComponent);
CWForce cwForce = default(CWForce);
CWLoadsAndRestraintsManager LBCMgr = default(CWLoadsAndRestraintsManager);
CWResults CWFeatobj = default(CWResults);
bool isSelected;
float maxDisp = 0.0f;
float maxStress = 0.0f;
int intStatus = 0;
int errors = 0;
int errCode = 0;
int warnings = 0;
int ran_m, ran_width, ran_thickness, ran_l1, ran_l2, ran_l3, ran_F1, ran_F2, ran_F3;
Excel.Application exlApp;
Excel.Workbook exlBook;
Excel.Worksheet exlSheet;
#endregion
#region read data from UI
try
{
m = int.Parse(textBox_hm_cells.Text);
width = double.Parse(textBox_width.Text) / 1000.0;
thickness = double.Parse(textBox_thickness.Text) / 1000.0;
l1 = double.Parse(textBox_l1.Text) / 1000.0;
l2 = double.Parse(textBox_l2.Text) / 1000.0;
l3 = double.Parse(textBox_l3.Text) / 1000.0;
F1 = double.Parse(textBox_F1.Text);
F2 = double.Parse(textBox_F2.Text);
F3 = double.Parse(textBox_F3.Text);
}
catch (Exception)
{
MessageBox.Show("input format error, check input data");
return;
}
l4 = l2;
l1_ = l1 - 2 * width;
l2_ = l2;
l3_ = l3 + 2 * width;
l4_ = l4;
cell_length = l1 + l3;
cell_height = l2 + width;
#endregion
#region geometrie
Console.WriteLine("Start creating new Germetrie");
swModel = swApp.NewPart();
swModel.Extension.SelectByID("前视基准面", "PLANE", 0, 0, 0, false, 1, null);
swModel.InsertSketch2(true);
#region sketch
for (int i = 0; i < m; i++)
{
if (i == 0)
{
swModel.SketchManager.CreateLine(0, 0, 0, l1 - width, 0, 0);
swModel.SketchManager.CreateLine(l1 - width, 0, 0, l1 - width, -l2, 0);
swModel.SketchManager.CreateLine(l1 - width, -l2, 0, l1 + l3 - width, -l2, 0);
swModel.SketchManager.CreateLine(l1 + l3 - width, -l2, 0, l1 + l3 - width, 0, 0);
swModel.SketchManager.CreateLine(0, 0, 0, 0, -width, 0);
swModel.SketchManager.CreateLine(0, -width, 0, l1_, -width, 0);
swModel.SketchManager.CreateLine(l1_, -width, 0, l1_, -(width + l2_), 0);
swModel.SketchManager.CreateLine(l1_, -(width + l2_), 0, l1_ + l3_, -(width + l2_), 0);
swModel.SketchManager.CreateLine(l1_ + l3_, -(width + l2_), 0, l1_ + l3_, -width, 0);
}
else
{
x_o = i * cell_length - width;
y_o = 0;
x_o_ = x_o + width;
y_o_ = y_o - width;
swModel.SketchManager.CreateLine(x_o, y_o, 0, x_o + l1, 0, 0);
swModel.SketchManager.CreateLine(x_o + l1, 0, 0, x_o + l1, -l2, 0);
swModel.SketchManager.CreateLine(x_o + l1, -l2, 0, x_o + l1 + l3, -l2, 0);
swModel.SketchManager.CreateLine(x_o + l1 + l3, -l2, 0, x_o + l1 + l3, 0, 0);
swModel.SketchManager.CreateLine(x_o_, y_o_, 0, x_o_ + l1_, y_o_, 0);
swModel.SketchManager.CreateLine(x_o_ + l1_, y_o_, 0, x_o_ + l1_, y_o_ - l2_, 0);
swModel.SketchManager.CreateLine(x_o_ + l1_, y_o_ - l2_, 0, x_o_ + l1_ + l3_, y_o_ - l2_, 0);
swModel.SketchManager.CreateLine(x_o_ + l1_ + l3_, y_o_ - l2_, 0, x_o_ + l1_ + l3_, y_o_, 0);
}
}
swModel.SketchManager.CreateLine(x_o_ + l1_ + l3_, y_o_, 0, x_o_ + l1_ + l3_, 0, 0);
swModel.SketchManager.CreateLine(x_o_ + l1_ + l3_, 0, 0, x_o + l1 + l3, 0, 0);
#endregion
swModel.FeatureManager.FeatureExtrusion2(
true, false, false, 0, 0, thickness, 0, false, false, false, false, 0, 0, false, false, false, false, true, true, true, 0, 0, true
);
swModel.SaveAsSilent("D:\\TUD\\7.Semeter\\SA\\SA_code\\c#\\W_Form_analyse_FEM_NN_GUI\\Geometrie.sldprt", true);
swApp.CloseAllDocuments(true);
swApp.OpenDoc("D:\\TUD\\7.Semeter\\SA\\SA_code\\c#\\W_Form_analyse_FEM_NN_GUI\\Geometrie.sldprt", (int)swOpenDocOptions_e.swOpenDocOptions_Silent);
Console.WriteLine("Geometrie success");
#endregion
#region simulaiton
string path_to_cosworks_dll = @"C:\Program Files\SOLIDWORKS Corp\SOLIDWORKS\Simulation\cosworks.dll";
errors = swApp.LoadAddIn(path_to_cosworks_dll);
COSMOSObject = (CwAddincallback)swApp.GetAddInObject("SldWorks.Simulation");
try
{
COSMOSWORKS = (CosmosWorks)COSMOSObject.CosmosWorks;
}
catch (Exception)
{
Console.WriteLine("something wrong in Simulaiton Add In, start a new one");
swApp.CloseAllDocuments(true);
return;
}
COSMOSWORKS = COSMOSObject.CosmosWorks;
//Get active document
ActDoc = (CWModelDoc)COSMOSWORKS.ActiveDoc;
//Create new static study
StudyMngr = (CWStudyManager)ActDoc.StudyManager;
Study = (CWStudy)StudyMngr.CreateNewStudy("static study", (int)swsAnalysisStudyType_e.swsAnalysisStudyTypeStatic, 0, out errCode);
//Add materials
//get MaterialLib
strMaterialLib = swApp.GetExecutablePath() + "\\lang\\english\\sldmaterials\\solidworks materials.sldmat";
SolidMgr = Study.SolidManager;
SolidComp = SolidMgr.GetComponentAt(0, out errCode);
SolidBody = SolidComp.GetSolidBodyAt(0, out errCode);
intStatus = SolidBody.SetLibraryMaterial(strMaterialLib, "AISI 1020");
//fixed restraints
LBCMgr = Study.LoadsAndRestraintsManager;
swModel = (ModelDoc2)swApp.ActiveDoc;
swModel.ShowNamedView2("", (int)swStandardViews_e.swIsometricView);
selectionMgr = (SelectionMgr)swModel.SelectionManager;
isSelected = swModel.Extension.SelectByID2("", "FACE", 0, -width / 2.0, thickness / 2.0, false, 0, null, 0);
if (isSelected)
{
object selectedFace = (object)selectionMgr.GetSelectedObject6(1, -1);
object[] fixedFaces = { selectedFace };
CWRestraint restraint = (CWRestraint)LBCMgr.AddRestraint((int)swsRestraintType_e.swsRestraintTypeFixed, fixedFaces, null, out errCode);
}
swModel.ClearSelection2(true);
//add force
selectionMgr = (SelectionMgr)swModel.SelectionManager;
isSelected = swModel.Extension.SelectByID2("", "FACE", x_o + l1 + l3 + (width / 2.0), 0, thickness / 2.0, false, 0, null, 0);
if (isSelected)
{
object selectedFace = (object)selectionMgr.GetSelectedObject6(1, -1);
object[] forceAdd = { selectedFace };
selectionMgr = (SelectionMgr)swModel.SelectionManager;
swModel.Extension.SelectByID2("", "FACE", x_o + l1 + l3 + width, -(l4 + width) / 2.0, thickness / 2.0, false, 0, null, 0);
object selectedFaceToForceDir = (object)selectionMgr.GetSelectedObject6(1, -1);
double[] distValue = null;
double[] forceValue = null;
double[] Force = { F2, F3, F1 };
cwForce = (CWForce)LBCMgr.AddForce3((int)swsForceType_e.swsForceTypeForceOrMoment, (int)swsSelectionType_e.swsSelectionFaceEdgeVertexPoint,
2, 0, 0, 0,
(distValue),
(forceValue),
false, true,
(int)swsBeamNonUniformLoadDef_e.swsTotalLoad,
0, 7, 0.0,
Force,
false, false,
(forceAdd),
(selectedFaceToForceDir),
false, out errCode); //i have tried to figure out these arguments for one day, keep them and dont't change them.
//the way to check cwForce : cwForce.GetForceComponentValues()
//ForceComponet: [int b1, // 1 if x-direction hat Force-komponent, else 0
// int b2, // 1 if y-direction hat Force-komponent, else 0
// int b3, // 1 if z-direction hat Force-komponent, else 0
// double d1, // Force-komponent in x
// double d2, // Force-komponent in y
// double d3 // Force-komponent in z
// ]
//PS: the definition of xyz seems like not the same as the global XYZ system in SW.
swModel.ClearSelection2(true);
//meshing
CWMesh CWMeshObj = default(CWMesh);
CWMeshObj = Study.Mesh;
CWMeshObj.MesherType = (int)swsMesherType_e.swsMesherTypeStandard;
CWMeshObj.Quality = (int)swsMeshQuality_e.swsMeshQualityDraft;
errCode = Study.CreateMesh(0, MeshEleSize, MeshTol);
CWMeshObj = null;
//run analysis
errCode = Study.RunAnalysis();
if (errCode != 0)
{
Console.WriteLine(string.Format("RunAnalysis errCode = {0}", errCode));
Console.WriteLine(string.Format("RunAnalysis failed"));
errors = swApp.UnloadAddIn(path_to_cosworks_dll);
swApp.CloseAllDocuments(true);
Console.WriteLine(string.Format("please start a new one"));
return;
}
Console.WriteLine("RunAnalysis successed, ready to get results");
//get results
CWFeatobj = Study.Results;
//get max von Mieses stress
Stress = (object[])CWFeatobj.GetMinMaxStress((int)swsStressComponent_e.swsStressComponentVON,
0, 0, null,
(int)swsStrengthUnit_e.swsStrengthUnitNewtonPerSquareMillimeter,
out errCode);
maxStress = (float)Stress[3]; //Stress: {node_with_minimum_stress, minimum_stress, node_with_maximum_stress, maximum_stress}
Console.WriteLine(maxStress);
if (maxStress >= 351.6)
{
Console.WriteLine("out of yield stress, start a new example");
errors = swApp.UnloadAddIn(path_to_cosworks_dll);
swApp.CloseAllDocuments(true);
textBox_FEM_result.Text = "out of yield stress";
return;
}
//get max URES displacement
Disp = (object[])CWFeatobj.GetMinMaxDisplacement((int)swsDisplacementComponent_e.swsDisplacementComponentURES,
0, null,
(int)swsLinearUnit_e.swsLinearUnitMillimeters,
out errCode);
maxDisp = (float)Disp[3]; //Disp: {node_with_minimum_displacement, minimum_displacement, node_with_maximum_displacement, maximum_displacement}
CWFeatobj = null;
Console.WriteLine(string.Format("max Displacement: {0:f4} mm", maxDisp));
textBox_FEM_result.Text = maxDisp.ToString();
}
else
{
Console.WriteLine("not selected");
}
#endregion
errors = swApp.UnloadAddIn(path_to_cosworks_dll);
swApp.CloseAllDocuments(true);
}