Software: FLAC3D

Updates

4.00 32

4.00 64

3.10 32

3.10 64

FLAC3D 4.00.82

Download Update

Current Version 32-bit
Date: 1/06/2012
Size: 10.3 MB

These updates are to be used only by licensed FLAC3D users. They are not "trial" versions of the code, nor are they fully functional installations — before downloading the update you must install FLAC3D from the Itasca software installation CD-ROM.

FLAC3D 4.0 (Current Version)

Prerequisites Update Instructions Supplementary Files
Revision History

Starting with subversion 68, an update to the Microsoft (R) Visual C++ 2005 Runtime is REQUIRED. (The current subversion of the code update available here appears at the end of the Revision History — see link above.) If not present, on startup the code will return an error message similar to "The application has failed to start because the application configuration is incorrect. Reinstalling the application may fix this problem" (on Windows XP or later) or it will complain of missing DLLs (Windows 2000).

Please download and install the runtime that corresponds to your 32-bit or 64-bit operating system. The installation process need only be performed once; subsequent updates to FLAC3D will not require these steps to be repeated (although future updates to the runtime released by Microsoft may need to be installed).

Version	Description	Date
82	1/6/12	If a shell is created on a zone face, and an opposing zone face exists, the target zone for the shell's node-zone link was arbitrarily chosen. Target zones are now chosen based on the zone range specified when shells are created. E.g., "sel shell id 200 range z 1 tol 0.1 group granite" would create shells that are linked to granite zone faces.
81	12/2/11	Fix memory corruption if you CONFIG ZGROUP from one >1 value to another >1 value.
80	11/3/11	Default tolerance on orientation range set to 2 degrees as described in manuals. Added group slot support to IMPGRID and EXPGRID.
79	10/25/11	Fixed crash that occasionally occurred during plotting. Zero stiffness in gridpoint error might be encountered on cycling if a GEN SEPARATE or GEN MERGE command was executed and a MODEL command is not encountered (after GEN SEPARATE or MERGE) prior to cycling. This has been fixed. Generate error if non-null zone encountered without overlays during cycling.
78	10/14/11	Added checking to the command processor to make certain that a FISH symbol on the command line is actually a function and takes the number of arguments provided in tokens following it. Structural elements no longer obtain their default properties from structures with the lowest id number. Structural element properties, if not initialized are zero.
77	9/29/11	Hide null fluid and thermal plot item attributes did not work. FISHHALT function was being called twice per cycle. It is now called once at the end of a cycle. Added optional timeout period to SYSTEM command. By default, system commands timeout after 30 seconds, after which control is returned to FLAC3D. An optional interger time, in seconds, can be added to system commands to change this timeout period. E.g., SYSTEM "notepad.exe" 120 would time out after 2 minutes.
76	9/20/11	Added PLOT SHOW {VIEW viewname} command. Regenerates (refreshes) all plot items in the specified view and makes it the current view. If view is not specified PLOT SHOW refreshes the current view. Fix picking and selection bug in History Location plot.
75	8/30/11	Yield state of liner was plotted incorrectly. Also, sl_ryield FISH function did not return the correct yield state.
74	8/24/11	Duplicate zones were reported in list of zones connected to gridpoint (when outputting gridpoint information as well as in FISH function gp_linkzone list). This only occurred with non-hexahedral zones. FISH function sl_ryield(pls,ni) returns yield state of shear coupling spring at node ni={1,2,3}. Changed from ni={0,1,2} to make function consistent with other functions. FISH functions sl_rdisp(pls,ns,ni), sl_rstr(pls,ns,ni), sl_rstrdir(pls,ns,ni,d), sl_vrstrdir(pls,ns,ni,d) did not have a complete parameter list in the documentation. Parameter ns={1,2} refers to the normal and shear coupling spring respectively. Parameter ni={1,2,3} is the node index.
73	8/18/11	Bug fix, HISTORY SEL NODE command generated node histories not corresponding to specified keyword. Bug fix, FISH dydt returned incorrect value when structural elements were present. Bug fix to ZGroup plot item (coordinates of gridpoints and zones were not correct in information window). Plot Query mode information boxes will now scroll vertically, and a Copy To Clipboard button has been added.
72	6/24/11	Plots did not save/restore model range elements. Bug fix to SOLVE FOS when only interfaces are included, FLAC3D did not detect unstable state properly.
71	6/20/11	Typographic correction to plot commands with keyword containing the word minimum (e.g., pl zcontour sminimum). Minimum was improperly spelled as mininum. Bug fix to parsing SOLVE FOS command (keywords INCLUDE and EXCLUDE were sometimes treated as material property strings). Added a check for floating point status during cycling (checks for overflows, divide by zero, etc. are made). Further bug fix to edge attach, FLAC3D would generate an error message while trying to create edge attach conditions even though they were valid. Bug fix to SEL geometry plot, markers and system type options were not being copied, saved, restored.
70	5/24/11	Bug fix to grid generator. Newly generated zones with triangle faces were not always joined to the existing grid.
69	5/19/11	Bug fix to group logic, preventing an occasional crash on new command. Bug fix to edge attach, FLAC3D would generate an error message while trying to create edge attach conditions even though they were valid. Improved interface element area distribution calculation for the special case when two elements fully cover a quad zone face (more uniform stress distribution results when quad faces are irregularly shaped). Tetrahedral zones were not joined to other tetrahedral zones when a GEN ZONE REFLECT or COPY command was issued. This resulted in incorrect boundary plots as well as interfaces being created on shared triangle zone faces when an interface command was issued involving the shared faces.
68	5/6/11	Bug fix to fluid flow. On multi-core machines FLAC3D would output an error message that all fluid zones are null when they were not. Updated IMPGRID to work with unicode strings and files. Bug fix to gp_near FISH function. Bug fix to zgroup, fixed crash when config zgroup 1 was done after config zgroup with a larger number. Fixed crash if history of crtime, fltime, dytime, thtime are taken and FLAC3D is not configured for that option. All structural element links were not appearing in SEL link marker plot. Relaxed tolerance for stiffness matrix consistency check for line-type structural elements. Added readonly zone FISH functions z_nummechprops(zptr) and z_mechpropname(zptr, i). z_nummechprops(zptr) returns the number of mechanical properties a zone's constitutive model contains. z_mechpropname(zptr, i) returns a string corresponding to the zone's ith property name (values of i range from 1 to the value returned by FISH function z_nummechprops). Fixed crash if optional array argument was not provided in FISH function z_sonplane.
67	2/8/11	Fixed bug in cell space. Changed zone search implementation, fixing bug when finding zones located at points in space.
66	2/1/11	Added pretension command for piles (similar to that available for cables). Usage is "sel pile pretension value range". It is assumed that the pile segment that gets pretensioned is in an ungrouted portion of the pile and that only axial forces have developed in the segment. This command only updates axial forces in the pile, all others are assumed to be zero. Added readonly FISH function m=gp_gravmass(p) which returns gravitational mass stored at a gridpoint.
65	1/25/11	Changed 60 second GUI interlock timeout to 10 minutes. Fixed error message that falsely indicates that the view was deactivated instead of just that specific plotitem. Fixed bug in zone isosurface contours that could cause an Unknown Exception. Fixed problem when deleting interfaces when FISH pointers to them existed.
64	1/19/11	z_findface FISH function was returning a 0-based face index. It now returns a 1-based face index (e.g., faces on a hexahedral zone are numbered 1-6). When interactively picking, faces indices are now 1-based.
63	1/17/11	Slot number was being misread in embedded liner command (e.g., "sel liner embedded group clay slot 2 range x=3.9,4.1").
62	1/14/11	Output error message and abort GEN MERGE or GEN SEPARATE process if gridpoints are associated with attachment conditions.
61	1/13/11	Fixed crash during cycling with multistepping and attach logic.
60	1/7/11	Sel liner property plots of cs_sk, cs_ncut, cs_nk ,cs_scohres, slide ,slide_tol were not correct. Range now works with INTERFACE MAXEDGE command. GEN SEPARATE command now requires a range (e.g., gen separate range group rock slot 2, gen separate range x 4 6 y 3 7)
59	12/29/10	The gp_massadd() FISH function was changed, the inertial mass adjustment is now only made during X-component equations of motion.
58	11/30/10	Text percentage in progress bar at bottom of main window now shows completely. Improved user interactions to X,Y and Z range elements in the Range Dialog. The DXF Plotitem now properly saves,restores,and copies assigned filename. DXF data is now also only loaded on demand (when the plot view is first shown after a change). The Attach logic no longer automatically deletes attach conditions when zones they are associated with are Nulled, instead the Attach conditions are ignored. The Attach conditions will be restored if the zones are subsequently set to a non-null material model. The Attach logic checks to see if the master has been Nulled during each large strain update, in case FISH has been used to Null a zone. Some optimizations have been made to internal spatial searching, which should improve performance for GEN ZONE, ATTACH FACE, and other similar operations. A small optimization has been made to speed up the face sorting phase of the ATTACH FACE command. Removing the key during cycling should no longer set your model to the Demo state. The z_model() FISH function should perform much faster if moving to or from a Null model.
57	11/24/10	Interface slip plot did not include no slip case.
56	11/23/10	Solve FOS optimization. Plots of deformed grids were incorrect (Def Factor > 0 in the plot attributes). Accumulated small displacements used in interface logic calculations were erroneously included in the deformation.
55	11/12/10	Bug fix to constitutive models with property tables. Assigning a table property of zero (e.g., PROP TTABLE 0), essentially unassigning a softening table, had no effect until model was saved/restored.
54	11/5/10	Bug fix to APPLY XBODYFORCE, YBODYFORCE, ZBODYFORCE command. The applied magnitudes were incorrect if more than one component was applied (e.g., if xbodyforce and ybodyforce were both applied then the value of the ybodyforce would be the sum of the values of xbodyforce and ybodyforce).
53	11/3/10	Liner FISH function sl_nforce did not work in write mode.
52	10/15/10	Bug fix to find/replace dialog in editor that caused FLAC3D to crash. Bug fix with gridpoint cell-space. In some gp corner-cases it was possible for rounding error to place gp on a boundary outside of extent and be missed (affects gp search functions such as gp_near, gp_nearall, gp_nearlive, and finding a gp for gp histories). Generate an error message if non-existent tables are used in APPLY HISTORY TABLE command.
51	10/8/10	FISH function thtime was readonly. This has been fixed. Added thermal flux vector plot item (Zones->Vectors->Thermal Flux)
50	9/28/10	Fixed translational and rotational velocities were not considered when a node local coordinate system rotated during a large strain update. The following changes have been made to fix this: The SEL NODE INIT command now has a LOCAL keyword allowing you to initialize nodal velocities (xvel, yvel, zvel, xrvel, yrvel, zrvel) in local node coordinates as opposed to global coordinates. e.g., SEL NODE INIT XVEL LOCAL 0.0. Note that the SEL NODE FIX command fixes velocity in the specified directions in terms of the node-local system for all nodes in the range, thereby preventing these velocities from changing. If any translational/rotational velocity component is fixed, then all three components of the translational/rotational velocity should also be specified in the local coordinate system via the SEL NODE INIT {XVEL, XRVEL, YVEL, ...} LOCAL command." Bug fix. MAIL SET commands generated an error message when parsed. FISH function crtime output an error message that model was not configured for thermal instead of creep calculations. This has been fixed. Fix potential crash if read value of local FISH variable before that variable was initialized.
49	8/19/10	Added read/writeable gridpoint FISH function gp_massadd(gp). This function adds intertial mass to a gridpoint during the equations of motion in either static or dynamic modes.
48	8/11/10	Updated text output system to increase GUI responsiveness during LISTings, etc. Stopped appending a '.' character to filenames called with the execute button that had no extension. Extensions are now added automatically to filenames created with most file creation dialogs.
47	7/29/10	Correction to Drucker-Prager and Hydration-Drucker-Prager tension, shear flow rule domain bisector. FISH function substr did not allow a starting character index of 1. FISH function char did not allow character index of 1. Zone data interpolation algorithms (average, IDW, polynomial) were not aware of the Hide Null - Mechanical, Thermal, Fluid attributes of the plot item when doing the interpolation (data from null zones was always included when interpolating the field data). The hidden null zones are now excluded from the zone data interpolation routines (if checked in the plot item's Hide Null attributes). This bug fix affects plot items in the ZoneContour category as well as data extracted using the zone field data FISH functions. Zone field data now includes the following variables: ii = zfd_hidemechnull = ii, ii = zfd_hidethermnull = ii, ii = zfd_hidefluidnull = ii ; where ii is an integer and the variables are read/writeable. These variables indicate which zones should be excluded from the zone field data. Setting a value of 1 causes FLAC3D to exclude those null zones from the data set (e.g., zfd_hidemechnull = 1 would cause FLAC3D to exclude null mechanical zones from the data set). All of the above default to zero if not set (all zones are included in the data set even if they are null).
46	7/15/10	grand FISH function for Guassian distribution was never called (uniform distribution was called in error).
45	7/9/10	Bug fix, restoring save files with group hierarchies caused a crash. Bug fix, RANGE specified with INTERFACE command was only used with the keyword prior to the RANGE keyword in the command list. Now a range applies to all keywords given with the INTERFACE command. Beam axial force and axial stress plots were incorrect (PLOT SEL BCONTOUR {AXIALFORCE\|AXIALSTRESS}).
44	6/29/10	Bug fix to GENERATE ZONE {COPY\|REFLECT} command. The zone join list (zones connected to the faces of a zone) was not being updated properly for zones with less than 6 faces (tet, pyramid, wedge zones). Bug fix. Gridpoint vectors can be restricted to being plotted on model boundaries by selecting the "Boundary" attribute in the plot item. This feature did not work with non-hexahedral grids.
43	6/14/10	Interactive plot elements are now aware of plot magnification, and so rescale themselves properly. Sped up ZONE command by a very great deal. Changed default front and back view clipping planes to much farther away. PLOT commands now have a 5 minute timeout detect, instead of a 30 second one. Correctly use base-1 offset in FISH-plot commands set_color and set_dash_pattern. PlotItems with multi-keyword identifiers (like SEL LINCONT MX) will report command processing errors and help properly. Fix bug in plot of SEL Shell Maximum Fiber Stress.
42	6/10/10	Bug fix to NMD calculations when a mixed grid (tets and other zone types) was used and deterministic mode was on. Check normal vectors before joining faces during grid generation. Add support for unicode input data files, including filenames and strings. Stereonet plots now restore all plot settings correctly. The first "set log on" file now creates a "flac3d.log" file in the correct input directory.
41	6/1/10	Fixed 32/64 bit save file incompatibility problem in interface logic.
40	5/28/10	internal release
39	5/27/10	Bug fix to solving with NMD. Constitutive models accessing a tet in an overlay caused a crash (models affected: Burgers-Creep, Cam-Clay, CYSoil, Cviscous, C.W.I.P.P., Double-Yield, Finn, Hoek-Brown, P.W.I.P.P, Strain Softening, Soft-Ubiquitous-Joint, W.I.P.P.) Movie frames were not generated with the solve command.
38	5/10/10	Fix to FISH functions zfd_getdata(V3+) and nd_ssysx(ND_PTR, V3+). Parser did not read V3 properly when 3 separate parameters (representing x, y, z) were used for V3. Allow two overlays in tet zones. Optimization to stress-strain equations when a mixed-grid (tets and other zone types) is used and deterministic mode is on.
37	4/27/10	Stress-strength ratio zone contour plot renamed to strength-stress ratio plot. Bug fix to NMD calculations when a mixed grid (tets and other zone types) was used and deterministic mode was on. Bug fix to event queue. FLAC3D would sometimes stall in a wait state while solving long problems.
36	4/5/10	Bug fix to NMD calculations on mixed grids (grids containing tets and other zone types). FLAC3D would crash when NMD was on when cycling with a mixed grid. Bug fix to attach logic in dynamic problems. Gravitational mass was not accounted for. Added "About Itasca" help menu entry. Added smiddle, simiddle, srmiddle keywords to stress, strain plot commands. These are aliases for sintermediate, siintermediate, srintermediate keywords. Bug fix to INTERFACE WRAP command. Slot keyword did not work. Bug fix to PLOT SEL BCONTOUR {fish, xdisplacement, ydisplacement, zdisplacement, xvelocity, yvelocity, zvelocity,,axialforce, axialstress} The keywords in parentheses did not correspond to what was plotted. Bug fix to PLOT SEL PCONTOUR {fish, xdisplacement, ydisplacement, zdisplacement, xvelocity, yvelocity, zvelocity, ncstress, scstress} The keywords in parentheses did not correspond to what was plotted. Updated to Qt 4.6.2.
35	2/10/10	Bug fix to attach logic. Attaches were missed. Attach command now runs multiple passes up to a maximum of 5 passes or until no more attachments created.
34	1/29/10	Range logic now works with stereonet plots. Bug fix to attach logic. Attach command would attach more than once if attach face command was run more than once with the same range.
33	1/27/10	Bug fix to impgrid when reading group names of zones. Added check box to gridpoint vector plot item to only plot vectors on model boundaries Impgrid is now case insensitive when interpreting tokens from a grid file. Impgrid now issues a warning message when bad zone geometry is encountered. When solve command is issued, output warning if creep is on, age not specified and local damping is not on. When solve command is issued, output warning if creep is on, age is specified and combined damping is not on. Fixed bug in Polynomial extrapolation of zone values to gridpoints. Fixed error in texture coordinate calculation, contour interval boundaries should now fall in the correct place. Fish functions gp_xfapp, gp_yfapp, gp_zfapp, and gp_fapp did not include body forces, now they do. The DXF plotitem should now update automatically when colors are changed. The list of available Tables and Histories should now update automatically when a new table or history is added. 3 sided zone faces will no longer join to 4 sided zone faces. Vertex Buffer Object OpenGL extension was being turned off in error - this should default to on again. Allowed attempted generation of bitmaps up to 10240 pixels wide - an error message is given if the system fails to create one because it is too large. Picking now works for contour plotitems if vertex arrays are turned off. Fixed bug in zone interpolation using the averaging method -- affected zone contour plots. SEL applied force/moment plots would continue to display values even after they were removed. Linear SEL YDirection was being incorrectly saved/restored.
32	11/20/09	Some display settings weren't being stored in the registry, and weren't being saved to project files. Fixed memory allocation error, causing large amounts of wasted memory.
31	11/11/09	Updated the hydration constitutive model. Added automatic check for available updates - this shows as icon next to Help entry in menu bar. Added "Check Modifications" Help menu entry, loads modifications page into browser. Table read now treat commas as whitespace. Fix to DXF import, now reads DXF files exported from plotting logic correctly. The COLORLIST keyword can now be used to apply colors to plotitems using LABEL based color selection. Update Gp flags and ATTACH conditions when zones changes to/from model NULL. APPLY FACE conditions have the gradient applied on creation, not during cycling. LIST APPLY adds the current value being applied by history modifications. Fix bug in z_linkindex and gp_linkindex functions, were using base 0 instead of base 1. Many changes to example data files. Fixed debug warning output in console version during a SOLVE command. Eliminated possible small negative weighting when distributing forces to zone faces. Fixed error restoring save files with interfaces across 32/64 bit versions -- old save files are not 32/64 compatible. Fixed several sources of crashing if an error occurs during restore. Fixed minimal plot caption drawn when legend is off if the job and view titles are present. Model is no longer marked as modified after project clear or restore.
30	---	Fixed multithreading bug in the table logic (again?). NEW does not reset SET PROCESSORS. Fixed bug in fluid flow logic under SET DETERMINISTIC ON. Fixed Attach bug under SET DETERMINISTIC ON.
29	---	Fixed bug in OpenGL font rendering on some systems. Control panel tool buttons initialized properly on startup. Clipbox attribute controls initialized properly after copy. Polygon attributes were not being copied in zone plot item. PLOT ZONE with colorby uniform now includes color square. Zone block isosurface now includes interactive data (mouseover). Fixed corner case for isosurface generation, when entire model and surface all have the exact same value. Fixed bug in deterministic mode if null zones were present. Add check for dynamic multistepping during deterministic mode calculations. SET DETERMINISTIC ON is now the default.
28	---	Fixed save/restore bug in elastic properties of constitutive models. Added "set deterministic on" option. Eliminates varying accumulated roundoff error in multithreaded calculations, at the expense of 10-20% less speed and a little more memory. Require three characters for MACro command, and two for MAil command. Fix bug in ATTACH logic (gp to edge attaches), introduced in 27. Face-based APPLY commands should execute much faster on large problems.
27	---	Reduced the size of block memory, hopefully resulting in less wasted memory. Implemented alternate multiprocess distribution method, that autmatically load balances. This should result in significant performance improvements on real-world non-uniform problems. On an Intel Core i7 under Vista 64, runtime speeds using a simple timing test are now 1109 kZones/sec on 32bit and 1239 kZones/sec on 64bit. This is 48% faster than version 3.1. Performance should increase even further on non-uniform problems. Please note that on a 64bit system, the 64 bit version is significantly faster. Multithreaded fluid logic, on a Core i7 runs 148% faster. Experimentally multithreaded the mechanical interface cycling logic. Experimentally multithreaded the mechanical attach cycling logic. Increased the speed of the INTERFACE n FACE and INTERFACE n PROPERTY commands by several orders of magnitude on large problems.
26	---	Fixed internal array error - could cause lock-up or crash in plotting. Attempted to speed up view creation when restoring project files. Error message during rendering notifies user that the view has been deactivated. Retrieving the current step number no longer requires a kernel interrupt, speeding rendering during cycling. GUI thread no longer uses 100% of CPU resources while restoring a model during a project restore. The thread also times out only after 2 minutes of no activity reported from the kernel, instead of 30 minutes total restore time. The error dialog gives users the option of 1) Waiting forever, 2) Waiting two more minutes, or 3) aborting. The 60 second timeout while the GUI waits for access to kernel data has been moved back down to 30 seconds, and a more descriptive error message has been added. Fixed problem if error or interrupt occurs during FISH command processing. Attach logic now polls GUI periodically. Fixed 32/64 bit save file incompatibility if using Shear-Yield SEL links. Fixed bug in WATER TABLE FACE logic, PP weren't being calculated correctly. Increased default size of plot legend fonts. This will only be visible if you delete your current registry information in HKEY_CURRENT_USER\Software\Itasca\flac3d400 using regedit.exe. Updated OS reporting during a crash in errorlog.txt. Prevented interface access to group lists during save/restore. Fixed mispelling of "Displacement" in several places.
25	---	Updated example user-defined constitutive model project to reflect changed distribution layout. Updated example C++ FISH plugin project to reflect changed distribution layout. Display properties were not being saved/restored into plot files. SPACE can once again be used to skip cycling without exiting the data file. The plot mode (model, perspective, etc) wasn't being initialized properly. Fix tolerance bug on edge in ATTACH command. Fix memory bug in GEN MERGE, and increased performance slightly. gp_group(gp,num,slot) can have slot=0, as documented. Changing to-and-from NULL mechanical models causes plots to update. Fixed memory allocation bug that could cause kernel lock-up.
24	---	Implemented the "Help/Manual (Draft)" menu item. Gravitational mass update now triggered locally if saturation changes over a tolerance, which should speed up unsaturated fluid/mechanical problems. The (undocumented) command "SET FLUID TOL_SAT " can be used to change this, it defaults to 1e-3. Command processing and FISH settings no longer go into the registry, since they are now considered part of the model state. FIRST PRE-RELEASE
23	---	Fixed update bug in FISH Plot Items. Minor changes the FISH library plug-in interface. The interface definition is in the process of being documented with DOxygen. Water table generation should no longer hang the GUI. Increase the speed of adding water table faces when a grid is present. PLOT ZONE adds zone colorby info to element information info. Fix bug in gp_near and z_near. Many dynamic example problem updates. FLAC3D now pays attention to a true return of modulus_changed_ in the UDM interface. The constitutive model interface has been updated with new functions. The Demonstration mode limits have been changed to 600 zones and 20 structural elements. .DMP files should now be created on program crash in the 64 bit executable. Changed icons for Comment, goto, and syntax verification in the editor. Case sensitivity, safe mode, FISH case sensitivity, FISH safe mode, and FISH autocreate are reset on NEW, and saved and restored. In other words these are considered part of the model state. Fixed clipping bug in chart axes titles if job or view title was present. Changed timeout error waiting for file to restore during project restore, from 2 minutes to 30 minutes. RANGE command processing now gives error if unrecognized trailing tokens are found. Clipbox settings are now properly copied if the view is copied. Bug fix in clipbox command processing. The User-Defined-Model option has been renamed the C++ Plugin Option. Fixed restore bug when using hysteretic damping. Fixed crash if you had a zone plot up when resetting or restoring the model. Made an optimization that should radically increase the speed of the z_prop FISH intrinsic. Fixed save/restore bug in ZONE plotitem, colors were being lost. Fixed command processing for GP based plotitems. Fixed error when parsing multiple-words group names in impgrid.
22	---	All non-fixed velocities are set to zero if dynamic is turned on or off. Non-fixed velocities in a static mechanical model are scaled to conserve momentum if the inertial mass changes as a result of timestep changes. Merged 3.1 bug-fixes and additions made since 4.0 work began -- this includes the addition of several FISH functions (dipfromnorm, ddfromnorm, normfromdip, normfromdipdd, toupper, tolower, z_apply, z_iecreate, z_qualitytest). Scale attributes now disable "Target" option if auto is turned off. Axes once again allow interactive modification of font sizes. Fixed bug in SEL cyield plot and FISH access. Further changes to verification and example problems, mostly in SEL. Progress dialog now displays while loading a project. "Point" polygon mode (shown if both fill and wireframe are off) has been removed, it was just confusing people. Chart axes labels should render faster and cleaner. Fixed a crash if cut-plane faces were passed to the ZGROUP plotitem.
21	---	Fixed misspelling of "asterisk" in keywords. C++ FISH plugins now require a CONFIG command to execute, and a key option to cycle. LIST INFORMATION now lists current CONFIG settings. LIST FISH INTRINSIC now indicates if the intrinsic is a plugin. Editor: Comment/Uncomment toolbuttons should now be inactive if no text is selected. Models cycled under an academic license or a demonstration license are now indicated in the plot legend. Models cycled under a demonstration license have a "Demonstration Model" watermark added to plots. Movie and Legend options are once again present in the "Display" tab of the Tools/Options dialog. Axes direction labels are once again depth-tested. Suppressed page title-bar context menu. Repaired constitutive model example, the setProperty() method had not been updated. Fixed crash if a PRANGE range element is used to modify the model -- error is properly reported instead. Fixed SEL node stiffness listing output. Added HISTORY SEL SFORCE, which works for any SEL element type. Added "Null Faces Only" option to zone-based plots, to more easily allow plotting on the surface of excavations. Fixed bug in HISTORY SEL RECOVER.
20	---	Updated PRINT FISH MEMITEM (a little). The global shortcuts CTRL-ALT-C will Call and data file, and CTRL-ALT-R will Restore a save file directly without using the general Open Items dialog. Fixed problem where CTRL-S would not save the item in some circumstances. Improved "Edits made outside interface" check. Improved "changes made to edit file" system, should get less false change states. Added "demonstration" mode if registry entry present. Many updates to example and verification files, in preparation for manual. Fixed crash if you attempt to load the same constitutive model twice.
19	---	The Group range element can now refer to more than one group, both on the command line and in the GUI. GEN MERGE should operate much faster now. Plot "Fail" (Stress-Strength ratio) can now specify a model and properties to use, as well as use the actual model and properties assigned to the zone. PLOT ZGROUP - plots zones colored by groups showing internal group boundaries - defaults transparent. Fixed examples dialog so it works on installation layout. Added prototype HELP file. Separate command documentation into CommandReference.pdf and PlotCommandReference.pdf Misc. Bug fixes. INTERNAL BETA MARK
18	---	"Fail" added to zone plots. Added default color Listing controls to List Options tab. Memory allocation system goes all the way down to a single OS page in search of memory, other bugs fixed. Added "Zone Field Data" FISH functions, allowing generic data to be efficiently gained at arbitrary points in the model. Updated to Qt 4.5.2. Project page name in File menu context section changed. "Open" option in Project page context menu made bold to indicate it is the double click action. "Open in explorer" option added to Project page context menu. Further Editor changes: - CTRL-N is a hotkey for "Create New File". - Added Menu option to edit the "predfined data file header". - comment/uncomment block options disabled if there is no selection. - highlight zero-argument FISH intrinsics. - Merge "Define" and "Fish Variable" syntax highlighting color selection. - Test for external file updates when editor gains focus. - General fixes to syntax highlighting. Added Example Project Dialog to the Help menu. Added simple "Line Profile" plotitems (as chart and as line). Improved "Body" PlotItem. Fixed GEN POINT using assigned IDs. Bug fixes and behavior improvements to stereonet plots. Added zone-based "Condition number" plots - a measure of how badly deformed a zone is.
17	---	Fixed SaveAs icon in File menu. Added hotkeys to some editor commands. Improvements to behavior of dock window maximization. Body logic plotitems now update when the body data changes. Stereonet plotitem - principal stresses, strains, and zone extra. Zone and Gp spatial search optimization, affects gp_near, z_near, etc. FISH functions. Fixed bug in new memory system causing much wasted memory in plots. Fixed bug - display not automatically updating if regenerate pressed. Fixed occasional shutdown crash.
16	---	Fixed several bugs in new memory system. Fixed named range error. Temporary files created during SOLVE FOS no longer show up in the project. Added anti-aliasing check to startup. Many new icons. New features added to built-in editor, this work is ongoing, including - Search improved, Replace added. - More standardized shortcut keys. - Comment/Uncomment selection tool. - User-specified predined header file for new data files. - Can tab/untab selections. - Smart-indenting. - FISH intrinsic names highlighted in BOLD in FISH functions. - Testing "macro" expansions of common structures "def..", "com..", etc. Line dash patterns now more distinct on plots. Added SKIP and MAXNUMBER controls to vector plots. Updated NMD system, now more reliable & faster. 3x speedup on 4 processors. GEN ZONE COPY and GEN ZONE REFLECT much faster, GEN ZONE in general sped up. Fix RANGE VOLUME. Fix property assignment bug in cysoil. Many other minor bug fixes.
15	---	Updated and optimized memory use. Save, Restore, and especially New are faster. Models now use 30% less memory. Fix crash when closing docking window. Fixed bug when removing GP apply items.
14	---	Fixed some 32/64 bit save/project file bugs. Old save files still won't work.
12	---	CTRL-4 will now hide project page if it is already active. Fix bug reading values from vertex contours. Fix volume-weighted averaging for zone contours (the default method).
11	---	Merged "Execute" and "Stop" buttons on toolbar. Added new interactive options to move a simple cut plane normal to its surface. PlotItem list "modify range" icon now changes if a range has been specified. Added option to suppress tool tips. Fixed problem with "60 second timeout" error, it should no longer hang the GUI. Added "let the kernel do something" polling to zone interpolate calculations. Fixed APPLY boundary detection, was taking far too long. Added polling to model reset operations, hopefully preventing a timeout. Fixed bug if property missing from model in zone colorscale plot.
10	---	Removed LEGEND lines from DXF export, and suppressed chart plots as well. Made mouse mode tool buttons mutually exclusive again. "Last used" directories for Open Item, Export Bitmap, and Export DXF dialogs now reset if the current directory chagnes. Fixed bugs in how the range logic interacted with zone group slots. Interface properties should now recognize FISH replacements.
9	---	Added Zone Group slots, added GEN SURFACE. Fix bug in File/Grid/Import if data file selected already exists. Added Export to Bitmap tool button to plot context tool bar. Sped up saving. Fix bug in ATTACH FACE, some edges were being missed. Fix memory leak in GEN SURFACE. Added polling during APPLY REMOVE.
8	---	Stage 6 internal release.
7	---	Fixed FISH I/O bug in text mode, only affected by very large files.
6	---	Stage 5 internal release.
5	---	Changed version number to 4.00.
4	---	Fixed attribute bug in Scale Box plotitem. Fixed dropped polygons along corners or Wedge and Octant cut-planes (work still needed here). Fixed potential crash while plotting histories or tables. Fix bug in FISH error message processing. The ubiquitous joint model was corrupting joint orientation information in it's large-strain update. Copy Plot wasn't copying clip-boxes, cut-planes, or polygon attributes of plotitems. Contour-coloring of tensors, vectors, and color-scale plots could have been incorrect if a manual maximum was used. Further improvements to how consistent orientation information is set and updated in constitutive models. Constitutive model access to table logic was not thread safe, producing incorrect results. This affects all softening models, for example.
1	---	Created

