dda-cvs Mailing List for Discontinuous Deformation Analysis (Page 4)
Status: Beta
Brought to you by:
doolin
Menu
▾
▴
dda-cvs — cvs mailing list for dda project
You can subscribe to this list here.
| 2003 |
Jan
|
Feb
|
Mar
|
Apr
|
May
(7) |
Jun
(6) |
Jul
(15) |
Aug
|
Sep
|
Oct
|
Nov
|
Dec
(5) |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 2004 |
Jan
(10) |
Feb
(3) |
Mar
|
Apr
|
May
|
Jun
|
Jul
|
Aug
(7) |
Sep
|
Oct
|
Nov
|
Dec
(5) |
| 2006 |
Jan
|
Feb
|
Mar
(88) |
Apr
(32) |
May
(9) |
Jun
(63) |
Jul
(55) |
Aug
(7) |
Sep
(1) |
Oct
(3) |
Nov
(10) |
Dec
(3) |
| 2007 |
Jan
(2) |
Feb
(2) |
Mar
|
Apr
|
May
|
Jun
|
Jul
|
Aug
|
Sep
|
Oct
|
Nov
(8) |
Dec
|
| 2008 |
Jan
|
Feb
|
Mar
|
Apr
(12) |
May
|
Jun
|
Jul
|
Aug
|
Sep
|
Oct
|
Nov
|
Dec
|
| 2010 |
Jan
|
Feb
(2) |
Mar
|
Apr
|
May
|
Jun
|
Jul
|
Aug
|
Sep
|
Oct
|
Nov
|
Dec
|
|
From: David M. D. <do...@us...> - 2006-07-08 16:26:47
|
Update of /cvsroot/dda/htdocs/webref In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv24777/webref Added Files: README Log Message: Added Makefile to with development chores. --- NEW FILE: README --- This directory will eventually contain all the reference files. |
|
From: David M. D. <do...@us...> - 2006-07-08 16:11:22
|
Update of /cvsroot/dda/htdocs In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv19742 Added Files: research.php Log Message: Moved research page into php script for easier future modification. --- NEW FILE: research.php --- <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"> <html xmlns="http://www.w3.org/1999/xhtml"> <head> <?php include("header.php"); ?> <!- <link type="text/css" rel="stylesheet" href="dda.css" /> <link rel="SHORTCUT ICON" href="./Images/favicon.ico" /> <meta name="GENERATOR" content="Microsoft FrontPage 6.0" /> <meta name="ProgId" content="FrontPage.Editor.Document" /> <meta http-equiv="Content-Type" content= "text/html; charset=us-ascii" /> --> <title> Research </title> </head> <body> <table border="4" width="100%" id="table1" cellpadding="10"> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> </td> <td width="76%" bgcolor="#CCCCFF"> <p class="c2"> <span class="c1">DDA Research</span> </p> </td> </tr> <tr> <td width="17%" valign="top" bgcolor="#CCCCFF"> <?php include("menu.php"); ?> </td> <td width="78%"> <div id="maincol"> <p class="c9"> This web site has been developed for sharing work of all members for Discontinuous Deformation Analysis software (DDA). Currently available <span class="c8">researches area contains of:</span> </p> </div> <ol> <li> <p class="c9"> Increasing the graphical ability to create the geometry of model easily. </p> </li> <li> <p class="c9"> <span class="c8">Adding the material models that simulate viscoelastic and viscoplastic (creep) and elastoplastic behavior.</span> </p> </li> <li> <p class="c9"> Implementing the reinforcement elements such as fully grouted cable bolt, lining. </p> </li> <li> <p class="c9"> Graphical output in a wide range of industry-standard formats such as such as Excel, Perl, TCL, Matlab, Octave. </p> </li> <li> <p class="c9"> Creating HTML-based Help accessed from the Help menu describes use of the user-interface; and provides access to the full <span class= "nmc">DDA</span> Manual. </p> </li> </ol> <p> </p> </td> </tr> </table> <p class="c10"> </p> </body> </html> |
|
From: David M. D. <do...@us...> - 2006-07-08 16:10:35
|
Update of /cvsroot/dda/htdocs In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv19385 Added Files: download.php header.php link.php members.php menu.php photos.php references.php Log Message: Moved research page into php script for easier future modification. --- NEW FILE: menu.php --- <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> <html> <head> <title>menu</title> </head> <body> <ul> <li> <p class="c4"> <a class="c12" href="index.php">Home</a> </p> </li> <li> <p class="c4"> <a class="c12" href= "research.php">Research</a> </p> </li> <li> <p class="c4"> <a class="c12" href= "Pages/Download_Page.htm">Download</a> </p> </li> <li> <p class="c4"> <a class="c12" href= "Pages/Members_Page.htm">Members</a> </p> </li> <li> <p class="c4"> <a class="c12" href= "Pages/References_Page.htm">References</a> </p> </li> <li> <p class="c4"> <a class="c12" href= "Pages/Links_Page.htm">Links</a> </p> </li> <li> <p class="c4"> <a class="c12" href= "Pages/Photos_Page.htm">Photos</a> </p> </li> </ul> <hr> <address><a href="mailto:"></a></address> <!-- Created: Sat Jul 08 08:20:07 Pacific Daylight Time 2006 --> <!-- hhmts start --> Last modified: Sat Jul 08 08:20:09 Pacific Daylight Time 2006 <!-- hhmts end --> </body> </html> --- NEW FILE: header.php --- <link type="text/css" rel="stylesheet" href="dda.css" /> <link rel="SHORTCUT ICON" href="./Images/favicon.ico" /> <meta name="GENERATOR" content="Microsoft FrontPage 6.0" /> <meta name="ProgId" content="FrontPage.Editor.Document" /> <meta name="author" content="Roozbeh Grayeli"> <meta name="author" content="David M. Doolin"> <meta http-equiv="Content-Type" content= "text/html; charset=us-ascii" /> --- NEW FILE: photos.php --- <html> <head> <meta name="GENERATOR" content="Microsoft FrontPage 6.0"> <meta name="ProgId" content="FrontPage.Editor.Document"> <meta http-equiv="Content-Type" content="text/html; charset=windows-1252"> <title>Photos</title> </head> <body> <font FACE="Arial" SIZE="4"><i><b> <table border="4" width="100%" id="table1" cellpadding="10"> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> </td> <td width="78%" bgcolor="#CACAFF"><font FACE="Times New Roman"> <p align="center"><font size="7">Photos</font></td> </tr> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> <b> <font FACE="Arial" SIZE="4"> <ul> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="../index.htm"><font color="#000000" size="4">Home page</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="Research_Page.htm"><font color="#000000" size="4">Research</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font color="#0000FF" FACE="Times New Roman"> <a href="Download_Page.htm"><font color="#000000" size="4">Download</font></a><font color="#000000" size="4">s</font></font></li> <li> <p align="justify" style="line-height: 150%"> <font face="Times New Roman"> <a href="Members_Page.htm"><font color="#000000" size="4">Members</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="References_Page.htm"><font color="#000000" size="4">References</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Links_Page.htm"><font color="#000000" size="4">Links</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Photos_Page.htm"><font color="#000000" size="4">Photos</font></a></font></li> </ul> </font></b> <p> </td> <td width="78%"><font FACE="Times New Roman" SIZE="3"> <font FACE="Arial" SIZE="4"><i><b> <p align="center"><font size="6">Under Construction</font></p> </b></i></font> <p ALIGN="justify"> </p></font> <p> </td> </tr> </table> </b></i></font> </body> </html> --- NEW FILE: members.php --- <html> <head> <meta http-equiv="Content-Language" content="en-us"> <meta name="GENERATOR" content="Microsoft FrontPage 6.0"> <meta name="ProgId" content="FrontPage.Editor.Document"> <meta http-equiv="Content-Type" content="text/html; charset=windows-1252"> <title>Members</title> </head> <body> <font FACE="Arial" SIZE="4"><i><b> <table border="4" width="100%" id="table1" cellpadding="10"> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> </td> <td width="78%" bgcolor="#CACAFF"><font FACE="Times New Roman"> <p align="center"><font size="7">Members</font></td> </tr> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> <b> <font FACE="Arial" SIZE="4"> <ul> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="../index.htm"><font color="#000000" size="4">Home page</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="Research_Page.htm"><font color="#000000" size="4">Research</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font color="#0000FF" FACE="Times New Roman"> <a href="Download_Page.htm"><font color="#000000" size="4">Download</font></a><font color="#000000" size="4">s</font></font></li> <li> <p align="justify" style="line-height: 150%"> <font face="Times New Roman"> <a href="Members_Page.htm"><font color="#000000" size="4">Members</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="References_Page.htm"><font color="#000000" size="4">References</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Links_Page.htm"><font color="#000000" size="4">Links</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Photos_Page.htm"><font color="#000000" size="4">Photos</font></a></font></li> </ul> </font></b> <p> </td> <td width="78%"><font FACE="Times New Roman" SIZE="3"> <font FACE="Arial" SIZE="4"><i><b> <p align="center"><font size="6">Under Construction</font></p> </b></i></font> <p ALIGN="justify"> </p></font> <p> </td> </tr> </table> </b></i></font> </body> </html> --- NEW FILE: references.php --- <html> <head> <meta name="GENERATOR" content="Microsoft FrontPage 6.0"> <meta name="ProgId" content="FrontPage.Editor.Document"> <meta http-equiv="Content-Type" content="text/html; charset=windows-1256"> <title>References</title> </head> <body> <i><b> <table border="4" width="100%" id="table2" cellpadding="10"> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> </td> <font FACE="Arial" SIZE="4"> <td width="78%" bgcolor="#CACAFF"><font FACE="Times New Roman"> <p align="center"><font size="7">DDA's Reference list</font></td> </tr> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> <b> <font FACE="Arial" SIZE="4"> <ul> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="../index.htm"><font color="#000000" size="4">Home page</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="Research_Page.htm"><font color="#000000" size="4">Research</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font color="#0000FF" FACE="Times New Roman"> <a href="Download_Page.htm"><font color="#000000" size="4">Download</font></a><font color="#000000" size="4">s</font></font></li> <li> <p align="justify" style="line-height: 150%"> <font face="Times New Roman"> <a href="Members_Page.htm"><font color="#000000" size="4">Members</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="References_Page.htm"><font color="#000000" size="4">References</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Links_Page.htm"><font color="#000000" size="4">Links</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Photos_Page.htm"><font color="#000000" size="4">Photos</font></a></font></li> </ul> </font></b> <p> </td> <td width="78%"><font FACE="Times New Roman"> N. Barton. <B>The Shear Strength of Rock and Rock Joints.</B> 1976. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 13:255-279. <P> N. Bi'cani'c and D. Ponniah and J. Robinson. 2001. Chapter in <I>Computational Modelling of Masonry, Brickwork and Blockwork Structures.</I> Saxe-Coburg Publications. <P> Y. Cai. <B>Numerical simulation od dynamic process of the Tangshan earthquake by a new method --- LDDA.</B> 2000. <I>Pure. appl. geophys.,</I> 157:2083-2104. <P> C. S. Chang. <B>Nonlinear dynamic discontinuous deformation analysis with finite element meshed block system.</B> 1994. Ph. D. thesis, University of California, Berkeley. <P> G. Chen and Y. Ohnishi. <B>Slope Stability Analysis using Discontinuous Deformation Analysis Method.</B> 1999. in <I>Rock Mechanics for Industry,</I> pages 535-541. Proceedings of the 37th U. S. Rock Mechanics Symposium. <P> W.-S. Chen and X.-G. Gu and X.-R. Ge. <B>A new discrete element method and its application in engineering.</B> 1994. in <I>Computer Methods and Advances in Geomechanics,</I> pages 883-887. Proceedings of the Eighth International Conference on Computer Methods and Advances in Geomechanics. <P> S. Chen and J. C. Chern and C. Y. Koo. <B>Performance prediction of tunnel excavation in clean cobble-gravel deposits by DDA method.</B> 1996. in <I>Rock Mechanics: Tools and Techniques---Proceedings of the 2nd North American Rock Mechanics Symposium,</I> pages 2017. <P> Y. M. Cheng and Y. H. Zhang. <B>Rigid body rotation and block internal discretization in DDA analysis.</B> 2000. <I>International Journal for Numerical and Analytical Methods in Geomechanics,</I> 24:567-578. <P> D. M. Doolin and N. Sitar. <B>Time integration in discontinuous deformation analysis DDA.</B> 2004. <I>ASCE Journal of Engineering Mechanics,</I> 130:249-258. <P> D. M. Doolin and N. Sitar. <B>Unified framework for discontinuous deformation analysis and the distinct element method.</B> 2003. <I>IJNME.</I> <P> D. M. Doolin. <B>Preliminary results unifying discontinuous deformation analysis (DDA) and the distinct element method (DEM).</B> 2003. in <I>Development and Application ,</I> pages 119-128. Balkema, Rotterdam. <P> D. M. Doolin. <B>Mathmetical structure and numerical accuracy of discontinuous deformation analysis.</B> 2002. Ph. D. thesis, University of California, Berkeley. <P> D. M. Doolin and N. Sitar. <B>Displacement accuracy of discontinuous deformation analysis method applied to sliding block.</B> 2002. <I>ASCE Journal of Engineering Mechanics,</I> 128:1158-1168. <P> D. M. Doolin and N. Sitar. <B>Accuracy of the DDA method for frictionless impact.</B> 2002. in <I>Mining and Tunnelling Innovation and Opportunity,</I> pages 1289-1296. North American Rock Mechanics Symposium--Tunneling Association of Canada. Balkema, Rotterdam, Toronto, Canada. <P> D. M. Doolin and N. Sitar. <B>DDAML --- Discontinuous deformation analysis markup language.</B> 2001. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 38:467-474. <P> D. M. Doolin and N. Sitar. <B>Accuracy of the DDA method with respect to a single sliding block.</B> 2001. in <I>Rock Mechanics in the National Interest, proceedings of the 38th U.S. Rock Mechanics Symposium,</I> pages 1429-1435. American Rock Mechanics Association. Balkema, Rotterdam, Washington D.C.. <P align="justify"> <a name="Moosavi"></a>M. Moosavi, R. Grayeli<font FACE="Times New Roman" SIZE="3">. <b>A model for cable bolt-rock mass interaction: Integration with discontinuous deformation analysis (DDA) algorithm</b>. 2006. <i>Int. J. Rock Mech. & Min. Sci.</i>, 43:661-670. </font> <P> V. K. Garga and B-L. Wang. <B>A Numerical method for modelling large displacements of jointed rocks. II. Modeling of rock bolts and groundwater and applications.</B> 1993. <I>Canadian Geotechnical Journal,</I> 30:109-123. <P> R. Hamajima and T. Suzuki and Y. Nakai and S. Han. <B>Analysis for discontinuous medium considering elemental deformation.</B> 1994. in <I>Computer Methods and Advances in Geomechanics,</I> pages 821-826. Proceedings of the Eighth International Conference on Computer Methods and Advances in Geomechanics. <P> Y. H. Hatzor and Y. Zaslavsky and A. A. Arzi and A. Shapira. <B>Fully dynamic stability analysis of jointed rock slopes using the DDA method: king Herod's palace, Masada.</B> 2004. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 41:813-832. <P> Y. H. Hatzor and A. Feintuch. <B>The validity of dynamic block displacement prediction using DDA.</B> 2001. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 38:599-606. <P> Y. H. Hatzor. <B>Dynamic rock slope stability analysis at Masada National Monument using Block Theory and DDA.</B> 1999. in <I>Proceedings of 37th U.S. Rock Mechanics Symposium,</I> pages 63-70. Balkema, Rotterdam. <P> Y. H. Hatzor and R. Benary. <B>The Stability of a Laminated Voussoir Beam: Back Analysis of a Historic Roof Collpse Using DDA.</B> 1998. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 35:165-181. <P> L. B. Hilbert, Jr. and W. Yi and N. G. W. Cook and Y. Cai and G.-P. Liang. <B>A new discontinuous finite element method for interaction of many deformable bodies in geomechanics.</B> 1994. in <I>Computer Methods and Advances in Geomechanics,</I> pages 831-836. Proceedings of the Eighth International Conference on Computer Methods and Advances in Geomechanics. <P> L. B. Hilbert. <B>Substructuring and Domain decomposition Techniques for Analysis of Discontinuous Media.</B> 1995. Ph. D. thesis, University of California, Berkeley. <P> S. M. Hsiung and G.-H. Shi. <B>Simulation of earthquake effects on underground excavations using Discontinuous Deformation Analysis (DDA).</B> 2001. in <I>Rock Mechanics in the National Interest, proceedings of the 38th U.S. Rock Mechanics Symposium,</I> pages 1413-1420. American Rock Mechanics Association. Balkema, Rotterdam, Washington D.C.. <P> S. M. Hsiung. <B>Discontinuous Deformation Analysis (DDA) with nth order polynomial displacement functions.</B> 2001. in <I>Rock Mechanics in the National Interest, proceedings of the 38th U.S. Rock Mechanics Symposium,</I> pages 1437-1444. American Rock Mechanics Association. Balkema, Rotterdam, Washington D.C.. <P> T.-K. Huang. <B>Mechanical behavior of interconnected concrete-block retaining wall.</B> 1997. <I>ASCE Journal of Geotechnical and Geoenvironmental Engineering,</I> 123:197-203. <P> Q. H. Jiang and M. R. Yeung. <B>A Model of Point-to-Face Contact for Three-Dimensional Discontinuous Deformation Analysis.</B> 2003. <I>Rock Mechanics and Rock Engineering.</I> <P> L. Jing and Y. Ma and Z. Fang. <B>Modeling of fluid flow and solid deformation for fractured rocks with discontinuous deformation analysis (DDA) method.</B> 2001. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 38:343-355. <P> L. Jing. <B>Block System construction for three-dimensional discrete element models of fractured rocks.</B> 2000. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 37:645-659. <P> L. Jing. <B>Formulation of discontinuous deformation analysis (DDA)---an implicit discrete element model for block systems.</B> 1998. <I>Engineering Geology,</I> 49:371-381. <P> D. A. Horner and J. F. Peters and A. Carillo. <B>Large Scale Discrete Element Modeling of Vehicle-Soil Interaction.</B> 2001. <I>ASCE Journal of Engineering Mechanics,</I> pages 1027-1032. <P> M. Jurasius. <B>Rock slope kinematics of the Mission Peak landslide, Fremont, California.</B> 2002. Master's thesis, San Jose State University. <P> T.-C. Ke and J. Bray. <B>Modeling of Particulate Media Using Discontinuous Deformation Analysis.</B> 1995. <I>Journal of Engineering Mechanics,</I> 121:1234-1243. <P> T-C. Ke and J-H. Tang. <B>Modeling of solid-fluid interactions using the manifold method.</B> 1996. in <I>Rock Mechanics Tools and Techniques,</I> pages 1815-1822. 2nd North American Rock Mechanics Symposium. <P> T. C. Ke. <B>Simulated testing of two-dimensional heterogeneous and discontinuous rock masses using discontinuous deformation analysis.</B> 1993. Ph. D. thesis, University of California, Berkeley. <P> Y.-I. Kim and B. Amadei and E. Pan. <B>Modeling the effect of water, excavation sequence and rock reinforcement with discontinuous deformation analysis.</B> 1999. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 36:949-970. <P> Y. I. Kim and B. Amadei and E. Pan. <B>Modeling the Effect of geology on uplift in concrete gravity dam foundations.</B> 1999. in <I>Rock Mechanics for Industry,</I> pages 527-534. Proceedings of the 37th U. S. Rock Mechanics Symposium. <P> D. S. Kieffer. <B>Rock Slumping: a compound failure mode of jointed hard rock slopes.</B> 1998. Ph. D. thesis, University of California, Berkeley, California. <P> C. Y. Koo and J. C. Chern. <B>Modification of the DDA Method for Rigid Block Problems.</B> 1998. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 35:683-693. <P> H. K. Law and I. P. Lam. <B>Evaluation of seismic perfomance for tunnel retrofit project.</B> 2003. <I>Journal of Geotechnical and Geoenvironmental Engineering,</I> 129:575-589. <P> C. T. Lin and B. Amadei and S. Sture. <B>Using an Augmented Lagrangian Method and block fracturing in the DDA method.</B> 1994. in <I>Computer Methods and Advances in Geomechanics,</I> pages 837-842. Proceedings of the Eighth International Conference on Computer Methods and Advances in Geomechanics. <P> C. T. Lin and B. Amadei and J. Jung and J. Dwyer. <B>Extensions of Discontinuous Deformation Analysis for Jointed Rock Masses.</B> 1996. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 33:671-694. <P> Y. Ohnishi et al. <B>Modifications of the DDA for Elastoplastic Analysis with Illustrative Generic Problems.</B> 1995. in <I>Proceedings of the 35th U.S. Symposium on Rock Mechanics, Reno, Nevada,</I> pages 45-50. Balkema, Rotterdam. <P> M. M. MacLaughlin and N. Sitar and D. M. Doolin and T. Abbot. <B>Investigation of slope stability kinematics using discontinuous deformation analysis.</B> 2001. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 38:753-762. <P> M. M. MacLaughlin. <B>Discontinuous Deformation Analysis of the Kinematics of Landslides.</B> 1997. Ph. D. thesis, University of California, Berkeley. <P> M. M. MacLaughlin and N. Sitar. <B>Kinematics of sliding on a hinged failure surface.</B> 1996. in <I>Rock Mechanics Tools and Techniques,</I> pages 2025-2030. 2nd North American Rock Mechanics Symposium. <P> M. M. MacLaughlin and N. Sitar and T. S. Dobson. <B>Slope Stability Evaluation Using Discontinuous Deformation Analysis.</B> 1995. in <I>Proceedings of the 31st Symposium on Engineering Geology and Geotechnical Engineering,</I> pages 144-156. Logan, Utah. <P> M. M. MacLaughlin and N. Sitar. <B>DDA for Windows: Discontinuous Deformation Analysis for the Windows PC Environment.</B> 1995. Geotechnical Engineering Report No. UCB$/$GT$/$95-04, University of California, Berkeley, California. (urlhttp://www.ce.berkeley.edu/Programs/Geotech/DDA) <P> A. T. McBride and F. Scheele. <B>Validation of Discontinuous Deformation Analysis using a Physical Model.</B> 2001. in <I>Proceedings of Structural Engineering and Mechanics Conference,</I> pages 719-726. Cape Town, South Africa. <P> A. T. McBride. <B>An investigation of Discontinuous De- formation Analysis.</B> 2000. Master's thesis, University of Cape Town. <P> I. H. P. Mamaghani and S. Baba and "O. Aydan and Y. Shimizu. <B>Discrete finite element method for blocky systems.</B> 1994. in <I>Computer Methods and Advances in Geomechanics,</I> pages 843-850. Proceedings of the Eighth International Conference on Computer Methods and Advances in Geomechanics. <P> Y. Ohnishi and K. Yamamukai. <B>Application of DDA in Rockfall Analysis.</B> 1996. in <I>Rock Mechanics: Tools and Techniques---Proceedings of the 2nd North American Rock Mechanics Symposium,</I> pages 2031. <P> G. Rein Soto-Yarritu. <B>Simulaci'on num'erica de material granular mediante el M'etodo de las Deformaciones Discontinuas.</B> 2001. Ph. D. thesis, Universidad Pontificia Comillas, Madrid. <P> T. Sasaki and D. Ishii and Y. Ohnishi and R. Yoshinaka. <B>Stability analysis of jointed rock foundations by DDA.</B> 1995. in <I>Rock Foundation Workshop, ISRM 8th International Conference,</I> pages 337-342. <P> T. Sasaki and Y. Ohnishi and R. Yoshinaka. <B>Discontinuous Deformation Analysis and its Application to Rock Mechanics Problems.</B> 1994. in <I>Proc. JSCE.,</I> pages 493/III-27: 11-20. <P> G.-H. Shi. <B>Three dimensional Discontinuous Deformation Analysis.</B> 2001. in <I>Rock Mechanics in the National Interest, proceedings of the 38th U.S. Rock Mechanics Symposium,</I> pages 1421-1428. American Rock Mechanics Association. Balkema, Rotterdam, Washington D.C.. <P> G.-H. Shi. <B>Discontinuous deformation analysis of block systems.</B> 1996. Technical Report, US Army Engineer Waterways Experiment Station, 3909 Halls Ferry Road, Vicksburg, MS 39180. <P> G-H. Shi. <B>Working Forum on the Manifold Method of Material Analysis.</B> 1995. Technical Report, US Army Corps of Engineers, US Army Corps of Engineers Waterways Experiment Station, Vicksburg, MS. <P> G-H. Shi. <B>Block System Modeling by Discontinuous Deformation Analysis.</B> 1993. Computation Mechanics Publications, Southampton, UK. <P> G-H. Shi. <B>Modeling Rock Blocks by Manifold Method.</B> 1992. in <I>Rock Mechanics,</I> pages 639-648. Proceedings of the 33rd U. S. Symposium. <P> G.-H. Shi. <B>Forward and Backward discontinuous deformation analysis of rock block systems.</B> 1990. in <I>Proceedings of the International Conference on Rock Joints, Loen, Norway,</I> pages 731-743. ISRM. Balkema, Rotterdam. <P> G-H. Shi. <B>Discontinuous Deformation Analysis---A New Model for the Statics and Dynamics of Block Systems.</B> 1988. Ph. D. thesis, University of California, Berkeley. <P> G.-H. Shi and R. E. Goodman. <B>Discontinuous Deformation Analysis---A new method for computing stress, strain and sliding of block systems.</B> 1988. in <I>Key Questions in Rock Mechanics: Proceeding of the 29th U. S. Symposium,</I> pages 381-393. <P> G. H. Shi and R. Goodman and J. Tinucci. <B>The kinematics of block inter-penetrations.</B> 1985. in <I>Research & Engineering Applications in Rock Masses,</I> pages 121-130. 26th U.S. Symposium on Rock Mechanics. <P> G-H. Shi and R. E. Goodman. <B>Two Dimensional Discotinuous Deformation Analysis.</B> 1985. <I>International Journal for Numerical and Analytical Methods in Geomechanics,</I> 9:541-556. <P> G.-H. Shi and R. E. Goodman. <B>Stability Analysis of Infinite Block Systems Using Block Theory.</B> 1985. Unpublished. (Presented July 23, 1985 at the John Bray Colloquium, Imperial College, UK) <P> G.-H. Shi and R. E. Goodman and J. P. Tinucci. <B>Application of block theory to simulated joint trace maps.</B> 1985. in <I>Proceedings of the International symposium on Fundamentals of Rock Joints,</I> pages 367-?. ???. <P> G-H. Shi and R. E. Goodman. <B>Discontinuous deformation analysis.</B> 1984. in <I>Proceedings, 25th U. S. Symposium on Rock Mechanics,</I> pages 269-277. <P> G. H. Shi. <B>Keyblock bolting.</B> 1983. in <I>Proceedings of the International Conference on Rock Bolting</I>. Abisko. <P> G.-H. Shi. <B>Direct-Iterative Solution of Ill-Conditioned Finite Element Matrices.</B> 1982. <I>International Journal for Numerical Methods in Engineering,</I> 18:181-194. <P> G.-H. Shi. <B>On the Least Number of Fixed Points for Infinite Complexes.</B> 1982. <I>Pacific Journal of Mathematics,</I> 103:377-387. <P> G.-H. Shi. <B>A geometric method of stability analysis of rock mass.</B> 1978. Technical Report, Water Conservancy and Hydroelectric Power Scientific Research Institute, Ministry of Water Conservancy and Electric Power, Academica Sinica, The People's Republic of China. <P> N. Sitar and M. MacLaughlin. <B>Kinematics and Discontinuous Deformation Analysis of Landslide Movement (Invited keynote lecture in Session on Mass Movement).</B> 1997. in <I>Proceedings of the 2nd Panamerican Symposium on Landslides,</I> pages 65-73. Rio de Janiero, Brazil. <P> P. A. Thomas and J. D. Bray. <B>Capturing Nonspherical shape of Granular Media with Disk Clusters.</B> 1999. <I>Journal of Geotechnical and Geoenvironmental Engineering,</I> 125:169-178. <P> R.-J. Tsay and Y.-J. Chiou and W.-L. Chuang. <B>Crack Growth Prediction by Manifold Method.</B> 1999. <I>Journal of Engineering Mechanics,</I> 125:884-890. <P> M. Tsesarsky and Y. H. Hatzor and Sitar N. <B>Dynamic displacement of a block on an inclined plane: analytical, experimental, and DDA results.</B> In press. <I>Rock Mechanics and Rock Engineering.</I> <P> B-L. Wang. <B>Deformable block modelling of jointed rocks with the Block-Spring model.</B> 1996. in <I>Rock Mechanics Tools and Techniques,</I> pages 1845-1850. 2nd North American Rock Mechanics Symposium. <P> B-L. Wang and V. K. Garga. <B>A Numerical method for modelling large displacements of jointed rocks. I. Fundamentals.</B> 1993. <I>Canadian Geotechnical Journal,</I> 30:96-108. <P> J. Wartman. <B>Physical model studies of seisically induced deformation in slopes.</B> 1999. Ph. D. thesis, University of California, Berkeley, California. <P> T. Yamamoto and H. L. Konig and H. Sellmejjer and E. V. Hijum. <B>On the response of a poro-elastic bed to water waves.</B> 1978. <I>Journal of fluid mechanics,</I> 87:193-206. <P> M. R. Yeung and Q. H. Jiang and N. Sun. <B>Validation of block theory and three-dimensional discontinuous deformation analysis as wedge stability analysis methods.</B> 2003. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr.,</I> 40:265-275. <P> M. R. Yeung and N. Sun and Q. H. Jiang. <B>Physical modeling and three-dimensional discontinuous deformation analysis of wedge failure.</B> 2003. in <I>Soil and Rock America</I>. <P> M. R. Yeung and S. C. Blair. <B>DDA back analysis of large block test data.</B> 2000. in <I>Proceedings of Pacific Rocks, the 4th North American Rock Mechanics Symposium, Seattle, WA,</I> pages 937-944. Balkema, Rotterdam. <P> M. R. Yeung. <B>Analysis of a Mine Roof using DDA.</B> 1993. <I>Int. J. Rock Mech. Min. Sci. \& Geomech. Abstr..</I> <P> M. R. Yeung and R. E. Goodman. <B>Multi-block stability analysis using the DDA method.</B> 1992. in <I>Proceedings of the Conference on Fractured and Jointed Rock Masses, Lake Tahoe, CA,</I> pages 701-707. Balkema, Rotterdam. <P> M-C. R. Yeung. <B>Application of Shi's Discontinuous Deformation Analysis to the Study of Rock Behavior.</B> 1991. Ph. D. thesis, University of California, Berkeley. <P> </font></td> </tr> </table> </font></b></i> <address> </address> <!-- Created: Tue Aug 31 08:53:22 Pacific Daylight Time 2004 --> <!-- hhmts start --> <p>Last modified: Aprill 18 P 2006 </p> </body> </html> --- NEW FILE: link.php --- <html> <head> <meta name="GENERATOR" content="Microsoft FrontPage 6.0"> <meta name="ProgId" content="FrontPage.Editor.Document"> <meta http-equiv="Content-Type" content="text/html; charset=windows-1252"> <title>Link</title> </head> <body> <font FACE="Arial" SIZE="4"><i><b> <table border="4" width="100%" id="table1" cellpadding="10"> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> </td> <td width="78%" bgcolor="#CACAFF"><font FACE="Times New Roman"> <p align="center"><font size="7">Links</font></td> </tr> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> <b> <font FACE="Arial" SIZE="4"> <ul> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="../index.htm"><font color="#000000" size="4">Home page</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="Research_Page.htm"><font color="#000000" size="4">Research</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font color="#0000FF" FACE="Times New Roman"> <a href="Download_Page.htm"><font color="#000000" size="4">Download</font></a><font color="#000000" size="4">s</font></font></li> <li> <p align="justify" style="line-height: 150%"> <font face="Times New Roman"> <a href="Members_Page.htm"><font color="#000000" size="4">Members</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="References_Page.htm"><font color="#000000" size="4">References</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Links_Page.htm"><font color="#000000" size="4">Links</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Photos_Page.htm"><font color="#000000" size="4">Photos</font></a></font></li> </ul> </font></b> <p> </td> </font> <font SIZE="4"> <td width="78%"><font FACE="Times New Roman" SIZE="3"> <p align="justify"><b>Useful Links:</b></p> <ul> <li> <p style="margin-bottom: 16px"> <span style="font-family: Times New Roman; font-size: 12.0pt; color: black"> <a style="color: blue; text-decoration: underline; text-underline: single" href="http://sourceforge.net/projects/dda/"> SourceForge.net: Discontinuous Deformation Analysis</a> </span> <span style="font-size: 12pt; color: black">The world's largest development and download repository of Open Source code and applications.</span></li> <li> <p style="margin-bottom: 16px"> <span style="font-size: 12.0pt; color: black"> <a target="_blank" style="color: blue; text-decoration: underline; text-underline: single" href="http://www.ce.berkeley.edu/geo/research/DDA/index.html"> DDA for Windows 95/NT (Discontinous Deformation Analysis)</a> Dave Doolin and Nicholas Sitar.</span></li> <li> <p align="justify" style="margin-bottom: 16px"> <span style="font-size: 12.0pt; font-family: Arial; color: black"> <a style="color: blue; text-decoration: underline; text-underline: single" href="http://www.itasca.ca/"> <span style="font-family:"Times New Roman"">Itasca</span><span style="font-family: Times New Roman"> Consulting Canada Inc.</span></a><span style="font-family: Times New Roman"> </span></span><span style="font-size: 12pt; color: black">mine engineering, design, sequencing, rock mechanics, ground control, instrumentation, consulting, modeling.</span></li> <li> <p align="justify" style="margin-bottom: 16px"> <span style="font-family:Arial;color:black"> <a style="color: blue; text-decoration: underline; text-underline: single" href="http://www.rocscience.com/"> <span style="font-family:"Times New Roman"">Rocscience - Slope Stability Software, Geotechnical software</span></a></span><span style="color: black"> </span><span style="font-size: 12pt; color: black">slope stability software, rock engineering, rocscience, rockscience, examine3D, examine2D, examinetab, dips, rocfall, rocdata, unwedge, swedge, slide, ...</span></li> <li> <p align="justify" style="margin-bottom: 16px"> <span style="font-size: 12.0pt; font-family: Times New Roman"> <a href="http://www.rockfield.co.uk">Rockfield Software - World leaders in Advanced Finite Element</a> </span> <span style="font-size: 12pt; color: black">Rockfield Software - Developers & Suppliers of Finite Element Discrete Element Software Systems and Numerical Analysis Consultancy Service.</span></li> <li> <p align="justify" style="margin-bottom: 16px"> <span style="font-family: Times New Roman; font-size: 12.0pt; color: black"> <a style="color: blue; text-decoration: underline; text-underline: single" href="http://www.ansys.com/"> Welcome to ANSYS, Inc. - Corporate Homepage</a> </span> <span style="font-size: 12pt; color: black">Welcome to ANSYS, Inc. - computer-aided engineering technology and engineering design analysis software products and services.<br> <br> </span></li> </ul> <p align="justify"> </p> </font> </font> <font FACE="Arial" SIZE="4"> <p> </td> </tr> </table> </font> </b></i> </body> </html> --- NEW FILE: download.php --- <html> <head> <meta http-equiv="Content-Language" content="en-us"> <meta name="GENERATOR" content="Microsoft FrontPage 6.0"> <meta name="ProgId" content="FrontPage.Editor.Document"> <meta http-equiv="Content-Type" content="text/html; charset=windows-1252"> <title>Download</title> <style> <!-- li.MsoNormal {mso-style-parent:""; margin-bottom:.0001pt; font-size:12.0pt; font-family:"Times New Roman"; margin-left:0in; margin-right:0in; margin-top:0in} span.grame {} span.spelle {} --> </style> </head> <body> <i><b> <table border="4" width="100%" id="table1" cellpadding="10"> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> </td> <font FACE="Arial" SIZE="4"> <td width="78%" bgcolor="#CACAFF"><font FACE="Times New Roman" SIZE="7"> <p align="center"><font color="#0000FF" FACE="Times New Roman"> <a href="Download_Page.htm" style="text-decoration: none"><font color="#000000"> Download</font></a><font color="#000000">s</font></font></td> </tr> <tr> <td width="17%" valign="top" bgcolor="#CACAFF"> <b> <font FACE="Arial" SIZE="4"> <ul> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="../index.htm"><font color="#000000" size="4">Home page</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="Research_Page.htm"><font color="#000000" size="4">Research</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font color="#0000FF" FACE="Times New Roman"> <a href="Download_Page.htm"><font color="#000000" size="4">Download</font></a><font color="#000000" size="4">s</font></font></li> <li> <p align="justify" style="line-height: 150%"> <font face="Times New Roman"> <a href="Members_Page.htm"><font color="#000000" size="4">Members</font></a></font></li> <li> <p align="justify" style="line-height: 150%"> <font FACE="Times New Roman"> <a href="References_Page.htm"><font color="#000000" size="4">References</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Links_Page.htm"><font color="#000000" size="4">Links</font></a></font></li> <li> <p align="justify" style="line-height: 150%"><font face="Times New Roman"> <a href="Photos_Page.htm"><font color="#000000" size="4">Photos</font></a></font></li> </ul> </font></b> <p> </td> <td width="78%"><font FACE="Times New Roman" SIZE="3"> </font> <p align="justify"><b><font face="Times New Roman">DDA's Download link:</font></b><ul> <li> <p align="justify"><font face="Times New Roman">Latest version of DDA created by members of this groups (<a href="http://prdownloads.sourceforge.net/dda/InstallDDA.zip?download">Download</a>)</font></li> </ul> <p align="justify"><b><font face="Times New Roman">Useful papers:</font></b></p> <ul> <li> <p align="justify" style="margin-bottom:12px"> <font face="Times New Roman">Sitar, N., and MacLaughlin, M.M., Kinematics and Discontinuous Deformation Analysis of Landslide Movement, Invited Keynote Lecture, II Panamerican Symposium on Landslides, Rio de Janeiro, Nov. 10-14th, 1997.</font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.ce.berkeley.edu/geo/research/DDA/pubs/"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <span style="font-size: 12.0pt; font-family: Times New Roman">Sitar, N., and MacLaughlin, M.M., Doolin, D.M. and Abbot, T. Investigation of slope stability kinematics using discontinuous deformation analysis. International Journal of Rock Mechanics and Mining Sciences, v. 38, pp. 753-762, 2001.</span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.ce.berkeley.edu/geo/research/DDA/pubs"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <font face="Times New Roman">Doolin, D.M. and Sitar, N., Accuracy of the DDA method with respect to a single sliding block. In: Rock Mechanics in the National Interest, proceedings of the 38th U.S. Rock Mechanics Symposium, Washington D.C., July 5-7, 2001. Balkema, Rotterdam, pp. 1429-1435.</font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.ce.berkeley.edu/geo/research/DDA/pubs"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <font face="Times New Roman">Doolin, D.M. and Sitar, N., DDAML --- discontinuous deformation dnalysis markup language.Int. J. of Rock Mechanics and Mining Sciences, V38, N3, April, 2001, 467-474.</font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.ce.berkeley.edu/geo/research/DDA/pubs"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <font face="Times New Roman">Sitar, N., and MacLaughlin, M.M., Investigation of slope stability kinematics using discontinuous deformation analysis. Invited Keynote Lecture, II Panamerican Symposium on Landslides, Rio de Janeiro, Nov. 10-14th, 1997.</font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.ce.berkeley.edu/geo/research/DDA/pubs"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <font face="Times New Roman">Tsesarsky, M. <span class="grame">and</span> Hatzor, Y. H., 2005. Tunnel roof deflection as a function of joint spacing and friction in blocky rock masses a parametric study using Discontinuous Deformation Analysis (DDA). Tunneling and Underground Space Technology.</font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <span style="font-size: 12.0pt; font-family: Times New Roman"> Tsesarsky M., Hatzor, Y. H. and Sitar N.2005.Dynamic Displacement of a Block on an Inclined Plane: Analytical, Experimental and DDA Results. Technical Note. Rock Mechanics and Rock Engineering. Vol. 38, No. 2, pp. 153-167.</span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <span style="font-size: 12.0pt; font-family: Times New Roman">Hatzor, Y. H., Y. Zaslavsky, A. A. Arzi, and A. Shapira. 2004. Dynamic stability analysis of jointed rock slopes using the DDA method: King Herods Palace, Masada, Israel. International Journal of Rock Mechanics and Mining Sciences Vol.41, pp. 813832.</span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <span style="font-size: 12.0pt; font-family: Times New Roman"> Tsesarsky, M., Hatzor, Y. H. and Sitar, N. 2002. Dynamic block displacement prediction validation of DDA using analytical solutions and shaking table experiments. Stability of Rock Structures: Proceedings of the 5<sup>th</sup> International Conference on Analysis of Discontinuous Deformation (ed. Y. H. Hatzor), Balkema Publishers, Lisse, pp. 195 203.</span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <span style="font-size: 12.0pt; font-family: Times New Roman">Hatzor, Y. H. and A. Feintuch, 2001. The validity of dynamic displacement prediction using DDA. International Journal of Rock Mechanics and Mining Sciences. Vol. 38, No. 4, pp. 599-606.</span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <span style="font-size: 12.0pt; font-family: Times New Roman">Hatzor, Y. H. and R. Benary, 1998. Stability of a laminated Voussoir beam: back analysis of a historic roof collapse using DDA. International Journal of Rock Mechanics and Mining Sciences., Vol. 35, No. 2, pp. 165-181.</span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <font face="Times New Roman"> <span style="font-size: 12pt; font-family: TimesNewRomanPSMT">Kamai, R. and Hatzor, Y.H. 2005. Dynamic back analysis of structural failures in archeological sites to obtain paleo-seismic parameters using DDA. Stability Analyses using Discontinuous Methods: Proceedings of the 7<sup>th</sup> International Conference on Analysis of Discontinuous Deformation (ed. M. MacLaughlin and N. Sitar), Honolulu, Hawaii, </span> <span style="font-size: 12pt">pp. 121-136.</span></font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <span lang="IT" style="font-size: 12pt; font-family: Times New Roman"> Hatzor, Y. H., Tsesarsky, M. And Eimermacher, R. C. 2006. Structural stability of historic underground openings in rocks: two case studies from Israel. In: </span><span class="spelle"> <span style="font-size: 12pt; font-family: Times New Roman">Stavros</span></span><span style="font-size: 12pt; font-family: Times New Roman"> K. Ed. Fracture and Failure of Natural Building Stones (<span class="spelle">ed.S</span>. <span class="grame">K. </span><span class="spelle">Kourkoulis</span><span class="grame">)</span></span><span class="grame"><span lang="IT" style="font-size: 12pt; font-family: Times New Roman">, Kluwer Academic Publishers.</span></span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%"> <font face="Times New Roman"><span style="font-size: 12pt">Tsesarsky, M. <span class="grame">and</span> Hatzor, Y. H., 2005. Tunnel roof deflection as a function of joint spacing and friction in blocky rock masses a parametric study using Discontinuous Deformation Analysis (DDA). Tunneling and Underground Space Technology.</span></font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p style="line-height: 150%; margin-bottom:12px"> <span style="font-size: 12pt; font-family: Times New Roman"> Tsesarsky, M. Hatzor, Y. H., Leviathan, Saltzman, U., Sokolowsky, M. 2005. Structural control on the stability of overhanging, discontinuous rock slopes. 40th U.S. Symposium on Rock Mechanics (USRMS): Rock Mechanics for Energy, Mineral and Infrastructure Development in the Northern Regions. Anchorage, Alaska, June 25-29, 2005. CD-ROM paper ARMA/USRMS 05-771.</span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <font face="Times New Roman"><span style="font-size: 12pt">Tsesarsky, M. and Hatzor, Y. H. 2003. Deformation and kinematics of vertically jointed rock layers in underground openings. In: Development and Application of Discontinuous Modelling for Rock Engineering:</span><span style="font-size: 12pt; color: blue"> </span><span style="font-size: 12pt">ICADD 6, Proceedings of the 6<sup>th</sup> International Conference on Analysis of Discontinuous Deformation (ed. M. Lu), Balkema Publishers, Lisse, 93-101.</span></font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <font face="Times New Roman"><span style="font-size: 12pt">Tsesarsky, M., Hatzor, Y. H. and Sitar, N. 2002. Dynamic block displacement prediction validation of DDA using analytical solutions and shaking table experiments. Stability of Rock Structures: Proceedings of the 5<sup>th</sup> International Conference on Analysis of Discontinuous Deformation (ed. Y. H. Hatzor), Balkema Publishers, Lisse, pp. 195 203.</span></font><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> <li> <p align="justify" style="line-height: 150%; margin-bottom:12px"> <span style="font-size: 12pt; font-family: Times New Roman">Hatzor, Y, H., 1999. The Vousoir beam reaction curve. ICADD-3, Proceedings of the 3rd Intl. Con. On Analysis of Discontinuous Deformation (ed. B. Amadei) ARMA, Alexandria, pp. 117-126.</span><i><b><font face="Times New Roman" SIZE="4"><a href="http://www.bgu.ac.il/geol/hatzor"><img border="0" src="../Images/list.gif" width="20" height="21"></a></font></b></i></li> </ul> <p align="justify"> </p> <p align="justify"> </p> <p align="justify"> </p> <p align="justify"> </p> <p align="justify"> </p> <p align="justify"> </p> <p align="justify"><i><br> </i> </p> <p align="justify"> </p> <p align="justify"> </p> <p align="justify"> </p> <p class="MsoNormal" style="margin-left: .25in"> </p> <p> </p> <p> </p> <p> </p> <p> </td> </tr> </table> </font></b></i> </body> </html> |
|
From: David M. D. <do...@us...> - 2006-07-08 16:10:32
|
Update of /cvsroot/dda/htdocs/Images In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv19385/Images Added Files: favicon.ico Log Message: Moved research page into php script for easier future modification. --- NEW FILE: favicon.ico --- (This appears to be a binary file; contents omitted.) |
|
From: David M. D. <do...@us...> - 2006-07-08 16:09:11
|
Update of /cvsroot/dda/htdocs In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv18972 Added Files: dda.css Log Message: Refactored style information from index.htm into dda.css for global style changes. --- NEW FILE: dda.css --- p.c2 {text-align: center} p.c4 {font-weight: bold; line-height: 150%; text-align: justify} p.c6 {line-height: 150%; text-align: justify} p.c9 {color: #0000FF; font-family: Times New Roman; line-height: 150%; text-align: justify} p.c10 {text-align: justify} p.c11 {font-weight: bold} span.c1 {font-size:36.0pt} span.c8 {color: #000000; font-size: 120%} a.c5 {color: #000000; font-family: Times New Roman; font-size: 120%} a.c3 {color: #000000; font-family: Arial; font-size: 120%} a.c7 {color: #000000; font-size: 120%} a.c12 {color: blue; text-decoration: underline; text-underline: single; font-family: Arial; font-size: 120%; } |
|
From: David M. D. <do...@us...> - 2006-07-08 16:08:08
|
Update of /cvsroot/dda/htdocs In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv18439 Modified Files: index.php Log Message: Rewrote index.php page to allow dynamic generation of dda web pages, with modularity. Index: index.php =================================================================== RCS file: /cvsroot/dda/htdocs/index.php,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** index.php 21 Apr 2006 12:05:56 -0000 1.6 --- index.php 8 Jul 2006 16:08:04 -0000 1.7 *************** *** 1,10 **** ! <?php ! header("Location: http://dda.sourceforge.net/index.htm"); ! exit; ! ?> ! <html> ! <head> ! </head> ! <body> ! If you have not been automatically redirected, you should click <a href="http://dda.sourceforge.net/index.htm/"> http://dda.sourceforge.net/index.htm </body> ! </html> \ No newline at end of file --- 1,89 ---- ! <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" ! "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"> ! <html xmlns="http://www.w3.org/1999/xhtml"> ! <head> ! ! <?php ! include("header.php"); ! ?> ! ! <title> ! Discontinuous Deformation Analysis ! </title> ! ! </head> ! <body> ! <table border="1" width="100%" id="AutoNumber1" cellpadding="10"> ! <tr> ! <td width="17%" bgcolor="#CCCCFF"> ! <img border="0" src="Images/DDA1.gif" width="200" ! alt="DDA Analyses" height="150" /> ! </td> ! ! <td width="76%" bgcolor="#CCCCFF"> ! <p class="c2"> ! <span class="c1">Discontinuous Deformation ! Analysis</span> ! </p> ! </td> ! ! </tr> ! ! <tr> ! ! <td width="17%" valign="top" bgcolor="#CCCCFF"> ! <?php ! include("menu.php"); ! ?> ! </td> ! ! ! <td width="78%"> ! <h2>What is DDA?</h2> ! ! <p class="c10"> ! <?php ! include("whatisdda.txt"); ! ?> ! </p> ! ! ! <h2 class="c11">Using DDA</h2> ! ! <p class="c10"> ! <?php ! include("usingdda.txt"); ! ?> ! </p> ! ! <!-- TODO: Move this to a file, use php include ! to allow other things to get swapped out. ! --> ! ! <p class="c2" dir="ltr"> ! <img border="0" src="Images/cableu.jpg" width="240" ! height="316" /><img border="0" src= ! "Images/cablew.jpg" width="240" height="315" /> ! </p> ! <p class="c2" dir="ltr"> ! (a) ! (b) ! </p> ! <p class="c2"> ! Axial load along the cable bolt that installed at ! roof of<br /> ! underground excavation: (a) without plate and (b) ! with plate (<a class="c12" href= ! "Pages/References_Page.htm#Moosavi">M. Moosavi, R. ! Grayeli</a>)<b>.</b> ! </p> ! <p class="c2" dir="ltr"> ! ! </p> ! </td> ! </tr> ! ! ! </table> ! </body> ! </html> |
|
From: David M. D. <do...@us...> - 2006-07-08 16:07:17
|
Update of /cvsroot/dda/htdocs In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv18340 Added Files: whatisdda.txt usingdda.txt Log Message: Refactored text from DDA index page. --- NEW FILE: whatisdda.txt --- Discontinuous Deformation Analysis (DDA) is a numerical method, similar to the finite element method, for solving stress-displacement problems, but for a jointed rock mass. DDA uses an implicit algorithm for simultaneous solution of the equations of equilibrium by minimizing the total potential energy of the blocky rock mass system. The relation between adjacent blocks is governed by equations of contact interpenetration and accounts for friction. DDA adopts a stepwise approach to solve for the large displacements that accompany discontinuous movements between blocks. The blocks are said to be "simply deformable". Since the method accounts for the inertial forces of the blocks' mass, it can be used to solve the full dynamic problem of block motion. --- NEW FILE: usingdda.txt --- DDA was originally presented as a complete methodology for numerically modeling discontinuous rock mass behavior. Because the behavior of fractured rock masses is largely controlled by the geometry of the fracture system, fracture geometry pre-processing is integral part of an analysis. DDA for Windows currently accepts lists of fractures and reduces the fractures to a block system. All this geometry adds complexity for the user as well as the programmer. The end result is that even in the Windows versions, DDA should be considered at least as complex as the average finite element code, with additional considerations for fracture geometry, plus all the difficulties common to any particle code (i.e., contact mechanics, etc.). |
|
From: David M. D. <do...@us...> - 2006-07-06 15:52:32
|
Update of /cvsroot/dda/htdocs In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv15518 Modified Files: index.htm Log Message: Cleaned up index file home page. Index: index.htm =================================================================== RCS file: /cvsroot/dda/htdocs/index.htm,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** index.htm 1 May 2006 17:44:53 -0000 1.4 --- index.htm 6 Jul 2006 15:52:28 -0000 1.5 *************** *** 1,95 **** ! <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" ! "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"> ! ! <html xmlns="http://www.w3.org/1999/xhtml"> ! ! <head> ! <meta name="generator" content= ! "Microsoft FrontPage 6.0" /> ! <title>Discontinuous Deformation Analysis</title> ! </head> ! ! <body> ! ! ! <table border="1" width="100%" id="AutoNumber1" cellpadding="10"> ! <tr> ! <td width="17%" bgcolor="#CCCCFF"> ! <img border="0" src="Images/DDA1.gif" width="200" height="150"></td> ! <td width="76%" bgcolor="#CCCCFF"> ! <p class="MsoNormal" align="center"><span style="font-size:36.0pt"> ! Discontinuous Deformation Analysis</span></td> ! </tr> ! <tr> ! <td width="17%" valign="top" bgcolor="#CCCCFF"> ! <ul> ! <b> ! <font FACE="Arial" SIZE="4"> ! <li> ! <p align="justify" style="line-height: 150%"> ! <font FACE="Times New Roman"> ! <a href="index.dwt"><font color="#000000" size="4">Home page</font></a></font></li> ! <li> ! <p align="justify" style="line-height: 150%"> ! <font FACE="Times New Roman"> ! <a href="Pages/Research_Page.htm"><font color="#000000" size="4">Research</font></a></font></li> ! <li> ! <p align="justify" style="line-height: 150%"> ! <font color="#0000FF" FACE="Times New Roman"> ! <a href="Pages/Download_Page.htm"><font color="#000000" size="4">Download</font></a><font color="#000000" size="4">s</font></font></li> ! <li> ! <p align="justify" style="line-height: 150%"> ! <font face="Times New Roman"> ! <a href="Pages/Members_Page.htm"><font color="#000000" size="4">Members</font></a></font></li> ! <li> ! <p align="justify" style="line-height: 150%"> ! <font FACE="Times New Roman"> ! <a href="Pages/References_Page.htm"><font color="#000000" size="4">References</font></a></font></li> ! <li> ! <p align="justify" style="line-height: 150%"><font face="Times New Roman"> ! <a href="Pages/Links_Page.htm"><font color="#000000" size="4">Links</font></a></font></li> ! <li> ! <p align="justify" style="line-height: 150%"><font face="Times New Roman"> ! <a href="Pages/Photos_Page.htm"><font color="#000000" size="4">Photos</font></a></font></li> ! </ul> ! </font> ! </b></td> ! <td width="78%"><b>What is DDA?</b><p align="justify">DDA, Discontinuous Deformation ! Analysis, is a numerical method, similar to the finite element method, for ! solving stress-displacement problems, but for a jointed rock mass. DDA uses ! an implicit algorithm for simultaneous solution of the equations of ! equilibrium by minimizing the total potential energy of the blocky rock mass ! system. The relation between adjacent blocks is governed by equations of ! contact interpenetration and accounts for friction. DDA adopts a stepwise ! approach to solve for the large displacements that accompany discontinuous ! movements between blocks. The blocks are said to be "simply deformable". ! Since the method accounts for the inertial forces of the blocks' mass, it ! can be used to solve the full dynamic problem of block motion.</p> ! <p><b>Using DDA</b></p> ! <p align="justify">DDA was originally presented as a complete methodology for numerically ! modeling discontinuous rock mass behavior. Because the behavior of fractured ! rock masses is largely controlled by the geometry of the fracture system, ! fracture geometry pre-processing is integral part of an analysis. In the ! current incarnation, DDA accepts lists of fractures and reduces the ! fractures to a block system. All this geometry adds complexity for the user ! as well as the programmer. The end result is that even in the Windows ! versions, DDA should be considered at least as complex as the average finite ! element code, with additional considerations for fracture geometry, plus all ! the difficulties common to any particle code (i.e., contact mechanics, ! etc.).</p> ! <p> </p> ! <p align="center" dir="ltr"> ! <img border="0" src="Images/cableu.jpg" width="240" height="316"><img border="0" src="Images/cablew.jpg" width="240" height="315"></p> ! <p align="center" dir="ltr"> ! (a) ! (b)</p> ! <p align="center" style="text-align:center">Axial load along the cable bolt ! that installed at roof of<br> ! underground excavation: (a) without plate and (b) with plate (<a style="color: blue; text-decoration: underline; text-underline: single" href="Pages/References_Page.htm#Moosavi">M. ! Moosavi, R. Grayeli</a>)<b>.</b></p> ! <p align="center" dir="ltr"> </td> ! </tr> ! </table> ! ! </body> ! </html> \ No newline at end of file --- 1,156 ---- ! <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" ! "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"> ! <html xmlns="http://www.w3.org/1999/xhtml"> ! <head> ! <meta name="generator" content="Microsoft FrontPage 6.0" /> ! <meta name="author" content="Roozbeh Grayeli"> ! <meta name="author" content="David M. Doolin"> ! <title> ! Discontinuous Deformation Analysis ! </title> ! <style type="text/css"> ! /*<![CDATA[*/ ! a.c12 {color: blue; text-decoration: underline; text-underline: single} ! p.c11 {font-weight: bold} ! p.c10 {text-align: justify} ! p.c9 {color: #0000FF; font-family: Times New Roman; line-height: 150%; text-align: justify} ! span.c8 {color: #000000; font-size: 120%} ! a.c7 {color: #000000; font-size: 120%} ! p.c6 {line-height: 150%; text-align: justify} ! a.c5 {color: #000000; font-family: Times New Roman; font-size: 120%} ! p.c4 {font-weight: bold; line-height: 150%; text-align: justify} ! a.c3 {color: #000000; font-family: Arial; font-size: 120%} ! p.c2 {text-align: center} ! span.c1 {font-size:36.0pt} ! /*]]>*/ ! </style> ! </head> ! <body> ! <table border="1" width="100%" id="AutoNumber1" cellpadding= ! "10"> ! <tr> ! <td width="17%" bgcolor="#CCCCFF"> ! <img border="0" src="Images/DDA1.gif" width="200" ! height="150" /> ! </td> ! <td width="76%" bgcolor="#CCCCFF"> ! <p class="MsoNormal c2"> ! <span class="c1">Discontinuous Deformation ! Analysis</span> ! </p> ! </td> ! </tr> ! <tr> ! <td width="17%" valign="top" bgcolor="#CCCCFF"> ! <ul> ! <li> ! <p class="c4"> ! <a class="c3" href="index.dwt">Home page</a> ! </p> ! </li> ! <li> ! <p class="c6"> ! <a class="c5" href= ! "Pages/Research_Page.htm">Research</a> ! </p> ! </li> ! <li> ! <p class="c9"> ! <a class="c7" href= ! "Pages/Download_Page.htm">Download</a><span class="c8"> ! s</span> ! </p> ! </li> ! <li> ! <p class="c6"> ! <a class="c5" href= ! "Pages/Members_Page.htm">Members</a> ! </p> ! </li> ! <li> ! <p class="c6"> ! <a class="c5" href= ! "Pages/References_Page.htm">References</a> ! </p> ! </li> ! <li> ! <p class="c6"> ! <a class="c5" href= ! "Pages/Links_Page.htm">Links</a> ! </p> ! </li> ! <li> ! <p class="c6"> ! <a class="c5" href= ! "Pages/Photos_Page.htm">Photos</a> ! </p> ! </li> ! </ul> ! </td> ! <td width="78%"> ! <b>What is DDA?</b> ! <p class="c10"> ! DDA, Discontinuous Deformation Analysis, is a ! numerical method, similar to the finite element ! method, for solving stress-displacement problems, but ! for a jointed rock mass. DDA uses an implicit ! algorithm for simultaneous solution of the equations ! of equilibrium by minimizing the total potential ! energy of the blocky rock mass system. The relation ! between adjacent blocks is governed by equations of ! contact interpenetration and accounts for friction. ! DDA adopts a stepwise approach to solve for the large ! displacements that accompany discontinuous movements ! between blocks. The blocks are said to be "simply ! deformable". Since the method accounts for the ! inertial forces of the blocks' mass, it can be used ! to solve the full dynamic problem of block motion. ! </p> ! <p class="c11"> ! Using DDA ! </p> ! <p class="c10"> ! DDA was originally presented as a complete ! methodology for numerically modeling discontinuous ! rock mass behavior. Because the behavior of fractured ! rock masses is largely controlled by the geometry of ! the fracture system, fracture geometry pre-processing ! is integral part of an analysis. In the current ! incarnation, DDA accepts lists of fractures and ! reduces the fractures to a block system. All this ! geometry adds complexity for the user as well as the ! programmer. The end result is that even in the ! Windows versions, DDA should be considered at least ! as complex as the average finite element code, with ! additional considerations for fracture geometry, plus ! all the difficulties common to any particle code ! (i.e., contact mechanics, etc.). ! </p> ! <p> ! ! </p> ! <p class="c2" dir="ltr"> ! <img border="0" src="Images/cableu.jpg" width="240" ! height="316" /><img border="0" src= ! "Images/cablew.jpg" width="240" height="315" /> ! </p> ! <p class="c2" dir="ltr"> ! (a) ! (b) ! </p> ! <p class="c2"> ! Axial load along the cable bolt that installed at ! roof of<br /> ! underground excavation: (a) without plate and (b) ! with plate (<a class="c12" href= ! "Pages/References_Page.htm#Moosavi">M. Moosavi, R. ! Grayeli</a>)<b>.</b> ! </p> ! <p class="c2" dir="ltr"> ! ! </p> ! </td> ! </tr> ! </table> ! </body> ! </html> |
|
From: David M. D. <do...@us...> - 2006-07-06 02:28:41
|
Update of /cvsroot/dda/ntdda/include In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv30991 Modified Files: analysisdata.h contacts.h datalog.h dda.h ddamemory.h functions.h geometrydata.h graphics.h gravity.h statusbar.h toolbar.h Log Message: Cleaned up code in unittests to ensure compile under cygwin. Index: dda.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/dda.h,v retrieving revision 1.19 retrieving revision 1.20 diff -C2 -d -r1.19 -r1.20 *** dda.h 7 Mar 2006 05:14:59 -0000 1.19 --- dda.h 6 Jul 2006 02:28:37 -0000 1.20 *************** *** 127,131 **** ! void dda_set_output_directory(const char * dirname, size_t dirnamesize); --- 127,131 ---- ! void dda_set_output_directory(const char * dirname, size_t dirnamesize); *************** *** 141,143 **** #endif ! #endif /* __DDA_H__ */ \ No newline at end of file --- 141,144 ---- #endif ! #endif /* __DDA_H__ */ ! Index: statusbar.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/statusbar.h,v retrieving revision 1.2 retrieving revision 1.3 diff -C2 -d -r1.2 -r1.3 *** statusbar.h 26 May 2002 23:47:23 -0000 1.2 --- statusbar.h 6 Jul 2006 02:28:37 -0000 1.3 *************** *** 50,52 **** #endif ! #endif /* __STATUSBAR_H__ */ \ No newline at end of file --- 50,53 ---- #endif ! #endif /* __STATUSBAR_H__ */ ! Index: contacts.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/contacts.h,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** contacts.h 28 Oct 2002 13:46:57 -0000 1.6 --- contacts.h 6 Jul 2006 02:28:37 -0000 1.7 *************** *** 71,73 **** ! #endif /* __CONTACTS_H__ */ \ No newline at end of file --- 71,74 ---- ! #endif /* __CONTACTS_H__ */ ! Index: ddamemory.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/ddamemory.h,v retrieving revision 1.3 retrieving revision 1.4 diff -C2 -d -r1.3 -r1.4 *** ddamemory.h 14 Oct 2002 16:02:46 -0000 1.3 --- ddamemory.h 6 Jul 2006 02:28:37 -0000 1.4 *************** *** 5,8 **** --- 5,12 ---- * $Date$ * $Log$ + * Revision 1.4 2006/07/06 02:28:37 doolin + * Cleaned up code in unittests to ensure compile + * under cygwin. + * * Revision 1.3 2002/10/14 16:02:46 doolin * Lots of small clean up done. No major changes in this commit. *************** *** 57,59 **** #endif ! #endif /* __DDAMEMORY_H__ */ \ No newline at end of file --- 61,64 ---- #endif ! #endif /* __DDAMEMORY_H__ */ ! Index: functions.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/functions.h,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** functions.h 26 Nov 2002 13:13:32 -0000 1.4 --- functions.h 6 Jul 2006 02:28:37 -0000 1.5 *************** *** 139,141 **** #endif ! #endif /* __FUNCTIONS_H__ */ \ No newline at end of file --- 139,143 ---- #endif ! #endif /* __FUNCTIONS_H__ */ ! ! Index: datalog.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/datalog.h,v retrieving revision 1.8 retrieving revision 1.9 diff -C2 -d -r1.8 -r1.9 *** datalog.h 27 Oct 2002 20:53:16 -0000 1.8 --- datalog.h 6 Jul 2006 02:28:37 -0000 1.9 *************** *** 154,156 **** #endif ! #endif /* __DATALOG_H__ */ \ No newline at end of file --- 154,157 ---- #endif ! #endif /* __DATALOG_H__ */ ! Index: gravity.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/gravity.h,v retrieving revision 1.3 retrieving revision 1.4 diff -C2 -d -r1.3 -r1.4 *** gravity.h 27 Oct 2002 20:53:16 -0000 1.3 --- gravity.h 6 Jul 2006 02:28:37 -0000 1.4 *************** *** 91,93 **** ! #endif /* __GRAVITY_H__ */ \ No newline at end of file --- 91,94 ---- ! #endif /* __GRAVITY_H__ */ ! Index: geometrydata.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/geometrydata.h,v retrieving revision 1.17 retrieving revision 1.18 diff -C2 -d -r1.17 -r1.18 *** geometrydata.h 17 Dec 2003 23:36:36 -0000 1.17 --- geometrydata.h 6 Jul 2006 02:28:37 -0000 1.18 *************** *** 305,307 **** ! #endif /* __GEOMETRYDATA_H__ */ \ No newline at end of file --- 305,308 ---- ! #endif /* __GEOMETRYDATA_H__ */ ! Index: toolbar.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/toolbar.h,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** toolbar.h 26 May 2002 23:47:23 -0000 1.4 --- toolbar.h 6 Jul 2006 02:28:37 -0000 1.5 *************** *** 35,37 **** #endif ! #endif /* __TOOLBAR_H__ */ \ No newline at end of file --- 35,38 ---- #endif ! #endif /* __TOOLBAR_H__ */ ! Index: analysisdata.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/analysisdata.h,v retrieving revision 1.25 retrieving revision 1.26 diff -C2 -d -r1.25 -r1.26 *** analysisdata.h 20 Dec 2002 05:20:01 -0000 1.25 --- analysisdata.h 6 Jul 2006 02:28:37 -0000 1.26 *************** *** 486,487 **** --- 486,488 ---- #endif /* __ANALYSISDATA_H__ */ + Index: graphics.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/graphics.h,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** graphics.h 28 May 2002 12:18:31 -0000 1.5 --- graphics.h 6 Jul 2006 02:28:37 -0000 1.6 *************** *** 47,49 **** int freeGraphicStruct(GRAPHICS *); ! #endif /* __GRAPHICS_H__ */ \ No newline at end of file --- 47,50 ---- int freeGraphicStruct(GRAPHICS *); ! #endif /* __GRAPHICS_H__ */ ! |
|
From: David M. D. <do...@us...> - 2006-07-06 02:28:30
|
Update of /cvsroot/dda/ntdda/src In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv30965 Modified Files: analysisdata.c constants.c dda.c ddafile.c dxf.c ghssolver.c loadpoint.c Log Message: Cleaned up code in unittests to ensure compile under cygwin. Index: dxf.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/dxf.c,v retrieving revision 1.2 retrieving revision 1.3 diff -C2 -d -r1.2 -r1.3 *** dxf.c 1 Jul 2006 14:35:44 -0000 1.2 --- dxf.c 6 Jul 2006 02:28:27 -0000 1.3 *************** *** 220,224 **** jy2[count2] = y2; type[count2] = ctype; ! } if (strcmp (str[i], "MTEXT") == 0) { if (strcmp (str[i + 9], "62") == 0) { --- 220,226 ---- jy2[count2] = y2; type[count2] = ctype; ! } ! ! if (strcmp (str[i], "MTEXT") == 0) { if (strcmp (str[i + 9], "62") == 0) { Index: analysisdata.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/analysisdata.c,v retrieving revision 1.27 retrieving revision 1.28 diff -C2 -d -r1.27 -r1.28 *** analysisdata.c 7 Mar 2006 05:39:53 -0000 1.27 --- analysisdata.c 6 Jul 2006 02:28:27 -0000 1.28 *************** *** 554,556 **** adata_get_grav_accel(Analysisdata * ad){ return ad->gravaccel; ! } \ No newline at end of file --- 554,558 ---- adata_get_grav_accel(Analysisdata * ad){ return ad->gravaccel; ! } ! ! Index: dda.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/dda.c,v retrieving revision 1.13 retrieving revision 1.14 diff -C2 -d -r1.13 -r1.14 *** dda.c 7 Mar 2006 05:47:17 -0000 1.13 --- dda.c 6 Jul 2006 02:28:27 -0000 1.14 *************** *** 11,14 **** --- 11,18 ---- #include <direct.h> #define getcwd _getcwd + #define mkdir _mkdir + #define putenv _putenv + #else + #include <unistd.h> #endif *************** *** 188,192 **** * done with mkdir. Smells like a kludge. */ ! checkval = _mkdir(outdir); --- 192,196 ---- * done with mkdir. Smells like a kludge. */ ! checkval = mkdir(outdir); *************** *** 227,230 **** --- 231,235 ---- + #if 0 #ifdef STANDALONE int *************** *** 239,243 **** } #endif /* STANDALONE */ ! --- 244,248 ---- } #endif /* STANDALONE */ ! #endif Index: ddafile.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/ddafile.c,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** ddafile.c 2 Jul 2006 23:07:27 -0000 1.6 --- ddafile.c 6 Jul 2006 02:28:27 -0000 1.7 *************** *** 9,13 **** --- 9,18 ---- #include <string.h> #include <assert.h> + + #ifdef WIN32 #include <direct.h> + #else + #include <unistd.h> + #endif #include "dda.h" *************** *** 389,391 **** } ! #endif \ No newline at end of file --- 394,397 ---- } ! #endif ! Index: constants.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/constants.c,v retrieving revision 1.9 retrieving revision 1.10 diff -C2 -d -r1.9 -r1.10 *** constants.c 7 Mar 2006 05:39:53 -0000 1.9 --- constants.c 6 Jul 2006 02:28:27 -0000 1.10 *************** *** 368,376 **** printer(stream,I1"</AConstants>\n\n"); /* ! openclose = .0002; /* s0 = .0002 ! opencriteria = .0000002; /* f0 = .0000002; ! norm_spring_pen = .0004;//g3 = .0004; /* g3 = .0004; ! angle_olap = 3; /* h1 = 3; ! shear_norm_ratio = 2.5; /* h2 = 2.5; */ } --- 368,376 ---- printer(stream,I1"</AConstants>\n\n"); /* ! openclose = .0002; // s0 = .0002 ! opencriteria = .0000002; // f0 = .0000002; ! norm_spring_pen = .0004;//g3 = .0004; // g3 = .0004; ! angle_olap = 3; // h1 = 3; ! shear_norm_ratio = 2.5; // h2 = 2.5; */ } Index: ghssolver.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/ghssolver.c,v retrieving revision 1.14 retrieving revision 1.15 diff -C2 -d -r1.14 -r1.15 *** ghssolver.c 17 Dec 2003 23:36:36 -0000 1.14 --- ghssolver.c 6 Jul 2006 02:28:27 -0000 1.15 *************** *** 991,993 **** fprintf(fp.dfile,"\n"); ! } /* close writeSolutionVector() */ \ No newline at end of file --- 991,994 ---- fprintf(fp.dfile,"\n"); ! } /* close writeSolutionVector() */ ! Index: loadpoint.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/loadpoint.c,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** loadpoint.c 17 Dec 2003 23:36:36 -0000 1.6 --- loadpoint.c 6 Jul 2006 02:28:27 -0000 1.7 *************** *** 211,213 **** #ifdef __cplusplus } ! #endif \ No newline at end of file --- 211,214 ---- #ifdef __cplusplus } ! #endif ! |
|
From: David M. D. <do...@us...> - 2006-07-06 02:28:09
|
Update of /cvsroot/dda/ntdda/src/unittests In directory sc8-pr-cvs12.sourceforge.net:/tmp/cvs-serv30714 Modified Files: Makefile Log Message: Cleaned up code in unittests to ensure compile under cygwin. Index: Makefile =================================================================== RCS file: /cvsroot/dda/ntdda/src/unittests/Makefile,v retrieving revision 1.17 retrieving revision 1.18 diff -C2 -d -r1.17 -r1.18 *** Makefile 31 Oct 2002 15:35:04 -0000 1.17 --- Makefile 6 Jul 2006 02:28:05 -0000 1.18 *************** *** 25,28 **** --- 25,29 ---- ../ddamemory.c \ ../loadpoint.c \ + ../dda.c \ ../analysisreader.c \ ../ddafile.c *************** *** 42,47 **** material_test \ bolttest \ - analysisdata_test \ - geometrydata_test \ inpoly_test \ ddadlist_test \ --- 43,46 ---- |
|
From: David M. D. <do...@us...> - 2006-07-02 23:37:50
|
Update of /cvsroot/dda/ntdda/include In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv9444/include Modified Files: dxf.h mohrcoulomb.h Log Message: Roozbeh's Mohr-Coulomb code added to project. Index: mohrcoulomb.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/mohrcoulomb.h,v retrieving revision 1.2 retrieving revision 1.3 diff -C2 -d -r1.2 -r1.3 *** mohrcoulomb.h 2 Jul 2006 23:22:45 -0000 1.2 --- mohrcoulomb.h 2 Jul 2006 23:37:45 -0000 1.3 *************** *** 50,53 **** --- 50,56 ---- double sigma[]); + #ifdef __cplusplus + } + #endif #endif /* __MOHRCOULOMB_H__ */ Index: dxf.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/dxf.h,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** dxf.h 1 Jul 2006 05:25:53 -0000 1.1 --- dxf.h 2 Jul 2006 23:37:45 -0000 1.2 *************** *** 20,23 **** --- 20,26 ---- + #ifdef __cplusplus + } + #endif #endif /* __DXF_H__ */ |
|
From: David M. D. <do...@us...> - 2006-07-02 23:37:50
|
Update of /cvsroot/dda/ntdda In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv9444 Modified Files: ntdda.dsp Log Message: Roozbeh's Mohr-Coulomb code added to project. Index: ntdda.dsp =================================================================== RCS file: /cvsroot/dda/ntdda/ntdda.dsp,v retrieving revision 1.40 retrieving revision 1.41 diff -C2 -d -r1.40 -r1.41 *** ntdda.dsp 1 Jul 2006 05:27:09 -0000 1.40 --- ntdda.dsp 2 Jul 2006 23:37:45 -0000 1.41 *************** *** 170,173 **** --- 170,177 ---- # End Source File # End Group + # Begin Source File + + SOURCE=.\src\mohrcoulomb.c + # End Source File # End Group # Begin Group "gui" |
|
From: David M. D. <do...@us...> - 2006-07-02 23:31:41
|
Update of /cvsroot/dda/ntdda/src In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv6957 Modified Files: mohrcoulomb.c Log Message: File level variables declared static so namespace doesn't get polluted. THIS IS IMPORTANT Index: mohrcoulomb.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/mohrcoulomb.c,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** mohrcoulomb.c 2 Jul 2006 23:29:30 -0000 1.6 --- mohrcoulomb.c 2 Jul 2006 23:31:37 -0000 1.7 *************** *** 25,31 **** ! double sigm, dsbar, theta, cs[5], pl[5][5]; ! int i, j; /* --- 25,31 ---- ! static double sigm, dsbar, theta, cs[5], pl[5][5]; ! static int i, j; /* |
|
From: David M. D. <do...@us...> - 2006-07-02 23:29:33
|
Update of /cvsroot/dda/ntdda/src In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv6082 Modified Files: mohrcoulomb.c Log Message: Added directives to allow mohrcoulomb to be compiled into a free standing executable for testing purposes. Index: mohrcoulomb.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/mohrcoulomb.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** mohrcoulomb.c 2 Jul 2006 23:18:40 -0000 1.5 --- mohrcoulomb.c 2 Jul 2006 23:29:30 -0000 1.6 *************** *** 17,20 **** --- 17,28 ---- */ + #ifdef __cplusplus + extern "C" { + #endif + #if 0 + } + #endif + + double sigm, dsbar, theta, cs[5], pl[5][5]; *************** *** 201,202 **** --- 209,230 ---- } + + #ifdef __cplusplus + } + #endif + + #ifdef STANDALONE + /** + * Free standing program for testing Mohr-Coulomb + * computations. + * + * Compile with: + * gcc -Wall -o mc mohrcoulomb.c -I../include -DSTANDALONE + * + */ + int + main(int argc, char ** argv) { + + return 0; + } + #endif /* STANDALONE */ |
|
From: David M. D. <do...@us...> - 2006-07-02 23:22:49
|
Update of /cvsroot/dda/ntdda/include In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv3231 Modified Files: mohrcoulomb.h Log Message: reformatted header to bring into modec standard. Index: mohrcoulomb.h =================================================================== RCS file: /cvsroot/dda/ntdda/include/mohrcoulomb.h,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** mohrcoulomb.h 2 Jul 2006 23:07:00 -0000 1.1 --- mohrcoulomb.h 2 Jul 2006 23:22:45 -0000 1.2 *************** *** 16,24 **** ! void invar(); ! double mocouf(double c,double phi); ! void mocopl(double phi,double psi,double nu, double ymod); ! void mohrcoulomb(double c, double phi, double psi,double nu,double ymod, ! double eps[],double * stress,double sigma[]); --- 16,52 ---- ! /** ! * Add explanatory comment here. ! */ ! void invar (void); ! ! ! /** ! * Add explanatory comment here. ! */ ! double mocouf (double c, ! double phi); ! ! ! /** ! * Add explanatory comment here. ! */ ! void mocopl (double phi, ! double psi, ! double nu, ! double ymod); ! ! ! /** ! * Add explanatory comment here. ! */ ! void mohrcoulomb (double c, ! double phi, ! double psi, ! double nu, ! double ymod, ! double eps[], ! double * stress, ! double sigma[]); |
|
From: David M. D. <do...@us...> - 2006-07-02 23:18:43
|
Update of /cvsroot/dda/ntdda/src In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv1478 Modified Files: mohrcoulomb.c Log Message: Added more white space for readability, fully braced around if anf for statements. Index: mohrcoulomb.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/mohrcoulomb.c,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** mohrcoulomb.c 2 Jul 2006 23:16:02 -0000 1.4 --- mohrcoulomb.c 2 Jul 2006 23:18:40 -0000 1.5 *************** *** 46,53 **** xj3 = dx * dy * dz - dz * pow (txy, 2); sine = -13.5 * xj3 / pow (dsbar, 3); ! if (sine > 1) sine = 1; ! if (sine < -1) sine = -1; theta = asin (sine) / 3; } --- 46,55 ---- xj3 = dx * dy * dz - dz * pow (txy, 2); sine = -13.5 * xj3 / pow (dsbar, 3); ! if (sine > 1) { sine = 1; ! } ! if (sine < -1) { sine = -1; + } theta = asin (sine) / 3; } *************** *** 103,116 **** d3 = dx * dy * dz - dz * txy * txy; th = -3 * sq3 * d3 / (2 * pow (d2, 3)); ! if (th > 1) th = 1; ! if (th < -1) th = -1; th = asin (th) / 3; snth = sin (th); if (fabs (snth) > 0.49) { sig = -1; ! if (snth < 0) sig = 1; rph = snph * (1 + nu) / 3; rps = snps * (1. + nu) / 3; --- 105,121 ---- d3 = dx * dy * dz - dz * txy * txy; th = -3 * sq3 * d3 / (2 * pow (d2, 3)); ! if (th > 1) { th = 1; ! } ! if (th < -1) { th = -1; + } th = asin (th) / 3; snth = sin (th); if (fabs (snth) > 0.49) { sig = -1; ! if (snth < 0) { sig = 1; + } rph = snph * (1 + nu) / 3; rps = snps * (1. + nu) / 3; *************** *** 171,178 **** double fnew; ! for (i = 1; i <= 4; i++) elso[i] = 0; ! for (i = 1; i <= 4; i++) cs[i] = stress[i + 3] + sigma[i]; invar (); --- 176,186 ---- double fnew; ! for (i = 1; i <= 4; i++) { elso[i] = 0; ! } ! ! for (i = 1; i <= 4; i++) { cs[i] = stress[i + 3] + sigma[i]; + } invar (); *************** *** 188,192 **** } ! for (i = 1; i <= 4; i++) stress[i + 3] = stress[i + 3] + sigma[i] - elso[i]; } --- 196,202 ---- } ! for (i = 1; i <= 4; i++) { stress[i + 3] = stress[i + 3] + sigma[i] - elso[i]; + } + } |
|
From: David M. D. <do...@us...> - 2006-07-02 23:16:06
|
Update of /cvsroot/dda/ntdda/src In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv32731 Modified Files: mohrcoulomb.c Log Message: Code formatting brought more inline with existing DDA code using indent -br -nut -ce mohrcoulomb.c. Index: mohrcoulomb.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/mohrcoulomb.c,v retrieving revision 1.3 retrieving revision 1.4 diff -C2 -d -r1.3 -r1.4 *** mohrcoulomb.c 2 Jul 2006 23:12:25 -0000 1.3 --- mohrcoulomb.c 2 Jul 2006 23:16:02 -0000 1.4 *************** *** 16,31 **** * theta: lode angle */ - - double sigm,dsbar,theta,cs[5],pl[5][5]; ! int i,j; /* * first invariant "sigm", second invariant "dsbar" * and lode angle "theta" for a stress tensor "stress" ! */ ! void invar() { ! double sx,sy,sz,txy,dx,dy,dz,xj3,sine; sx = cs[1]; --- 16,32 ---- * theta: lode angle */ ! double sigm, dsbar, theta, cs[5], pl[5][5]; ! ! int i, j; /* * first invariant "sigm", second invariant "dsbar" * and lode angle "theta" for a stress tensor "stress" ! */ ! void ! invar () { ! double sx, sy, sz, txy, dx, dy, dz, xj3, sine; sx = cs[1]; *************** *** 33,52 **** txy = cs[3]; sz = cs[4]; ! sigm = (sx+sy+sz)/3; ! dsbar = sqrt(pow((sx-sy),2)+pow((sy-sz),2)+pow((sz-sx),2)+6.*pow(txy,2))/sqrt(2.); ! if(dsbar < 1.e-10) ! { ! theta = .0; ! } ! else ! { ! dx = (2*sx-sy-sz)/3; ! dy = (2*sy-sz-sx)/3; ! dz = (2*sz-sx-sy)/3; ! xj3 = dx*dy*dz-dz*pow(txy,2); ! sine = -13.5*xj3/pow(dsbar,3); ! if(sine > 1) sine = 1; ! if(sine < -1) sine = -1; ! theta = asin(sine)/3; } --- 34,54 ---- txy = cs[3]; sz = cs[4]; ! sigm = (sx + sy + sz) / 3; ! dsbar = ! sqrt (pow ((sx - sy), 2) + pow ((sy - sz), 2) + pow ((sz - sx), 2) + ! 6. * pow (txy, 2)) / sqrt (2.); ! if (dsbar < 1.e-10) { ! theta = .0; ! } else { ! dx = (2 * sx - sy - sz) / 3; ! dy = (2 * sy - sz - sx) / 3; ! dz = (2 * sz - sx - sy) / 3; ! xj3 = dx * dy * dz - dz * pow (txy, 2); ! sine = -13.5 * xj3 / pow (dsbar, 3); ! if (sine > 1) ! sine = 1; ! if (sine < -1) ! sine = -1; ! theta = asin (sine) / 3; } *************** *** 58,70 **** * and cohesion "c" */ ! double mocouf(double c,double phi) ! { ! double phir,snph,csph,csth,snth,f; ! phir = phi*4*atan(1)/180; ! snph = sin(phir); ! csph = cos(phir); ! csth = cos(theta); ! snth = sin(theta); ! f = snph*sigm+dsbar*(csth/sqrt(3.)-snth*snph/3.)-c*csph; return f; } --- 60,73 ---- * and cohesion "c" */ ! double ! mocouf (double c, double phi) ! { ! double phir, snph, csph, csth, snth, f; ! phir = phi * 4 * atan (1) / 180; ! snph = sin (phir); ! csph = cos (phir); ! csth = cos (theta); ! snth = sin (theta); ! f = snph * sigm + dsbar * (csth / sqrt (3.) - snth * snph / 3.) - c * csph; return f; } *************** *** 75,155 **** * ,"nu" for a Mohr-Coulomb material */ ! void mocopl(double phi,double psi,double nu, double ymod) { ! double row[5],col[5],sx,sy,txy,sz,pi,phir,psir; ! double dx,dy,dz,d2,d3,th,snth,sig,rph,rps,cps,snps,sq3; ! double cc,cph,alp,ca,sa,dd,snph,ee,s1,s2; ! sx = cs[1]; ! sy = cs[2]; ! txy = cs[3]; ! sz = cs[4]; ! pi = 4*atan(1); ! phir = phi*pi/180; ! psir = psi*pi/180; ! snph = sin(phir); ! snps = sin(psir); ! sq3 = sqrt(3); ! cc = 1-2*nu; ! dx = (2*sx-sy-sz)/3; ! dy = (2*sy-sz-sx)/3; ! dz = (2*sz-sx-sy)/3; ! d2 = sqrt(-dx*dy-dy*dz-dz*dx+txy*txy); ! d3 = dx*dy*dz-dz*txy*txy; ! th = -3*sq3*d3/(2*pow(d2,3)); ! if(th > 1)th = 1; ! if(th < -1)th = -1; ! th = asin(th)/3; ! snth = sin(th); ! if(fabs(snth) > 0.49) ! { ! sig = -1; ! if(snth < 0)sig = 1; ! rph = snph*(1+nu)/3; ! rps = snps*(1.+nu)/3; ! cps = 0.25*sq3/d2*(1+sig*snps/3); ! cph = 0.25*sq3/d2*(1+sig*snph/3); ! col[1] = rph+cph*((1-nu)*dx+nu*(dy+dz)); ! col[2] = rph+cph*((1-nu)*dy+nu*(dz+dx)); ! col[3] = cph*cc*txy; ! col[4] = rph+cph*((1-nu)*dz+nu*(dx+dy)); ! row[1] = rps+cps*((1-nu)*dx+nu*(dy+dz)); ! row[2] = rps+cps*((1-nu)*dy+nu*(dz+dx)); ! row[3] = cps*cc*txy; ! row[4] = rps+cps*((1-nu)*dz+nu*(dx+dy)); ! ee = ymod/((1+nu)*cc*(rph*snps+2*cph*cps*d2*d2*cc)); ! } ! else ! { ! alp = atan(fabs((sx-sy)/(2*txy))); ! ca = cos(alp); ! sa = sin(alp); ! dd = cc*sa; ! s1=1; ! s2=1; ! if((sx-sy) < 0)s1 = -1; ! if(txy < 0)s2 = -1; ! col[1] = snph+s1*dd; ! col[2] = snph-s1*dd; ! col[3] = s2*cc*ca; ! col[4] = 2*nu*snph; ! row[1] = snps+s1*dd; ! row[2] = snps-s1*dd; ! row[3] = s2*cc*ca; ! row[4] = 2*nu*snps; ! ee = ymod/(2*(1.+nu)*cc*(snph*snps+cc)); } ! for(i=1; i<=4; i++) ! { ! for(j=1; j<=4; j++) ! { ! pl[i][j] = ee*row[i]*col[j]; ! } } } --- 78,159 ---- * ,"nu" for a Mohr-Coulomb material */ ! void ! mocopl (double phi, double psi, double nu, double ymod) { ! double row[5], col[5], sx, sy, txy, sz, pi, phir, psir; ! double dx, dy, dz, d2, d3, th, snth, sig, rph, rps, cps, snps, sq3; ! double cc, cph, alp, ca, sa, dd, snph, ee, s1, s2; ! sx = cs[1]; ! sy = cs[2]; ! txy = cs[3]; ! sz = cs[4]; ! pi = 4 * atan (1); ! phir = phi * pi / 180; ! psir = psi * pi / 180; ! snph = sin (phir); ! snps = sin (psir); ! sq3 = sqrt (3); ! cc = 1 - 2 * nu; ! dx = (2 * sx - sy - sz) / 3; ! dy = (2 * sy - sz - sx) / 3; ! dz = (2 * sz - sx - sy) / 3; ! d2 = sqrt (-dx * dy - dy * dz - dz * dx + txy * txy); ! d3 = dx * dy * dz - dz * txy * txy; ! th = -3 * sq3 * d3 / (2 * pow (d2, 3)); ! if (th > 1) ! th = 1; ! if (th < -1) ! th = -1; ! th = asin (th) / 3; ! snth = sin (th); ! if (fabs (snth) > 0.49) { ! sig = -1; ! if (snth < 0) ! sig = 1; ! rph = snph * (1 + nu) / 3; ! rps = snps * (1. + nu) / 3; ! cps = 0.25 * sq3 / d2 * (1 + sig * snps / 3); ! cph = 0.25 * sq3 / d2 * (1 + sig * snph / 3); ! col[1] = rph + cph * ((1 - nu) * dx + nu * (dy + dz)); ! col[2] = rph + cph * ((1 - nu) * dy + nu * (dz + dx)); ! col[3] = cph * cc * txy; ! col[4] = rph + cph * ((1 - nu) * dz + nu * (dx + dy)); ! row[1] = rps + cps * ((1 - nu) * dx + nu * (dy + dz)); ! row[2] = rps + cps * ((1 - nu) * dy + nu * (dz + dx)); ! row[3] = cps * cc * txy; ! row[4] = rps + cps * ((1 - nu) * dz + nu * (dx + dy)); ! ee = ymod / ((1 + nu) * cc * (rph * snps + 2 * cph * cps * d2 * d2 * cc)); ! } else { ! alp = atan (fabs ((sx - sy) / (2 * txy))); ! ca = cos (alp); ! sa = sin (alp); ! dd = cc * sa; ! s1 = 1; ! s2 = 1; ! if ((sx - sy) < 0) ! s1 = -1; ! if (txy < 0) ! s2 = -1; ! col[1] = snph + s1 * dd; ! col[2] = snph - s1 * dd; ! col[3] = s2 * cc * ca; ! col[4] = 2 * nu * snph; ! row[1] = snps + s1 * dd; ! row[2] = snps - s1 * dd; ! row[3] = s2 * cc * ca; ! row[4] = 2 * nu * snps; ! ee = ymod / (2 * (1. + nu) * cc * (snph * snps + cc)); } ! for (i = 1; i <= 4; i++) { ! for (j = 1; j <= 4; j++) { ! pl[i][j] = ee * row[i] * col[j]; ! } } } *************** *** 159,185 **** * correct stress using Mohr-Coulomb criteria */ ! void mohrcoulomb(double c,double phi,double psi,double nu,double ymod, ! double eps[4],double * stress,double sigma[4]) { ! double elso[5]; ! double fnew; ! for(i = 1; i <= 4; i++) elso[i] = 0; ! for(i = 1; i <= 4; i++) cs[i] = stress[i+3] + sigma[i]; ! invar(); ! fnew = mocouf(c,phi); ! if(fnew > 0) ! { ! mocopl(phi,psi,nu,ymod); ! elso[1] = pl[1][1] * eps[1] + pl[1][2] * eps[2] + pl[1][3] * eps[3]; ! elso[2] = pl[2][1] * eps[1] + pl[2][2] * eps[2] + pl[2][3] * eps[3]; ! elso[3] = pl[3][1] * eps[1] + pl[3][2] * eps[2] + pl[3][3] * eps[3]; ! elso[4] = pl[4][1] * eps[1] + pl[4][2] * eps[2] + pl[4][3] * eps[3]; ! } ! ! for(i = 1; i <= 4; i++) stress[i+3] = stress[i+3] + sigma[i] - elso[i]; } --- 163,192 ---- * correct stress using Mohr-Coulomb criteria */ ! void ! mohrcoulomb (double c, double phi, double psi, double nu, double ymod, ! double eps[4], double *stress, double sigma[4]) { ! double elso[5]; ! double fnew; ! for (i = 1; i <= 4; i++) ! elso[i] = 0; ! for (i = 1; i <= 4; i++) ! cs[i] = stress[i + 3] + sigma[i]; ! invar (); ! fnew = mocouf (c, phi); ! if (fnew > 0) { ! mocopl (phi, psi, nu, ymod); ! elso[1] = pl[1][1] * eps[1] + pl[1][2] * eps[2] + pl[1][3] * eps[3]; ! elso[2] = pl[2][1] * eps[1] + pl[2][2] * eps[2] + pl[2][3] * eps[3]; ! elso[3] = pl[3][1] * eps[1] + pl[3][2] * eps[2] + pl[3][3] * eps[3]; ! elso[4] = pl[4][1] * eps[1] + pl[4][2] * eps[2] + pl[4][3] * eps[3]; ! } ! ! for (i = 1; i <= 4; i++) ! stress[i + 3] = stress[i + 3] + sigma[i] - elso[i]; } |
|
From: David M. D. <do...@us...> - 2006-07-02 23:12:29
|
Update of /cvsroot/dda/ntdda/src In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv31446 Modified Files: mohrcoulomb.c Log Message: Reformat comments to match existing DDA code. Index: mohrcoulomb.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/mohrcoulomb.c,v retrieving revision 1.2 retrieving revision 1.3 diff -C2 -d -r1.2 -r1.3 *** mohrcoulomb.c 2 Jul 2006 23:08:21 -0000 1.2 --- mohrcoulomb.c 2 Jul 2006 23:12:25 -0000 1.3 *************** *** 6,19 **** #include "mohrcoulomb.h" ! /* sigm: mean stress invariant ! /* dsbar: deviatoric stress invariant ! /* fnew or f: value of failure function ! /* phi: friction angle ! /* c: cohesion ! /* psi: dilation angle ! /* tnph: tan(phi) ! /* snph: sin(phi) ! /* theta: lode angle ! */ double sigm,dsbar,theta,cs[5],pl[5][5]; --- 6,19 ---- #include "mohrcoulomb.h" ! /* sigm: mean stress invariant ! * dsbar: deviatoric stress invariant ! * fnew or f: value of failure function ! * phi: friction angle ! * c: cohesion ! * psi: dilation angle ! * tnph: tan(phi) ! * snph: sin(phi) ! * theta: lode angle ! */ double sigm,dsbar,theta,cs[5],pl[5][5]; *************** *** 22,29 **** /* ! /* first invariant "sigm", second invariant "dsbar" ! /* and lode angle "theta" for a stress tensor "stress" ! */ ! void invar() { --- 22,28 ---- /* ! * first invariant "sigm", second invariant "dsbar" ! * and lode angle "theta" for a stress tensor "stress" ! */ void invar() { *************** *** 55,63 **** /* ! /* Mohr-Coulomb yield function F from invariants "sigm" ! /* and "dsbar" lode angle "theta (radian) of friction "phi" ! /* and cohesion "c" ! */ ! double mocouf(double c,double phi) { --- 54,61 ---- /* ! * Mohr-Coulomb yield function F from invariants "sigm" ! * and "dsbar" lode angle "theta (radian) of friction "phi" ! * and cohesion "c" ! */ double mocouf(double c,double phi) { *************** *** 73,81 **** /* ! /* plastic stress-strain matrix "pl" from stresses "stress" ! /* angle of friction "phi", dilation "psi" and elastic stiffness "E" ! /* ,"nu" for a Mohr-Coulomb material ! */ ! void mocopl(double phi,double psi,double nu, double ymod) { --- 71,78 ---- /* ! * plastic stress-strain matrix "pl" from stresses "stress" ! * angle of friction "phi", dilation "psi" and elastic stiffness "E" ! * ,"nu" for a Mohr-Coulomb material ! */ void mocopl(double phi,double psi,double nu, double ymod) { *************** *** 159,165 **** /* ! /* main function in order to obtain the ! /* correct stress using Mohr-Coulomb criteria ! */ void mohrcoulomb(double c,double phi,double psi,double nu,double ymod, double eps[4],double * stress,double sigma[4]) --- 156,162 ---- /* ! * main function in order to obtain the ! * correct stress using Mohr-Coulomb criteria ! */ void mohrcoulomb(double c,double phi,double psi,double nu,double ymod, double eps[4],double * stress,double sigma[4]) *************** *** 186,188 **** for(i = 1; i <= 4; i++) stress[i+3] = stress[i+3] + sigma[i] - elso[i]; ! } \ No newline at end of file --- 183,185 ---- for(i = 1; i <= 4; i++) stress[i+3] = stress[i+3] + sigma[i] - elso[i]; ! } |
|
From: David M. D. <do...@us...> - 2006-07-02 23:08:25
|
Update of /cvsroot/dda/ntdda/src In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv29725 Modified Files: mohrcoulomb.c Log Message: Line ending fixes. Index: mohrcoulomb.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/mohrcoulomb.c,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** mohrcoulomb.c 2 Jul 2006 23:07:28 -0000 1.1 --- mohrcoulomb.c 2 Jul 2006 23:08:21 -0000 1.2 *************** *** 1,188 **** ! ! #include <stdio.h> ! #include <stdlib.h> ! #include <math.h> ! ! #include "mohrcoulomb.h" ! ! /* sigm: mean stress invariant ! /* dsbar: deviatoric stress invariant ! /* fnew or f: value of failure function ! /* phi: friction angle ! /* c: cohesion ! /* psi: dilation angle ! /* tnph: tan(phi) ! /* snph: sin(phi) ! /* theta: lode angle ! */ ! ! double sigm,dsbar,theta,cs[5],pl[5][5]; ! ! int i,j; ! ! /* ! /* first invariant "sigm", second invariant "dsbar" ! /* and lode angle "theta" for a stress tensor "stress" ! */ ! ! void invar() ! { ! double sx,sy,sz,txy,dx,dy,dz,xj3,sine; ! ! sx = cs[1]; ! sy = cs[2]; ! txy = cs[3]; ! sz = cs[4]; ! sigm = (sx+sy+sz)/3; ! dsbar = sqrt(pow((sx-sy),2)+pow((sy-sz),2)+pow((sz-sx),2)+6.*pow(txy,2))/sqrt(2.); ! if(dsbar < 1.e-10) ! { ! theta = .0; ! } ! else ! { ! dx = (2*sx-sy-sz)/3; ! dy = (2*sy-sz-sx)/3; ! dz = (2*sz-sx-sy)/3; ! xj3 = dx*dy*dz-dz*pow(txy,2); ! sine = -13.5*xj3/pow(dsbar,3); ! if(sine > 1) sine = 1; ! if(sine < -1) sine = -1; ! theta = asin(sine)/3; ! } ! ! } ! ! /* ! /* Mohr-Coulomb yield function F from invariants "sigm" ! /* and "dsbar" lode angle "theta (radian) of friction "phi" ! /* and cohesion "c" ! */ ! ! double mocouf(double c,double phi) ! { ! double phir,snph,csph,csth,snth,f; ! phir = phi*4*atan(1)/180; ! snph = sin(phir); ! csph = cos(phir); ! csth = cos(theta); ! snth = sin(theta); ! f = snph*sigm+dsbar*(csth/sqrt(3.)-snth*snph/3.)-c*csph; ! return f; ! } ! ! /* ! /* plastic stress-strain matrix "pl" from stresses "stress" ! /* angle of friction "phi", dilation "psi" and elastic stiffness "E" ! /* ,"nu" for a Mohr-Coulomb material ! */ ! ! void mocopl(double phi,double psi,double nu, double ymod) ! { ! ! double row[5],col[5],sx,sy,txy,sz,pi,phir,psir; ! double dx,dy,dz,d2,d3,th,snth,sig,rph,rps,cps,snps,sq3; ! double cc,cph,alp,ca,sa,dd,snph,ee,s1,s2; ! ! sx = cs[1]; ! sy = cs[2]; ! txy = cs[3]; ! sz = cs[4]; ! pi = 4*atan(1); ! phir = phi*pi/180; ! psir = psi*pi/180; ! snph = sin(phir); ! snps = sin(psir); ! sq3 = sqrt(3); ! cc = 1-2*nu; ! dx = (2*sx-sy-sz)/3; ! dy = (2*sy-sz-sx)/3; ! dz = (2*sz-sx-sy)/3; ! d2 = sqrt(-dx*dy-dy*dz-dz*dx+txy*txy); ! d3 = dx*dy*dz-dz*txy*txy; ! th = -3*sq3*d3/(2*pow(d2,3)); ! if(th > 1)th = 1; ! if(th < -1)th = -1; ! th = asin(th)/3; ! snth = sin(th); ! if(fabs(snth) > 0.49) ! { ! sig = -1; ! if(snth < 0)sig = 1; ! rph = snph*(1+nu)/3; ! rps = snps*(1.+nu)/3; ! cps = 0.25*sq3/d2*(1+sig*snps/3); ! cph = 0.25*sq3/d2*(1+sig*snph/3); ! ! col[1] = rph+cph*((1-nu)*dx+nu*(dy+dz)); ! col[2] = rph+cph*((1-nu)*dy+nu*(dz+dx)); ! col[3] = cph*cc*txy; ! col[4] = rph+cph*((1-nu)*dz+nu*(dx+dy)); ! ! row[1] = rps+cps*((1-nu)*dx+nu*(dy+dz)); ! row[2] = rps+cps*((1-nu)*dy+nu*(dz+dx)); ! row[3] = cps*cc*txy; ! row[4] = rps+cps*((1-nu)*dz+nu*(dx+dy)); ! ! ee = ymod/((1+nu)*cc*(rph*snps+2*cph*cps*d2*d2*cc)); ! } ! else ! { ! alp = atan(fabs((sx-sy)/(2*txy))); ! ca = cos(alp); ! sa = sin(alp); ! dd = cc*sa; ! s1=1; ! s2=1; ! if((sx-sy) < 0)s1 = -1; ! if(txy < 0)s2 = -1; ! col[1] = snph+s1*dd; ! col[2] = snph-s1*dd; ! col[3] = s2*cc*ca; ! col[4] = 2*nu*snph; ! ! row[1] = snps+s1*dd; ! row[2] = snps-s1*dd; ! row[3] = s2*cc*ca; ! row[4] = 2*nu*snps; ! ! ee = ymod/(2*(1.+nu)*cc*(snph*snps+cc)); ! } ! for(i=1; i<=4; i++) ! { ! for(j=1; j<=4; j++) ! { ! pl[i][j] = ee*row[i]*col[j]; ! } ! } ! } ! ! /* ! /* main function in order to obtain the ! /* correct stress using Mohr-Coulomb criteria ! */ ! void mohrcoulomb(double c,double phi,double psi,double nu,double ymod, ! double eps[4],double * stress,double sigma[4]) ! { ! ! double elso[5]; ! double fnew; ! ! for(i = 1; i <= 4; i++) elso[i] = 0; ! for(i = 1; i <= 4; i++) cs[i] = stress[i+3] + sigma[i]; ! ! invar(); ! fnew = mocouf(c,phi); ! ! if(fnew > 0) ! { ! mocopl(phi,psi,nu,ymod); ! ! elso[1] = pl[1][1] * eps[1] + pl[1][2] * eps[2] + pl[1][3] * eps[3]; ! elso[2] = pl[2][1] * eps[1] + pl[2][2] * eps[2] + pl[2][3] * eps[3]; ! elso[3] = pl[3][1] * eps[1] + pl[3][2] * eps[2] + pl[3][3] * eps[3]; ! elso[4] = pl[4][1] * eps[1] + pl[4][2] * eps[2] + pl[4][3] * eps[3]; ! } ! ! for(i = 1; i <= 4; i++) stress[i+3] = stress[i+3] + sigma[i] - elso[i]; } \ No newline at end of file --- 1,188 ---- ! ! #include <stdio.h> ! #include <stdlib.h> ! #include <math.h> ! ! #include "mohrcoulomb.h" ! ! /* sigm: mean stress invariant ! /* dsbar: deviatoric stress invariant ! /* fnew or f: value of failure function ! /* phi: friction angle ! /* c: cohesion ! /* psi: dilation angle ! /* tnph: tan(phi) ! /* snph: sin(phi) ! /* theta: lode angle ! */ ! ! double sigm,dsbar,theta,cs[5],pl[5][5]; ! ! int i,j; ! ! /* ! /* first invariant "sigm", second invariant "dsbar" ! /* and lode angle "theta" for a stress tensor "stress" ! */ ! ! void invar() ! { ! double sx,sy,sz,txy,dx,dy,dz,xj3,sine; ! ! sx = cs[1]; ! sy = cs[2]; ! txy = cs[3]; ! sz = cs[4]; ! sigm = (sx+sy+sz)/3; ! dsbar = sqrt(pow((sx-sy),2)+pow((sy-sz),2)+pow((sz-sx),2)+6.*pow(txy,2))/sqrt(2.); ! if(dsbar < 1.e-10) ! { ! theta = .0; ! } ! else ! { ! dx = (2*sx-sy-sz)/3; ! dy = (2*sy-sz-sx)/3; ! dz = (2*sz-sx-sy)/3; ! xj3 = dx*dy*dz-dz*pow(txy,2); ! sine = -13.5*xj3/pow(dsbar,3); ! if(sine > 1) sine = 1; ! if(sine < -1) sine = -1; ! theta = asin(sine)/3; ! } ! ! } ! ! /* ! /* Mohr-Coulomb yield function F from invariants "sigm" ! /* and "dsbar" lode angle "theta (radian) of friction "phi" ! /* and cohesion "c" ! */ ! ! double mocouf(double c,double phi) ! { ! double phir,snph,csph,csth,snth,f; ! phir = phi*4*atan(1)/180; ! snph = sin(phir); ! csph = cos(phir); ! csth = cos(theta); ! snth = sin(theta); ! f = snph*sigm+dsbar*(csth/sqrt(3.)-snth*snph/3.)-c*csph; ! return f; ! } ! ! /* ! /* plastic stress-strain matrix "pl" from stresses "stress" ! /* angle of friction "phi", dilation "psi" and elastic stiffness "E" ! /* ,"nu" for a Mohr-Coulomb material ! */ ! ! void mocopl(double phi,double psi,double nu, double ymod) ! { ! ! double row[5],col[5],sx,sy,txy,sz,pi,phir,psir; ! double dx,dy,dz,d2,d3,th,snth,sig,rph,rps,cps,snps,sq3; ! double cc,cph,alp,ca,sa,dd,snph,ee,s1,s2; ! ! sx = cs[1]; ! sy = cs[2]; ! txy = cs[3]; ! sz = cs[4]; ! pi = 4*atan(1); ! phir = phi*pi/180; ! psir = psi*pi/180; ! snph = sin(phir); ! snps = sin(psir); ! sq3 = sqrt(3); ! cc = 1-2*nu; ! dx = (2*sx-sy-sz)/3; ! dy = (2*sy-sz-sx)/3; ! dz = (2*sz-sx-sy)/3; ! d2 = sqrt(-dx*dy-dy*dz-dz*dx+txy*txy); ! d3 = dx*dy*dz-dz*txy*txy; ! th = -3*sq3*d3/(2*pow(d2,3)); ! if(th > 1)th = 1; ! if(th < -1)th = -1; ! th = asin(th)/3; ! snth = sin(th); ! if(fabs(snth) > 0.49) ! { ! sig = -1; ! if(snth < 0)sig = 1; ! rph = snph*(1+nu)/3; ! rps = snps*(1.+nu)/3; ! cps = 0.25*sq3/d2*(1+sig*snps/3); ! cph = 0.25*sq3/d2*(1+sig*snph/3); ! ! col[1] = rph+cph*((1-nu)*dx+nu*(dy+dz)); ! col[2] = rph+cph*((1-nu)*dy+nu*(dz+dx)); ! col[3] = cph*cc*txy; ! col[4] = rph+cph*((1-nu)*dz+nu*(dx+dy)); ! ! row[1] = rps+cps*((1-nu)*dx+nu*(dy+dz)); ! row[2] = rps+cps*((1-nu)*dy+nu*(dz+dx)); ! row[3] = cps*cc*txy; ! row[4] = rps+cps*((1-nu)*dz+nu*(dx+dy)); ! ! ee = ymod/((1+nu)*cc*(rph*snps+2*cph*cps*d2*d2*cc)); ! } ! else ! { ! alp = atan(fabs((sx-sy)/(2*txy))); ! ca = cos(alp); ! sa = sin(alp); ! dd = cc*sa; ! s1=1; ! s2=1; ! if((sx-sy) < 0)s1 = -1; ! if(txy < 0)s2 = -1; ! col[1] = snph+s1*dd; ! col[2] = snph-s1*dd; ! col[3] = s2*cc*ca; ! col[4] = 2*nu*snph; ! ! row[1] = snps+s1*dd; ! row[2] = snps-s1*dd; ! row[3] = s2*cc*ca; ! row[4] = 2*nu*snps; ! ! ee = ymod/(2*(1.+nu)*cc*(snph*snps+cc)); ! } ! for(i=1; i<=4; i++) ! { ! for(j=1; j<=4; j++) ! { ! pl[i][j] = ee*row[i]*col[j]; ! } ! } ! } ! ! /* ! /* main function in order to obtain the ! /* correct stress using Mohr-Coulomb criteria ! */ ! void mohrcoulomb(double c,double phi,double psi,double nu,double ymod, ! double eps[4],double * stress,double sigma[4]) ! { ! ! double elso[5]; ! double fnew; ! ! for(i = 1; i <= 4; i++) elso[i] = 0; ! for(i = 1; i <= 4; i++) cs[i] = stress[i+3] + sigma[i]; ! ! invar(); ! fnew = mocouf(c,phi); ! ! if(fnew > 0) ! { ! mocopl(phi,psi,nu,ymod); ! ! elso[1] = pl[1][1] * eps[1] + pl[1][2] * eps[2] + pl[1][3] * eps[3]; ! elso[2] = pl[2][1] * eps[1] + pl[2][2] * eps[2] + pl[2][3] * eps[3]; ! elso[3] = pl[3][1] * eps[1] + pl[3][2] * eps[2] + pl[3][3] * eps[3]; ! elso[4] = pl[4][1] * eps[1] + pl[4][2] * eps[2] + pl[4][3] * eps[3]; ! } ! ! for(i = 1; i <= 4; i++) stress[i+3] = stress[i+3] + sigma[i] - elso[i]; } \ No newline at end of file |
|
From: David M. D. <do...@us...> - 2006-07-02 23:07:31
|
Update of /cvsroot/dda/ntdda/src In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv29287 Modified Files: ddafile.c Added Files: mohrcoulomb.c Log Message: Initial mohrcoulomb file from roozbeh. Index: ddafile.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/ddafile.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** ddafile.c 8 Mar 2006 16:39:20 -0000 1.5 --- ddafile.c 2 Jul 2006 23:07:27 -0000 1.6 *************** *** 137,147 **** // output folder is created when programe is ran. #define BUFSIZE 1024 ! char wdbuf[BUFSIZE]; char * outdir; outdir = malloc(BUFSIZE); ! getcwd(wdbuf,BUFSIZE); dda_set_output_directory("output",sizeof("output")); ! strncpy(outdir,wdbuf,sizeof(wdbuf)); strncat(outdir,"\\",sizeof("\\")); strncat(outdir,"output",sizeof("output")); --- 137,147 ---- // output folder is created when programe is ran. #define BUFSIZE 1024 ! char cwdbuf[BUFSIZE]; char * outdir; outdir = malloc(BUFSIZE); ! getcwd(cwdbuf,BUFSIZE); dda_set_output_directory("output",sizeof("output")); ! strncpy(outdir,cwdbuf,sizeof(cwdbuf)); strncat(outdir,"\\",sizeof("\\")); strncat(outdir,"output",sizeof("output")); --- NEW FILE: mohrcoulomb.c --- #include <stdio.h> #include <stdlib.h> #include <math.h> #include "mohrcoulomb.h" /* sigm: mean stress invariant /* dsbar: deviatoric stress invariant /* fnew or f: value of failure function /* phi: friction angle /* c: cohesion /* psi: dilation angle /* tnph: tan(phi) /* snph: sin(phi) /* theta: lode angle */ double sigm,dsbar,theta,cs[5],pl[5][5]; int i,j; /* /* first invariant "sigm", second invariant "dsbar" /* and lode angle "theta" for a stress tensor "stress" */ void invar() { double sx,sy,sz,txy,dx,dy,dz,xj3,sine; sx = cs[1]; sy = cs[2]; txy = cs[3]; sz = cs[4]; sigm = (sx+sy+sz)/3; dsbar = sqrt(pow((sx-sy),2)+pow((sy-sz),2)+pow((sz-sx),2)+6.*pow(txy,2))/sqrt(2.); if(dsbar < 1.e-10) { theta = .0; } else { dx = (2*sx-sy-sz)/3; dy = (2*sy-sz-sx)/3; dz = (2*sz-sx-sy)/3; xj3 = dx*dy*dz-dz*pow(txy,2); sine = -13.5*xj3/pow(dsbar,3); if(sine > 1) sine = 1; if(sine < -1) sine = -1; theta = asin(sine)/3; } } /* /* Mohr-Coulomb yield function F from invariants "sigm" /* and "dsbar" lode angle "theta (radian) of friction "phi" /* and cohesion "c" */ double mocouf(double c,double phi) { double phir,snph,csph,csth,snth,f; phir = phi*4*atan(1)/180; snph = sin(phir); csph = cos(phir); csth = cos(theta); snth = sin(theta); f = snph*sigm+dsbar*(csth/sqrt(3.)-snth*snph/3.)-c*csph; return f; } /* /* plastic stress-strain matrix "pl" from stresses "stress" /* angle of friction "phi", dilation "psi" and elastic stiffness "E" /* ,"nu" for a Mohr-Coulomb material */ void mocopl(double phi,double psi,double nu, double ymod) { double row[5],col[5],sx,sy,txy,sz,pi,phir,psir; double dx,dy,dz,d2,d3,th,snth,sig,rph,rps,cps,snps,sq3; double cc,cph,alp,ca,sa,dd,snph,ee,s1,s2; sx = cs[1]; sy = cs[2]; txy = cs[3]; sz = cs[4]; pi = 4*atan(1); phir = phi*pi/180; psir = psi*pi/180; snph = sin(phir); snps = sin(psir); sq3 = sqrt(3); cc = 1-2*nu; dx = (2*sx-sy-sz)/3; dy = (2*sy-sz-sx)/3; dz = (2*sz-sx-sy)/3; d2 = sqrt(-dx*dy-dy*dz-dz*dx+txy*txy); d3 = dx*dy*dz-dz*txy*txy; th = -3*sq3*d3/(2*pow(d2,3)); if(th > 1)th = 1; if(th < -1)th = -1; th = asin(th)/3; snth = sin(th); if(fabs(snth) > 0.49) { sig = -1; if(snth < 0)sig = 1; rph = snph*(1+nu)/3; rps = snps*(1.+nu)/3; cps = 0.25*sq3/d2*(1+sig*snps/3); cph = 0.25*sq3/d2*(1+sig*snph/3); col[1] = rph+cph*((1-nu)*dx+nu*(dy+dz)); col[2] = rph+cph*((1-nu)*dy+nu*(dz+dx)); col[3] = cph*cc*txy; col[4] = rph+cph*((1-nu)*dz+nu*(dx+dy)); row[1] = rps+cps*((1-nu)*dx+nu*(dy+dz)); row[2] = rps+cps*((1-nu)*dy+nu*(dz+dx)); row[3] = cps*cc*txy; row[4] = rps+cps*((1-nu)*dz+nu*(dx+dy)); ee = ymod/((1+nu)*cc*(rph*snps+2*cph*cps*d2*d2*cc)); } else { alp = atan(fabs((sx-sy)/(2*txy))); ca = cos(alp); sa = sin(alp); dd = cc*sa; s1=1; s2=1; if((sx-sy) < 0)s1 = -1; if(txy < 0)s2 = -1; col[1] = snph+s1*dd; col[2] = snph-s1*dd; col[3] = s2*cc*ca; col[4] = 2*nu*snph; row[1] = snps+s1*dd; row[2] = snps-s1*dd; row[3] = s2*cc*ca; row[4] = 2*nu*snps; ee = ymod/(2*(1.+nu)*cc*(snph*snps+cc)); } for(i=1; i<=4; i++) { for(j=1; j<=4; j++) { pl[i][j] = ee*row[i]*col[j]; } } } /* /* main function in order to obtain the /* correct stress using Mohr-Coulomb criteria */ void mohrcoulomb(double c,double phi,double psi,double nu,double ymod, double eps[4],double * stress,double sigma[4]) { double elso[5]; double fnew; for(i = 1; i <= 4; i++) elso[i] = 0; for(i = 1; i <= 4; i++) cs[i] = stress[i+3] + sigma[i]; invar(); fnew = mocouf(c,phi); if(fnew > 0) { mocopl(phi,psi,nu,ymod); elso[1] = pl[1][1] * eps[1] + pl[1][2] * eps[2] + pl[1][3] * eps[3]; elso[2] = pl[2][1] * eps[1] + pl[2][2] * eps[2] + pl[2][3] * eps[3]; elso[3] = pl[3][1] * eps[1] + pl[3][2] * eps[2] + pl[3][3] * eps[3]; elso[4] = pl[4][1] * eps[1] + pl[4][2] * eps[2] + pl[4][3] * eps[3]; } for(i = 1; i <= 4; i++) stress[i+3] = stress[i+3] + sigma[i] - elso[i]; } |
|
From: David M. D. <do...@us...> - 2006-07-02 23:07:09
|
Update of /cvsroot/dda/ntdda/include In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv28889 Added Files: mohrcoulomb.h Log Message: Initial mohrcoulomb file from roozbeh. --- NEW FILE: mohrcoulomb.h --- /* * Mohr-Coulomb criteria */ #ifndef __MOHRCOULOMB_H__ #define __MOHRCOULOMB_H__ #ifdef __cplusplus extern "C" { #endif #if 0 } #endif void invar(); double mocouf(double c,double phi); void mocopl(double phi,double psi,double nu, double ymod); void mohrcoulomb(double c, double phi, double psi,double nu,double ymod, double eps[],double * stress,double sigma[]); #endif /* __MOHRCOULOMB_H__ */ |
|
From: David M. D. <do...@us...> - 2006-07-01 17:48:28
|
Update of /cvsroot/dda/ntdda/src/win32gui In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv8627 Modified Files: winmain.c Log Message: Cleaned up title bar handling for main window using a function call with arguments for setting text. Index: winmain.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/win32gui/winmain.c,v retrieving revision 1.43 retrieving revision 1.44 diff -C2 -d -r1.43 -r1.44 *** winmain.c 1 Jul 2006 14:58:07 -0000 1.43 --- winmain.c 1 Jul 2006 15:20:58 -0000 1.44 *************** *** 57,67 **** /* Whole bunch of global variables that need to disappear fast. */ ! char szAppName[] = "Berkeley DDA for Windows"; char mess[80]; /* mess appears to be a temporary variable. Might be able to get rid of it. */ // Get rid of this as well. - char mainWinTitle[120]; - #define ABOUT "UC Berkeley DDA for Windows 95/NT,\n", \ --- 57,65 ---- /* Whole bunch of global variables that need to disappear fast. */ ! char szAppName[] = "DDA for Windows"; char mess[80]; /* mess appears to be a temporary variable. Might be able to get rid of it. */ // Get rid of this as well. #define ABOUT "UC Berkeley DDA for Windows 95/NT,\n", \ *************** *** 144,147 **** --- 142,154 ---- + void + set_mainwindow_titlebar(HWND hwMain, char * filetype, char * filename) { + + char mainWinTitle[120]; + + sprintf(mainWinTitle, "%s --- %s: %s",szAppName,filetype,filepath.gfile); + SetWindowText(hwMain, (LPCTSTR) mainWinTitle); + } + void *************** *** 625,630 **** whatToDraw = NOTHING; filepath.gfile[0] = '\0'; ! sprintf(mainWinTitle, "%s for Windows 95/NT", (LPSTR) szAppName); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); dda_set_menu_state(dda,GEOM_STATE | ABORTED); } --- 632,639 ---- whatToDraw = NOTHING; filepath.gfile[0] = '\0'; ! //sprintf(mainWinTitle, "%s for Windows 95/NT", (LPSTR) szAppName); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! set_mainwindow_titlebar(hwMain,NULL,NULL); ! dda_set_menu_state(dda,GEOM_STATE | ABORTED); } *************** *** 645,650 **** strcpy(filepath.oldfile, filepath.gfile); filepath.gfile[0] = '\0'; ! sprintf(mainWinTitle, "%s for Windows 95/NT", (LPSTR) szAppName); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); hInst = (HINSTANCE) GetWindowLong(hwMain, GWL_HINSTANCE); --- 654,660 ---- strcpy(filepath.oldfile, filepath.gfile); filepath.gfile[0] = '\0'; ! //sprintf(mainWinTitle, "%s for Windows 95/NT", (LPSTR) szAppName); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! set_mainwindow_titlebar(hwMain,NULL,NULL); hInst = (HINSTANCE) GetWindowLong(hwMain, GWL_HINSTANCE); *************** *** 655,660 **** if(DialogBoxParam(hInst, "DRAWDLG", hwMain, (DLGPROC)DrawDlgProc, (LPARAM)12 ) ) { ! sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", (LPSTR) szAppName, filepath.gfile); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); } --- 665,672 ---- if(DialogBoxParam(hInst, "DRAWDLG", hwMain, (DLGPROC)DrawDlgProc, (LPARAM)12 ) ) { ! //sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", (LPSTR) szAppName, filepath.gfile); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! set_mainwindow_titlebar(hwMain, "Geometry",filepath.gfile); ! } *************** *** 680,685 **** if(filepath.gfile[0] != '\0') { ! sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", (LPSTR) szAppName, filepath.gfile); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); if (strcmp(filepath.gfile, filepath.oldfile) != 0 && MessageBox(hwMain, "Apply geometry?", "GEOMETRY", MB_YESNO) == IDYES) { --- 692,699 ---- if(filepath.gfile[0] != '\0') { ! //sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", (LPSTR) szAppName, filepath.gfile); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! set_mainwindow_titlebar(hwMain, "Geometry",filepath.gfile); ! if (strcmp(filepath.gfile, filepath.oldfile) != 0 && MessageBox(hwMain, "Apply geometry?", "GEOMETRY", MB_YESNO) == IDYES) { *************** *** 707,713 **** else { ! sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", (LPSTR) szAppName, (LPSTR) filepath.gfile); ! ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); dda_set_menu_state(dda,GEOM_STATE); --- 721,727 ---- else { ! //sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", (LPSTR) szAppName, (LPSTR) filepath.gfile); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! set_mainwindow_titlebar(hwMain, "Geometry",filepath.gfile); dda_set_menu_state(dda,GEOM_STATE); *************** *** 805,813 **** else { ! sprintf(mainWinTitle, ! "%s for Windows 95/NT --- Geometry = %s, Analysis File = %s", ! (LPSTR) szAppName, (LPSTR) filepath.gfile, (LPSTR) filepath.afile); ! ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); //updateMainMenu(hwMain, anastate*readystate); //dda->menustate = ANA_STATE | READY_STATE; --- 819,830 ---- else { ! ! // TODO: Find a way to handle both Geometry and Analysis files ! // with a titlebar call. ! //sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s, Analysis File = %s", (LPSTR) szAppName, (LPSTR) filepath.gfile, (LPSTR) filepath.afile); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! ! set_mainwindow_titlebar(hwMain, "Analysis",filepath.afile); ! //updateMainMenu(hwMain, anastate*readystate); //dda->menustate = ANA_STATE | READY_STATE; *************** *** 926,936 **** Geometrydata * geomdata = dda_get_geometrydata(dda); strcpy(filepath.oldfile, filepath.afile); ! filepath.afile[0] = '\0'; ! sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", ! (LPSTR) szAppName, filepath.gfile); ! ! hInst = (HINSTANCE) GetWindowLong(hwMain, GWL_HINSTANCE); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); if (!DialogBoxParam(hInst, "ANALDLG2", hwMain, (DLGPROC)AnalDlgProc,(LPARAM)geomdata )) --- 943,952 ---- Geometrydata * geomdata = dda_get_geometrydata(dda); strcpy(filepath.oldfile, filepath.afile); ! filepath.afile[0] = '\0'; ! hInst = (HINSTANCE) GetWindowLong(hwMain, GWL_HINSTANCE); ! //sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", LPSTR) szAppName, filepath.gfile); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! set_mainwindow_titlebar(hwMain, "Geometry",filepath.gfile); if (!DialogBoxParam(hInst, "ANALDLG2", hwMain, (DLGPROC)AnalDlgProc,(LPARAM)geomdata )) *************** *** 939,949 **** } ! if(filepath.afile[0] != '\0') ! { ! sprintf(mainWinTitle, ! "%s for Windows 95/NT --- Geometry = %s, Analysis File = %s", ! (LPSTR) szAppName, filepath.gfile, filepath.afile); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); if (MessageBox(hwMain, "Run analysis?", "ANALYSIS", MB_YESNO) == IDYES) --- 955,963 ---- } ! if(filepath.afile[0] != '\0') { ! //sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s, Analysis File = %s", (LPSTR) szAppName, filepath.gfile, filepath.afile); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! set_mainwindow_titlebar(hwMain, "Analysis",filepath.afile); if (MessageBox(hwMain, "Run analysis?", "ANALYSIS", MB_YESNO) == IDYES) *************** *** 974,982 **** //iface->setdisplay((unsigned int)hwMain); ! if(filepath.afile[0] != '\0') ! { ! sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s, Analysis File = %s", ! (LPSTR) szAppName, filepath.gfile, filepath.afile); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); if (MessageBox(hwMain, "Run analysis?", "ANALYSIS", MB_YESNO) == IDYES) { --- 988,997 ---- //iface->setdisplay((unsigned int)hwMain); ! if(filepath.afile[0] != '\0') { ! ! //sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s, Analysis File = %s", LPSTR) szAppName, filepath.gfile, filepath.afile); ! //SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! set_mainwindow_titlebar(hwMain, "Analysis",filepath.afile); ! if (MessageBox(hwMain, "Run analysis?", "ANALYSIS", MB_YESNO) == IDYES) { *************** *** 1713,1717 **** - // This subroutine is wrote by Roozbeh to read the .dxf // format and to create the .geo format --- 1728,1731 ---- *************** *** 1724,1728 **** DDA * dda = (DDA *)GetWindowLong(hwMain,GWL_USERDATA); ! LPCTSTR szFilter[] = {"Autocad Dxf files (*.dxf)\0*.dxf\0All files (*.*)\0*.*\0\0"}; fileBrowse(hwMain, &ofn, szFilter, filepath.gpath, filepath.gfile, "dxf"); --- 1738,1744 ---- DDA * dda = (DDA *)GetWindowLong(hwMain,GWL_USERDATA); ! LPCTSTR szFilter[] = {"Autocad Dxf files (*.dxf)\0*.dxf\0All files (*.*)\0*.*\0\0"}; ! // TODO: Add a dxf file type for this, instead of shoving into the geometry file. ! // That way, the dxf file can be modified independently at a future date. fileBrowse(hwMain, &ofn, szFilter, filepath.gpath, filepath.gfile, "dxf"); *************** *** 1731,1737 **** return 0; /* user pressed cancel */ } ! ! sprintf(mainWinTitle, "%s for Windows --- Geometry = %s", (LPSTR) szAppName, (LPSTR) filepath.gfile); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); dda_set_menu_state(dda,GEOM_STATE); --- 1747,1752 ---- return 0; /* user pressed cancel */ } ! ! set_mainwindow_titlebar(hwMain,"DXF",filepath.gfile); dda_set_menu_state(dda,GEOM_STATE); *************** *** 2112,2117 **** hddawin = CreateWindowEx(WS_EX_ACCEPTFILES, ! szAppName, ! "Berkeley DDA for Windows NT", WS_OVERLAPPEDWINDOW | WS_CLIPCHILDREN, --- 2127,2132 ---- hddawin = CreateWindowEx(WS_EX_ACCEPTFILES, ! szAppName, ! szAppName, WS_OVERLAPPEDWINDOW | WS_CLIPCHILDREN, |
|
From: David M. D. <do...@us...> - 2006-07-01 14:58:10
|
Update of /cvsroot/dda/ntdda/src/win32gui In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv32509 Modified Files: winmain.c Log Message: Geometry file is now automatically loaded after being read from dxf import. Index: winmain.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/win32gui/winmain.c,v retrieving revision 1.42 retrieving revision 1.43 diff -C2 -d -r1.42 -r1.43 *** winmain.c 1 Jul 2006 14:51:25 -0000 1.42 --- winmain.c 1 Jul 2006 14:58:07 -0000 1.43 *************** *** 1716,1720 **** // This subroutine is wrote by Roozbeh to read the .dxf // format and to create the .geo format ! static int handleDxfBrowse(HWND hwMain, LPARAM lParam) //Added by Roozbeh { --- 1716,1720 ---- // This subroutine is wrote by Roozbeh to read the .dxf // format and to create the .geo format ! int handleDxfBrowse(HWND hwMain, LPARAM lParam) //Added by Roozbeh { *************** *** 1727,1764 **** fileBrowse(hwMain, &ofn, szFilter, filepath.gpath, filepath.gfile, "dxf"); ! if( !GetOpenFileName(&ofn) ) { ! strcpy(filepath.gpath, filepath.oldpath); return 0; /* user pressed cancel */ - } ! ! //else { ! ! ! sprintf(mainWinTitle, "%s for Windows --- Geometry = %s", (LPSTR) szAppName, (LPSTR) filepath.gfile); ! ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! ! dda_set_menu_state(dda,GEOM_STATE); ! /* Initialize the path. Note that the .dxf file is already loaded into ! * the struct, so it doesn't have to be reloaded. ! */ ! strcpy(temp, filepath.gfile); ! strcpy (filepath.rootname, strtok(temp, ".")); ! // dda_set_output_directory("output",sizeof("output")); //Disabled By Roozbeh ! // The dxf file will be read from this line and lines, polylines ! // and arcs will be recognized from dxf file ! fp1 = fopen(filepath.gfile,"r"); ! strcpy (filepath.gfile, strcat (filepath.rootname, ".geo")); ! // Pass everything in to readDXF. ! dxf_read_file(fp1,filepath.gfile); ! // Send message to GeomApply to have geometry automatically loaded. ! ! //} return 1; --- 1727,1756 ---- fileBrowse(hwMain, &ofn, szFilter, filepath.gpath, filepath.gfile, "dxf"); ! if(!GetOpenFileName(&ofn)) { strcpy(filepath.gpath, filepath.oldpath); return 0; /* user pressed cancel */ } ! ! sprintf(mainWinTitle, "%s for Windows --- Geometry = %s", (LPSTR) szAppName, (LPSTR) filepath.gfile); ! SetWindowText(hwMain, (LPCTSTR) mainWinTitle); ! dda_set_menu_state(dda,GEOM_STATE); ! /* Initialize the path. Note that the .dxf file is already loaded into ! * the struct, so it doesn't have to be reloaded. ! */ ! strcpy(temp, filepath.gfile); ! strcpy (filepath.rootname, strtok(temp, ".")); ! //dda_set_output_directory("output",sizeof("output")); //Disabled By Roozbeh ! // The dxf file will be read from this line and lines, polylines ! // and arcs will be recognized from dxf file ! fp1 = fopen(filepath.gfile,"r"); ! // WARNING: This overwrites the name of the currently open file fp1! ! strcpy (filepath.gfile, strcat (filepath.rootname, ".geo")); ! dxf_read_file(fp1,filepath.gfile); ! // Send message to GeomApply to have geometry automatically loaded. ! SendMessage(hwMain, WM_COMMAND, GEOM_APPLY, lParam); return 1; |
|
From: David M. D. <do...@us...> - 2006-07-01 14:51:32
|
Update of /cvsroot/dda/ntdda/src/win32gui In directory sc8-pr-cvs4.sourceforge.net:/tmp/cvs-serv29802 Modified Files: winmain.c Log Message: Deleted #if'ed out code in winmain for handling dxf which is now in it's own file. Index: winmain.c =================================================================== RCS file: /cvsroot/dda/ntdda/src/win32gui/winmain.c,v retrieving revision 1.41 retrieving revision 1.42 diff -C2 -d -r1.41 -r1.42 *** winmain.c 1 Jul 2006 05:25:26 -0000 1.41 --- winmain.c 1 Jul 2006 14:51:25 -0000 1.42 *************** *** 59,63 **** char szAppName[] = "Berkeley DDA for Windows"; char mess[80]; /* mess appears to be a temporary variable. Might be able to get rid of it. */ ! char mainWinTitle[120]; --- 59,65 ---- char szAppName[] = "Berkeley DDA for Windows"; char mess[80]; /* mess appears to be a temporary variable. Might be able to get rid of it. */ ! ! ! // Get rid of this as well. char mainWinTitle[120]; *************** *** 1710,2019 **** - #if 0 - void - readDXF(FILE * fp1) { - - double *jx1,*jy1,*jx2,*jy2; - char str[5000][50]; - FILE *fp2; - int count1=0,count2=0,i,j; - long *type,pn; - int cabcou=0,dec,ii,pltype,ctype,n1; - int fixn=0,measn=0,holen=0,loadn=0; - double r0,se,u,v,x1,x2,x3,y1,y2,y3,x0,y0,deg; - double fx[25],fy[25],lx[25],ly[25]; - double mx[25],my[25],hx[25],hy[25]; - int nline=0,npoly=0,narc=0,ncir=0,ntext=0; - //OPENFILENAME ofn; - //char temp[200]; - - - while(!feof(fp1)) { - count1 +=1; - fscanf(fp1,"%s",str[count1]); - - if(strcmp(str[count1],"LINE")==0) - nline+=1; - - if(strcmp(str[count1],"LWPOLYLINE")==0) - npoly+=1; - - if(strcmp(str[count1],"ARC")==0) - narc+=1; - - if(strcmp(str[count1],"CIRCLE")==0) - ncir+=1; - } - - // Array Allocation - n1=nline+40*npoly+40*narc+40*ncir+1; - jx1 =(double *)calloc(n1,sizeof(double)); - jy1 =(double *)calloc(n1,sizeof(double)); - jx2 =(double *)calloc(n1,sizeof(double)); - jy2 =(double *)calloc(n1,sizeof(double)); - type =(long *)calloc(n1,sizeof(long)); - - for(i=1;i<=count1;i++) { - - // The dxf file will be read from this line and lines, polylines - // and arcs will be recognized from dxf file - - if(strcmp(str[i],"LINE")==0) { - - count2 +=1; - if(strcmp(str[i+9],"62")==0) { - - type[count2]=atoi(str[i+10])+1; - dec=2; - } else { - type[count2]=1; - dec=0; - } - - jx1[count2]=atof(str[i+12+dec]); - jy1[count2]=atof(str[i+14+dec]); - jx2[count2]=atof(str[i+18+dec]); - jy2[count2]=atof(str[i+20+dec]); - - } - - if(strcmp(str[i],"LWPOLYLINE")==0) { - - if(strcmp(str[i+9],"62")==0) { - - pltype=atoi(str[i+10])+1; - dec=2; - } else { - - pltype=1; - dec=0; - } - pn=atoi(str[i+12+dec]); - for(j=0;j<pn-1;j++) - { - ii=4*j; - count2 +=1; - jx1[count2]=atof(str[i+18+ii+dec]); - jy1[count2]=atof(str[i+18+ii+2+dec]); - jx2[count2]=atof(str[i+18+ii+4+dec]); - jy2[count2]=atof(str[i+18+ii+6+dec]); - type[count2]=pltype; - } - count2 +=1; - if(strcmp(str[i+9],"62")==0) - { - type[count2]=atoi(str[i+10])+1; - dec=2; - } - else - { - type[count2]=1; - dec=0; - } - ii=4*3; - jx1[count2]=atof(str[i+18+ii+dec]); - jy1[count2]=atof(str[i+18+ii+2+dec]); - jx2[count2]=atof(str[i+18+dec]); - jy2[count2]=atof(str[i+20+dec]); - type[count2]=pltype; - } - // From this line first, center and second points of an arc will be extracted from dxf - // and the lines of arc are created with a simple mathematical source - - if(strcmp(str[i],"ARC")==0) - { - if(strcmp(str[i+9],"62")==0) - { - ctype=atoi(str[i+10])+1; - dec=2; - } - else - { - ctype=1; - dec=0; - } - se=30; - deg=atof(str[i+24+dec]); - r0=((3.1415926535)*deg)/(se*180); - x1=atof(str[i+12+dec])+atof(str[i+18+dec]); - y1=atof(str[i+14+dec]); - x2=atof(str[i+12+dec]); - y2=atof(str[i+14+dec]); - x0=x2; - y0=y2; - x3=atof(str[i+12+dec])-atof(str[i+18+dec]); - y3=atof(str[i+14+dec]); - - for(j=1;j<se;j++) - { - u=((x1-x0)*(cos(r0)-1))-((y1-y0)*sin(r0)); - v=((x1-x0)*sin(r0))+((y1-y0)*(cos(r0)-1)); - - x2=x1+u; - y2=y1+v; - count2 +=1; - jx1[count2]=x1; - jy1[count2]=y1; - jx2[count2]=x2; - jy2[count2]=y2; - type[count2]=ctype; - x1=x2; - y1=y2; - - } - - x2=x3; - y2=y3; - count2 +=1; - jx1[count2]=x1; - jy1[count2]=y1; - jx2[count2]=x2; - jy2[count2]=y2; - type[count2]=ctype; - } - // From this line first, center and second points of an circle will be extracted from dxf - // and the lines of arc are created with a simple mathematical source - if(strcmp(str[i],"CIRCLE")==0) - { - if(strcmp(str[i+9],"62")==0) - { - ctype=atoi(str[i+10])+1; - dec=2; - } - else - { - ctype=1; - dec=0; - } - se=30; - deg=360; - r0=((3.1415926535)*deg)/(se*180); - x1=atof(str[i+12+dec])+atof(str[i+18+dec]); - y1=atof(str[i+14+dec]); - x2=atof(str[i+12+dec]); - y2=atof(str[i+14+dec]); - x0=x2; - y0=y2; - x3=x1; - y3=y1; - - for(j=1;j<se;j++) - { - u=((x1-x0)*(cos(r0)-1))-((y1-y0)*sin(r0)); - v=((x1-x0)*sin(r0))+((y1-y0)*(cos(r0)-1)); - - x2=x1+u; - y2=y1+v; - count2 +=1; - jx1[count2]=x1; - jy1[count2]=y1; - jx2[count2]=x2; - jy2[count2]=y2; - type[count2]=ctype; - x1=x2; - y1=y2; - - } - - x2=x3; - y2=y3; - count2 +=1; - jx1[count2]=x1; - jy1[count2]=y1; - jx2[count2]=x2; - jy2[count2]=y2; - type[count2]=ctype; - } - if(strcmp(str[i],"MTEXT")==0) - { - if(strcmp(str[i+9],"62")==0) - { - dec=2; - } - else - { - dec=0; - } - if(strcmp(str[i+26+dec],"F")==0) - { - fixn +=1; - fx[fixn]=atof(str[i+12+dec]); - fy[fixn]=atof(str[i+14+dec]); - } - if(strcmp(str[i+26+dec],"L")==0) - { - loadn +=1; - lx[loadn]=atof(str[i+12+dec]); - ly[loadn]=atof(str[i+14+dec]); - } - if(strcmp(str[i+26+dec],"M")==0) - { - measn +=1; - mx[measn]=atof(str[i+12+dec]); - my[measn]=atof(str[i+14+dec]); - } - if(strcmp(str[i+26+dec],"H")==0) - { - holen +=1; - hx[holen]=atof(str[i+12+dec]); - hy[holen]=atof(str[i+14+dec]); - } - } - } - fclose(fp1); - // From this line the data that extracted from dxf file will be sorted - // and will be wrote into geo file. - strcpy (filepath.gfile, strcat(filepath.rootname, ".geo")); - fp2 = fopen(filepath.gfile,"w"); - - fprintf(fp2,"<?xml version=\"1.0\" standalone=\"no\"?>\n"); - fprintf(fp2,"<!DOCTYPE DDA SYSTEM \"geometry.dtd\">\n"); - fprintf(fp2,"<Berkeley:DDA xmlns:Berkeley=\"http://www.tsoft.com/~bdoolin/dda\">\n"); - fprintf(fp2,"<!-- Bogus comment to keep ddaml tree-stripping\n"); - fprintf(fp2,"from seg faulting on bad child node. -->\n\n"); - fprintf(fp2,"<Geometry>\n"); - fprintf(fp2," <Edgenodedist distance=\"%lf\"/>\n",0.001); - fprintf(fp2," <Jointlist>\n"); - for(i=1;i<=count2;i++) - { - fprintf(fp2," <Joint type=\"%d\">",type[i]); - fprintf(fp2," %lf %lf %lf %lf",jx1[i],jy1[i],jx2[i],jy2[i]); - fprintf(fp2," </Joint>\n"); - } - fprintf(fp2,"</Jointlist>\n"); - - fprintf(fp2,"<Fixedpointlist>\n"); - for(i=1;i<=fixn;i++)fprintf(fp2,"<Line> %lf %lf %lf %lf </Line>\n",fx[i],fy[i],fx[i],fy[i]); - fprintf(fp2,"</Fixedpointlist>\n"); - - fprintf(fp2,"<Loadpointlist>\n"); - for(i=1;i<=loadn;i++)fprintf(fp2,"<Point> %lf %lf </Point>\n",lx[i],ly[i]); - fprintf(fp2,"</Loadpointlist>\n"); - - fprintf(fp2,"<Measuredpointlist>\n"); - for(i=1;i<=measn;i++)fprintf(fp2,"<Point> %lf %lf </Point>\n",mx[i],my[i]); - fprintf(fp2,"</Measuredpointlist>\n"); - - fprintf(fp2,"<Holepointlist>\n"); - for(i=1;i<=holen;i++)fprintf(fp2,"<Point> %lf %lf </Point>\n",hx[i],hy[i]); - fprintf(fp2,"</Holepointlist>\n"); - - fprintf(fp2,"</Geometry>\n"); - fprintf(fp2,"</Berkeley:DDA>\n"); - - fclose(fp2); - - - - free(jx1); - free(jy1); - free(jx2); - free(jy2); - free(type); - - return; - - } - #endif --- 1712,1715 ---- *************** *** 2023,2036 **** handleDxfBrowse(HWND hwMain, LPARAM lParam) //Added by Roozbeh { ! FILE *fp1;//,*fp2; ! //int count1=0,count2=0,i,j; ! //char str[5000][50]; ! //long *type,pn; ! //int cabcou=0,dec,ii,pltype,ctype,n1; ! int fixn=0,measn=0,holen=0,loadn=0; ! //double r0,se,u,v,x1,x2,x3,y1,y2,y3,x0,y0,deg; ! //double fx[25],fy[25],lx[25],ly[25]; ! //double mx[25],my[25],hx[25],hy[25]; ! int nline=0,npoly=0,narc=0,ncir=0,ntext=0; OPENFILENAME ofn; char temp[200]; --- 1719,1723 ---- handleDxfBrowse(HWND hwMain, LPARAM lParam) //Added by Roozbeh { ! FILE *fp1; OPENFILENAME ofn; char temp[200]; *************** *** 2045,2052 **** return 0; /* user pressed cancel */ ! } else { ! sprintf(mainWinTitle, "%s for Windows 95/NT --- Geometry = %s", (LPSTR) szAppName, (LPSTR) filepath.gfile); SetWindowText(hwMain, (LPCTSTR) mainWinTitle); --- 1732,1741 ---- return 0; /* user pressed cancel */ ! } ! ! //else { ! sprintf(mainWinTitle, "%s for Windows --- Geometry = %s", (LPSTR) szAppName, (LPSTR) filepath.gfile); SetWindowText(hwMain, (LPCTSTR) mainWinTitle); *************** *** 2068,2073 **** // Pass everything in to readDXF. dxf_read_file(fp1,filepath.gfile); ! ! } return 1; --- 1757,1764 ---- // Pass everything in to readDXF. dxf_read_file(fp1,filepath.gfile); ! ! // Send message to GeomApply to have geometry automatically loaded. ! ! //} return 1; |
