[ GiDlist ] visualization of gauss point values

[Aldo Ghisi]

Hello,

I implemented a GiD model with two different types of 2D elements,
triangles (one gauss point) and quadrilaterals (four GP).

Now, I'm trying to visualize the stresses and strains at gauss points in my
example (see attachments), but in GiD 6.1.5b I can see at once only one of
the two different domains (tria or quad). Is it possible to visualize the
results in the gauss points for all the two domains (quadrilaterals _and_
triangles)?

Note that for the quadrilaterals the gauss points values (from the second
row to the fourth) are not realistic in es03b.flavia.res, but this should
not matter.

Thank you in advance and regards.
Aldo

ALDO GHISI
Department of Structural Engineering
Politecnico di Milano
Piazza Leonardo da Vinci 32, 20133 Milano, Italy
Tel.: +39 02 2399 4218
Fax.: +39 02 2399 4220
E_mail: aldo.ghisi at polimi.it
-------------- next part --------------
GiD Post Results File 1.0
GaussPoints "Tria3_1_GP" ElemType Triangle "Tria3_1_GP"
Number of Gauss Points: 1
Natural Coordinates: Internal
End GaussPoints
GaussPoints "Quad4_4_GP" ElemType Quadrilateral "Quad4_4_GP"
Number of Gauss Points: 4
Natural Coordinates: Internal
End GaussPoints
Result "Displacements" "Titolo" 1 Vector OnNodes
ComponentNames "X-DISPL", "Y-DISPL", "Z-DISPL"
Values
1 0.00000E+00 0.00000E+00
2 -6.86845E-04 -5.53398E-04
3 3.92894E-06 1.78447E-05
4 -6.78444E-04 2.13458E-05
5 0.00000E+00 0.00000E+00
6 -1.75166E-03 -9.27109E-04
7 -8.52555E-04 6.17414E-04
8 -1.73433E-03 2.70467E-05
9 -1.75616E-03 1.02228E-03
10 -2.85734E-03 -9.86065E-04
11 -2.90828E-03 2.83479E-05
12 -3.12001E-03 1.16769E-03
End Values
Result "Stress" "Titolo" 1 Matrix OnNodes
ComponentNames "XX-STR", "YY-STR", "XY-STR", "ZZ-STR"
Values
1 -1.8839E+00 -4.4562E+00 -1.9020E+00 -1.5438E+00
2 -7.5555E-01 -2.4206E+00 -8.1245E-01 -1.1560E+00
3 -6.0784E-01 -1.4614E+00 -6.2076E-01 -1.1360E+00
4 -1.6667E-01 7.1489E-01 1.9710E-01 -1.0694E+00
5 3.5565E-01 2.1994E+00 7.6651E-01 -1.1039E+00
6 -3.3196E-01 -1.9817E+00 -4.8138E-01 -9.2480E-01
7 -7.6968E-02 3.3601E+00 8.6123E-01 -7.9987E-01
8 -3.5712E-01 -2.0002E-01 -1.3230E-01 -8.4685E-01
9 -9.0847E-01 1.5384E+00 9.8198E-02 -8.0979E-01
10 -4.4375E-01 -5.7841E-01 -2.5554E-01 -8.7097E-01
11 -1.3169E+00 2.9780E-01 -2.5478E-01 -1.0505E+00
12 -2.7428E+00 1.2572E+00 -3.7138E-01 -1.2572E+00
End Values
Result "Strains" "Titolo" 1 Matrix OnNodes
ComponentNames "XX-EPS", "YY-EPS", "XY-EPS", "ZZ-EPS"
Values
1 7.8579E-07 -1.1068E-04 0.0000E+00 -1.3380E-04
2 2.0284E-06 -5.6583E-05 0.0000E+00 -6.0575E-05
3 5.6010E-07 -3.6426E-05 0.0000E+00 -9.8456E-05
4 -7.5831E-06 1.0640E-05 0.0000E+00 -5.0567E-05
5 -1.7804E-05 6.2092E-05 0.0000E+00 -9.5670E-05
6 3.2399E-06 -3.9940E-05 0.0000E+00 -3.4796E-05
7 -2.0411E-05 6.0200E-05 0.0000E+00 -3.4058E-05
8 -3.7766E-06 -1.2052E-05 0.0000E+00 -3.1415E-05
9 -1.7976E-05 1.6103E-05 0.0000E+00 -3.2713E-05
10 -3.3609E-06 -5.7655E-06 0.0000E+00 -3.1106E-05
11 -1.8967E-05 9.8676E-06 0.0000E+00 -3.7518E-05
12 -4.2346E-05 2.9082E-05 0.0000E+00 -4.4901E-05
End Values
Result "Stress_GP" "Titolo" 1 Matrix OnGaussPoints "Tria3_1_GP"
ComponentNames "XX-STR", "YY-STR", "ZZ-STR", "XY-STR"
Values
1 -1.961E-02 1.468E-02 -1.481E-03 -1.616E+00
2 -1.884E+00 -4.456E+00 -1.902E+00 -1.544E+00
3 7.998E-02 5.736E-02 4.120E-02 -2.484E-01
4 7.309E-01 4.384E+00 1.535E+00 -5.919E-01
7 -3.898E-02 -8.934E-01 -2.331E-01 -8.485E-01
8 -8.485E-01 -2.634E-01 -2.780E-01 -8.934E-01
9 -3.595E-01 -1.004E-01 -1.150E-01 -1.001E+00
10 -2.743E+00 1.257E+00 -3.714E-01 -1.257E+00
End Values
Result "Stress_GP" "Titolo" 1 Matrix OnGaussPoints "Quad4_4_GP"
ComponentNames "XX-STR", "YY-STR", "ZZ-STR", "XY-STR"
Values
5 -3.847E-01 -2.013E+00 -4.795E-01 -1.496E+00
-3.647E-01 -2.113E+00 -4.695E-01 -1.396E+00
-3.747E-01 -2.213E+00 -4.795E-01 -1.596E+00
-3.447E-01 -2.313E+00 -4.595E-01 -1.496E+00
6 -8.004E-01 1.938E+00 2.276E-01 -1.044E+00
-8.104E-01 1.438E+00 2.176E-01 -1.144E+00
-8.204E-01 1.338E+00 2.376E-01 -1.044E+00
-8.304E-01 1.538E+00 2.476E-01 -1.144E+00
End Values
Result "Strains_GP" "Titolo" 1 Matrix OnGaussPoints "Tria3_1_GP"
ComponentNames "XX-EPS", "YY-EPS", "XY-EPS", "ZZ-EPS"
Values
1 -7.858E-07 7.002E-07 0.000E+00 -1.400E-04
2 7.858E-07 -1.107E-04 0.000E+00 -1.338E-04
3 1.680E-06 7.002E-07 0.000E+00 -2.153E-05
4 -3.482E-05 1.235E-04 0.000E+00 -5.130E-05
7 3.466E-06 -1.179E-05 0.000E+00 -3.030E-05
8 -1.019E-05 2.602E-07 0.000E+00 -3.191E-05
9 -4.367E-06 2.602E-07 0.000E+00 -3.574E-05
10 -4.235E-05 2.908E-05 0.000E+00 -4.490E-05
End Values
Result "Strains_GP" "Titolo" 1 Matrix OnGaussPoints "Quad4_4_GP"
ComponentNames "XX-EPS", "YY-EPS", "XY-EPS", "ZZ-EPS"
Values
5 2.276E-06 -3.680E-05 0.000E+00 -7.183E-05
2.376E-06 -3.780E-05 0.000E+00 -7.283E-05
2.476E-06 -3.880E-05 0.000E+00 -7.383E-05
2.176E-06 -3.980E-05 0.000E+00 -7.483E-05
6 -2.467E-05 4.106E-05 0.000E+00 -5.013E-05
-2.367E-05 4.206E-05 0.000E+00 -5.113E-05
-2.567E-05 4.306E-05 0.000E+00 -5.213E-05
-2.667E-05 4.406E-05 0.000E+00 -5.313E-05
End Values
-------------- next part --------------
# FEAP Es03b