FLAC3D 3.1 (Previous Version)

Pre-Requisites Update Instructions Revision History

User-Defined Models Help File

Starting with subversion 424, an update to the Microsoft (R) Visual C++ 2005 Runtime is REQUIRED. (The current subversion of the code update available here appears at the end of the Revision History — see link above.) If not present, on startup the code will return an error message similar to "The application has failed to start because the application configuration is incorrect. Reinstalling the application may fix this problem" (on Windows XP or later) or it will complain of missing DLLs (Windows 2000).

Please use the link above that corresponds to the 32-bit or 64-bit version of FLAC3D. The installation process need only be performed once; subsequent updates to FLAC3D will not require these steps to be repeated.

Close This PanelX

FLAC3D 3.10 Modification List  (Type PRINT VERSION to check version number.)
Changes listed are included in  F3310_32.EXE and F3310_64.EXE.

Sub-
Version Date
Number  MM/DD/YY

================================================================================
420     12/19/06
================================================================================

Official release version of FLAC3D Version 3.10.

================================================================================
421     12/28/06
================================================================================

Addressed 64-bit portability problems.

================================================================================
422     01/02/07
================================================================================

Hoek-Brown model no longer forces a stress point back to Hoek-Brown failure-
surface if the iterative algorithm fails to converge.  An error message is
output instead.

