[baezacaljo]

Good afternoon dear all


i am quite new using openFoam and i am facing the next error when i tried to compile my code


error: no match for ‘operator/’ (operand types are ‘Foam::SphericalTensor<double>’ and ‘Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >’)


i am following some tutorials and info about plasticSolver in openfoam. However, they only apply for the previous version of openfoam.


(BTW let me tell you that i took the solidDisplacementFoam solver as a base case and i started doing my modifications)



when i am trying to compile my file called calculateStress.H the previous mentioned error appears.


i am guessing that such error appears when i am trying to make a multiplication by a variable. see the code next.

        forAll(fac, celli)
        {
            if (fac[celli] > 1.0)    //Plasticity occurs
            {
                sigma[celli] = 1.0/3.0*I*tr(sigma[celli]) + dev(sigma[celli])/fac[celli];
                symmTensor dsigma = sigma[celli] - sigma_old[celli];
                tensor deps_e = (1.0/3.0*I*tr(dsigma))/(3.0*lambda+2.0*mu)
                              + dev(dsigma)/(2.0*mu);
                tensor deps = 1.0/2.0*(graddD[celli] + graddD[celli].T());
                deps_p[celli] = deps - deps_e;
            }
            else    // only elasticity
            {
                deps_p[celli] = tensor::zero;
            }
        }

the error relies exactly on the next line

                tensor deps_e = (1.0/3.0*I*tr(dsigma))/(3.0*lambda+2.0*mu)
                              + dev(dsigma)/(2.0*mu);

when i remove the variable volScalarField mu and volScalarField lambda the code compiles.


any suggestions?


what can i do?


Thanks in advance
Alberto

[Yang Zhe]

I am also implementing the same code in openFOAM4.0 now.

but I encountered the following error:

calculateStress.H:43:49: error: no match for ‘operator=’ (operand types are ‘Foam::SymmTensor<double>’ and ‘const Foam::Tensor<double>’)
deps_p[celli] = tensor::zero;

Is there anybody could solve this problem?

[baezacaljo]

Originally Posted by Yang Zhe

I am also implementing the same code in openFOAM4.0 now.

but I encountered the following error:

calculateStress.H:43:49: error: no match for ‘operator=’ (operand types are ‘Foam::SymmTensor<double>’ and ‘const Foam::Tensor<double>’)
deps_p[celli] = tensor::zero;

Is there anybody could solve this problem?

Hi Yang Zhe. Look. Did you already tried to comment out the line

deps_p[celli] = tensor::zero;

if it compiles, change tensor::zero for another value and you will se what's the real problem.

I haven't solve mine, but it seems to me that your code compiles the lines on which I am struggling.

We have to keep doing trial and error.

Best regards
Alberto Campuzano

[Yang Zhe]

Thanks a lot !
It passed by your suggestion.
like this:
deps_p[celli] = I;


As for your problem, I am also facing with it now and I cannot solve it.


I think there are some missing rules controlling the calculation between tensors and scalars in openFOAM 4.0 (at least for 4.0 ). And I have referred to this book before:



http://foam.sourceforge.net/docs/Gui...mmersGuide.pdf


But it did not work. What a pity !


I would be very appreciated if you can recommend some reference books to me !

[baezacaljo]

Originally Posted by Yang Zhe

Thanks a lot !
It passed by your suggestion.
like this:
deps_p[celli] = I;


As for your problem, I am also facing with it now and I cannot solve it.


I think there are some missing rules controlling the calculation between tensors and scalars in openFOAM 4.0 (at least for 4.0 ). And I have referred to this book before:



http://foam.sourceforge.net/docs/Gui...mmersGuide.pdf


But it did not work. What a pity !


I would be very appreciated if you can recommend some reference books to me !

Hi Yang Zhe


I ve been searching for answers. However, nothing comes out.


Take a look, at the next link (Awesome slides from Philip Cardiff)



https://www.researchgate.net/profile...am-toolbox.pdf


Check the slide 109. Look, there are solving the system which we are looking for.


I think that applying the code from old OF versions (lets say you make your own code about solids taking as a references codes from old OF versions) and compiling it in new OF versions is a bad idea, usually compatibility problems are within the path.


what i did is this, compile from source foam extended 4.0 in Ubuntu 20.04 after i compiled solids4Foam, in there already exist the solver that we are looking for.


Notice, that now we have the advantage to take a look of the original code, copy it and try to compile it in OF8, and now we will see if our problem comes with compatibility or was a silly mistake in the code.


Best Regards
Alberto