版本为 foam-extend3.1
fluidStructureInterface 库下的 fluidStructureInterface.C 文件
void Foam::fluidStructureInterface::calcCurrentSolidZonePoints() const
{
// Find global face zones
if (currentSolidZonePointsPtr_)
{
FatalErrorIn
(
"void fluidStructureInterface::"
"calcCurrentSolidZonePoints() const"
)
<< "Current solid zone points alarady exist"
<< abort(FatalError);
}
currentSolidZonePointsPtr_ =
new vectorField(stress().currentFaceZonePoints(solidZoneIndex()));
}
该函数首先判断固体网格节点是否存在,如果存在则报错,如果不存在则进行计算赋值。
最后的赋值语句不太清楚意思,globe face zones是全部固体网格的意思把。