================================================================================
423     01/08/07
================================================================================

Addressed 64-bit window creation problem.

Bug fix in recording more than one history for recovered stresses in structural
elements.  e.g., HISTORY SEL RECOVER SRES SURFX ...
If more than one history was created, then in large strain mode it was possible
to get a "surface system not set at node" error.

Bug fix: File names that contained ASCII character codes > 127
(characters with accents, foreign character sets, ...)were not recognized.

================================================================================
424     01/16/07
================================================================================

Hysteretic damping bug fix.  If hysteretic damping was used with a model that
was run on more than one processor then FLAC3D could crash.

================================================================================
425     01/18/07
================================================================================

PRINT DYNAMIC command now outputs damping type information based on what is
set at zones and gridpoints.

Added a PRINT ZONE DYNAMIC command to output hysteretic damping information for
zones.

Plotting bug fix.  When block plots were created (e.g., PLOT BLOCK STATE), then
deactivated in the plot window, the model run for a period of time, then that
plot re-activated, FLAC3D would sometimes issue an error message.  This has
been corrected.

================================================================================
426     01/19/07
================================================================================

Put a lower cap on the hysteretic damping multiplier of 1e-3 to avoid shear
modulus from becoming zero.

================================================================================
427     01/25/07
================================================================================

The "versus step" radio button did not work in the modify history plot item
dialog.

