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)参数,用以输入泊松比;
Attachments:
Download this file (return-mapping-iteration-closest-point.pdf.pdf)return-mapping-iteration-closest-point.pdf.pdf[图5.4 返回映射法的迭代过程]张雄0 kB

用户登录