Hello everyone,
I have a UMAT subroutine that I'd like to apply to ABAQUS beam elements.
The inp file supplied with it uses an S4R shell element to test this subroutine (works perfectly).
I'd like to test it with a B31 beam element, for example.
The original Part was
*Part, name=Part-1
*Node
1, 0., 0., 0.
2, 1., 0., 0.
3, 0., 1., 0.
4, 1., 1., 0.
*Element, type=S4R
1, 1, 2, 4, 3
*Nset, nset=Set-1, generate
1, 4, 1
*Elset, elset=Set-1
1,
** Section: Section-UMAT
*Shell Section, elset=Set-1, material=TR1180_HAH20_U2
1., 3
*Transverse Shear Stiffness
62656.6, 62656.6, 0.
*End Part
I then suggest replacing this part with the following:
*Part, name=Part-1
*Node
1, 0., 0., 0.
2, 1., 0., 0.
*Element, type=B31
1, 1, 2
*Nset, nset=Set-1, internal, generate
1, 2, 1
*Elset, elset=Set-1, internal, generate
1,
** Section: Section-UMAT Profile: Profile-1
*Beam Section, elset=Set-1, material=TR1180_HAH20_U2, poisson = 0.29, temperature=GRADIENTS, section=RECT
0.1, 0.1
0.,0.,-1.
*Transverse Shear Stiffness
630.67,630.67
*End Part
From what I understand from the abaqus guide, I need to specify the transverse shear stiffness, K_α3 = kGA
Using the material's E and ν values and k=0.85 for a 0.1x0.1 rectangular section, I find K_α3 = 630.7, hence the values I put.
So do you see any other points to modify? (in addition of course to the fact that there are fewer nodes and therefore the boundary conditions will have to be corrected etc).
I saw that Deformation Gradient is unavailable for beams, but I don't see DFGRD being called up in my UMAT.
For reasons of shared licenses, I can not run calculations very frequently, so I would like to be sure that my inp is correct until I can use a license
Thank you very much
