| ²é¿´: 1496 | »Ø¸´: 1 | |||
zhangyanzhiľ³æ (СÓÐÃûÆø)
|
[½»Á÷]
lammpsÇóÆøÒº½çÃæµÄ±íÃæÕÅÁ¦ ÒÑÓÐ1È˲ÎÓë
|
|
ÇóÖú´óÉñ£¬ÈçºÎ¼ÆËãÆøÒº½çÃæµÄ±íÃæÕÅÁ¦¼°Æä·Ö²¼¡£ÏÂÃæÊÇÎÒinÎļþÖеÄÏàÓ¦²¿·Ö£¬µ«×ÜÊǼÆËã²»³öÎÄÏ×ÖбíÃæÕÅÁ¦µÄ·Ö²¼¹æÂÉ¡£ #--------------------calculate surface tension----------------- variable N equal 100 variable box equal (xhi-xlo)*(yhi-ylo)*(zhi-zlo) # volume of the box #--------------------surface tension(x-direction)-------------- compute p all stress/atom NULL compute k all reduce sum c_p[1] c_p[2] c_p[3] variable pn atom -c_p[1] variable pt atom -(c_p[2]+c_p[3])/2.0 variable rk atom ((v_pn-v_pt)/(v_box/v_N))*((xhi-xlo)/v_N)*1.0e-2 # unit mN/m #--------------------surface tension(global)------------------- variable tension equal 0.5*(xhi-xlo)*(c_k[1]/v_box-(c_k[2]+c_k[3])/(2.0*v_box))*1.0e-2 # unit mN/m fix mean all ave/time 100 5 1000 v_tension mode scalar ave one file tension.dat СµÜÔÚ´ËÆÚ´ý´óÉñµÄ°ïÖú! ·¢×ÔСľ³æAndroid¿Í»§¶Ë |
» ²ÂÄãϲ»¶
277Çóµ÷¼Á ÊýÒ»104·Ö
ÒѾÓÐ7È˻ظ´
288»·¾³×¨Ë¶,Çóµ÷²ÄÁÏ·½Ïò
ÒѾÓÐ16È˻ظ´
277¹¤¿ÆÇóµ÷¼Á
ÒѾÓÐ8È˻ظ´
266Çóµ÷¼Á
ÒѾÓÐ19È˻ظ´
»·¾³×¨Ë¶µ÷¼Á
ÒѾÓÐ3È˻ظ´
081200-11408-276ѧ˶Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
²ÄÁϵ÷¼Á
ÒѾÓÐ13È˻ظ´
¹¤¿Æ 22408 267ÇóÍÆ¼ö
ÒѾÓÐ4È˻ظ´
305Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
270·ÖÇóµ÷¼Á
ÒѾÓÐ5È˻ظ´

awuliÕÅ
гæ (СÓÐÃûÆø)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 253
- É¢½ð: 30
- Ìû×Ó: 97
- ÔÚÏß: 22.5Сʱ
- ³æºÅ: 15793324
- ×¢²á: 2019-07-31
- רҵ: »¯Ñ§¹¤³Ì¼°¹¤Òµ»¯Ñ§
2Â¥2019-10-18 18:42:12














»Ø¸´´ËÂ¥