[Yinghao Zhou]

Hello everyone!
I am using solidDisplacementFoam to solve thermal stress. I built a cube model and used the subsetMesh command to dig out some of the mesh to form a porous structure. But I got an error when I ran it.

I tried to merge the oldInternalFaces formed by the grid into other patches or define it as a separate patch, but it still reported an error. But I can run it without subsetMesh.
could anyone help me?
The following is the error message

Create time

Create mesh for time = 0

Reading mechanical properties

Normalising E : E/rho

Calculating Lame's coefficients

Plane Strain

Reading thermal properties

Normalising k : k/rho

Calculating thermal coefficients

Reading field D

Reading field T

Calculating stress field sigmaD

Calculating explicit part of div(sigma) divSigmaExp


Calculating displacement field

Iteration: 1

DICPCG:  Solving for T, Initial residual = 0.625951, Final residual = 0.0112757, No Iterations 17
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib64/libc.so.6"
#3  Foam::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) at ??:?
#4  Foam::DICPreconditioner::DICPreconditioner(Foam::lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#5  Foam::lduMatrix::preconditioner::addsymMatrixConstructorToTable<Foam::DICPreconditioner>::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#6  Foam::lduMatrix::preconditioner::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#7  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8  ? at ??:?
#9  ? at ??:?
#10  ? at ??:?
#11  ? at ??:?
#12  __libc_start_main in "/lib64/libc.so.6"
#13  ? at ??:?

[IMG]