search

somdattagoswami / IGAPack-PhaseField

Second and fourth-order adaptive phase field modeling of fracture using PHT-splines in the framework of IGA.
43 stars 16 forks source link

Reference to non-existent field 'constk'. #1

Closed asdasdqwdqwfsdf closed 4 years ago

asdasdqwdqwfsdf commented 4 years ago

Hello,

I may find an error, i.e.,


CantileverBeam
Deciding the phase field model: 
 Second-order model : 1 
 Fourth-order model : 2 
 Your choice: 2
Deciding the stress degradation function: 
 1: Quadratic stress-degradation function 
 2: Cubic stress-degradation function 
 Your choice: 2
Initializing boundary conditions on the initial geometry.
Precomputing shape functions and derivatives.
History function and phase field initialization.

istep =

     1

Assembling the stiffness matrix.
Reference to non-existent field 'constk'.

Error in gStiffnessUUcubic (line 38)
                    gphi = Fract.s*((1.0-phigp).^3 -(1.0-phigp).^2) + 3*(1.0-phigp).^2
                    - 2*(1.0-phigp).^3 +Fract.constk;

Error in
CantileverBeam>@(PHTelem,sizeBasis,numberElements,dgdx,shape,Fract,Mater,volume,tdisp,geometry)gStiffnessUUcubic(PHTelem,sizeBasis,numberElements,dgdx,shape,Fract,Mater,volume,tdisp,geometry)

Error in solver_4th (line 48)
        [stiffUU,elemRef] =
        PF.StiffUU(PHTelem,sizeBasis,numberElements,dgdx,shape,Fract,Mater,volume,tdisp,geometry);

Error in CantileverBeam (line 68)
        solver_4th(Input,PF,file_name);

Best regards!

somdattagoswami commented 4 years ago

We have fixed the error. Thank you for letting us know.