ABAQUS C++接口子程序及二次开发全攻略

call "C:\Program Files (x86)\Microsoft Visual Studio\2019\Community\VC\Auxiliary\Build\vcvars64.bat"
call "C:\Program Files (x86)\IntelSWTools\compilers_and_libraries_2020.4.311\windows\bin\ifortvars.bat" intel64 vs2019

@call ifortvars.bat intel64 vs2019

@call "C:\Program Files (x86)\IntelSWTools\compilers_and_libraries_2020.4.311\windows\bin\ifortvars.bat" intel64 vs2019

-----------------------------------------------------------------------------
Abaqus make utility with Fortran


         ...PASS

         Continuing...

-----------------------------------------------------------------------------
Abaqus make utility with C++


         ...PASS

         Continuing...

-----------------------------------------------------------------------------

#include <odb_API.h>
#include <iostream>
using namespace std;

int ABQmain(int argc, char **argv)
{
    //Insert user code here
    cout << "This is compiled by abaqus 2019!!!" << endl;

    odb_initializeAPI();
    odb_String odbPath;
    bool ifOdbName = false;
    odb_String elsetName;
    cin.get();
    return 0;
}

Abaqus JOB cpp_test
Begin Compiling User Post-Processing Program
2/10/2022 4:41:38 PM
用于 x64 的 Microsoft (R) C/C++ 优化编译器 19.16.27045 版
版权所有(C) Microsoft Corporation。保留所有权利。

cpp_test.cpp
End Compiling User Post-Processing Program
2/10/2022 4:41:39 PM
Begin Compiling User Post-Processing Program
2/10/2022 4:41:39 PM
用于 x64 的 Microsoft (R) C/C++ 优化编译器 19.16.27045 版
版权所有(C) Microsoft Corporation。保留所有权利。

main_7188.C
End Compiling User Post-Processing Program
2/10/2022 4:41:40 PM
Begin Linking User Post-Processing Program
2/10/2022 4:41:40 PM
2/10/2022 4:41:40 PM
End Linking User Post-Processing Program
Abaqus JOB cpp_test COMPLETED

/***************************************************************
Finding the maximum value of von Mises stress

odbMaxMises.cpp
Code to determine the location and value of the maximum
von-mises stress in an output database.
Usage: abaqus odbMaxMises -odb odbName -elset(optional)
	   elsetName
Requirements:
1. -odb   : Name of the output database.
2. -elset : Name of the assembly level element set.
			Search will be done only for element belonging
			to this set. If this parameter is not provided,
			search will be performed over the entire model.
3. -help  : Print usage
****************************************************************/
#if (defined(HP) && (! defined(HKS_HPUXI)))
#include <iostream.h>
#else
#include <iostream>
using namespace std;
#endif

#include <odb_API.h>
#include <sys/stat.h>
/*
***************
utility functions
***************
*/
bool fileExists(const odb_String  &string);
void rightTrim(odb_String  &string, const char* char_set);
void printExecutionSummary();
/***************************************************************/


