用于离散元计算的PFC软件

fname: pfc310pre.TXT PFC2D/PFC3D (Version 3.10, pre-release) PUBLIC MODIFICATIONS LIST Type PRINT VERSION DEVELOPMENT to display version information. ======================================================================== Version Date Description 3.10-158A 07/17/02 (1) Version 3.1 pre-released (ALPHA version) (-067-157-158) (2) Add periodic logic (v.3.0 also) 3.10-159A 08/01/02 (1) Add & modify periodic logic, more (v.3.0 also) (-072-157-159) 3.10-161A 08/20/02 (1) Add & modify periodic logic (v.3.0 also) (-074-159-161) (2) Add thermal logic (v.3.0 also) 3.10-170A 03/07/03 (1) BUG fixes (v.3.0 also): (-01_083-166-170) (a) parallel option: communication buffer overflow & interface region balls' detection (b) CPPUDM option (3D only): data "fb.ptt_dot" communication between PFC and DLLs (c) parallel option: data with EXTRA ARRAY (2) Change the way of Registry inq. (HKEY_CURRENT_USER to HKEY_LOCAL_MACHINE & CFG_INFO) (v.3.0 also) (3) Add FISH function c_gobj1(cp), c_gobj2(cp), c_go1clist(cp)and c_go2clist(cp), these are the same as c_ball1(cp), c_ball2(cp), c_b1clist(cp) and c_b2clist(cp) (v.3.0 also) (4) Add no-range ball generation in GEN command (GEN id 1 10 rad 1 1 x 0 10 y 5 5 z 0 10) (v.3.0 also) (5) Update pbrick logic for Linux PFC (Linux PFC only) 3.10-173A 04/10/03 (1) Add GUI for thermal option (temp., pipe, power and pob) (v.3.0 also) (_01-085-169-173) (2) Add functionality for WALL to have temperature (thermal option) (v.3.0 also) (3) BUG fix: (ver.3.0 also) (a) plotting problem with deleted line and flat walls (ver.3.0 also) (b) retrieve SAV file created a154 or before (4) Update pbrick logic 3.10-175A 04/24/03 (1) Change Command Syntax: SET DAMP VISCOUS NOTENSION ON/OFF: ON/OFF (_01-085-171-175) is needed after alpha version (2) BUG fix: (a) timestep calculation of Hertz contact model with viscous damping (b) 'print contact prop' Command for Hertz contact model 3.10-178A 06/09/03 (1) Add 'os' as a FISH Variable (_01-086-173-178) return 1: Windows 2: Linux 0: otherwise (2) Update pbrick & MPI logic (3) BUG fix: (a) 'wall ?' command (b) table logic & 'plot' command on Linux 3.10-179A 06/14/03 (1) BUG fix: (ver.3.0 also) (_01-086-174-179) (a) Motion calculation of volumetric center by angular velocity (general wall logic 2/3D) (b) Retrieve SAV file created 154A or before 3.10-180A 06/20/03 (1) Make another PFC2D/3D by intel compiler as well as standard. Check PRINT VER DEV. -c0 standard(VC++), -c1 intel 3.10-181A 06/27/03 (1) BUG fix: (a) the way of Registry inq. (HKEY_CURRENT_USER to HKEY_LOCAL_MACHINE & CFG_INFO) 3.10-182A 07/11/03 (1) Make PFC2D/3D by VC7 (.NET) (_01-089-175-182) (2) Update pbrick logic 3.10-183A 08/04/03 (1) Update pbrick logic (including update MPICH 1.2.5) (_01-090-175-183) (2) Add command "SET DAMP VISCOUS MODE 0(default)/1" and "DAMP VISCOUS MODE 0(default)/1" and the mode functionality in FD law (see below) --- Adding the viscous damping term in the PFC2D/3D contact model --- Putting the damping viscous term on contact models has two ways; mode 0 and 1.Users can set these by the command. The default is zero. SET DAMP VISCOUS MODE 0(default)/1 <range> The mode is set to 0 or 1 for all newly created contacts. DAMP VISCOUS MODE 0(default)/1 <range> The mode is set to 0 or 1 for all contacts in the range. Below is the mode description. Mode 0 (default): The viscous damping force is added AFTER checking contact bond breaking criteria and frictional slip criteria. Viscous damping in this case applies to the built-in contact models (the Linear contact model or Hertz-Mindlin contact model)or the user-defined contact models. Mode 1: The viscous damping force is added BEFORE checking contact bond breaking criteria and frictional slip criteria. Viscous damping in this case applies to the built-in contact models (the Linear contact model or Hertz-Mindlin contact model). NOTE that users who want to make a user-defined contact models with a viscous damping term within their code should include the viscous damping term in the calculation of contact shear force and set it to variable "cs_force" in FdBlock (in the updated contmodel.h). Also, mode 1 ("SET DAMP VISCOUS MODE 1" or "DAMP VISCOUS MODE 1") should be set o avoid extra damping forces added externally to the user-defined model. (3) Add solve ratio in cycle information during SOLVE. 3.10-185A 09/02/03 (1) Add information of Maxfric(on/off) from Print Wall Prop command (_01-092-177-185) (2) BUG fix: (a) Wrong message in PRINT CELL after SAVE RESTORE and BALL command (b) Retieve 0 for User-defined object pointer in extra array (CPPUWC option) (c) Fatal error with ductile model (user-defined contact model) (3) Add Mag(Udvect &) to get magnitude of Udvect class (CPPUWC option) (4) Add Udata::pointer_type() = 105 for User-defined object (CPPUWC option) (5) Add Udata::Uremap(Uobj *uop) function for User-defined object (CPPUWC option), see example below class Uobj_2 : public Uobj { }; class Uobj_1 : public Uobj { private: Uobj_1 *uobj1; Uobj_2 *uobj2; // .... }; void Uobj_1::read(Ufile &fi) {{ uojb1 = (Uobj_1 *)Udata::Uread(fi,uobj1); uobj2 = (Uobj_2 *)Udata::Uread(fi,uobj2); //.... } void Uobj_1::remap(void) { uojb1 = (Uobj_1 *)Udata::Uremap(uobj1); uobj2 = (Uobj_2 *)Udata::Uremap(uobj2); // .... } (6) Update pbrick (7) Add RESTORE fname <mpi_append=on/off> (8) Add "line3d", general Wall of line segment in 3D 3.10-187A 09/18/03 (1) Add FISH specification statements MATER_CONTROL (_01-094-177-187) and SHARE feature in MPI (2) Add new FISH function to install a parallel bond at a contact: pb = c_installpb(cp) ,where cp is the contact pointer and pb is the returned parallel bond pointer, or null if the bond cannot be created. A parallel bond will be installed in an existing contact even if the 2 balls are not touching or if the normal force is zero; it is up to the user to ensure that contact conditions are meaningful. The default bond radius is 0.25, and all other prop