# numnp,numel,nummat,ndm,ndf,nen,nen1
# 12 10 3 2 2 4 9


Mesh "Only_coordinates" Dimension 2 ElemType Triangle Nnode 3

Coordinates
1 0.0000 0.0000
2 0.0000 5.0000
3 5.0000 0.0000
4 5.0000 5.0000
5 10.0000 0.0000
6 0.0000 10.0000
7 10.0000 5.0000
8 5.0000 10.0000
9 10.0000 10.0000
10 0.0000 15.0000
11 5.0000 15.0000
12 10.0000 15.0000
End coordinates

Elements
# no elements entered, they will be in the successive meshes
End elements


Mesh "Tria3_1_GP" Dimension 2 ElemType Triangle Nnode 3

Coordinates
# no coordinates entered, they are already in the first mesh
End coordinates

Elements
1 3 5 4 1
2 2 1 3 1
3 3 4 2 1
4 7 4 5 1
7 6 8 10 3
8 8 11 10 3
9 8 9 11 3
10 9 12 11 3
End elements


Mesh "Quad4_4_GP" Dimension 2 ElemType Quadrilateral Nnode 4

Coordinates
# no coordinates entered, they are already in the first mesh
End coordinates

Elements
5 2 4 8 6 2
6 4 7 9 8 2
End elements