Bug fix which caused FLAC3D to occasionally crash when plotting.

================================================================================
428     01/29/07
================================================================================

Added USEGIVENIDS keyword to IMPGRID command.  With this keyword, the ID numbers
of gridpoints and zones specified in the *.FLAC3D file are used in the final
grid.  Otherwise gridpoints and zones are generated using the next available ID
numbers (the default).
Example usage:  IMPGRID yourfile.flac3d USEGIVENIDS

================================================================================
429     01/31/07
================================================================================

Optimizations to IMPGRID command when USEGIVENIDS keyword is specified.

================================================================================
430     02/05/07
================================================================================

Bug fix which caused FLAC3D to occasionally crash when plotting.

================================================================================
431     02/12/07
================================================================================

Gridpoints can now belong to an unlimited number of different groups.  Prior to
this change, gridpoints were only allowed to belong to a maximum of 12 different
groups.

================================================================================
432     02/13/07
================================================================================

Cylinder and cshell mesh primitives did not work properly if sweep angle was
greater than 90 degrees.

================================================================================
433     03/01/07
================================================================================

Set default plot background to white instead of light gray.

Plots were not always being refreshed when an existing plot window was reopened.

Dangling pointers result when a variable points to an object (e.g. a zone) which
no longer exists (e.g., the zone was deleted). Dangling pointers can cause
unpredictable behavior if used to access the deleted object.  FISH pointers are
now checked before writing to sav files to see if they point to real objects.    

SET RANDOM  command
sets the seed for the random-number generator (used by the PROPERTY command) to
the integer iseed. By default, iseed equals 10000.  If a seed value is not
given or is zero then the random number seed is reset to it default value of
10000.  The seed is not reset by the NEW command.

================================================================================
434     03/02/07
================================================================================

GEN ZONE REFLECT now works with degenerate brick (dbrick) type zones.

================================================================================
435     03/07/07
================================================================================

Added NO_PP and CALC_PP keywords to WATER TABLE command.

When a WATER TABLE FACE command is executed, the pore pressures are recalculated
for all gridpoints using all the previously stored water table faces.  For
models with a large number of water table faces the computation time becomes
excessive.  The NO_PP keyword suppresses the pore pressure calculation (only the
face is stored).  The CALC_PP keyword can be used without any face
specifications to force a recalculation of pore pressures using all stored water
table faces.  A WATER DENSITY command will also force a recalculation of pore
pressures using all stored water table faces.

Example usage:

; store a water table face only, no pore pressures are calculated
water no_pp table face  2.5490e+003, 1.0200e+004, 2.9591e+003 ...

; store a new face and recalculate pore pressures using previously stored faces
water table face  2.5490e+003, 1.0200e+004, 2.9591e+003 ...

; recalculate pore pressures using previously stored water table faces
water calc_pp

; recalculate pore pressures using previously stored water table faces and new
; water density 
water dens 1000

================================================================================
436     03/19/07
================================================================================

GEOM_TEST did not return proper values if run with a tetrahedral grid.

Cut plane dialog did not update the cut plane properly if user changed a dip or
dipdir value and selected OK without tabbing out of the field first.

SOLVE ELASTIC didn't work if a zone with a Null model was encountered (because
it wasn't a Mohr-Coulomb model).  Models with Null and Mohr-Coulomb zones now
work with SOLVE ELASTIC.

================================================================================
437     03/27/07
================================================================================

Bug fix to BLOCK STATE plot item.  Error would occur if there were no zones to
display.  Also user choices NOW, PAST, NOT were not saved 
(e.g., PLOT BLOCK STATE NOW).

================================================================================
438     04/03/07
================================================================================

Bug fix to contour legend.  Extra contour interval was sometimes shown in the
legend (e.g., a null interval from 0.00 to 0.00).

================================================================================
439     04/03/07
================================================================================

Added new FISH function for retrieving hysteretic damping values from zones.
value = z_hyst(p_z, vname)
vname is a string with the following possible values:
        'c0','c1','c2','c3','modfac','model'
        'e1','e2','e3','e4','e5','e6'
'c0', 'c1', 'c2', 'c3' = 1st through 4th parameter of hysteretic model
'model' = hysteretic model identification number, return values as below
     0 = no hysteretic damping
     1 = default model
     2 = off
     3 = sig3 model
     4 = sig4 model
     5 = hardin model
'modfac' = current modulus reduction factor for this zone
'e1',...'e6' - strain components (gamma1, ... gamma6 in equations 
     3.34-3.39 Optional Features manual)
If the zone does not have a hysteretic damping model then a value of zero is
returned.

Search time optimization to FISH functions p_z = find_zone(id) and
p_g = find_gp(id).  Search time reduced from O(n) to O(log n).

Faster importing of grids with IMPGRID command when USEGIVENIDS is used.

Contour plots of extra variables were not being refreshed during solving.

Bug fix to interfaces.  A math error (e.g., sqrt of a negative number) would
sometimes occur when interfaces were used.  This has been attributed to a bug
in the optimization compiler (optimizations are now turned off for interface
calculations).

If a zone is experiencing plastic flow and hysteretic damping is active, then
the strain path components (equations 3.34-3.39 Theory and Background) are
rolled back to their previous state.  The sequence of events is outlined below:
1) calculate hysteretic damping multiplier with accumulated strains
2) call the constitutive model using the multiplier
3) if all the tet overlays in a zone are yielding then subtract strain path
   components from the currently stored values in the hysteretic damping model.

================================================================================
440     04/27/07
================================================================================

In a dynamic analysis with Rayleigh or hysteretic damping, the constitutive 
model bViscous state flag is used to indicate whether or not plastic flow is
occurring in the sub-zones.  If plastic flow is occurring then bViscous is set
to false.  If the model is in the elastic range, bViscous is set to true.  
FLAC3D sets bViscous to false before calling the constitutive model for each
tetrahedron in the overlays.

In FLAC3D a value of false for bViscous inhibits the stiffness-proportional
Rayleigh damping term.  The bViscous flag was not set in several of the
constitutive models and as a result stiffness-proportional Rayleigh damping was
always inhibited in the following models: anisotropic, orthotropic, Drucker,
Ubiquitous-Joint, Strain-Softening, Soft-Ubiquitous-Joint, Cam-clay,
Double Yield, Finn, Hydration, Viscous, Power law, W.I.P.P., C.W.I.P.P.,
Burger-Creep, P.W.I.P.P., Burger creep-viscoplastic.  Stiffness-proportional
Rayleigh damping is now allowed when these materials are in the elastic range
(and inhibited otherwise).

The Hoek-Brown constitutive model always enabled stiffness-proportional Rayleigh
damping terms. These are now inhibited if the material has no strength
(i.e, Sigma_ci <= 0) or the material is yielding.

If hysteretic damping is active, FLAC3D also uses the bViscous flag to determine
whether or not to correct the hysteretic damping strains (see revision 439
above).  If all tetrahedral overlays in a zone have bViscous set to false (all
tets are plastically flowing) then the hysteretic damping strains are corrected.
The constitutive models mentioned in the previous two paragraphs will be affected
if they were used with hysteretic damping.

Nodal Mixed Discretization (NMD) is now enabled for all the constitutive models
except for P.W.I.P.P., Burger creep-viscoplastic and viscoplastic-power.
Prior to this change NMD was functioning in Mohr-Coulomb and Strain-Softening
models as well as the following models (which do not require a NMD plastic
stress correction): elastic, anisotropic, orthotropic, Viscous, Power law,
W.I.P.P., C.W.I.P.P., Burger-Creep.  This revision added the NMD plastic stress
correction to the following models: Drucker, Ubiquitous-Joint, 
Soft-Ubiquitous-Joint, Cam-Clay, Double Yield, Hoek-Brown, Finn, and Hydration.

These changes are also reflected in the user defined constitutive models source
code.

================================================================================
441     05/03/07
================================================================================

When a water table was created with the WATER TABLE FACE command, the algorithm
sometimes failed to assign pore pressures to some gridpoints.  These gridpoints
were located near the edges of water table faces.

================================================================================
442     05/24/07
================================================================================

Added a group name to tracks
e.g., track line 0 .1 0.5 0 .1 2.5 20 group AAA
and range logic for particle track plotting
e.g., plot add track cyan range group AAA
Groups can only be assigned to tracks when they are created.

An xyz type range will also work for track plotting, however, as the particle
travels with the fluid its position changes so it may move in and out of various
xyz ranges.
e.g., plot add track cyan range x -1 2
will plot all tracks currently within x -1 2 as cyan

Hoek-Brown constitutive model did not always converge, typically when the model
was in tension.  This is now fixed.

================================================================================
443     06/12/07
================================================================================

Restricted SOLVE FOS to small strain simulations.

In coupled fluid-mechanical simulations a saturation change causes the saturated
density contribution to body force to change.  This component was not being
updated.

When contouring zone related quantities, the zone values are first projected to
gridpoints using either a gradient (default) or average method. A least-squares
curve fitting process is used in the gradient method and degenerate matrices are
sometimes encountered (not-enough data to compute gradients or data does not
vary enough). In these degenerate cases, spurious contours sometimes appeared in
the contour plots.  This has been fixed. 

