24小时热门版块排行榜

返回列表

【奖励】本帖被评价15次，作者cmf0001增加金币 13 个

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cmf0001

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 126.1
帖子: 84
在线: 17.1小时
虫号: 302057

[资源] 【在线答疑】GULP 计算软件使用经验交流

我使用GULP软件做晶体性质计算工作两年多，并查阅了大量GULP软件相关计算文献，计算材料的各种物理和化学性质非常理想！希望与大家作进一步的交流，本人有相当的使用经验，基本掌握各种命令的使用方法！希望大家能就GULP软件的使用情况在这里讨论！
我已基本完成该软件说明书的翻译工作，希望以后对各位能有帮助！

[ Last edited by lei0736 on 2009-12-21 at 18:06 ]

回复此楼

» 收录本帖的淘帖专辑推荐

杂七杂八

» 猜你喜欢

交叉科学部支持青年基金，对三无青椒是个机会吗？已经有7人回复
青椒八年已不青，大家都被折磨成啥样了？已经有15人回复
多组分精馏求助已经有4人回复
免疫学博士有名额，速联系已经有4人回复
国家基金申请书模板内插入图片不可调整大小? 已经有6人回复
国家级人才课题组招收2026年入学博士已经有5人回复
Fe3O4@SiO2合成已经有6人回复
青年基金C终止已经有4人回复
26申博求博导推荐-遥感图像处理方向已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2009-01-20 16:08:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

swf1632001

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 155
帖子: 49
在线: 2.2小时
虫号: 648031

Material studio中如何输出某个原子在某刻的速度和力？

各位好！
我使用的是Materials Studio中的Gulp模块，在模拟结束后生成的output文件有一个“.xtd”文件，就是trajectory file，我想请教各位如何将某个原子在不同frame时对应的速度velocity、力force求出来？
我知道在Materials studio中有对应的Script，可以编写一个简单的Script进行连续输出，但是当我编写如下时：
#!perl

use strict;
use MaterialsScript qw(:all);

#define the input document and create the study table
my $doc=$Documents{"SiO2_Cylinder(10,100,110)SPY90_2.xtd"};
#create a new study table
my $newStudyTable=Documents->New("ForceXYZ.std"

;
my $calcSheet=$newStudyTable->ActiveSheet;

#creat the column heading in the study table
$calcSheet->ColumnHeading(0)="ForceAtomXYZ";

$calcSheet->ColumnHeading(1)="AtomicIndexNumber";
$calcSheet->ColumnHeading(2)="Mass";
$calcSheet->ColumnHeading(3)="MassNumber";
$calcSheet->ColumnHeading(4)="ElementName";
$calcSheet->ColumnHeading(5)="Force";
$calcSheet->ColumnHeading(6)="Velocity";

my $count=$doc->Atoms->Count;
print "$count";

for (my $i=0;$i<$count;++$i)
{
my $atom=$doc->DisplayRange->Atoms($i);
my $atomicNumber=$atom->AtomicNumber;
my $elementName=$atom->ElementName;
my $mass=$atom->Mass;
my $massNumber=$atom->MassNumber;
my $elementName=$atom->ElementName;

my $x=$atom->X;
my $y=$atom->Y;
my $z=$atom->Z;
my $point=Point(x=>$x,Y=>$y,Z=>$z);
my $vel=$atom->Velocity;
my $force=$atom->Force;
print "$vel";
print "$force";

#my $distance=$x*$x+$z*$z;
#place into study table

$calcSheet->Cell($i,1)=$atomicNumber;
$calcSheet->Cell($i,2)=$mass;
$calcSheet->Cell($i,3)=$massNumber;
$calcSheet->Cell($i,4)=$elementName;
$calcSheet->Cell($i,5)=$force;
$calcSheet->Cell($i,6)=$vel;
#$calcSheet->Cell($i,7)=$x;
#$calcSheet->Cell($i,8)=$y;
#$calcSheet->Cell($i,9)=$z;
#$calcSheet->Cell($i,8)=$distance;
}

$doc->Close;

运行不下去，老提示说，
Variant contains an object that is not a document (writing property "Cell"

at -e line 56.
5140ISDPoint3d [0x116b0c10]

像这种情况，如何将velocity和force通过某方法直接输出到study table中去？或者直接连续读出来？

望请各位赐教，急用。谢谢！