Notifications
Clear all

Creating user defined elements in ABAQUS Explicit

8 Posts
5 Users
0 Likes
1,163 Views
hayakuwu
Posts: 15
Topic starter
(@hayakuwu)
Active Member
Joined: 13 years ago

Hi there,

I am trying to develop a user defined element subroutine in abaqus explicit, VUEL. Im starting with a simple 4 node linear solid element. However still there are some questions regarding VUEL structure. Im wondering if anyone can introduce me a good source that explains how to write a VUEL or perhaps an example.

thanks

Topic Tags
7 Replies
ndaschakladar
Posts: 9
(@ndaschakladar)
Active Member
Joined: 13 years ago

Hi Fred,

I had to write a 3-D hexahedral VUEL as part of my Phd and apart from the ABAQUS documentation the only book I looked at was Finite Element Method: A Practical Course by Liu G.R. It doesnt give VUEL examples but it explains the formulation of each type of element (1-D, 2-D, 3-D) in detail and helped me a lot.

Topic Tags
7 Replies
hayakuwu
Posts: 15
Topic starter
(@hayakuwu)
Active Member
Joined: 13 years ago

Thanks Stephen,
Well, Im relatively familiar with finite element formulation although I definitely need to take a look at the book you introduced to finish the job,
Right now my biggest concern is that ABAQUS explicit uses logarithmic strain rate. In my previous experiences I used to use the relation e = B*U where e is strain, B is the derivatives of shape function and U is the deformations of the node to derive the strain at any point. This formula in face represent the fact that in my previous experiences strain, let&#8217,s say in x direction, is defined as e = dU/dx. However, I&#8217,m not sure what the approach is when the strain rate is logarithmic. Like I said right now my problems are more related to VUEL structure itself. Also if I want to write a nonlinear VUEL I need to know how to derive geometrically nonlinear logarithmic strain. Therefore, at this stage I really need to find some detailed explanation of VUEL.

Thanks again,

Reply
mahdibes
Posts: 1
(@mahdibes)
New Member
Joined: 13 years ago

StephenKiernan if u come across any other book that give VUEL examples in detail plz do tell......thnx

Reply
ndaschakladar
Posts: 9
(@ndaschakladar)
Active Member
Joined: 13 years ago

@fred1980
Maybe Im wrong but isnt the strain that is calculated from B*u non-linear already? This is the equation I use and the strain that comes out is the true logarithmic strain, not engineering.

@Issac44
If I see one Ill let you know πŸ™‚

Reply
stusapien
Posts: 22
(@stusapien)
Eminent Member
Joined: 13 years ago

I believe with VUEL you must handle the stress and strain definitions yourself in the coding. ABAQUS provides displacements and displacement increments, and expects the right-hand-side vector. You dont have to use logarithmic strains by the way just because ABAQUS uses them for its own formulations - thats part of the freedom of using your own elements! However, if you do use logarithmic or any other advanced strain measure then you have to derive the required B matrix relating those strains to displacements. In that case, I recommend reading books/papers on the topic of large deformations, such as the one recommended by StephenKiernan. If the element youre developing is 1D that process will be much easier than if it is 2D or 3D.

Reply
hayakuwu
Posts: 15
Topic starter
(@hayakuwu)
Active Member
Joined: 13 years ago

StephenKiernan:

well, I know the strain that ABAQUS provides for VUMAT is logarithmic, therefore I thought we should use the same logarithmic strain in VUEL but considering what FEguru said I think he is right and one do not neccessarily have to use logarithmic scale. However, in VUEL example provided by abaqus (the cases of trusses) logarithmic strain are used. Of course for a 3d brick element I think it should be much more difficult than truss elements to use logarithmic strain.

I think if you just simply calculate the strain from B*U it would be linear unless you consider the rotational terms in B matrix.

FEguru :

thanks for your reply, still I hope I can somehow implement logarithmic strain in VUEL. Ill look to see if I can find required formulations.

Reply
Page 1 / 2
Share: