ParaDiS Aniso Code Manuals: Difference between revisions

From Micro and Nano Mechanics Group
Jump to navigation Jump to search
 
(4 intermediate revisions by the same user not shown)
Line 1: Line 1:
Anisotropic Elasticity
=Aniso code in ParaDiS=
==Sylvie Aubry Dec 11 2009==
Sylvie Aubry Jan 28 2010.


The aniso module is similar to DDLab/aniso. For more information refer
===How to use the aniso code===
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
Before starting using the aniso code, the user decides if he/she wants to use the tables or not and if he/she wants to use the Voigt or the Reuss averages for mu and nu.
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:
To use the tables, turn on DEFS += D_TABLES in the makefile.
To create the table files, set QSBMETHOD = 1 in tabulate.c.
To read the table, make sure they are located in the inputs directory.


D_ANISOTROPIC : turns on or off anisotropic elasticity.
To set the mu and nu variables, set #define VOIGT and #define REUSS variables in ANI_Util.c


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:


===Flow chart of the aniso code===
'''ANI_Main'''


TableQ character Define the name of the file to read in tables.
defines the directory where files Qtable.dat, Stable.dat and Rtable.dat will be created or read.
Defaults is inputs/Qtable.dat
define the aniso0 structure
calls ANI_Init in ANI_Utils.


TableS character Define the name of the file to read in tables.
'''ANI_Init'''
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
initialize the aniso0 structure.
second one is by Steve Fitzgerald at Culham Centre for Fusion Energy
. define C, eps
+44 (0)1235 466350. Turn on _ANISO_Steve to run his method.


. define mu and nu. User can choose between the Voigt or the Reuss averages.


A test case isalso included in the aniso module, it is called
. calls the routine that defines Q, S, B. Uses file QSBMatrices.c
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
. user can choose to read in tables (QSBMETHOD = 1) in tabulate.c or not.
components of the implementation.

'''Param.c and Param.h'''

defines variables specific to aniso code: C12, C44, Cpr.
Also defines the tables files.

NodeForce.c

When the ANISOTROPIC flag is on, use anisotropic SelfForce.
The SeflForce routine calls QSBMatrices.c routines to get the QSB matrices.


'''LocalSegForces.c'''
When the ANISOTROPIC flag is on, use anisotropic SegSegForces.
This routine calls ComputeANIForces located in Stress.c.

'''Stress.c'''

ComputeANIForces.c calls the ANI_WillisSteeds function also located in Stress.c.
ANI_WillisSteeds decide whether to use the regular definition or the collinear one.

That's it!

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"