Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

  • Congratulations waross on being selected by the Tek-Tips community for having the most helpful posts in the forums last week. Way to Go!

Conversion to particles formulation (SPH) using user subroutine–based criterion

Status
Not open for further replies.

SuperSmile

Mechanical
Oct 31, 2014
2
0
0
GB
Hello to everybody,
I wrote a VUMAT subroutine for a linear elastic orthotropic material.
In it I defined a state variable controlling the element conversion. When the condition for the conversion is satisfied, the SDV variable is set to be zero.
The subroutine works fine until the conversion criterion is reached. When that happens I get the following error: Abaqus Error: The executable explicit_dp.exe aborted with system error code 1073741819.

I tried to fix it renaming the mkl_avx2.dll in (C:\SIMULIA\Abaqus\6.13-x\code\bin) to mkl_avx2.dll.11.0.0.1, but it does not work to me.
I reported the subroutine I wrote below hoping you can help me to understand where I am wrong:

do i=1,nblock
C STRESS UPDATE
stressNew(i,1)=stressOld(i,1)+strainInc(i,1)*C11+strainInc(i,2)*C12+strainInc(i,3)*C13
stressNew(i,2)=stressOld(i,2)+strainInc(i,1)*C21+strainInc(i,2)*C22+strainInc(i,3)*C23
stressNew(i,3)=stressOld(i,3)+strainInc(i,1)*C31+strainInc(i,2)*C32+strainInc(i,3)*C33
stressNew(i,4)=stressOld(i,4)+strainInc(i,4)*C44
stressNew(i,5)=stressOld(i,5)+strainInc(i,5)*C55
stressNew(i,6)=stressOld(i,6)+strainInc(i,6)*C66

if (stressNew(i,1).GT.1000) then
stateNew(i,8)=zero
endif
end do

I really appreciate any help and advice.
Thank you in advance.
 
Status
Not open for further replies.
Back
Top