================================================================================
444     06/13/07
================================================================================

Bug fix to restoring sav files.  A stream I/O error message would sometimes
appear when attempting to restore a sav file.

================================================================================
445     06/26/07
================================================================================

Improvements to contour plotting with gradient method.

The command "INTERFACE n WRAP GroupA GroupB" had the restriction that it would
only put an interface on open faces in GroupA that had gridpoints that
exclusively belonged to GroupA.  This restriction has been removed.

Added FISH read only access to gridpoint x- y- z-accelerations:
gp_xaccel(pg), gp_yaccel(pg), gp_zaccel(pg).
This is only available if the model is configured for dynamics.

Histories of gridpoint x- y- or z-accelerations were always available but are
not documented in the Command Reference.  These accelerations are only available
if the model is configured for dynamics.
e.g. 
history id n gp xacceleration {or yacceleration, zacceleration} x y z {or id m}

================================================================================
446     06/26/07
================================================================================

Bug fix to FISH library.  Any FISH code would crash.

================================================================================
447     06/26/07
================================================================================

Plotting bug fix.  If mechanical zones are made null their stress state is set
to zero.  Prior to this change, the zone stresses that were stored prior to
making the zone null were not zeroed, making it appear as if null zones could
carry stresses.

================================================================================
448     06/28/07
================================================================================

The command "IMPGRID yourfile.FLAC3D USEGIVENIDS" has been modified to allow
reading of partial grid definitions.  For example if you have the following
grid files:

    g1.FLAC3D - gridpoint definitions
    g2.FLAC3D - more gridpoint definitions
    g3.FLAC3D - more gridpoint definitions
   z13.FLAC3D - zone definitions using gridpoint ids from g1 and g3 
   z12.FLAC3D - zone definitions using gridpoint ids from g1 and g2
groups.FLAC3D - groups that zones fall into

... you can now import the grid piecemeal as follows:

IMPGRID     g1.FLAC3D USEGIVENIDS
IMPGRID     g3.FLAC3D USEGIVENIDS
IMPGRID    z13.FLAC3D USEGIVENIDS
IMPGRID     g2.FLAC3D USEGIVENIDS
IMPGRID    z12.FLAC3D USEGIVENIDS
IMPGRID groups.FLAC3D USEGIVENIDS

The only requirement for this to work is gridpoints must exist prior to creating
the zones that require them.  Also, zones must exist prior to placing them in
groups.  With the USEGIVENIDS keyword, FLAC3D will not renumber the gridpoint
or zone ids; it will use the ids as specified in the input files.

================================================================================
449     07/06/07
================================================================================

Added USEFREEGPS keyword to the GEN ZONE command.  Prior to this change the
GEN ZONE command would only reuse (merge with) gps that existed on boundary
faces of the existing grid.  If a user generated gps with the GP command and
these gps were not connected to any zones then these gps were not considered in
the GEN ZONE command when searching for gps to merge with.  These 'free' gps are
now used in the GEN ZONE command if the USEFREEGPS keyword is given.  This gives
users the potential to create gps with given ids at specified locations.

e.g.,
gp id 100 0 0 0
gp id 200 1 1 1
gp id 300 0.5 0 0.5
gen zone brick size 2 2 2 p0 0 0 0 p1 1 0 0 p2 0 1 0 p3 0 0 1 usefreegps
plot sur yell gpnum on

gps 100 200 300 exist in the grid when usefreeps is specified otherwise they
are not connected to any zones.

================================================================================
450     07/19/07
================================================================================

FLAC3D would occasionally crash if you tried printing to an inaccessible
printer.  This has been fixed.

================================================================================
451     08/01/07
================================================================================

A stress leak that could occur in partially saturated zones under volumetric
straining has been corrected.  The bug fix addresses coupled drained/undrained
fluid mechanical simulations under unsaturated conditions.  For those
simulations, stresses in the unsaturated part of the domain may have been
affected provided fluid bulk modulus was non-zero there.

================================================================================
452     08/20/07
================================================================================

Not all plot items were refreshed in a view during cycling.
E.g. in a small strain problem, plotting the grid geometry with a magnification
factor > 0 (plot grid mag 100) would result in the plot not getting refreshed.
All plot items are now refreshed during plotting.

Added readonly FISH function gp_ppfix(p_g).
Returns 1 if the gridpoint has a fixed pore pressure, returns 0 otherwise.

Added a new plot item which plots a marker (a circle) at gridpoints with a fixed
pore pressure.
e.g., PLOT PPFIX RED SCALE 0.01 RANGE X -0.1 0.1

================================================================================
453     08/23/07
================================================================================

The APPLY command would sometimes cause FLAC3D to crash when applied items were
removed.  This has been fixed.

FLAC3D would crash if a FISH scalar variable such as dytime was accessed before
a CONFIG DYNAMIC command was issued.  FLAC3D now returns 0 for FISH scalars
fltime, dytime, crtime, thtime, fldt, dydt, crdt, thdt, fluid_ratio, therm_ratio
if FLAC3D has not yet been configured for the associated model with command
CONFIG {DYNAMIC, FLUID, THERMAL}.

Modified functionality of SET VISIBLE {MECHANICAL, FLUID THERMAL} command.

The SET VISIBLE command specifies which zones to consider as non-null when
creating contour plots.  SET VISIBLE MECHANICAL, for example, will signify to
FLAC3D to show all zones with a non-null mechanical model in contour plots
(zones with a null mechanical model assigned will not be shown).  All zones with
the specified visible model (one of type mechanical, fluid, or thermal) will
contribute data to the contour plot independent of the variable being plotted.

For example, if SET VISIBLE FLUID is specified then all zones with a live fluid
model (fl_isotropic or fl_anisotropic) will contribute data to contour plot.
fl_null zones or zones that haven't been assigned a fluid model will not
contribute data to contour plots.  If a stress plot (e.g., PLOT CONT SZZ) is
created with SET VISIBLE FLUID specified then all zones containing a live fluid
model will contribute data to the average or gradient method used for contouring
SZZ whether or not that zone has a live mechanical model.  In this case if a
zone doesn't have a live mechanical model then that zone will contribute a zero
value to the average or gradient method used for contouring.  With
SET VISIBLE FLUID all zones with fl_null or unassigned fluid model will NOT
contribute data to the average or gradient calculations.

It is advisable to specify SET VISIBLE MECHANICAL when plotting any values that
are associated with a mechanical model (e.g., stress, strain, any zextra variables
that are derived from mechanical quantities, ...).  SET VISIBLE FLUID or THERMAL
should be specified when contouring data that is meaningful with respect to
a fluid or thermal model respectively.

Prior to this change, FLAC3D attempted to classify some of the plot quantities.
For example, a contour plot of SZZ would be classified as a mechanical model
type of plot and zones without a mechanical model or a null mechanical model
would not contribute to the average or gradient calculations regardless of what
SET VISIBLE was set to.  With the change in this sub-version, FLAC3D no longer
makes this assumption about plot quantities.

================================================================================
454     09/04/07
================================================================================

Fixed bug in uwedge gridding primitive.  FLAC3D was not able to generate a
uniform wedge in some cases where the aspect ratio of the uwedge was too high.

================================================================================
455     09/26/07
================================================================================

