Hi, thanks for the useful video. When we have geometry nonlinearity, we have to add this code or it can help convergency? I have a umat for an isotropic hyperelastic material which works very well under simple loading but it's hard to get converged for complex geometry and high deformation. It can be solved if I add this piece of code into my umat?

Please share the setting of your input file and your solver setting. Statev(2) = 0 alone cant delete elements.

Hey. I tried to implement it in a matlab code. I want to silmulate a uniaxial tensile Test. After the first plastic step i get a little stress in the 22 and 33 component of the stresstensor. this false components are then go into the new trial stress and over several loops add up to a very false result. what do you think i did wrong?

Thanks for your excellent tutorial, I need to assign material strength and modulus to different elements of three materials, if I go through with NOEL, the array will cross the boundary, how can I solve this problem, can you provide some suggestions

thank you for making this beautiful video

Sir, thank you for these exceptional videos on user subroutines in Abaqus. Would you please make a detailed video or a series of videos on writing and implementing VUEL for 2D and 3D uncoupled Cohesive Zone Models? It would be a great opportunity to learn from you about this topic. Please consider this request.

thank you very much for sharing the wonderful vedio! i understand field(1) controls the elastic modulus, but why statev(2) = 0 means the element is deleted can you explain that, thanks again!

The option deepVar is grey in my Abaqus Learning Edition 2023. Is this a limitation?

please help me what is this error

i get error found invalid intiger 0

hi，sir. I changed the geometry in the first model by modelchange. After , in the second model, I import the part by the odb file. The geometry wasn't same as the first one shown in the result. Please tell me if it were nomal.

Thank you very much, it seems very useful, but I don't know how to use this method. What I mean is that I now have a VUMAT subroutine with over 2000 lines, and I want to change it to a UMAT subroutine. What should I do using your method?

How do you set up the number of cpus using command?

thanks for your sharing, sir. I'm using your script but I'm having trouble getting it to run., can you help find where errors are? subroutine vdload ( C Read only (unmodifiable)variables - 1 nblock, ndim, stepTime, totalTime, 2 amplitude, curCoords, velocity, dirCos, jltyp, sname, C Write only (modifiable) variable - 1 value ) C include 'vaba_param.inc' C dimension curCoords(nblock,ndim), velocity(nblock,ndim), 1 dirCos(nblock,ndim,ndim), value(nblock) character*80 sname C real*8 a,b real*8 LA,LB real*8 p !real*8 slabW,v a=0.35 b=0.2 LA=5 LB=2.5 p=1000 slabW=10 v=20 disp=v*stepTime do km = 1, nblock if(curCoords(km,1)<(slabW-LB)/2. .and. curCoords(km,1)>(slabW-LB)/2.-b) then if(curCoords(km,3)<disp+a .and. curCoords(km,3)>disp) then value(km)=p else if(curCoords(km,3)<disp+a+LA .and. curCoords(km,3)>disp+LA) then value(km)=p else value(km)=0.0 end if else if(curCoords(km,1)>(slabW+LB)/2. .and. curCoords(km,1)<(slabW+LB)/2.+b) then if(curCoords(km,3)<disp+a .and. curCoords(km,3)>disp) then value(km)=p else if(curCoords(km,3)<disp+a+LA .and. curCoords(km,3)>disp+LA) then value(km)=p else value(km)=0.0 end if else value(km)=0.0 end if end do return end

@kianaghani I didn't see the concrete tension damage parameter in this plugin

Hi sir! Can I get your email?

Thank you for the video!!! I wish you all good health and happiness, as well as peace, love and harmony!!!

Do it uses stress strain diagram from code or test data?

I replicated your model, but the element were not deleted where the FIELD(1) =0 and SDV2=0, what could be the reason? thanks.

Please explain why did you assign filed(1) = 0 Later field(1) = 1 and finally 3

Suppose, if I wanted to do parameterization, the thickness will dimensions will change, in that case how can I get the corresponding dimensions dynamically in subroutine to compute E. please make a video

Hello bro. Please, how can I contact you. The vdload subroutine file is not compiling. Guess may be an error somewhere.

Is the plugin you have listed still valid for 2023 version? I can't seem to get it working on mine.

@kianaghani1709 How to contact you?

I want to contact you. Can you provide your email address?

The OPEN statement is passing the 72 column limit. Why isn't this causing an error? I thought ABAQUS used fortran 77, with fixed format.

thank you sir for this video

