论著
2014.11.6
第202页倒数第三行最右端项应该由sd(1)改为sd(2)。
2014.5.26
修正第145页的图5.4，新图见附件。
2014.4.4
修正程序的3处bug和4处注释，并进一步做了4处调整。更新后程序可在svn源中下载。（注：此次修改基于清华大学土木工程系孙玉进反馈的问题，在此表示感谢！）
- Bug 1：计算摩擦接触力的代码有误，对应的代码在 subroutine Lagr_NodContact()，书稿中为P98，第18行。正确的 代码为 cforce=nomforce*CP1%ndir + val_fslip*(fstick/val_fstick)
- Bug 2：计算背景网格结点质量截断值存在量纲问题，代码位于 DataIn module中的subroutine SetMaterial()，正确的 代码为CutOff=CutOff*Dcell**3*1e-5/matcount.
- Bug 3：在DataIn module中的 subroutine statinfo() 计算初始动能时丢掉了系数 0.5，正确的代码为 cke=cke+dot_product(pt%VXp,pt%VXp)*pt%mass*0.5d0
- 注释1：Particle.f90文件中的 type Body，第二行，nb_mat 改为 nb_component; （对应的书P80第11行）
- 注释2：在MPM3D.f90文件中Step 4的注释，删除 “and integrate momentum euqation”，对应书 P88
- 注释3：在update_step.f90文件中的subroutine ParticleStressUpdate() 程序内，计算应变率D和旋率张量W时，删除其后注释中的*DT，对应书P92的第2-9行。
- 注释4：在Grid.f90文件中的subroutine SetGridData()倒数第2行注释改为 !iJacobi4 = iJacobi/8
- 调整1：将计算步骤Step 2调整为指针对USF格式，即将其与Step 3 同时执行；
- 调整2：在MPM.f90文件中的倒数第3行中的“zone” 改为 particle
- 调整3：书P90，第13行， CellsNode(icell,1) 替换为 InflNode(1). 并修改相应的代码；
- 调整4：书P92，倒数第12行， grid_list(comdID,InflNode(n)) 替换为gd，并修改相应的代码。
2014.1.31
- 第252页：该页的输入文件范例存在几处笔误。正确的范例文件TNT3D.mpm可在MPM3D-F90的svn源里下载，地址为：http://svn.comdyn.cn/opensource/trunk/MPM3D-F90；
2014.1.6
- 第58页：将第3.5.1节第1段第3行至第4行的\varGamma_{c}改为\varGamma^{C}，共6处；
- 第182页：将式(5.259)中的V改为1/V；
- 第186页：将式(5.282)第一行右端最后一项中dp^3改为d^3p；
- 第187页：将式(5.288)第一行右端的s_1改为s；
- 第190页：删除式(5.300)第二行的第一个“(”；
2013.12.9
- 第250页：关键字point后面的两个关键字应改为 comID num(i)；
- 第250页：倒数第三行的关键字pmass(r)改为density(r)；相应的将倒数第1行中的“单个质点的质量pmass” 替换为“初始密度density”；
- 第251页：在第1行末尾处添加“因此该离散体各个质点均为density×dp3”；
- 第251页：将球体建模的关键字pmass，及文字介绍同上修改；
2013.11.13
- 第247页：关键字john、sjc和sjcf后面的参数列表中，在E(r)后面都需添加v(r)参数，用以输入泊松比；
- 第248页：关键字jcf后面的参数列表中，在E(r)后面需添加v(r)参数，用以输入泊松比；