UEL subroutine in ABAQUS part four

In this video, the UEL subroutine is coded and runned through ABAQUS.
Link to download the files:
github.com/KianAghani/UEL

Пікірлер: 31

  • @Valera197615
    @Valera1976154 ай бұрын

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

  • @qwm1007
    @qwm10072 жыл бұрын

    Thank you for such videos!

  • @ehsankk1636
    @ehsankk16367 ай бұрын

    Great job. Thanks a lot.

  • @hnrwagner
    @hnrwagner2 жыл бұрын

    Great video as always

  • @kianaghani

    @kianaghani

    2 жыл бұрын

    Thank you Dr. Wagner. Very heartwarming.

  • @maziartour
    @maziartour Жыл бұрын

    Hello. Thank you for your great videos. Do we have to define *PREPRINT, *UELPROPERTY, and *USERELEMENT in the inputfile? ABAQUS can not visualize the inputfile because it cannot define these commands

  • @kianaghani

    @kianaghani

    Жыл бұрын

    Hello. Yes. Watch the other complementary videos to get your answer.

  • @seifzeghichi
    @seifzeghichi2 жыл бұрын

    Thank you for the video, it was really helpful I have one question how can I edit the input file when I have two parts (a bimorph) with constraints between them thenks.

  • @kianaghani

    @kianaghani

    2 жыл бұрын

    Hi. Create the model using CAE, then write the INP. Next, edit the INP file.

  • @seifzeghichi

    @seifzeghichi

    2 жыл бұрын

    @@kianaghani I forgot to mention that I'm using UEL subroutine It's easy to edit the INP file for one part, but when I have a two part body, it becomes difficult to edit the INP file. The problem is in the surface created between the parts. I'm trying to find a model of a two part body, but I can't find any.

  • @kianaghani

    @kianaghani

    2 жыл бұрын

    @@seifzeghichi In that case, it’s better to generate the file using MATLAB.

  • @jackrider4510

    @jackrider4510

    5 ай бұрын

    Hello, did you solve the issue? I found an error if I'm trying to create a surface in UEL. Please suggest me to rectify this issue.

  • @seifzeghichi

    @seifzeghichi

    5 ай бұрын

    @@jackrider4510 I used "partition" feature which creates additional nodes to the part. My structure was not complicated though, just two layers of piezoelectric material, I can give you my WhatsApp if you want. If you need the surface to apply distributed load, there is a subroutine that allows you to do so in UEL.

  • @biswabhanupuhan7036
    @biswabhanupuhan70362 жыл бұрын

    Hello, thank you very much for these videos which are very very helpful to amateurs like me. I tried this code with 16 elements, then 64 elements, and 100 elements, it worked every time with minor modifications. But when I tried with 400 elements like 20*20 in x and y direction, it is running but diverging on the first step!! I am not sure why this is happening, not even 1 increment, kindly please help, it will help a lot of us.

  • @biswabhanupuhan7036

    @biswabhanupuhan7036

    2 жыл бұрын

    Okay so in the code the rhs has to be multiplied with determinant of jacobian

  • @altaflone3356
    @altaflone33567 ай бұрын

    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.

  • @kianaghani

    @kianaghani

    7 ай бұрын

    I think you are making a mistake in the process.

  • @kianaghani

    @kianaghani

    7 ай бұрын

    ​@@altaflone3356 At the end of the inp file, there is a line "*element output, elset=decoy, SDV ". With this line the data will be shown.

  • @142ginovarkey3
    @142ginovarkey3 Жыл бұрын

    Can anybody help me with the visualization part, I am not getting it. What is the need for saving the variables to be displayed using statevariable subroutine?

  • @kianaghani

    @kianaghani

    Жыл бұрын

    It is one way to utilize Visualization module's capabilities to show the results. Since, ABAQUS does not automatically show the results for a user element.

  • @mortazapirouz7482
    @mortazapirouz7482 Жыл бұрын

    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 !

  • @kianaghani

    @kianaghani

    11 ай бұрын

    There are some complications in your text. First, your A matrix has ONE column not two. Assuming that it is a typo error, you should write a code to make your integration points compatible with those of the CPE4R element. Statenew is a 2D array, while STATEV is a 1D array. Also, the way I used for visualization is not related to the size of STATEV.

  • @mortazapirouz7482

    @mortazapirouz7482

    11 ай бұрын

    @@kianaghani you are right, the size of A matrix is (80,2), including the nblock and the value for each block. Do you think if I implement the way you did it may work? I did once and no success, I may had a mistake in the input file somewhere. I will try again. Thank you !

  • @kianaghani

    @kianaghani

    11 ай бұрын

    @@mortazapirouz7482 Theoretically, it should work but it will be hard.

  • @mortazapirouz7482

    @mortazapirouz7482

    11 ай бұрын

    @@kianaghani Thank you for the input. I got it work, the thing I was missing is doing a simple stress calculation in VUMAT for dummy elements. This let me to visualize whatever is assigned to statenew. BTW, your videos really helped me to figure out how to get it done. Thank you and keep up the great work !!

  • @kianaghani

    @kianaghani

    11 ай бұрын

    @@mortazapirouz7482 Your welcome. If it is possible, share your subroutines ( or at least a part of it) with me for expanding the knowledge.

  • @compengaero4481
    @compengaero44814 ай бұрын

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

  • @kianaghani

    @kianaghani

    4 ай бұрын

    At the end of the command line type, for example, "cpus=4", you can use gpus too, "gpus=5"

  • @compengaero4481

    @compengaero4481

    4 ай бұрын

    Thanks for your quick reply. ​@@kianaghani

  • @compengaero4481

    @compengaero4481

    4 ай бұрын

    ​@@kianaghanithanks for your quick reply