int ABQmain(int argc, char **argv)
{
	odb_String odbPath;
	bool ifOdbName = false;
	odb_String elsetName;
	bool ifElset = false;
	odb_Set myElset;
	odb_String region = "over the entire model";
	char msg[256];
	char *abaCmd = argv[0];

	for (int arg = 0; arg < argc; arg++)
	{
		if (strncmp(argv[arg], "-o**", 2) == 0)
		{
			arg++;
			odbPath = argv[arg];
			rightTrim(odbPath, ".odb");
			if (!fileExists(odbPath))
			{
				cerr << "**ERROR** output database  " << odbPath.CStr()
					<< " does not exist\n" << endl;
				exit(1);
			}
			ifOdbName = true;
		}
		else if (strncmp(argv[arg], "-e**", 2) == 0)
		{
			arg++;
			elsetName = argv[arg];
			ifElset = true;
		}
		else if (strncmp(argv[arg], "-h**", 2) == 0)
		{
			printExecutionSummary();
			exit(0);
		}
	}
	if (!ifOdbName)
	{
		cerr << "**ERROR** output database name is not provided\n";
		printExecutionSummary();
		exit(1);
	}

	// Open the output database
	odb_Odb& myOdb = openOdb(odbPath);
	odb_Assembly& myAssembly = myOdb.rootAssembly();

	if (ifElset) {
		if (myAssembly.elementSets().isMember(elsetName)) {
			myElset = myAssembly.elementSets()[elsetName];
			region = " in the element set : " + elsetName;
		}
		else
		{
			cerr << "An assembly level elset " << elsetName.CStr()
				<< " does not exist in the output database :"
				<< myOdb.name().CStr() << endl;
			myOdb.close();
			exit(0);
		}
	}
	
	// Initialize maximum values.
	float maxMises = -0.1;
	int numFV = 0;
	int maxElem = 0;
	odb_String maxStep = "__None__";
	int maxFrame = -1;
	static const odb_String Stress = "S";
	bool isStressPresent = false;
	int numBD = 0, numElems = 0, numIP = 0, numComp = 0, position = 0;

	// Iterate over all available steps
	odb_StepRepository& sRep1 = myOdb.steps();
	odb_StepRepositoryIT sIter1(sRep1);
	for (sIter1.first(); !sIter1.isDone(); sIter1.next())
	{
		odb_Step& step = sRep1[sIter1.currentKey()];
		cout << "Processing Step: " << step.name().CStr() << endl;
		odb_SequenceFrame& frameSequence = step.frames();
		int numFrames = frameSequence.size();
		for (int f = 0; f < numFrames; f++)
		{
			odb_Frame& frame = frameSequence[f];
			odb_FieldOutputRepository& fieldRep = frame.fieldOutputs();
			if (fieldRep.isMember(Stress))
			{
				isStressPresent = true;
				odb_FieldOutput field = fieldRep.get(Stress);
				if (ifElset)
					field = field.getSubset(myElset);
				const odb_SequenceFieldBulkData& seqVal =
					field.bulkDataBlocks();
				int numBlocks = seqVal.size();
				for (int iblock = 0; iblock < numBlocks; iblock++)
				{
					const odb_FieldBulkData& bulkData = seqVal[iblock];
					numBD = bulkData.length();
					numElems = bulkData.numberOfElements();
					numIP = numBD / numElems;
					numComp = bulkData.width();
					float* mises = bulkData.mises();
					int* elementLabels = bulkData.elementLabels();
					int* integrationPoints = bulkData.integrationPoints();
					for (int elem = 0; elem < numElems; elem++)
					{
						for (int ip = 0; ip < numIP; ip++)
						{
							position = elem * numIP + ip;
							float misesData = mises[position];
							if (misesData > maxMises)
							{
								maxMises = misesData;
								maxElem = elementLabels[elem];
								maxStep = step.name();
								maxFrame = frame.incrementNumber();
							}
						}
					}
				}

			}
		}
	}
	cout << "Num of Elements:"<< numElems << endl;
	cout << "Num of ElementNodes:"<< numIP << endl;
	cout << "Num of Comp:"<< numComp << endl;
	if (isStressPresent)
	{
		cout << "Maximum von Mises stress " << region.CStr()
			<< " is " << maxMises << "  in element "
			<< maxElem << endl;
		cout << "Location: frame # " << maxFrame << " step:  "
			<< maxStep.CStr() << endl;
	}
	else
	{
		cout << " Stress output is not available in  the "
			<< "output database : " << myOdb.name().CStr() << endl;
	}
	// close the output database before exiting the program
	myOdb.close();
	return(0);
}

bool fileExists(const odb_String  &string)
{
	bool exists = false;
	struct  stat  buf;
	if (stat(string.CStr(), &buf) == 0)
		exists = true;
	return exists;
}

void rightTrim(odb_String  &string, const char* char_set)
{
	int length = string.Length();
	if (string.Find(char_set) == length)
		string.append(odb_String(char_set));
}

void printExecutionSummary()
{
	cout << " Code to determine the location and value of the\n"
		<< " maximum von-mises stress in an output database.\n"
		<< " Usage: abaqus odbMaxMises -odb odbName \n"
		<< " -elset(optional), -elsetName\n"
		<< " Requirements:\n"
		<< " 1. -odb   : Name of the output database.\n"
		<< " 2. -elset : Name of the assembly level element set.\n"
		<< "             Search will be done only for element \n"
		<< "             belonging to this set.\n"
		<< "             If this parameter is not provided, search \n"
		<< "             will be performed over the entire model.\n"
		<< " 3. -help  : Print usage\n";
}

https://www.gofarlic.com\ProjectData\ABAQUS_DB\202202-cpp-fortran>odbMaxMises.exe -odb Job-TM-E1500-1MPa.odb
Processing Step: Step-1
Num of Elements:417739
Num of ElementNodes:4
Num of Comp:6
Maximum von Mises stress over the entire model is 8.19462  in element 35625
Location: frame # 10 step:  Step-1