Close this search box.
Clear all

transversal isotropy for beam elements

4 Posts
3 Users
0 Reactions
Posts: 15
Topic starter
Active Member
Joined: 8 years ago

Hi all,
I am doing some simulations for UD pultruded composite beams in Abaqus. In my case the considered profiles are symmetric and solid. I wrote the following simple VUMAT to specify the material behaviour (1 being the fibre direction):

[FONT=courier new]e11 = props(1)
g12 = props(2)
v12 = props(3)

do k = 1, nblock

C Update stresses (direct & torsion shear)
stressNew(k,1) = stressOld(k,1) + e11 * strainInc(k,1)
stressNew(k,4) = stressOld(k,4) + g12 * strainInc(k,4)
C Define transverse strain increments (direct)
strainInc(k,2) = -v12 * strainInc(k,1)
strainInc(k,3) = strainInc(k,2)

end do[/FONT]

This seems to work fine, but I still have two questions:
1. The subroutine manual states that it is necessary to define the transverse strain increments. However, I notice the results are the same when this step is omitted. Does Abaqus actually use these strain increments, and how?
2. Using a VUMAT, the transverse shear stiffness needs to be defined separately. When K12 = k G12 A, can we use the same shear correction factor k as for isotropic material? In some first checks against analytical calculation of a cantilever beam and simulation of a solid model this seems to hold fairly well, though there exists some error.

Thank you.

3 Replies
Posts: 80
Trusted Member
Joined: 15 years ago


subscribe to and seek assistance from the
ABAQUS mailing list
[url] [/url]

Search the archive of the list before posting in it.

The list does not accept attachments.

[url] [/url]
or from the ResearchGate discussion forum.

For advice on mechanics in general, including modeling,
you can consult
[url] [/url]

If you need subroutines: as an introduction get the file
[url] [/url] User Subroutines with ABAQUS.pdf

Good luck


3 Replies
Posts: 15
Topic starter
Active Member
Joined: 8 years ago

Hi Frank,

Thanks for the links.
There is also the Simulia Learning Community (though its not the most user-friendly interface)

But in the past I have always found the most interesting information on this forum, which is why I post it here.

Posts: 3998
Joined: 5 years ago

The transverse strain increment is only used for certain element types, and only becomes important for certain loading modes. Yes, you can estimate the transverse shear stiffness from the isotropic terms. I have also found that the transverse shear stiffness is often not very important.