thanks for your sharing, sir. I write the subroutine script, I met a problem that I want introduce an initial stress field with a legendre polynomil. how can I do it. my script like this, can you help find where errors are? subroutine sigini(sigma,coords,ntens,ncrds,noel,npt,layer, $ kspt,lrebar,rebarn) C include 'aba_param.inc' C dimension sigma(ntens),coords(ncrds) character*80 rebarn C integer i double precision x,y, Pnx x=coords(1) y=coords(2) z=coords(3) Pnx=legendreP(3,y/38) C if (n.eq.3.and.abs(x).lt.60.0d0) then sigma(1)=Pnx end if return end C function legendreP(n,y) result(Pn) integer, intent(in) n double precision, intent(in) y double precision Pn if (n==0) then Pn=1.0d0 else if (n==1) then Pn=y else Pn=((2*n-1)*y*legendreP(n-1,y) $ -(n-1) * legendreP(n-2,y)/n end if end function legendreP

thank you sir for this video

Excellent work 👍❤

Great job. Thanks a lot.

Can run the NEW model from the GUI?

Sir, i am using C3D8 elements with 8 nodes, 8 integration points, 5 degrees of freedom at each integration point, 18 variables at each integration point, 13 properties. i politely request if i am using the correct modifications below ************************************************************* *User element, nodes=8, type=U1, properties=13, coordinates=3, var=144 1,2,3 ************************************************* *ELEMENT, TYPE=U1, ELSET=SOLID ************************************************************* also, sir after 1,2.3 do i need to add anything more? politely request for help sir.

Hi Kian, Thank you for your excellent tutorial. I have a question about assigning different Young's moduli to elements, where each element has a unique Young's modulus, in order to simulate a heterogeneous model. I am curious about how we can implement this solution when using VUMAT, considering that Abaqus calls elements block by block. Appreciate your guidance.

Sir is there any MATLAB resource to modify any .inp File for UEL Compatibility. Politely request for help sir.

When i run the input and fortran scripts through terminal, i get the odb as shown in video. However i don't get any results. can you please tell what is wrong. i am using input file provided by you without modifying it.

Hi，sir. I CANT RUN the command, it's warning that"Abaqus Error: The following file(s) could not be located: Job-1.mdl". May you have ever encountered such a problem, how to solve it?

Hi sir i can't run the model with oldjob data. Can you help me please

Can we further add any parts after getting the deformed geometry? or change the dimesions?

Hi, nice video! Do you happen to have umat version instead of vumat?

Sir, I am using combined hardning rule for steel damper subjected to cycle loading. I am trying to validate the results but I am getting high initial stiffness. Is there a way we can use reduced value of youngs modules for initial cycles and higher values for remaining cycles

Thank sir for shearing this video

Why have u used state variables for stress input and not stress values are available at the start of the subroutine???

Great work sir, In damage modeling using UMAT, the following error is coming.. ''Time increment required is less than minimum provided''. My model consist of more than one material. Once the damage reaches maximum value in one element this error is coming and other elements are not even reaching the peak stress.

Hello is the code supplied by E.M. Paneda valid for any 3d geometry (thick bodies) or what are the limitations? Thanks in advance

Thank you very much for your tutorial. Btw, can we run the new model (step 2) directly on the software interface, basically press the button submit in job module?

Hello sir kian ,can you please help me i have just small problem with abaqus

thank you in advaance

I have a question about VUEL-VUMAT, I guess since they may be similar to the case if UEL-VUMAT, I post it here, however I did not see any other video or forum to post my question. I am working on a VUEL-VUMAT for fracture analysis, where calculations are performed in VUEL, and the final results need to be exported to VUMAT for visualization. Here are the key points: 1. I have 80 nodes in total, and I defined a user-defined element with 5 to 13 nodes. After calculating the extension between nodes, I assign a value to each node representing the number of broken bonds. So, after the calculation, I have a matrix 'A' of size (80,1) for one time increment. The first column contains node numbers, and the second column represents the degree of damage. I assume that my nodes and integration points are the same, but corrections are welcome. 2. I created a module to share data files and access the matrix 'A' of size (80, 1) in VUMAT. Using a simple do loop, I assigned the values of matrix 'A' to 'statenew'. However, I couldn't see the results in Abaqus through SDV1. I assign one to one because I am using CPE4R element in the vumat to visualize: do k = 1, 80 statenew(k,1) = A(k) end do 3. I suspect the issue is related to using element CPE4R, which has reduced integration points (only 1 point) located at the middle of the element(?). The values assigned to 'statenew' are now at the corner nodes, which might be causing the problem(?) I am unsure whether I should directly assign 'A' to 'statenew' or create a new node list with different coordinates, forming 4 new nodes to represent the corner nodes around the statenew nodes. 4. At the end of step 2, I am able to see displacement, stress, strain, etc. interpolated using the node values, but I am struggling with how to assign the specific variable values to be shown through SDVs (State Dependent Variables). 5. Regarding the structure of 'statenew' (NBlock, NStatev), I have a question. Is it a simple 2D array? In comparison with UEL-UMAT structure, where we assign a 3D array to 'statev' for visualization in UMAT, I am confused. In the attachment file you will see my node and element distribution. Cross signs are VUEL elements, and in the background is VUMAT elements. As you see the value of SDV1 is 0 everwhere. MANY THANKS IN ADVANCE !

I want to use VDFLUX and VUSDFLD subroutine in one file. how to do this?