<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:12pt"><DIV>Dear All, </DIV>
<DIV> </DIV>
<DIV style="FONT-FAMILY: times new roman, new york, times, serif; FONT-SIZE: 12pt">
<DIV style="FONT-FAMILY: times new roman, new york, times, serif; FONT-SIZE: 12pt">
<DIV style="FONT-FAMILY: times new roman, new york, times, serif; FONT-SIZE: 12pt">
<DIV>I'm new in GEOQUS and I'm in a similar problem as Mr. Bassam Saad had. Suggestion from Prof. Kaspr works fine. Now I want to try with S. Ellis tricks but I do not access into the link you mentioned below. </DIV>
<DIV> </DIV>
<DIV>Pls. send me the paper of S. Ellis.</DIV>
<DIV> </DIV>
<DIV>To recall your memory, I pasted below the response to Basssam Saad. If you have any other suggestion pls. let me know.</DIV>
<DIV> </DIV>
<DIV>Look forward to hear from you and best regards.</DIV>
<DIV> </DIV>
<DIV>Anup</DIV>
<DIV>+1 514 935 5431 </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV><STRONG><U>Prof.'s reply to Bassam Saad request</U></STRONG></DIV>
<DIV>Hi,<BR><BR>I assume you have allready tried to do a *geosatic calculation and this still results in the 10 cm displacement. I would suggest that you try also to constrain all nodes to have zero displacement in the *geostatic step and then relax all constrains but the one you have to leave in your model in a subsequent step. If that still doesn't help try to calculate the stress by abaqus itself an import them in your calculations as initial conditions. This approach was suggested by S. Ellis and you may have to take a look at her poster, which is available at<BR><A href="http://www.geophysik.ruhr-uni-bochum.de/geoqus2007" rel=nofollow target=_blank><FONT color=#810081>http://www.geophysik.ruhr-uni-bochum.de/geoqus2007</FONT></A><BR><BR>Best regards<BR><BR>Kasper<BR><BR>PS: Please subscribe to the geoqus mailing list. Otherwise I have to<BR>mannually mark all of your e-mails as accepted. Instructions can be<BR>found at <A
href="http://www.ruhr-uni-bochum.de/geoqus" rel=nofollow target=_blank><FONT color=#810081>http://www.ruhr-uni-bochum.de/geoqus</FONT></A><BR><BR>+-----------------------------------------------------------------------+<BR>|<I> Dr. Kasper D. Fischer <<A href="http://lists.ruhr-uni-bochum.de/mailman/listinfo/geoqus" rel=nofollow target=_blank><FONT color=#810081>kasper.fischer at ruhr-uni-bochum.de</FONT></A>> |<BR></I>|<I> Ruhr-Universität Bochum |<BR></I>|<I> Institut für Geologie, Mineralogie und
Geophysik |<BR></I>|<I> NA 3 / 174 Tel: +49 234 3227574 |<BR></I>|<I> D-44780 Bochum Fax: +49 234 3214181 |<BR></I>|<I> Germany <A href="http://www.geophysik.ruhr-uni-bochum.de/" target=_blank>www.geophysik.ruhr-uni-bochum.de</A>
|<BR></I>+-----------------------------------------------------------------------+<BR><BR><BR><BR><A href="http://lists.ruhr-uni-bochum.de/mailman/listinfo/geoqus" rel=nofollow target=_blank><FONT color=#810081>bsaad at po-box.mcgill.ca</FONT></A> schrieb:<BR><BR>><I>Does any one have some remarks on how can I establish equilibrium and <BR></I>><I>get zero displacement out of Geostatic Step. I have avery simple plain <BR></I>><I>strain (2D) model and I applied the following initial stress that <BR></I>><I>balances the gravity load of g=-10 m/s^2 operating in z direction <BR></I>><I>throughout a zone that extends vertically from z=0 to z=100 m<BR></I>><I><BR></I>><I>I am estimating the intial vertical stress as follows. I have the zone <BR></I>><I>density of 2066 kg/m^3 and therefore bulk unit weight is (2066 * g=-10) <BR></I>><I>=-20060 N/m^3<BR></I>><I>the stress at depth of 100 is 100*-20060= -2006000
N/M^2<BR></I>><I><BR></I>><I>Therefore Initial Geostatic Stress is<BR></I>><I><BR></I>><I>*initial conditions, type=STRESS, geostatic<BR></I>><I>Ground, -2006000, 0, 0, 100, 0.38, 0<BR></I>><I>Where 0.38 is the coefficient of earth pressure calculated from <BR></I>><I>K0=v/1-v (v is poission ratio).However I am getting ridicules <BR></I>><I>displacement : in the order of 0.10 m<BR></I>><I><BR></I>><I>Background<BR></I>><I>In fact (in plain strain) the exact geostatic stresses in the body is <BR></I>><I>done through Calculating stress and displacement fields resulted from <BR></I>><I>gravity (assuming elastic media) is done through elasticity relations <BR></I>><I>Inserting initial stresses based on<BR></I>><I>\sigma_v = \gamma*h, \sigma_h = K*\sigma_v, where h and K are distance <BR></I>><I>of the calculated point from the surface and Coef earth pressure, <BR></I>><I>respectively, is just a
very approximate way for the existing stress <BR></I>><I>field and hence displacement field. One should play with these <BR></I>><I>approximate field until it becomes very close to the above exact stress <BR></I>><I>field (calculated based on FE elasticity). This game may regretfully <BR></I>><I>last hours and hours.<BR></I>><I><BR></I>><I><BR></I></DIV></DIV><BR>
<META content=on http-equiv=x-dns-prefetch-control></DIV></DIV><!-- cg20.c2.mail.re1.yahoo.com compressed/chunked Wed Nov 4 16:29:15 PST 2009 --></div><br>
</body></html>