PARADISCYL:How-To-Get-Nodal-Force - Revision history

Cweinber at 17:04, 27 October 2008

2008-10-27T17:04:13Z

← Older revision		Revision as of 17:04, 27 October 2008
Line 4:		Line 4:
	<DIV>		<DIV>

	This document describes the extension of the ParaDiS program that implements the cylindrical free surface boundary condition. This extension is based on Weinberger and Cai, J. Mech. Phys. Solids, 55, 2027 (2007)([http://micro.stanford.edu/~caiwei/papers/Weinberger07jmps-image.pdf PDF]). The code is written by Chris Weinberger. It is ~~owned by Chris~~ and is not included in the standard distribution of ParaDiS. Keonwook Kang wrote the first draft of this document, which was later revised by Chris Weinberger and Wei Cai. Questions on this document should be addressed to [[Chris Weinberger]]. Please read the manual of ParaDiS before reading this document.		This document describes the extension of the ParaDiS program that implements the cylindrical free surface boundary condition. This extension is based on Weinberger and Cai, J. Mech. Phys. Solids, 55, 2027 (2007)([http://micro.stanford.edu/~caiwei/papers/Weinberger07jmps-image.pdf PDF]). The code is written by Chris Weinberger and Wei Cai and is not included in the standard distribution of ParaDiS. Keonwook Kang wrote the first draft of this document, which was later revised by Chris Weinberger and Wei Cai. Questions on this document should be addressed to [[Chris Weinberger]]. Please read the manual of ParaDiS before reading this document.

	== Motivation ==		== Motivation ==

Caiwei at 08:51, 25 October 2008

2008-10-25T08:51:29Z

← Older revision		Revision as of 08:51, 25 October 2008
Line 4:		Line 4:
	<DIV>		<DIV>

	This document describes the extension of the ParaDiS program that implements the cylindrical free surface boundary condition. This extension is based on Weinberger and Cai, J. Mech. Phys. Solids, 55, 2027 (2007)([http://micro.stanford.edu/~caiwei/papers/Weinberger07jmps-image.pdf PDF]). The code is written by Chris Weinberger. It is owned by Chris and is not included in the standard distribution of ParaDiS. Keonwook Kang wrote the first draft of this document, which was later revised by Chris Weinberger and Wei Cai. Please read the manual of ParaDiS before reading this document.		This document describes the extension of the ParaDiS program that implements the cylindrical free surface boundary condition. This extension is based on Weinberger and Cai, J. Mech. Phys. Solids, 55, 2027 (2007)([http://micro.stanford.edu/~caiwei/papers/Weinberger07jmps-image.pdf PDF]). The code is written by Chris Weinberger. It is owned by Chris and is not included in the standard distribution of ParaDiS. Keonwook Kang wrote the first draft of this document, which was later revised by Chris Weinberger and Wei Cai. Questions on this document should be addressed to [[Chris Weinberger]]. Please read the manual of ParaDiS before reading this document.

	== Motivation ==		== Motivation ==

Caiwei at 08:23, 25 October 2008

2008-10-25T08:23:36Z

← Older revision		Revision as of 08:23, 25 October 2008
Line 3:		Line 3:
	How to get nodal force</H1>		How to get nodal force</H1>
	<DIV>		<DIV>

			This document describes the extension of the ParaDiS program that implements the cylindrical free surface boundary condition. This extension is based on Weinberger and Cai, J. Mech. Phys. Solids, 55, 2027 (2007)([http://micro.stanford.edu/~caiwei/papers/Weinberger07jmps-image.pdf PDF]). The code is written by Chris Weinberger. It is owned by Chris and is not included in the standard distribution of ParaDiS. Keonwook Kang wrote the first draft of this document, which was later revised by Chris Weinberger and Wei Cai. Please read the manual of ParaDiS before reading this document.

			== Motivation ==

	In ParaDiS, a dislocation line is represented by piecewise straight segments connecting a set of discrete nodal points. For diagnostic purposes, we may want to find out the driving force on every node. Here we explain how to extract this information from the ParaDiS Cylinder code.		In ParaDiS, a dislocation line is represented by piecewise straight segments connecting a set of discrete nodal points. For diagnostic purposes, we may want to find out the driving force on every node. Here we explain how to extract this information from the ParaDiS Cylinder code.

Caiwei: /* Computing nodal force */

2008-10-25T08:14:05Z

Computing nodal force

← Older revision		Revision as of 08:14, 25 October 2008
Line 12:		Line 12:
	Then compile ParaDiS cylinder codes again.		Then compile ParaDiS cylinder codes again.

	Next, run the executable with the input file, such as that in [[PARADISCYL:Test-Run]],		Next, run the executable with the input file, such as that in [[PARADISCYL:Test-Run \|M02 Test Run]],

	$ bin/paradiscyl runs/concentric_loop_test.ctrl.		$ bin/paradiscyl runs/concentric_loop_test.ctrl.

Caiwei: /* Computing nodal force */

2008-10-25T08:13:21Z

Computing nodal force

← Older revision		Revision as of 08:13, 25 October 2008
Line 12:		Line 12:
	Then compile ParaDiS cylinder codes again.		Then compile ParaDiS cylinder codes again.

	Next, run the executable with the input file, such as that in [PARADISCYL:Test-Run]],		Next, run the executable with the input file, such as that in [[PARADISCYL:Test-Run]],

	$ bin/paradiscyl runs/concentric_loop_test.ctrl.		$ bin/paradiscyl runs/concentric_loop_test.ctrl.

Caiwei: /* Computing nodal force */

2008-10-25T08:12:50Z

Computing nodal force

← Older revision		Revision as of 08:12, 25 October 2008
Line 8:		Line 8:
	== Computing nodal force ==		== Computing nodal force ==

	Modify '''makefile''' in ParaDiS/cylinder directory ~~such~~ that ~~'''DEFS'''~~ ~~macro~~ ~~include~~ ~~definition~~ ~~of '''WRITENODEFORCE''' or~~		Modify '''makefile''' in ParaDiS/cylinder directory so that the following line is active.
	DEFS += -D_WRITENODEFORCE		DEFS += -D_WRITENODEFORCE
	Then compile ParaDiS cylinder codes again.		Then compile ParaDiS cylinder codes again.

			Next, run the executable with the input file, such as that in [PARADISCYL:Test-Run]],
	Let's say we want to know nodal forces of a dislocation loop inside a cylinder used in test simulation run.
	Run the control file, '''runs/concentric_loop_test.ctrl'''.

	$ bin/paradiscyl runs/concentric_loop_test.ctrl.		$ bin/paradiscyl runs/concentric_loop_test.ctrl.

	Note that, even if the total simulation step is set as 100 in the control file, the ParaDiS simulation does one time force evaluation and will quit ~~after~~ ~~then~~.		Note that, even if the total simulation step is set as 100 in the control file, the ParaDiS simulation does one time force evaluation and will quit immediately afterwards.

	~~You will see that the~~ force data is dumped out as '''force.out''' in '''outputs/concentric_loop_test/''' directory.		The force data is dumped out as '''force.out''' in the '''outputs/concentric_loop_test/''' directory. Each line of this file has 6 values for nodal point coordinate (''x'', ''y'', ''z'') and nodal force components (''f''<sub>x</sub>, ''f''<sub>y</sub>, ''f''<sub>z</sub>), respectively.
	Each line of the force data file has 6 numbers of nodal point coordinate (''x'', ''y'', ''z'') and nodal force components (''f''<sub>x</sub>, ''f''<sub>y</sub>, ''f''<sub>z</sub>).

			Because the default stress unit is Pa and length unit is '''burgMag''', the unit of nodal forces is Pa.burgMag^2. Of course, this is the case for cylinder radius R = 1. If we want to interprete the results for the case of R not equal to one, then the force will have to be scaled according to [[PARADISCYL:Scale-Rule \|M03 How units are scaled]].
	('''What is the unit of force???''')
	----		----

Caiwei at 08:04, 25 October 2008

2008-10-25T08:04:55Z

← Older revision		Revision as of 08:04, 25 October 2008
Line 3:		Line 3:
	How to get nodal force</H1>		How to get nodal force</H1>
	<DIV>		<DIV>

	== Overview ==

	In ParaDiS, a dislocation line is represented by piecewise straight segments connecting a set of discrete nodal points. For diagnostic purposes, we may want to find out the driving force on every node. Here we explain how to extract this information from the ParaDiS Cylinder code.		In ParaDiS, a dislocation line is represented by piecewise straight segments connecting a set of discrete nodal points. For diagnostic purposes, we may want to find out the driving force on every node. Here we explain how to extract this information from the ParaDiS Cylinder code.

	----

	== Computing nodal force ==		== Computing nodal force ==

Caiwei at 08:04, 25 October 2008

2008-10-25T08:04:25Z

← Older revision		Revision as of 08:04, 25 October 2008
Line 1:		Line 1:
	<H1 ALIGN="CENTER"><FONT SIZE="-1">Manual 04 for ParaDiS Cylinder Codes</FONT>		<H1 ALIGN="CENTER"><FONT SIZE="-1">Manual 04 for ParaDiS Cylinder Codes</FONT>
	<BR>		<br><br>
	How to get nodal force</H1>		How to get nodal force</H1>
	<DIV>		<DIV>

	<P ALIGN="CENTER"><STRONG>Keonwook Kang and Wei Cai</STRONG></P>
	</DIV>
	<P ALIGN="CENTER">Original date : [[ Oct 2 ]], [[2008]]</P>
	<P ALIGN="CENTER">Latest update on [[ Oct 2 ]], [[2008]]</P>
	<P>
	<BR><HR>
	<!--Table of Child-Links-->

	== Overview ==		== Overview ==

	In ParaDiS, a dislocation is ~~idealized~~ as a ~~piecewise~~ ~~line object:~~ discrete nodal points ~~and~~ ~~lines~~ ~~connecting~~ ~~them. Sometimes~~, we ~~need~~ to ~~know~~ ~~how~~ ~~much~~ force ~~is acting~~ on ~~each~~ node ~~under given dislocation geometry~~. In ~~this~~ ~~manual,~~ ~~it will be explained~~ how we ~~export~~ ~~nodal~~ ~~force data~~ from ParaDiS ~~simulation~~.		In ParaDiS, a dislocation line is represented by piecewise straight segments connecting a set of discrete nodal points. For diagnostic purposes, we may want to find out the driving force on every node. Here we explain how to extract this information from the ParaDiS Cylinder code.

	----		----

	== ~~Procedure~~ ==		== Computing nodal force ==

	Modify '''makefile''' in ParaDiS/cylinder directory such that '''DEFS''' macro include definition of '''WRITENODEFORCE''' or		Modify '''makefile''' in ParaDiS/cylinder directory such that '''DEFS''' macro include definition of '''WRITENODEFORCE''' or

Caiwei: How-To-Get-Nodal-Force moved to PARADISCYL:How-To-Get-Nodal-Force

2008-10-25T07:09:38Z

How-To-Get-Nodal-Force moved to PARADISCYL:How-To-Get-Nodal-Force

← Older revision	Revision as of 07:09, 25 October 2008
(No difference)

Kwkang at 22:01, 23 October 2008

2008-10-23T22:01:53Z

← Older revision		Revision as of 22:01, 23 October 2008
Line 1:		Line 1:
	<H1 ALIGN="CENTER"><FONT SIZE="-1">Manual 03 for ParaDiS Cylinder Codes</FONT>		<H1 ALIGN="CENTER"><FONT SIZE="-1">Manual 04 for ParaDiS Cylinder Codes</FONT>
	<BR>		<BR>
	How to get nodal force</H1>		How to get nodal force</H1>