The line "*** Add xxxx to all labels ***" would appear above the contour label
caption if the max and min contour limits were almost identical.  This allowed
small variations in contours to be discerned (e.g., if the highest contour value
was 2.0000001e10 and the lowest contour value was 2e10 then the caption would
display "Add 2e10 to all labels" with the contour intervals between 0 and 
0.0000001e10.  This feature has been removed.

Bug fix with PLOT PRINT VIEW command.  If a plot window was not open then
PLOT PRINT VIEW would cause a Windows error.

Increased the search angle tolerance when FLAC3D searches for free faces to
connect liner nodes to zones with node-zone links.  Prior to this, liner nodes
may not have had a node-zone link created.  This problem appeared if zone faces
were not planar.

================================================================================
456     11/01/07
================================================================================

Increased the search angle tolerance when FLAC3D searches for free faces to
connect liner nodes to zones with node-zone links.  Prior to this, liner nodes
may not have had a node-zone link created.  This problem appeared if zone faces
were not planar.

Bug fix to Finn constitutive model.  Not all Finn state variables were saved.
Results would differ if a run was continued from a restored state rather than 
running the model from start to finish completely (without an intermediate 
save restore).

If gravity is non-zero and a mechanical calculation is performed, FLAC3D now
checks for positive densities in zones.

================================================================================
457     11/16/07
================================================================================

Bug fix for thermal calculations with the th_null model.  The null thermal model
(th_null) did not behave as an adiabatic boundary in the following scenario:
A model was created with a non-null thermal model (e.g., th_iso).
This model was run for one or more cycles.
After cycling some of the non-null thermal zones were replaced by th_null zones
and cycling was continued.

Proper adiabatic behaviour was modeled if th_null zones were installed prior to
any cycling.

The Hoek-Brown constitutive model was unable to bring the stress point back to
the Hoek-Brown surface in some tensile stress states and an error message was
displayed.  This has been fixed.

================================================================================
458     11/16/07
================================================================================

Further enhancement to Hoek-Brown model.  With unrealistic initial conditions,
such as having stresses outside the yield surface without any accumulated
strain, the Hoek-Brown model would use a flow paramater, gamma, equal to zero.
The flow parameter is now calculated from the current stress state whether or
not incoming strain is non-zero.

================================================================================
459     11/29/07
================================================================================

Stress tensors with extremely small components could not always be resolved into
principal components and an error message would be displayed.  This has been
fixed.

================================================================================
460     12/13/07
================================================================================

Added a COLBYORDER (color by order) keyword to the STENSOR plot item.  When
COLBYORDER is ON, principal stress components are plotted with three colors.
SMin (most negative, most compressive) principal stress is plotted in color 1.
SMax (most positive, most tensile) principal stress is plotted in color 2.
SMid (the intermediate) principal stress is plotted in color 3.

Prior to this change the STENSOR plot item displayed the principal stress triad
in two colors:  color 1 for compressive stresses and color 2 for tensile
stresses.

Example usage:
PLOT STENSOR COLBYORDER=ON
PLOT STENSOR CO=ON YELLOW GREEN ORANGE

When assigning floating point values to FISH variables, FISH only recognized
a lowercase 'e' when an exponential format was used (e.g., value = 1e8).  FISH
now also recognizes an uppercase 'E' (e.g., value = 1E8).

================================================================================
461     01/04/08
================================================================================

Added FISH function pz = z_containing(x,y,z).  This function returns a pointer
to a zone that contains an arbitrary point with coordinates x,y,z.  Null is
returned if x,y,z does not reside in any zone.

Bug fix to structural element creation.  When structural elements were created
with nodes near the faces of zones, FLAC3D would sometimes have difficulty
finding a zone to create new node-zone links with.

================================================================================
462     01/15/08
================================================================================

FLAC3D would crash during solving if a structural element node had a node-node
link pointing to a deleted node.  FLAC3D now outputs an error message.

Hoek-Brown model would not coverge when s1 = s3 and both principal stresses were
tensile and below the apex of the Hoek-Brown strength envelope.  These stress
points are now moved to the apex.

================================================================================
463     02/14/08
================================================================================

The GEN SEPARATE groupname command now works with zones that are partially
connected to other zones.  Prior to this change, GEN SEPARATE would only
separate gridpoints on a zone face if that zone face shared all its gridpoints
with an adjacent zone.  The separation criterion is now based on the groups a
gridpoint belongs to.  If a gridpoint is shared by a zone belonging to groupname
and a zone not belonging to groupname then that gridpoint is duplicated with
the new gridpoint replacing the old gridpoint in groupname.

Added a plane-strain free-field.  This free-field is invoked with the 
keyword Y_PLANE or X_PLANE on the APPLY FF command.  If the Y_PLANE keyword is
given then displacements in the y-direction are automatically fixed, initial
y-velocities are set to zero and an appropriate 1D free-field is created on the
x-min and x-max extents of the grid.  Similarly if the X_PLANE keyword is given
then displacements in the x-direction are fixed etc.

FLAC3D 64-bit now works with remote desktop (FLAC3D 32-bit already worked with
remote desktop).

================================================================================
464     02/19/08
================================================================================

FLAC3D would output a divide by zero error in a dynamic analysis if structural
elements did not have any free degrees of freedom (all nodes were fixed).  The
stable timestep is inversely related to the maximum frequency of a structual
element node (which was zero in this case).  If there aren't any free DOFs for 
structural elements then the stable timestep is governed by other components 
in the model (e.g., zone stiffness).

Added z_vol_deformed(zp) FISH function which returns the volume of a zone using
gridpoint positions that take into account the accumulated displacement at the
gridpoints (intended to be used in small strain calculations).

If a model contained zones and interfaces and zones were deleted followed by
saving the model, FLAC3D would sometimes crash.  Interface nodes point to their
host zone (the zone it is connected to) as well as to a target zone (the zone it
interacts with).  When a zone was deleted the interface node pointer to the
deleted target was not updated.  This has been fixed.  Stepping the model before
saving also resolves this problem because target pointers are refreshed at the
start of stepping.

================================================================================
465     02/27/08
================================================================================

Added the capability to delete interfaces completely.  Prior to this change
an INTERFACE N DELETE command would only delete the interfaces and nodes and
elements (the interface with 0 nodes and 0 elements would still exist).

In FISH, the interface pointers i_node_head and i_elem_head now return null
if the list is empty.  Prior to this change a node or element with id number = 0
was returned.

Added FISH functions 
numnodes = s_get_nforce(s_ptr, ndforcearray, ndptrarray) and
nwritten = s_put_nforce(s_ptr, ndforcearray)

Function s_get_nforce gets the nodal forces from the structural element pointed
to by s_ptr.  On return, the nodal forces are placed into array ndforcearray
which must be size of at least 6 * number of nodes the SEL contains.  This
function returns the number of nodes the sel has (2 or 3) and also places the
pointers to these nodes in array ndptrarry which must be size of at least the
number of nodes the SEL has.  The forces in ndforcearray are in the global
coordinate system and are arranged as follows:

ndforcearray(1) =  force_x for node 1
ndforcearray(2) =  force_y for node 1
ndforcearray(3) =  force_z for node 1
ndforcearray(4) = moment_x for node 1
ndforcearray(5) = moment_y for node 1
ndforcearray(6) = moment_z for node 1
ndforcearray(7) =  force_x for node 2
ndforcearray(8) =  force_y for node 2
...
ndforcearray(18) = moment_z for node 3

Function s_put_nforce puts the nodal forces that are stored in array
ndforcearray into the structual element's nodes.  On return nwritten holds
the number of values that were written to the SEL nodes (either 2*6 or 3*6
depending on the type of SEL).  ndforcearray is assumed to have nodal forces
in the global coordinate system.

================================================================================
466     04/08/08
================================================================================

Bug fix to strings.  Strings are reference counted.  If the number of instances
of a string exceeded 2^16 then FLAC3D would crash.  The reference count limit
has been increased to 2^32.

If a PLOT EXTRACT point fell precisely on a zone surface this point would 
sometimes be treated as falling outside the zone in the PLOT EXTRACT logic.
This has been fixed.

When creating a new node-node link on side2 of an embedded liner with the
SEL LINK sid SIDE2 TARGET NODE TGT_NUM tid command
the side2 keyword was ignored (the side2 keyword was not ignored when creating
a node-zone link).

================================================================================
467     04/10/08
================================================================================

Plots containing text labels that were output to bitmap based files (bmp, pcx,
dcx, avi, Windows clipboard) had a white background surrounding the text.  Text
now has a transparent background in these file types and matches the image
displayed in the plot window.

Modified friction scaling in Solve FOS routine.
If friction is scaled during bracketing (non-associated flow rule) and dilation
angle is greater than the scaled friction angle then dilation is set equal to
friction.  Friction angle is now clamped to fall within the 0 to 89 degree range
when bracketing.

When plotting multiple graphs on the same plot (e.g. PLOT TABLE 1 LINE 2 LINE)
extraneous lines sometimes appeared on the plot.

Screen was not refreshed after issuing a PLOT SET DISTANCE command.

================================================================================
468     05/09/08
================================================================================

When plotting applied nodal forces (PLOT SEL FAPP) the forces are plotted as a
mix of global and local coordinate systems and it was impossible to decipher
which applied force was global and which was local.  Also, forces applied in a
local coordinate system were plotted in the global coordinate system without an
appropriate rotation (e.g., a force in the local x direction was plotted in the
global x direction).  Applied nodal forces are now plotted in the global 
coordinate system.

Added NOW and PAST switches to PLOT INTERFACE SLIP SHEAR command.  These
switches are only available on the command line.  
PLOT INTERFACE SLIP SHEAR PAST will only plot points that have slipped in the
past (points that are currently slipping are not plotted).
PLOT INTERFACE SLIP SHEAR NOW will only plot points that are currently slipping
(points that have slipped in the past are not plotted).  If NOW and PAST are
not specified then slipping now and slipped in the past points are plotted.

================================================================================
469     05/16/08
================================================================================

The FLACEXTRUDE command was unable to locate the FLAC 6.0 executable.  This has
been fixed.

================================================================================
470     05/28/08
================================================================================

Optimizations made to structural element calculations.

FLAC3D would return with an access denied error when trying to print to some
network printers.  This has been fixed.

================================================================================
471     07/03/08
================================================================================

The Hoek-Brown model makes three attempts to bring a stress point to the yield
surface.  
1) a fast Newton solver is tried.
2) if 1 did not converge the stress point is checked to see if it falls below
the apex of the Hoek-Brown envelope (crossed the sigma1 = sigma3 line).  If this
is the case then the stress point is set to the apex.
3) if 1 and 2 did not work a bisection method is used to find the stress point
on the yield surface.

The hb_ind property reflects the attempts made to bring a stress point to the
yield surface.  Each zone contains a number of tetrahedral subzones which are 
sent in sequence to the constitutive model.  When the first subzone is received
by the Hoek-Brown model, the Hoek-Brown model sets hb_ind to 0.  If case 1 works
hb_ind is set to max(hb_ind, number of iterations required by Newton solver).
If case 2 is encountered hb_ind is set to max(hb_ind, 1000)
If case 3 is encountered hb_ind is set to 
max(hb_ind, 1000 + iterations required by the bisection algorithm)
If none of the cases work then hb_ind is set to 9999.

hb_ind contains the highest value (of hb_ind) encountered for all the subzones
in a zone encountered during a cycle.  Possible values for hb_ind are summarized
below:
hb_ind = 0 - all subzones had stresses in the elastic range
0 < hb_ind < 1000 - Newton solver was used on one or more of the subzones
         (hb_ind = max number of iterations required by Newton solver)
hb_ind = 1000 - stress point forced back to apex on one or more of the subzones
hb_ind > 1000 - bisection solver was used on one or more of the subzones
        (number of iterations of bisection solver = hb_ind - 1000)
hb_ind = 9999 - Hoek-Brown model did not converge

================================================================================
472     08/13/08
================================================================================

In the FISH CASE statement, the maximum range has been increased to 1000, from
255.
E.g., the following is now supported
CASEOF expression
  CASE 0
  ...
  CASE 999
  ; some code
  CASE 1000
  ; some code
ENDCASE

The ubiquitous joint model when run in large strain mode updates the joint dip
and dip direction properties (rigid body rotation correction).  However, the
updated dip and dipdir were stored in radians instead of degrees.

================================================================================
473     09/02/08
================================================================================

The command PRINT INTERFACE n TARGET outputs information about an interface
node's target face.  The last column that is output by this command 
(labelled update) should reflect whether or not the interface's update option
has been set.  The update keyword on the interface command allows (update on) or 
prevents (update off) the search for new contacts after movement occurs on an
interface.  The status of the update switch was displayed incorrectly (No was
displayed instead of Yes).  This has been fixed and changed to display On/Off 
instead of Yes/No.  On/Off is now displayed in the update field to indicate
update on/off respectively.

System commands run from FLAC3D would sometimes produce error messages when none
were actually generated.

================================================================================
474     09/04/08
================================================================================

FLAC3D would occasionally crash if a plot containing zones and gridpoints was
displayed and zones or gridpoints were deleted.

Added FISH function z_apply(zptr, iface, sAppItem, sOption)
zptr = pointer to a zone
iface = integer index to a face on the zone (1-6)
sAppItem = string, one of 'sxx','syy','szz','sxy','syz','sxz',
                          'dstress','sstress','nstress'
sOption = string, one of '', 'add', 'mul', 'remove'
z_apply can be used to apply stresses to a zone face.  It must be used on the
left side of an equals sign.  
Example, 

z_apply(zptr, 4, 'szz', '') = -2.0 
applies an szz stress of -2.0 to face 4 of zone pointed to by zptr

z_apply(zptr, 4, 'szz', 'add') = -7.0
adds -7.0 to szz on the same zone face as above

z_apply(zptr, 4, 'szz', 'mul') = 2.0
multiplies the currently stored applied value of szz by 2.0

z_apply(zptr, 4, 'szz', 'remove') = 0.0
removes remove the szz apply item 

================================================================================
475     10/17/08
================================================================================

Fixed FLAC3D crash which occasionally occurred during making movies.

System commands would sometimes not wait for the process to complete before
returning control to FLAC3D.

Added logic to account for direct generation of pore pressure increment by a
UDM (using state variable dPPInc).

Hysteretic damping bug fixes when NMD is on.  
- If hysteretic damping was used with a model that was run on more than one
processor then FLAC3D could crash.
- Put a lower cap on the hysteretic damping multiplier of 1e-3 to avoid shear
modulus from becoming zero.

================================================================================
476     10/28/08
================================================================================

Added FISH function z_iecreate to create interface elements on a zone face.
n = z_iecreate(zptr, iface, id)
where
n     = number of interface elements created
zptr  = pointer to a zone
iface = face number on the zone (integer, 1-6)
id    = ID number of the interface to be created (integer > 0)

Modified FISH function z_facenorm which returns the area and normal to a 
zone face as described below:
array idarr[3]      ; contains gp id numbers or a face index
array facenormal[3] ; returns with x,y,z components of a face normal vector
facearea = z_facenorm(zptr, idarr, facenormal)
idarr can contain gridpoint id numbers of 3 unique gridpoints on the zone face
(this is the original implementation).  If idarr[2] = 0 and idarr[3] = 0 then
idarr[1] can contain an index to the face on the zone (an integer between 1 and
6 inclusive).  The latter is the updated functionality.