[Pablo Perez del Castillo]

Hi:
I have two questions:

1. It is possible create contact elements node to node (linear volumen) between
shells with 0 gap?.

2. When i create contact elements to nonlinear analisys, i use piramidal
elements, so every node has relation with four nodes, but how GID read
piramidal elements? to write calculation file.

3. I want create a load hidrostatic condition, so it has calculate a distance
from the coordenate system to CG of every plate element, it is possible
implement in GID or i have to do a calculate.exe.

Thanks, and congratulation to everybody with make possible GID.

Note
To me ACIIS sat import works perfect.

[Enrique Escolano]

----- Original Message -----
From: "Pablo Perez del Castillo" pablopdc at terra.es
To: gidlist at gatxan.cimne.upc.es
Sent: Thursday, October 11, 2001 4:43 AM
Subject: [ GiDlist ] contact between shells


Hi:
I have two questions:

1. It is possible create contact elements node to node (linear volumen) between
shells with 0 gap?.

Yes, see Create-Contact-Volume and select two surfaces in the same position.
(Contact-Separated Volume is for two surfaces with a rigid body movement)

Meshing are created prismas (if the surfaces are meshed with triangles) or hexahedral elements (for surfaces with quadrilaterals).

2. When i create contact elements to nonlinear analisys, i use piramidal
elements, so every node has relation with four nodes, but how GID read
piramidal elements? to write calculation file.


Gid cannot read piramidal elements, but this element can be split in two tetrahedral to read it


3. I want create a load hidrostatic condition, so it has calculate a distance
from the coordenate system to CG of every plate element, it is possible
implement in GID or i have to do a calculate.exe.

You can use a #FUNC field in the hidrostatic condition definition (file *.cnd), for example

NUMBER: 1 CONDITION: Hidrostatic
CONDTYPE: over surfaces
CONDMESHTYPE: over nodes
QUESTION: Z_Water:
VALUE: 0
QUESTION: Density:
VALUE: 1.0
QUESTION: Pressure#FUNC#(Cond(2,REAL)*(z-Cond(1,REAL)))
VALUE: Autocalculate
END CONDITION

After meshing, the Pressure field is calculated for each node, using z of the node (or the element center in codition over elements).

Best regards


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://listas.cimne.upc.edu/pipermail/gidlist/attachments/20011011/33361e14/attachment.htm