FEAP User Forum
FEAP => General questions => Topic started by: Gianluca Gabellini on March 08, 2013, 04:46:30 AM
-
Hello,
I'm using the mixed element Q1P0 and I am looking for a way to print the stresses computed at the Gauss points instead of the average value for the element. Is it possible to do that with standard commands? I'm using FEAP 8.2. Thank you.
Kind regards,
Gianluca
-
Which Q1P0 element? 2D/3D small or finite deformation?
Some elements do print all the results. For example, sld2d2.f (the small deformation mixed element does). All elements compute the stresses, however some then average to print a single point. Since the element you are using does not print all it will be necessary for you to modify the code yourself.
To do this:
Go to the directory where the element is located (e.g. /elements/solid3d contains the 3-d elements)
Edit the element routine (sld3d2.f for small fld3d2.f for finite - for 3d)
Locate the region where 'isw = 4' (may be a go to or an isw.eq.4 transfer)
Follow throught the computation of the stress (stresses are named sig... (maybe sigl).
Print values of the stress for eacy quadrature point (inside a do l = 1,lint loop).
You may need to add a format or you may be able to use the existing one.
-
Thank you for your kind reply.
I'm working with 3D elements and finite deformation so I will modify fld3d2.f following your instructions. Thank you again.
Kind regards,
Gianluca