POLYHEDRON
所属分类:matlab编程
开发工具:matlab
文件大小:76KB
下载次数:10
上传日期:2014-04-22 09:57:45
上 传 者:
伍海军
说明:
采用线积分均匀任意形状的多面体源的全部重力张量的分析计算
(Analytical computation of the full gravity tensor of a homogeneous arbitrarily shaped polyhedral source using line integrals)
文件列表:
POLYHEDRON (0, 2012-03-24)
POLYHEDRON\xyzposnew (136, 2010-06-11)
POLYHEDRON\polaut2_out (13467, 2010-06-11)
POLYHEDRON\wuerfel.f (13477, 2012-03-24)
POLYHEDRON\additionaltables.pdf (67756, 2011-06-21)
POLYHEDRON\topoaut (66, 2008-07-09)
POLYHEDRON\wuerfel_dat (280, 2010-06-11)
POLYHEDRON\dataut (16, 1998-12-10)
POLYHEDRON\wuerfel2_out (18586, 2010-06-11)
POLYHEDRON\polyhedron.f (26189, 2012-03-24)
POLYHEDRON\ptw_dat (250, 1999-07-31)
The appended electronic supplement contains following 9 files:
1. polyhedron.f
This file uses following input files
2. xyzposnew
3. topoaut
4. dataut
and writes to the output file
5. polaut2_out
Furthermore appended is the code for the computation of the
gravity tensor of a right rectangular prism according
to Mader (1951). The code is given by
6. wuerfel.f
which uses the input files
7. wuerfel_dat
8. ptw_dat
and writes to the output file
9. wuerfel2_out
Note that the appended-type files 'polaut2_out' and 'wuerfel2_out'
contain already 28 discrete cases as mentioned in the manuscript,
where for each case file 'xyzposnew' expressing the exact geometric
setting has also been added manually.
For running or testing code wuerfel.f for any rectangular prismatic
source you may need to remove the comments from lines 341, 391, 425
and 460 which refer to MSIMSL and read
USE MSIMSL
Module MSIMSL is specific to some Fortran compilers, e.g. Visual Compaq,
and is necessary for computing the hyperbolic arcsine function. Most
of the widely used compilers (gfortran, gfortran44 and ifort) already
include all of the mathematical functions used in wuerfel.f and for these
compilers the USE MSIMSL statement is redundant.
For running or testing code polyhedron.f for any polyhedral
source following steps are to be taken
A. Insertion or 'save as' in a file named 'xyzposnew' of all vertex
coordinates with respect to the computation point as X, Y, Z (the
origin is located at the computation point).
B. Definition/Computation of the topology matrix 'topoaut' of this
source where vertex numbers should comply with the order of the
vertices in file 'xyzposnew'.
C. Definition of the varying number of segments building each of the
polyhedral faces in file 'dataut'. File 'dataut' can be built simply
by counting the vertices of each line (record) of file 'topoaut'.
D. Manual editing of parameters 'nop' (number of planes) and 'nov'
(number of vertices) in Line 127 of polyhedron.f with the actual
total number of faces and vertices that build the source. Line 127
refers currently to the source body of Figure 2 of the paper and reads
parameter(nop=7,nov=8,noed=100)
Note that parameter 'noed'
has been set to a conveniently large number
that expresses the varying number of vertices building each face. If
any of the polyhedral faces consists of more than 100 vertices (or, equally,
segments) this parameter has to be changed accordingly as well.
E. Compile and run 'polyhedron.f'. The results will be appended to
the end of file 'polaut2_out'
If you do not change anything to the present input files and simply
compile and run 'polyhedron.f' and 'wuerfel.f' you should get an exact
copy of the last 10 lines of files 'polaut2_out' and 'wuerfel2_out'
appended at the end of these files.
Dimitris Tsoulis, August 21, 2011
近期下载者:
相关文件:
收藏者: