mirror of
https://github.com/OpenVSP/OpenVSP.git
synced 2026-03-13 10:13:08 +08:00
8d3f61159f
Update CFDMesh.vspscript and FESMesh.vspscript so that they demonstrate running an Analysis the old and the new way
247 lines
8.7 KiB
Plaintext
247 lines
8.7 KiB
Plaintext
|
|
void main()
|
|
{
|
|
Print( string( "Begin FEA Mesh Script\n" ) );
|
|
Print( string( "" ) );
|
|
|
|
string struct_id = GenerateStructure();
|
|
|
|
NewFEAMeshAnalysisExample( struct_id );
|
|
|
|
OldFEAMeshAnalysisExample( struct_id );
|
|
}
|
|
|
|
string GenerateStructure()
|
|
{
|
|
Print( string( "--> Adding Geometry and Creating Structure For Analysis\n" ) );
|
|
|
|
ClearVSPModel();
|
|
|
|
//==== Add Pod Geometry ====//
|
|
string pod_id = AddGeom( "POD" );
|
|
double length = 15.0;
|
|
|
|
SetParmValUpdate( pod_id, "X_Rel_Location", "XForm", 5.0 );
|
|
SetParmValUpdate( pod_id, "X_Rel_Rotation", "XForm", 90 );
|
|
SetParmValUpdate( pod_id, "Length", "Design", length );
|
|
|
|
Update();
|
|
|
|
//==== Set Structure Units ====//
|
|
string veh_id = FindContainer( "Vehicle", 0 );
|
|
SetParmVal( FindParm( veh_id, "StructUnit", "FeaStructure" ), BFT_UNIT );
|
|
|
|
//==== Add FeaStructure to Pod ====//
|
|
int struct_ind = AddFeaStruct( pod_id );
|
|
|
|
//==== Create FeaMaterial ====//
|
|
string mat_id = AddFeaMaterial();
|
|
|
|
SetParmVal( FindParm( mat_id, "MassDensity", "FeaMaterial" ), 0.016 );
|
|
SetParmVal( FindParm( mat_id, "ElasticModulus", "FeaMaterial" ), 2.0e6 );
|
|
SetParmVal( FindParm( mat_id, "PoissonRatio", "FeaMaterial" ), 0.4 );
|
|
|
|
//==== Create FeaProperty ====//
|
|
string prop_id = AddFeaProperty();
|
|
|
|
SetParmVal( FindParm( prop_id, "FeaMaterialIndex", "FeaProperty" ), 4 );
|
|
SetParmVal( FindParm( prop_id, "Thickness", "FeaProperty" ), 0.01 );
|
|
|
|
Update();
|
|
|
|
//==== Add Floor ====//
|
|
string floor_id = AddFeaPart( pod_id, struct_ind, FEA_SLICE );
|
|
|
|
SetParmVal( FindParm( floor_id, "IncludedElements", "FeaPart" ), FEA_SHELL_AND_BEAM );
|
|
SetParmVal( FindParm( floor_id, "RelCenterLocation", "FeaPart" ), 0.34 );
|
|
SetParmVal( FindParm( floor_id, "OrientationPlane", "FeaSlice" ), XZ_BODY );
|
|
|
|
SetParmVal( FindParm( floor_id, "FeaPropertyIndex", "FeaPart" ), 2 );
|
|
|
|
//==== Add Bulkead ====//
|
|
string bulkhead_id = AddFeaPart( pod_id, struct_ind, FEA_SLICE );
|
|
|
|
SetParmVal( FindParm( bulkhead_id, "IncludedElements", "FeaPart" ), FEA_SHELL_AND_BEAM );
|
|
SetParmVal( FindParm( bulkhead_id, "RelCenterLocation", "FeaPart" ), 0.15 );
|
|
SetParmVal( FindParm( bulkhead_id, "OrientationPlane", "FeaSlice" ), SPINE_NORMAL );
|
|
|
|
//==== Add Dome ====//
|
|
string dome_id = AddFeaPart( pod_id, struct_ind, FEA_DOME );
|
|
|
|
SetParmVal( FindParm( dome_id, "IncludedElements", "FeaPart" ), FEA_SHELL );
|
|
SetParmVal( FindParm( dome_id, "X_Location", "FeaDome" ), 0.7 * length );
|
|
SetParmVal( FindParm( dome_id, "A_Radius", "FeaDome" ), 1.5 );
|
|
|
|
//==== Add Stiffener ====//
|
|
string stiffener_id = AddFeaPart( pod_id, struct_ind, FEA_SLICE );
|
|
|
|
SetParmVal( FindParm( stiffener_id, "IncludedElements", "FeaPart" ), FEA_BEAM );
|
|
SetParmVal( FindParm( stiffener_id, "RelCenterLocation", "FeaPart" ), 0.45 );
|
|
SetParmVal( FindParm( stiffener_id, "OrientationPlane", "FeaSlice" ), SPINE_NORMAL );
|
|
|
|
//==== Add LineArray ====//
|
|
string line_array_id = AddFeaSubSurf( pod_id, struct_ind, SS_LINE_ARRAY );
|
|
|
|
SetParmVal( FindParm( line_array_id, "ConstLineType", "SS_LineArray" ), 1 ); // Constant W
|
|
SetParmVal( FindParm( line_array_id, "Spacing", "SS_LineArray" ), 0.25 );
|
|
SetParmVal( FindParm( line_array_id, "StartLocation", "SS_LineArray" ), 0.125 );
|
|
|
|
//==== Add Hole ====//
|
|
string hole_id = AddFeaSubSurf( pod_id, struct_ind, SS_RECTANGLE );
|
|
|
|
SetParmVal( FindParm( hole_id, "IncludedElements", "SubSurface" ), FEA_BEAM );
|
|
SetParmVal( FindParm( hole_id, "Center_U", "SS_Rectangle" ), 0.65 );
|
|
SetParmVal( FindParm( hole_id, "Center_W", "SS_Rectangle" ), 0.5 );
|
|
SetParmVal( FindParm( hole_id, "U_Length", "SS_Rectangle" ), 0.1 );
|
|
SetParmVal( FindParm( hole_id, "W_Length", "SS_Rectangle" ), 0.1 );
|
|
SetParmVal( FindParm( hole_id, "Test_Type", "SS_Rectangle" ), INSIDE );
|
|
|
|
Update();
|
|
|
|
string struct_id = GetFeaStructID( pod_id, struct_ind );
|
|
|
|
//==== Save Vehicle to File ====//
|
|
Print( string( "--> Saving Vehicle to FEAMeshTest.vsp3" ) );
|
|
string fname = "FEAMeshTest.vsp3";
|
|
WriteVSPFile( fname, 0 ); // SET_ALL
|
|
|
|
return struct_id;
|
|
}
|
|
|
|
void NewFEAMeshAnalysisExample( string struct_id )
|
|
{
|
|
Print( string( "This demonstrates how to use the FEA Mesh Analysis\n" ) );
|
|
|
|
int struct_ind = GetFeaStructIndex( struct_id );
|
|
SetFeaMeshStructIndex( struct_ind );
|
|
|
|
//Get analysis type
|
|
string analysis_name = "FeaMeshAnalysis";
|
|
Print( analysis_name );
|
|
|
|
// Set defaults
|
|
SetAnalysisInputDefaults( analysis_name );
|
|
|
|
//Set Inputs
|
|
array < double > baseNums;
|
|
baseNums.push_back( 0.5 );
|
|
SetDoubleAnalysisInput( analysis_name, "BaseLen", baseNums, 0 );
|
|
|
|
array < double > minNums;
|
|
minNums.push_back( 0.1 );
|
|
SetDoubleAnalysisInput( analysis_name, "MinLen", minNums, 0 );
|
|
|
|
array < double > maxGapNums;
|
|
maxGapNums.push_back( 0.005 );
|
|
SetDoubleAnalysisInput( analysis_name, "MaxGap", maxGapNums, 0 );
|
|
|
|
array < double > nCircSegNums;
|
|
nCircSegNums.push_back( 16.0 );
|
|
SetDoubleAnalysisInput( analysis_name, "NCircSeg", nCircSegNums, 0 );
|
|
|
|
array < double > growthRatioNums;
|
|
growthRatioNums.push_back( 1.3 );
|
|
SetDoubleAnalysisInput( analysis_name, "GrowthRatio", growthRatioNums, 0 );
|
|
|
|
array < double > relCurveTolNums;
|
|
relCurveTolNums.push_back( 0.01 );
|
|
SetDoubleAnalysisInput( analysis_name, "RelCurveTol", relCurveTolNums, 0 );
|
|
|
|
array < double > sTEPTolNums;
|
|
sTEPTolNums.push_back( 1e-06 );
|
|
SetDoubleAnalysisInput( analysis_name, "STEPTol", sTEPTolNums, 0 );
|
|
|
|
//Set Files to be Generated
|
|
array < string > stl_file_name;
|
|
stl_file_name.push_back( "Example_Test_STL.stl" );
|
|
SetStringAnalysisInput( analysis_name, "STLFileName", stl_file_name );
|
|
|
|
array < string > gmsh_file_name;
|
|
gmsh_file_name.push_back( "Example_Test_GMSH.gmsh" );
|
|
SetStringAnalysisInput( analysis_name, "GMSHFileName", gmsh_file_name );
|
|
|
|
array < string > mass_file_name;
|
|
mass_file_name.push_back( "Example_Test_MASS.txt" );
|
|
SetStringAnalysisInput( analysis_name, "MASSFileName", mass_file_name );
|
|
|
|
array < string > nastranf_file_name;
|
|
nastranf_file_name.push_back( "Example_Test_NASTRANF.dat" );
|
|
SetStringAnalysisInput( analysis_name, "NASTRANFileName", nastranf_file_name );
|
|
|
|
array < string > nkey_file_name;
|
|
nkey_file_name.push_back( "Example_Test_NKEY.nkey" );
|
|
SetStringAnalysisInput( analysis_name, "NKEYFileName", nkey_file_name );
|
|
|
|
array < string > calculix_file_name;
|
|
calculix_file_name.push_back( "Example_Test_CALCULIX.dat" );
|
|
SetStringAnalysisInput( analysis_name, "CALCULIXFileName", calculix_file_name );
|
|
|
|
array < string > curve_file_name;
|
|
curve_file_name.push_back( "Example_Test_CURVE.curv" );
|
|
SetStringAnalysisInput( analysis_name, "CURVFileName", curve_file_name );
|
|
|
|
array < string > p3d_file_name;
|
|
p3d_file_name.push_back( "Example_Test_P3D.p3d" );
|
|
SetStringAnalysisInput( analysis_name, "P3DFileName", p3d_file_name );
|
|
|
|
array < string > srf_file_name;
|
|
srf_file_name.push_back( "Example_Test_SRF.srf" );
|
|
SetStringAnalysisInput( analysis_name, "SRFFileName", srf_file_name );
|
|
|
|
array < int > iges_file_flag;
|
|
iges_file_flag.push_back( 0 );
|
|
SetIntAnalysisInput( analysis_name, "IGESFileFlag", iges_file_flag );
|
|
|
|
array < int > step_file_flag;
|
|
step_file_flag.push_back( 0 );
|
|
SetIntAnalysisInput( analysis_name, "STEPFileFlag", step_file_flag );
|
|
|
|
// list inputs, type, and current values
|
|
PrintAnalysisInputs( analysis_name );
|
|
|
|
Print( "--> Executing Analysis\n" );
|
|
string resid = ExecAnalysis( analysis_name );
|
|
|
|
// No results generated. See output files
|
|
|
|
//==== Check For API Errors ====//
|
|
while ( GetNumTotalErrors() > 0 )
|
|
{
|
|
ErrorObj err = PopLastError();
|
|
Print( err.GetErrorString() );
|
|
}
|
|
|
|
Print( string( "COMPLETE\n" ) );
|
|
}
|
|
|
|
void OldFEAMeshAnalysisExample( string struct_id )
|
|
{
|
|
Print( string( "This demonstrates how to use the FEA Mesh API functions\n" ) );
|
|
|
|
string parent_id = GetFeaStructParentGeomID( struct_id );
|
|
int struct_ind = GetFeaStructIndex( struct_id );
|
|
|
|
//==== Adjust FeaMeshSettings ====//
|
|
SetFeaMeshVal( parent_id, struct_ind, CFD_MAX_EDGE_LEN, 0.75 );
|
|
SetFeaMeshVal( parent_id, struct_ind, CFD_MIN_EDGE_LEN, 0.2 );
|
|
|
|
//=== Set Export File Name ===//
|
|
string export_name = "FEAMeshTest_calculix.dat";
|
|
SetFeaMeshFileName( parent_id, struct_ind, FEA_CALCULIX_FILE_NAME, export_name );
|
|
Print( string( "--> FEA Mesh Export File Name Set to FEAMeshTest_calculix.dat" ) );
|
|
|
|
//==== Generate FEA Mesh and Export ====//
|
|
Print( string( "--> Generating FeaMesh " ) );
|
|
|
|
ComputeFeaMesh( parent_id, struct_ind, FEA_CALCULIX_FILE_NAME );
|
|
// Could also call ComputeFeaMesh ( struct_id, FEA_CALCULIX_FILE_NAME );
|
|
|
|
//==== Check For API Errors ====//
|
|
while ( GetNumTotalErrors() > 0 )
|
|
{
|
|
ErrorObj err = PopLastError();
|
|
Print( err.GetErrorString() );
|
|
}
|
|
|
|
Print( string( "COMPLETE\n" ) );
|
|
} |