transversal isotropy for beam elements
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.
Hello
subscribe to and seek assistance from the
ABAQUS mailing list
[url] https://groups.yahoo.com/neo/groups/Abaqus/info [/url]
Search the archive of the list before posting in it.
The list does not accept attachments.
Or
[url] http://www.eng-tips.com/threadminder.cfm?pid=799 [/url]
or from the ResearchGate discussion forum.
For advice on mechanics in general, including modeling,
you can consult
[url] http://imechanica.org/ [/url]
If you need subroutines: as an introduction get the file
[url] http://imechanica.org/files/Writing [/url] User Subroutines with ABAQUS.pdf
Good luck
Frank
Hi Frank,
Thanks for the links.
There is also the Simulia Learning Community 3ds.com/slc (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.
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.
-Jorgen
- Need help validating hyperelastic VUMAT2 weeks ago
- question about vumat_nh3 years ago
- PTFE modelling with Abaqus/CAE 6.7-13 years ago
- Help regarding Johnson Cook VUHARD4 years ago
- Thermal residual stresses in orthotropic material using UMAT in ABAQUS5 years ago
Latest Post: BB Model for OptiStruct Nonlinear Static Our newest member: NursingAssignmentsHelp Recent Posts Unread Posts Tags
Forum Icons: Forum contains no unread posts Forum contains unread posts
Topic Icons: Not Replied Replied Active Hot Sticky Unapproved Solved Private Closed