24小时热门版块排行榜    

查看: 921  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

truemantwo

木虫 (正式写手)

[交流] 【求助】Discover模块里如何计算材料的弹性模量?

我想用Discover计算一下铁的弹性模量,把MS自带的模型文件Fe.msi导入,而后再选择Discover -> Analysis -> Mechnical -> Elastic properties (static),选择“Current structre in Fe.xsd”,计算完成后,在Fe.out文件里找不到与弹性模量有关的信息,怎么回事呢?当然,这个问题CASTEP可以做,但用Discover可否搞定?

还有一个,用MS的Discover模块做高分子链的MD模拟,想在链末端的碳原子上加外力,不知能否做到这一点?

期盼各位高手出手相助啊,谢谢!
回复此楼

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

仿真建模与计算 程序语言箱

» 猜你喜欢

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

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

abinitio

版主 (著名写手)

用discover和castep的dynamics计算有什么不同吗?
5楼2008-04-23 16:15:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

truemantwo

木虫 (正式写手)

我把.inp文件和.out文件贴出来,大家看看是哪里的问题。
.inp文件:

#BIOSYM btcl 3
#
#  Input File For Discover Generated By Materials Studio
#  Input Client Model Document: D:\1--Materials_Studio_Results\Exercise_Discover_Analysis\FeYoungsModulus\First Files\Documents\Fe.xsd
#  Job: [ZCH78] - Fe Disco  Elastic properties (static)
#
  autoEcho off
#
# Begin Forcefield Section
  begin forcefield = compass
#
# Nonbond section:
  forcefield nonbond \
      +separate_coulomb \
    vdw \
    summation_method = atom_based \
      cutoff = 9.50 \
      spline_width = 1.00 \
      buffer_width = 0.50 \
    coulomb \
    summation_method = ewald \
      ewald_accuracy = 1.e-002 \
      update_width = 1.00 \
    dielectric_value = 1.0000
#
# vdW long range energy correction
  vdWTailCorrection 9.50
# End Forcefield Section
#
analysis \
         method = Elastic_Static
3楼2008-04-23 15:37:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

truemantwo

木虫 (正式写手)

.out文件(省去了一些没用的东西):

randomSeed is set to 849775

Line   7:BTCL> autoEcho off

INPUT FILES
___________

File Type                Name
---------                ----

Forcefield                C:\PROGRA~1\Accelrys\MATERI~1.1\Discover\..\Data\Discover\res/compass.frc
Molecular data                Fe.mdf
Coordinate                Fe.car

Periodic Boundary Conditions
____________________________

                Length (A)                Angle (degrees)
                ----------                ---------------

                a  8.59920                alpha  90.00000
                b  8.59920                beta   90.00000
                c  8.59920                gamma  90.00000


MOLECULAR TOPOLOGY
__________________

Number of molecules:                     1
Number of residues:                     1
Number of atoms:                    54 (asymmetric unit: 54)
Number of atom types:                     2

FORCEFIELD OPTIONS
__________________

Filename                  : compass.frc
Definition name                  : cff91
Version                          : 2.7
Last modification date          : 7/28/2006
# of automatic parameters : 0

NONBOND ENERGY CUTOFFS
______________________

                Cutoff (A)        Spline Width (A)        Buffer Width (A)

vdW                  9.50                     1.00                     0.50
Coulomb                   --                       --                       --

NONBOND ENERGY EWALD PARAMETERS
_______________________________

                Accuracy (kcal/mol)        Update Width (A)

vdW                   --                               --
Coulomb                  0.01000000                     1.00

Summation method for vdW        : Atom based
Summation method for Coulomb        : Ewald
Dielectric                        : 1.00
NOTE: Pressure and energy tail corrections will be added to all subsequent energy evaluations


Total time used by DISCOVER: 0.22 secs

Completion date: Tue Apr 22 15:36:15 2008
Exiting Discover: status OK.


好像没有与应力、应变以及模量有关的信息啊,怎么回事呢?
期盼各位高手不吝赐教,谢谢!
4楼2008-04-23 15:40:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

truemantwo

木虫 (正式写手)

呵呵,当然没有什么不同。
我做的是高聚物的力学性能,模型里的原子比较多,用CASTEP算不过来,所以才想用discover来做。
现在的想法是先做一个简单的,就做铁,把路摸熟再搞高聚物。
6楼2008-04-23 16:43:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见