Added FISH functions for plane orientation
dipddfromnorm(normal, dip_dd) ; calculate dip and dipdir from normal vector
normfromdipdd(normal, dip_dd) ; calculate normal vector from dip and dipdir
array normal(3) ; normal vector (+x,+y,+z)=(East,North,Up)
array dip_dd(2) ; dip and dipdir in degrees
                ; dip measured positive down from horizontal (0 degrees)
                ; dipdir measured positive clockwise from North (0 degrees)

================================================================================
477     11/20/08
================================================================================

Interface creation/delete commands now outputs number of nodes and elements
created/deleted .

Geometric ranges used in interface plotting are now based on interface element
centroids and node positions rather than the average position of the entire
interface.

e.g., plot interface location range id 5 7 x 7.0 9.5
will plot interface elements that have their centroids within x 7.0 9.5
and also belong to interfaces 5 through 7

e.g, plot interface slip shear range id 1 x 7.0 9.5
will plot shear slip state at interface 1, with nodes located between x 7.0 9.5

Inertial density used in the equations of motion did not include fluid 
contributions.

================================================================================
478     12/29/08
================================================================================

FISH function in_tweight(p_in, int) was returning the host weight instead of
the target weight for an interface node.

================================================================================
479     01/21/09
================================================================================

Block contour plots within a user specified min max range would extend the range
to fall exactly on a contour interval.  Now block contour plots will fill
contours within the precise min max range.

The SOLVE FOS procedure always stores the last unstable state it encounters
during its bisection search for the FOS.  If an ESC key is pressed during 
the SOLVE FOS procedure the current state was stored.  This state could be
stable or unstable.  Now, when an ESC key is pressed the current state
is not saved (the last unstable state that was encountered is still saved).

================================================================================
480     01/29/09
================================================================================

Added string FISH functions
newstring = stolower(oldstring) ; return lower case string
newstring = stoupper(oldstring) ; return upper case string
NumChars = slength(string) ; return number of characters in string
newstring = substr(oldstring, iStart, iLen) ; return a substring of oldstring 
; the substring starts at character iStart in oldstring and is of length iLen
; 1 <= iStart <= length of oldstring
; the length of string returned will never exceed the length of oldstring
; if iStart <= 0 or iLen <= 0 a null string is returned

================================================================================
481     04/01/09
================================================================================

If a zero length normal vector is specified in the cutting plane dialog, an 
error message would be generated and FLAC3D would freeze.  Now, if a zero length
vector is entered, the z-component is set to 1.

min, max FISH functions would return an int if called with min(int,float)
or min(float,int).  Now these functions return an int if both arguments are
integers.  A float is returned otherwise.

Bug fix to uwedge gridding primitive.  Generating a grid with the uwedge
primitive would sometimes generate a gridding error message when valid
uwedge parameters were used.

================================================================================
482     04/22/09
================================================================================

If PLOT BLOCK STATE AVERAGE was executed before any cycles were executed or
zones did not have materials assigned then a math error message would appear.
This has been fixed.

For users writing UDMs, increased working array sizes in State data structure
in conmodel.h.  max_working was increased to 50 from 15 and max_iworking was
increased to 10 from 2.

================================================================================
483     04/29/09
================================================================================

When a heat source is specified with FISH, with gp_psource(p_gp), the thermal 
ratio was not updated correctly.  An absolute value of the specified heat
source is now used for the thermal ratio calculation.

================================================================================
484     06/03/09
================================================================================

Optimization to apply command for gridpoint quantities.  Applying gp quantities
to large grids sometimes required an extensive amount of time to complete.

The interface logic assumed that a real fluid model (not null) or a real thermal
model existed on both sides of an interface.  No provision was made for a null
fluid or null thermal model on one side of an interface.  This change allows a
null and non-null fluid (and thermal) model to exist on either side of the
interface as a no-flow or adiabatic boundary.

Added FISH array index boundary checking to highest (rightmost) index of the 
array.

================================================================================
485     07/23/09
================================================================================

Bug fix to GEN SEPARATE command.  Even though dbrick, wedge, pyramid, and tet
zones operate on 7, 6, 5, 4 gridpoints respectively, FLAC3D always stores
8 gridpoints per zone.  For these non-brick type zones, duplicate gridpoints are
stored in the upper indices (e.g., a wedge has 6 gridpoints with gridpoints 7
and 8 being duplicates of gridpoints 4 and 6).  These duplicate gridpoins were
not separated by the GEN SEPARATE command and would cause errors when the
equations of motion were solved if some of the non-brick zones contained a
null material model.

Bug fix to FISH array index boundary checking to highest (rightmost) index of
the array.

================================================================================
486     07/30/09
================================================================================

Fix to timestep in cases where mechanical model coupled with fluid or thermal
models and no master-slave relationship set up. 

================================================================================
487     10/20/09
================================================================================

If zones containing a non-null mechanical model were deleted and those zones
were joined to null zones then a "Zero stiffness in grid-point" error message
would appear.  This has been fixed.

If two zones are connected with the attach logic and one or both zones is set
to a null mechanical model, the attach is now deleted.  Prior to this change
an error message "Attach error, null master zone" would be output when cycling.
Attaches must be recreated if a non-null model is assigned to a null zone
that formerly had an attach.

When a PRINT APPLY command was issued the current value of history type apply
items was not output.  A column called Val*CurHist has been added to the output.
This column represents the value of the applied item using the model time when
the print command is issued.

================================================================================
488     11/05/09
================================================================================

Interface nodes use an area weighting method to distribute stiffness and force
information to the gridpoints on that face.  Small negative weights would
cause gridpoints to acquire small negative inertial stiffness and result in an
error message "Zero stiffness in gridpoint".  This has been fixed.

================================================================================
489     11/12/09
================================================================================

Fluid and thermal logic assumed full-face contact at interfaces.  A single host
face node contacting a target face would result in pore pressure and temperature
initialized to zero on the entire host face face instead of just at the node in
contact.  The code has been changed to only initialize and send information to
individual gps that are in contact on the host side (instead of to the entire
face).  Fluid and thermal logic is only active when the MAXEDGE keyword is NOT
specified when creating interfaces.

Host and target faces are now prohibited from sharing the same gridpoints.  The
host will now exclude these connected faces as potential target faces when
searching for a target.

Three sided host and target faces are now fully supported by the interface
logic.

================================================================================
490     12/18/09
================================================================================

Bug fix in FISH parser.  $ in FISH symbol names was treated as a blank.  Now it
is interpreted as a valid character.

Plotting bug fix.  Applied forces and moments at nodes were plotted even after
these applied forces and moments were removed.  This is a plotting bug only.
The mechanical calculations did not use removed applied forces.

Modification to Hoek-Brown model to address convergence problems when sending
a stress point that is near the apex to the yield surface.

================================================================================
491     02/01/10
================================================================================

Bug fix to attach logic.  Attach command would attach more than once if attach
face command was run more than once with the same range.

================================================================================
492     02/08/10
================================================================================

Copyright symbol changed to (c) in title.

================================================================================
493     02/10/10
================================================================================

Bug fix to attach logic.  Attaches were missed.  Attach command now runs
multiple passes up to a maximum of 5 passes or until no more attachments
created.

================================================================================
494     10/04/10
================================================================================

Fixed translational and rotational velocities were not considered when a node
local coordinate system rotated during a large strain update.  The following
changes have been made to fix this:
      
The SEL NODE INIT command now has a LOCAL keyword allowing you to initialize
nodal velocities (xvel, yvel, zvel, xrvel, yrvel, zrvel) in local node
coordinates as opposed to global coordinates.  
e.g., SEL NODE INIT XVEL LOCAL 0.0
      
Note that the SEL NODE FIX command fixes velocity in the specified directions in
terms of the node-local system for all nodes in the range, thereby preventing
these velocities from changing.  If any translational/rotational velocity
component is fixed, then all three components of the translational/rotational
velocity should also be specified in the local coordinate system via the
SEL NODE INIT ... LOCAL command.

Bug fix to apply items that had gradient keyword, when restoring save files
earlier than revision 487.

================================================================================
495     11/05/10
================================================================================
Bug fix to APPLY XBODYFORCE, YBODYFORCE, ZBODYFORCE command.  The applied
magnitudes were incorrect if more than one component was applied
(e.g., if xbodyforce and ybodyforce were both applied then the value of the
ybodyforce would be the sum of the values of xbodyforce and ybodyforce).

================================================================================

Close This PanelX

FLAC3D 3.00 Modification List  (Type PRINT VERSION to check version number.)
Changes listed are included in  F3300.EXE.

