Tek-Tips is the largest IT community on the Internet today!

Members share and learn making Tek-Tips Forums the best source of peer-reviewed technical information on the Internet!

  • Congratulations JAE on being selected by the Eng-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
Joined
Oct 31, 2014
Messages
2
Location
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.

Part and Inventory Search

Sponsor

Back
Top