| ²é¿´: 1497 | »Ø¸´: 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¿Í»§¶Ë |
» ²ÂÄãϲ»¶
ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£©
ÒѾÓÐ12È˻ظ´
µ÷¼ÁÇóÖú£¨ÉúÎïÓëÒ½Ò©£©
ÒѾÓÐ3È˻ظ´
Ò»Ö¾Ô¸¼ª´ó»¯Ñ§327Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
Ò»Ö¾Ô¸ÉúÎïÓëÒ½Ò©£¬296·Ö£¬Çóµ÷¼Á
ÒѾÓÐ13È˻ظ´
327Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
³õÊÔ301£¬´úÂë085701»·¾³¹¤³Ì£¬±¾Ë¶Ò»Ö£¬ËÄÁù¼¶Òѹý£¬ÓжþÇøÒ»×÷£¬¹²·¢±í5ƪÂÛÎÄ
ÒѾÓÐ8È˻ظ´
327¿¼Ñе÷¼ÁÍÆ¼ö
ÒѾÓÐ6È˻ظ´
328Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
304Çóµ÷¼Á
ÒѾÓÐ9È˻ظ´
»¯¹¤Çóµ÷¼Á£¡
ÒѾÓÐ16È˻ظ´

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














»Ø¸´´ËÂ¥