Sub-Version   Date        Description
Number

  251         10/1/04     

              Version 3.00 released as "pre-release" to users.

  252         11/08/04

              Coupled unsaturated fluid flow logic bug fix.  A slow build up
	      of pore pressure occurred at the water table level when the
	      water table level corresponded exactly to the position of a
	      gridpoint.

  254         12/01/04

              nx, ny, nz properties were ignored when they were being set in
	      the anisotropic and orthotropic models.  Fixed an error in the
	      orthotropic model rotation matrix.

  255         01/24/05

              When a component of velocity was removed at a gridpoint (e.g.,
	      "apply remove zvel") and other applied velocity components
	      existed at that gridpoint then a fatal error would result when
	      cycling.

  256         02/21/05

              With the radcylinder gridding primitive, when an odd number of
	      zones was specified for the 3rd size parameter and the fill
	      option was used (e.g., gen zone radcyl size 2 2 3 2 fill) then
	      the zones would not match between the radcylinder opening and
	      the filled cylinder.

              Added z_frr and z_fri for FISH access to the full rotation rate
	      and full rotation increment tensors (equation 1.50 of the v2.10
	      Theory and Background manual).

  257         03/14/05

              Pile elements did not do bending yielding calculations unless
	      the rockbolt property was turned on. Piles now always perform
	      the bending yielding calculations regardless of the rockbolt
	      property setting.

  258         03/31/05

              AVI file is closed properly if user quits FLAC3D or types a
	      "new" command while an AVI file is open.

              If a temperature or pore pressure history of a gridpoint
	      specified with x y z coordinates was taken, and that gridpoint
	      fell in a null zone (null mechanical constitutive model), then
	      that gridpoint would not be found even though that gridpoint
	      belonged to a non-null fluid or thermal zone.

              Cable yield plot items were not always created or saved
	      properly.  Users may have problems restoring sav files if they
	      contained cable yield plot items.

              When FISH I/O was used to read an ASCII file and the user
	      requested a single line to be read then the function would
	      return 0 if an EOF was encountered.  A value of 1 is now
	      returned if an EOF is encountered in this case.

  259         06/22/05

              Added "gen separate groupname" command to allow a zones to be
	      unmerged from the rest of the model.

              User specified rgb triplets in plotting commands were ignored.

              FISH function redefinition bug fixed.  E.g., when a user defined
	      FISH functions FA and FB, where FB calls FA and then FA is
	      subsequently redefined then the call to FA from FB would no
	      longer work.

              Added "plot extract viewname _xarr _yarr _zarr _datarr" command
	      to extract data from a contour plot into fish array _datarr.
	      _xarr _yarr _zarr are fish arrays containing the x, y, z
	      coordinates of the extraction points.

              A cut plane's move increment is now decoupled from the camera's
	      move increment.  A plane's move increment can be set through the
	      "plot set plane moveincr value" command or through the
	      Edit-->Cutting Plane dialog.

              When double clicking on a face, a message box is now displayed
	      with zone information and mouse coordinates instead of
	      displaying the data at the bottom of the plot.

              The plot item caption in the legend is usually set to the color
	      of the plot item.  The caption color is now set to the
	      foreground color if the caption color is the same as the
	      background color.

              Added color dialog for choosing outline color on surface plots.

              Attach face creation has been improved.

  260         07/08/05

              Added "remainder" keyword to group command.
              group remainder  
              Assigns group  to all NULL group zones (zones that don't
	      have a group assigned to them) that also fall with the range.

              group none  unassigns groups.

              Added z_facenorm() FISH function.
              E.g.
              array _gpid(3)
              array _norm(3)
              _area = z_facenorm(zp, _gpid, _norm)
              Returns the area and normal to a zone face that is specified by
	      3 gridpoint ids.

              Fixed bug with "plot block state {any, average}".  Keywords any
	      or average were not recognized.

              Added "delete volumes" command to delete all volumes, surfaces
	      and points.

              Optimizations to solving equations of motion.

  261         08/19/05

              Memory leaks fixed in fluid and thermal models.

              Bug fix when structural elements used with creep material
	      models.  Scaled masses for structural element nodes could end up
	      being zero (an error) depending on the creep timestep.

              Bug fix with grayscale history plots.  FLAC3D would sometimes
	      hang (go into an endless loop) when plotting histories in
	      grayscale.


  263         09/06/05

              Adoption of new security logic - provides access to normal and
	      network keys in a single executable.


  265         10/10/05

              Fixed problem with new security logic and older (Win 98, etc)
	      versions of windows.

  266         10/11/05

              Increased efficiency of z_model() fish function if neither old
	      or new model is null.


  267         10/21/05

              Plotting bug fix.  Plotting with user specified colors would
	      sometimes cause a fatal exception if the plot item's modify
	      dialog box was instantiated.

              Attach logic bug fix.  If two zones were attached and the master
              was made null mechanically, the attached gridpoints on the slave
              zone would not have equations of motion calculations done.
              Now, When the master is made null the slave gridpoints have
              their eoms restored.

  268         11/11/05

              Bug fix to allow printing FISH arrays containing strings.

  269         11/18/05

              In the Help->About dialog box, the user serial number displayed
	      as zeros. Customer titles were not being read from the registry.

  270         12/02/05

              Enabled ranges on vector plots on plane (zone based).

  271         12/02/05

              Plotting bug fix.  FLAC3D would sometimes freeze when rotating a
	      plot. Memory leak fix with histories.

  272         12/21/05

              Default creep timestep did not always default to the minimum
	      default value of 1e-20.  mindt and maxdt in the set creep
	      command now only apply when automatic timestepping is on (set
	      creep dt auto on).

              Blank spaces allowed in group names when importing a 3DShop
	      file.  When blanks are included in the group name then the name
	      must be surrounded by either single or double quotes.

              Bug fix to Hoek-Brown constitutive model.  If a stress point is
	      outside the yield surface, the Hoek-Brown model sometimes failed
	      to return the point to the yield surface (an error message would
	      result).

  273         01/11/06

              Fluid discharge vector calculations are now done during cycling
	      if particle tracking is on. Before this change it was assumed
	      that fluid flow was already in steady state and the discharge
	      vectors did not change during cycling (the discharge vectors
	      were only calculated once prior to cycling).

  278         01/24/06

              The name of the general concrete hydration model was changed
	      from th_hyd to th_hyd_general.  The general hydration model is
	      now used when specifed as "model th_hyd".  Before this change
	      the th_hyd_concrete1 model was used when "model th_hyd" was
	      specified.

              The loads that are applied with FISH gp_xload,gp_yload, and
	      gp_zload are now applied during each step and not at the
	      initialization phase.

              Bug fix in large strain update of particle positions (when
	      tracking particles in a flow field).  Before this, the large
	      strain update could not always calculate the new particle
	      position and an error message was displayed.


  279         02/06/06

              Bug fix in large strain update of particle positions (when
	      tracking particles in a flow field).

  280         02/22/06

              Bug fix to confined and bulk modulus calculation in anisotropic
	      model.  These values are only used in stability calculations and
	      not the stress calculations.  The orthotropic model is also
	      affected by this change.

              Reversal detection bug fix in hysteretic damping.  Strain
	      increments multiplied by a factor of two so 2D and 3D hysteretic
	      damping formulations would produce the same results.

              Added gp_fload FISH function to add a value to the gp out of
	      balance flow.

  281         02/23/06

              Bug fix to gen zone command.  Command would sometimes indicate
	      that the mesh primitive did not conform to the node numbering
	      convention when it actually did.

  282         02/24/06

              Enhancements to support orthotropic shells where the principal
	      directions of orthotropy do not coincide with the SEL local
	      coordinate systems.

  283         02/27/06

              Only a partial property name is required to match a constitutive
              model property in FLAC3D.  If more than one match is found then
	      FLAC3D will now return the shortest property name that matches.

  284         03/07/06

              For logarithmic plots, if values <= 0 are specified as plot
	      limits or as plot data then these are set to 1.0e-100 for
	      plotting (i.e., they will be plotted as value = -100 =
	      log10(1.0e-100)).

  285         03/08/06

              Bug fix in interface property FISH function in_prop.

  286         03/10/06

              Bug fix in property command, if more than one matching property
	      found. Bug fix in fluid flow if attach logic was used.

  287         03/11/06

              Bug fix in property command.  In some cases, matching properties
	      could not be found.

              03/13/06

              Added note about Windows Remote Desktop to security dialog.

  288         03/15/06

              FISH function z_thmodel would always return undefined even if a
	      thermal model was defined for a zone.

  289         03/15/06

              Bug fix in property command.

  290         03/28/06

              Bug fix in shell elements.  Default properties weren't being
	      initialized properly.

  291         03/31/06

              Plotting bug fix with fluid thermal model overlays.  If a fluid
	      or thermal plot was created prior to cycling, then the fluid or
	      thermal overlays were not always present and a divide by zero
	      error would result.

  292         04/26/06

              Extended orthotropic shells to allow specification of separate
	      bending and membrane material-stiffness coefficients.

  293         04/27/06

              Displacement contour plots of shell-type structural elements
	      were not being filtered by the specified range and as a result
	      displacements for all shells were displayed.

  294         05/04/06

              CSHELL primitive didn't generate its gridpoints properly unless
	      user specified control points formed an orthogonal system.

  295         05/04/06

              Added keyword "sweepaxis" to the generate zone command which
	      informs the grid generator to use control points P0 and P2 as an
	      axis of rotation for cylindrical types of primitives (cylinder,
	      cshell, radcylinder).  This keyword is ignored for all other
	      primitives.

  296         05/08/06

              Sweepaxis option (above) didn't work properly with elliptical
	      arcs.

  297         05/10/06

              Optimization to grid generation of cylindrical type primitives.
	      Allow expgrid to support group names with spaces (group names
	      are surrounded by single quotes).

  298         05/18/06

              Stress tensor plots were always drawn with respect to an
	      undistorted grid (even if the grid was deformed in large strain
	      mode or if a magnification factor was specified for the plot).

  299         06/09/06

              "plot flow" command now accepts a range.

              When free-fields are created, the gridpoints in these free-field
	      zones are merged if they are within a certain tolerance.  This
	      tolerance value has been lowered to accomodate finer
	      discretizations of the model geometry.  Also, a the free-field
	      zones are now subjected to a geometry test.  Before this it was
	      possible for the free-field to contain invalid zones.

   300        07/10/06

              Shear modulus was not being adjusted when hysteretic damping was
	      used with a double yield constitutive model.

              Bug fix to anisotropic constitutive model.  When the components
	      of the normal vector describing the orientation of the
	      anisotropy (dNx, dNy, dNz) are set in the setPropertey method a
	      new dip and dip direction were calculated but were not being
	      converted back to degrees.  Internally everything is done in
	      radians but dip and dd are stored in degrees and the user enters
	      degrees.  This would cause a problem if the user read the dip or
	      dip direction or set the dip or dip direction after setting one
	      of the normal vector component.

   301        07/19/06

              Plot extract command would generate an error if fluid or thermal
	      modules were not active.

              A FISH symbol, which is not an array, would not generate an
	      error message if used with a set of parenthesis on the rhs of an
	      equation.  E.g.  xxx = gp_pqrstuv(p) would interpret gp_pqrstuv
	      as a symbol (which has an initial value of zero).  gp_pqrstuv is
	      not an intrinsic FLAC3D FISH function.


   302        07/24/06

              Bug fix when describing a plane's orientation with a normal
	      direction when used for plotting or constitutive model
	      properties. The conversion to dip and dip direction was
	      incorrect in the following cases: When the normal has a
	      y-component == 0 then the dip direction should be +/- 90 degrees
	      depending on the sign of the x component. When the normal has a
	      z-component == 0 then the dip angle should be 90 degrees.

              When a free field grid is created from a grid that has a fluid
	      and thermal model the free field did not get fluid or thermal
	      overlays allocated.

   303        08/08/06

              Updated to new security logic, requires key driver version
	      7.3.0. Windows Remote Desktop should not interfere with direct
	      key detection anymore. Fix to PLOT EXTRACT command when null
	      zones are present.

   304        09/15/06

              Increased the ASCII line length supported by FISH I/0 (read and
	      write) from 100 characters to 255 characters.

   305        10/10/06

              Cut plane contour interpolation was based on the cut plane's
	      position in the undeformed grid when a magnification factor was
	      used in the plot.  Cut plane contours are now based on the
	      plane's position in the deformed grid when a magnification
	      factor is used.

   306        12/11/06

              Bug fix in sel property command for orthotropic and anisotropic
	      shells, liners, geogrids. If the default keywords are not
	      specified the parser generates an error message.

              Poor geometry encountered when importing a grid now causes
	      warnings instead of errors.

   307        02/01/07

              Hoek-Brown model no longer forces a stress point back to
	      Hoek-Brown failure-surface if the iterative algorithm fails to
	      converge.  An error message is output instead.

   308        08/01/07

              Bug fix in recording more than one history for recovered
	      stresses in structural elements. e.g., HISTORY SEL RECOVER SRES
	      SURFX ... If more than one history was created, then in large
	      strain mode it was possible to get a "surface system not set at
	      node" error.

   309        09/13/07

              Cylinder and cshell mesh primitives did not work properly if
	      sweep angle was greater than 90 degrees.

   310        09/28/07

              The error message
              "Fluid or thermal zone geometry must support two overlays."
	      would appear when trying to add a free-field to a grid
	      containing a non-null fluid model.

Software: FLAC3D

Updates

FLAC3D 4.00.82

Download Update

FLAC3D 4.00.82

Download Update

FLAC3D 3.10.495

Download Update

FLAC3D 3.10.495

Download Update