What does this stress output for quad element mean?
Moderator: Moderators
-
- Posts: 109
- Joined: Sat May 05, 2007 12:28 pm
- Location: Houston, TX
What does this stress output for quad element mean?
Hi,
While I was trying to get some sort of stress output in previous versions, I produced a stress output file for one quad element (v1.7.5Beta). My toy model consists of 1 quad element,1 pile,1 PyLiq1 Spring. The recorder command I used looks like:
recorder Element -ele 1 -time -file PressurePDMY.out material 3 stress
where material 3 is the pressure dependent multi yield material.
The weird thing with the output file is this: It has 6 columns, the first being the the time steps. Do you have any idea what the other 5 columns represent?
0.0025 -2.00236e-005 -0.00919024 -2.79097e-008 3.18891e-006 0.0152104
0.005 -0.000448881 -0.0401242 -0.000429611 2.95748e-006 0.0644439
0.0075 -0.00343411 -0.0890313 -0.00341712 2.40216e-006 0.134397
0.01 -0.00951481 -0.146497 -0.00950023 1.98453e-006 0.20637
........................
Thanks,
While I was trying to get some sort of stress output in previous versions, I produced a stress output file for one quad element (v1.7.5Beta). My toy model consists of 1 quad element,1 pile,1 PyLiq1 Spring. The recorder command I used looks like:
recorder Element -ele 1 -time -file PressurePDMY.out material 3 stress
where material 3 is the pressure dependent multi yield material.
The weird thing with the output file is this: It has 6 columns, the first being the the time steps. Do you have any idea what the other 5 columns represent?
0.0025 -2.00236e-005 -0.00919024 -2.79097e-008 3.18891e-006 0.0152104
0.005 -0.000448881 -0.0401242 -0.000429611 2.95748e-006 0.0644439
0.0075 -0.00343411 -0.0890313 -0.00341712 2.40216e-006 0.134397
0.01 -0.00951481 -0.146497 -0.00950023 1.98453e-006 0.20637
........................
Thanks,
Bayram Aygun
Graduate Student, Civil&Env. Eng.
Rice University
Graduate Student, Civil&Env. Eng.
Rice University
Hi Bayram,
Here is what I see in the source code for the stress recorder (found in the getresponse and then getcommittedstress functions)
static Vector temp5(5);
workV6 = currentStress.t2Vector();
temp5[0] = workV6[0];
temp5[1] = workV6[1];
temp5[2] = workV6[2];
temp5[3] = workV6[3];
temp5[4] = scale;
So these would be your 5 columns after the time. currentStress is a t2vector so you just need to determine what each of the components of the t2vector are. Perhaps one of the authors, Frank or Silvia could shed some light on this.
This should give you a start: http://opensees.berkeley.edu/OpenSees/a ... tml#l00132
Robbie
Here is what I see in the source code for the stress recorder (found in the getresponse and then getcommittedstress functions)
static Vector temp5(5);
workV6 = currentStress.t2Vector();
temp5[0] = workV6[0];
temp5[1] = workV6[1];
temp5[2] = workV6[2];
temp5[3] = workV6[3];
temp5[4] = scale;
So these would be your 5 columns after the time. currentStress is a t2vector so you just need to determine what each of the components of the t2vector are. Perhaps one of the authors, Frank or Silvia could shed some light on this.
This should give you a start: http://opensees.berkeley.edu/OpenSees/a ... tml#l00132
Robbie
-
- Posts: 109
- Joined: Sat May 05, 2007 12:28 pm
- Location: Houston, TX
Hi,
I'm trying to understand the examples in <http://cyclic.ucsd.edu/opensees >. Let's consider Example 15: Single quadrilateral element, subjected to sinusoidal base shaking.
I look at the outputs and the post-processing Matlab files but it's really confusing. What are this "%integration point 1 excess pore water pressure" or "%integration point 3 confinement p VS. deviatoric q relation" figures under the "MATLAB Plotting File" link? Where is this integration point 1?
It's probably an easy thing to understand for a geotech. student but I'm puzzled at the moment.
Anyone, out there?
I'm trying to understand the examples in <http://cyclic.ucsd.edu/opensees >. Let's consider Example 15: Single quadrilateral element, subjected to sinusoidal base shaking.
I look at the outputs and the post-processing Matlab files but it's really confusing. What are this "%integration point 1 excess pore water pressure" or "%integration point 3 confinement p VS. deviatoric q relation" figures under the "MATLAB Plotting File" link? Where is this integration point 1?
It's probably an easy thing to understand for a geotech. student but I'm puzzled at the moment.
Anyone, out there?
Bayram Aygun
Graduate Student, Civil&Env. Eng.
Rice University
Graduate Student, Civil&Env. Eng.
Rice University
Re: What does this stress output for quad element mean?
Hope this helps.
From the manual...
For 2D problems, the stress output follows this order: sxx, syy, szz, sxy, hr, where hr is the ratio
between the shear (deviatoric) stress and peak shear strength at the current confinement
(0<= hr<=1.0). The strain output follows this order: exx, eyy, gxy.
For 3D problems, the stress output follows this order: sxx, syy, szz, sxy, syz, szx, hr, and the strain output
follows this order: exx, eyy, ezz, gxy, gyz, gzx.
From the manual...
For 2D problems, the stress output follows this order: sxx, syy, szz, sxy, hr, where hr is the ratio
between the shear (deviatoric) stress and peak shear strength at the current confinement
(0<= hr<=1.0). The strain output follows this order: exx, eyy, gxy.
For 3D problems, the stress output follows this order: sxx, syy, szz, sxy, syz, szx, hr, and the strain output
follows this order: exx, eyy, ezz, gxy, gyz, gzx.
bayram_aygun wrote:Hi,
While I was trying to get some sort of stress output in previous versions, I produced a stress output file for one quad element (v1.7.5Beta). My toy model consists of 1 quad element,1 pile,1 PyLiq1 Spring. The recorder command I used looks like:
recorder Element -ele 1 -time -file PressurePDMY.out material 3 stress
where material 3 is the pressure dependent multi yield material.
The weird thing with the output file is this: It has 6 columns, the first being the the time steps. Do you have any idea what the other 5 columns represent?
0.0025 -2.00236e-005 -0.00919024 -2.79097e-008 3.18891e-006 0.0152104
0.005 -0.000448881 -0.0401242 -0.000429611 2.95748e-006 0.0644439
0.0075 -0.00343411 -0.0890313 -0.00341712 2.40216e-006 0.134397
0.01 -0.00951481 -0.146497 -0.00950023 1.98453e-006 0.20637
........................
Thanks,
-
- Posts: 109
- Joined: Sat May 05, 2007 12:28 pm
- Location: Houston, TX
-
- Posts: 109
- Joined: Sat May 05, 2007 12:28 pm
- Location: Houston, TX
jaugalde,
I checked the (.pdf) manual and the online manual, but I didn't come across a passage where the stress outputs are explained. Can you tell me where you found this information?
One more thing, as far as what you wrote, we can't monitor the vertical stresses in the quad elements, is that it? I need to double-check the pore pressure ratio output I get, that's why I need to get the vertical stresses before after the ground shaking.
Thanks for your time,
I checked the (.pdf) manual and the online manual, but I didn't come across a passage where the stress outputs are explained. Can you tell me where you found this information?
One more thing, as far as what you wrote, we can't monitor the vertical stresses in the quad elements, is that it? I need to double-check the pore pressure ratio output I get, that's why I need to get the vertical stresses before after the ground shaking.
Thanks for your time,
Bayram Aygun
Graduate Student, Civil&Env. Eng.
Rice University
Graduate Student, Civil&Env. Eng.
Rice University
pest control
[url=http://pest-control.yoyo.pl/control-home-lutz-pest.html]Control home lutz pest[/url]
[url=http://pest-control.yoyo.pl/control-hom ... pider.html]Control home pest remedy spider[/url]
[url=http://pest-control.yoyo.pl/control-hou ... rmite.html]Control houston pest termite[/url]
[url=http://pest-control.yoyo.pl/control-humphreys-pest.html]Control humphreys pest[/url]
[url=http://pest-control.yoyo.pl/control-it- ... rself.html]Control it pest termite yourself[/url]
[url=http://pest-control.yoyo.pl/control-hom ... pider.html]Control home pest remedy spider[/url]
[url=http://pest-control.yoyo.pl/control-hou ... rmite.html]Control houston pest termite[/url]
[url=http://pest-control.yoyo.pl/control-humphreys-pest.html]Control humphreys pest[/url]
[url=http://pest-control.yoyo.pl/control-it- ... rself.html]Control it pest termite yourself[/url]
Re:
bayram_aygun wrote:
> jaugalde,
>
> I checked the (.pdf) manual and the online manual, but I didn't come across
> a passage where the stress outputs are explained. Can you tell me where you
> found this information?
>
> One more thing, as far as what you wrote, we can't monitor the vertical
> stresses in the quad elements, is that it? I need to double-check the pore
> pressure ratio output I get, that's why I need to get the vertical stresses
> before after the ground shaking.
>
> Thanks for your time,
This may be outdated, but I was also wondering where you found the stress output information.
> jaugalde,
>
> I checked the (.pdf) manual and the online manual, but I didn't come across
> a passage where the stress outputs are explained. Can you tell me where you
> found this information?
>
> One more thing, as far as what you wrote, we can't monitor the vertical
> stresses in the quad elements, is that it? I need to double-check the pore
> pressure ratio output I get, that's why I need to get the vertical stresses
> before after the ground shaking.
>
> Thanks for your time,
This may be outdated, but I was also wondering where you found the stress output information.
-
- Posts: 1
- Joined: Fri Feb 03, 2012 4:26 pm
Re: What does this stress output for quad element mean?
thanks for the info...
Re: What does this stress output for quad element mean?
I am impressed by the quality of information on this website. There are a lot of good resources here.
I am sure I will visit this place again soon.
I am sure I will visit this place again soon.