ParaDiS Aniso Code Manuals: Difference between revisions

From Micro and Nano Mechanics Group
Jump to navigation Jump to search
 
 
(9 intermediate revisions by the same user not shown)
Line 1: Line 1:
Anisotropic Elasticity
=DDAniso=
===Sylvie Aubry (11/14/08)===
Sylvie Aubry Jan 28 2010.


The aniso module is similar to DDLab/aniso. For more information refer
This code is being developed in collaboration with Steve Fitzgerald at
to Jie Yin, David Barnett and Wei Cai's paper titled "Computing forces
Culham Science Centre, Abingdon, UK.
on dislocation segments in anisotopic elasticity".


The aniso module contains only the files that are modified from the
source code. Mainly, the aniso structure has to be defined, the
Seg-Seg forces and the Self-Force for anisotropic elasticity are
computed.


Several flags are available:


D_ANISOTROPIC : turns on or off anisotropic elasticity.
The aniso code is a part of ParaDiSv2.2.


You have the option of computing the Q, S, B matrices with or
ANI_Main.c is a copy of Main.c. On top of regular ParaDiS routines, it includes calls to specific anisotropic routines such as
without tables. If you want to use tables, then you have to turn on
the _USE_QSB_TABLES flag. Tables files have to be present in the
inputs directory. To create the tables, turn on the
_CREATE_QSB_TABLES flag (along with the _USE_QSB_TABLES flag). Tables
will be automatically created in the right directory. You can also
specify the name of the table files in the input file. The syntax is
as follows:


ANI_Util.c
This file contains ANI_Init, where specific anisotropic variables are defined such as C12, C44, C[i,j,k,l], Q, S, B.


TableQ character Define the name of the file to read in tables.
tabulate.c
Defaults is inputs/Qtable.dat
This file either create Q S B tables or read from files called Qtable.dat, Stable.dat and Btable.dat. To create the tables files, turn the variable
QSBMETHOD by setting it to 1, to read from existing tables files, set the variable to 0.


TableS character Define the name of the file to read in tables.
NodeForce.c which is called in ParadisStep.c
Defaults is inputs/Stable.dat
TableB character Define the name of the file to read in tables.
Defaults is inputs/Btable.dat


Two methods have been implemented for anisotropic elasticity. The
Forces and
second one is by Steve Fitzgerald at Culham Centre for Fusion Energy
<math>\sigma \cdot b </math>
+44 (0)1235 466350. Turn on _ANISO_Steve to run his method.
are computed using anisotropic routines and not regular ParaDiS ones.


This file calls ANI_Stress which in turns calls the Willis-Steed stress calculations and also SelfForceSPF, which are specific anisotropic self-forces.


A test case isalso included in the aniso module, it is called
aniso_fr_src.*. Forces from the Matlab code are also added for
Other files modifying the current ParaDiD version are:
comparison. The file is compare_aniso_matlab.txt.
LocalSegForces.c

Also several tests files are in the aniso module. They verify several
SegSelfForce.c
components of the implementation.
Param.c

Latest revision as of 19:38, 28 January 2010

Anisotropic Elasticity Sylvie Aubry Jan 28 2010.

The aniso module is similar to DDLab/aniso. For more information refer to Jie Yin, David Barnett and Wei Cai's paper titled "Computing forces on dislocation segments in anisotopic elasticity".

The aniso module contains only the files that are modified from the source code. Mainly, the aniso structure has to be defined, the Seg-Seg forces and the Self-Force for anisotropic elasticity are computed.

Several flags are available:

D_ANISOTROPIC : turns on or off anisotropic elasticity.

You have the option of computing the Q, S, B matrices with or without tables. If you want to use tables, then you have to turn on the _USE_QSB_TABLES flag. Tables files have to be present in the inputs directory. To create the tables, turn on the _CREATE_QSB_TABLES flag (along with the _USE_QSB_TABLES flag). Tables will be automatically created in the right directory. You can also specify the name of the table files in the input file. The syntax is as follows:


TableQ character Define the name of the file to read in tables. Defaults is inputs/Qtable.dat

TableS character Define the name of the file to read in tables. Defaults is inputs/Stable.dat TableB character Define the name of the file to read in tables. Defaults is inputs/Btable.dat

Two methods have been implemented for anisotropic elasticity. The second one is by Steve Fitzgerald at Culham Centre for Fusion Energy +44 (0)1235 466350. Turn on _ANISO_Steve to run his method.


A test case isalso included in the aniso module, it is called aniso_fr_src.*. Forces from the Matlab code are also added for comparison. The file is compare_aniso_matlab.txt.

Also several tests files are in the aniso module. They verify several components of the implementation.

Retrieved from "http://micro.stanford.edu/mediawiki/index.php?title=ParaDiS_Aniso_Code_Manuals&oldid=4368"