Glassvision/JuliaFEM.jl
1
0
Fork 0
mirror of https://github.com/JuliaFEM/JuliaFEM.jl.git synced 2026-08-06 04:21:33 +00:00
Code Issues Packages Projects Releases Wiki Activity

perseesta 2

Browse Source
This commit is contained in:
ovainola
2015-08-01 15:50:52 +03:00
parent 83552c4e2d
commit b3c4c1750f
4 changed files with 104 additions and 84 deletions
Download Patch File Download Diff File Expand all files Collapse all files
docs/api/JuliaFEM.abaqus_reader.rst
+3 -4
View File
@@ -4,9 +4,8 @@ JuliaFEM.abaqus_reader
Internal
--------
.. function:: add_handler(section, function_name)
.. jl:function:: add_handler(section, function_name)
source:
-------
[JuliaFEM/src/abaqus_reader.jl:18](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/abaqus_reader.jl#L18)
Register new handler for parser
docs/api/JuliaFEM.elasticity_solver.rst
+94 -66
View File
@@ -4,90 +4,118 @@ JuliaFEM.elasticity_solver
Internal
--------
.. function:: assemble!(fe, eldofs_, I, V)
.. jl:function:: assemble!(fe, eldofs_, I, V)
Assemble global RHS to I,V ready for sparse format
Parameters
----------
fe : local vector
eldofs_ : Array
degrees of freedom
I,V : Arrays for sparse matrix
Notes
-----
eldofs can also be node ids for convenience. In that case dimension
is calculated and eldofs are "extended" to problem dimension.
:param fe : local vector
:param eldofs_ : Array
:param I,V : Arrays for sparse matrix
:notes: eldofs can also be node ids for convenience. In that case dimension
is calculated and eldofs are "extended" to problem dimension.
source:
-------
[JuliaFEM/src/elasticity_solver.jl:171](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/elasticity_solver.jl#L171)
.. function:: assemble!(ke, eldofs_, I, J, V)
.. jl:function:: assemble!(ke, eldofs_, I, J, V)
Assemble global stiffness matrix to I,J,V ready for sparse format
Parameters
----------
ke : local matrix
eldofs_ : Array
degrees of freedom
I,J,V : Arrays for sparse matrix
Notes
-----
eldofs can also be node ids for convenience. In that case dimension
is calculated and eldofs are "extended" to problem dimension.
:param ke : local matrix
:param eldofs_ : Array
:param I,J,V : Arrays for sparse matrix
:notes: eldofs can also be node ids for convenience. In that case dimension
is calculated and eldofs are "extended" to problem dimension.
source:
-------
[JuliaFEM/src/elasticity_solver.jl:130](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/elasticity_solver.jl#L130)
.. jl:function:: calc_local_matrices!(X, u, R, Kt, N, dNdchi, lambda_, mu_, ipoints, iweights)
.. function:: calc_local_matrices!(X, u, R, Kt, N, dNdchi, lambda_, mu_, ipoints, iweights)
Calculate local tangent stiffness matrix and residual force vector R = T - F
source:
-------
[JuliaFEM/src/elasticity_solver.jl:68](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/elasticity_solver.jl#L68)
.. function:: eliminate_boundary_conditions(dirichletbc, I, J, V)
.. jl:function:: eliminate_boundary_conditions(dirichletbc, I, J, V)
Eliminate Dirichlet boundary conditions from matrix
Parameters
----------
dirichletbc : array [dim x nnodes]
I, J, V : sparse matrix arrays
Returns
-------
I, J, V : boundary conditions removed
Notes
-----
pros:
- matrix assembly remains positive definite
cons:
- maybe inefficient because of extra sparse matrix operations. (It's hard to remove stuff from sparse matrix.)
- if u != 0 in dirichlet boundary requires extra care
Raises
------
Exception, if displacement boundary conditions given, i.e.
DX=2 for some node, for example.
:param dirichletbc : array [dim x nnodes]
:param I, J, V : sparse matrix arrays
:returns: I, J, V : boundary conditions removed
:notes: pros:
- matrix assembly remains positive definite
cons:
- maybe inefficient because of extra sparse matrix operations. (It's hard to remove stuff from sparse matrix.)
- if u != 0 in dirichlet boundary requires extra care
:raises: Exception, if displacement boundary conditions given, i.e.
DX=2 for some node, for example.
source:
-------
[JuliaFEM/src/elasticity_solver.jl:218](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/elasticity_solver.jl#L218)
.. function:: eliminate_boundary_conditions(dirichletbc, I, V)
.. jl:function:: eliminate_boundary_conditions(dirichletbc, I, V)
Eliminate Dirichlet boundary conditions from vector
Parameters
----------
dirichletbc : array [dim x nnodes]
I, V : sparse vector arrays
Returns
-------
I, V : boundary conditions removed
Notes
-----
pros:
- matrix assembly remains positive definite
cons:
- maybe inefficient because of extra sparse matrix operations. (It's hard to remove stuff from sparse matrix.)
- if u != 0 in dirichlet boundary requires extra care
Raises
------
Exception, if displacement boundary conditions given, i.e.
DX=2 for some node, for example.
:param dirichletbc : array [dim x nnodes]
:param I, V : sparse vector arrays
:returns: I, V : boundary conditions removed
:notes: pros:
- matrix assembly remains positive definite
cons:
- maybe inefficient because of extra sparse matrix operations. (It's hard to remove stuff from sparse matrix.)
- if u != 0 in dirichlet boundary requires extra care
:raises: Exception, if displacement boundary conditions given, i.e.
DX=2 for some node, for example.
source:
-------
[JuliaFEM/src/elasticity_solver.jl:257](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/elasticity_solver.jl#L257)
.. function:: interpolate{T<:Real}(field::Array{T<:Real, 1}, basis::Function, ip)
.. jl:function:: interpolate{T<:Real}(field::Array{T<:Real, 1}, basis::Function, ip)
Interpolate field variable using basis functions f for point ip.
This function tries to be as general as possible and allows interpolating
lot of different fields.
Parameters
----------
field :: Array{Number, dim}
Field variable
basis :: Function
Basis functions
ip :: Array{Number, 1}
Point to interpolate
:param field : Array{Number, dim}
:param basis : Function
:param ip : Array{Number, 1}
source:
-------
[JuliaFEM/src/elasticity_solver.jl:30](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/elasticity_solver.jl#L30)
.. jl:function:: solve_elasticity_increment!(X, u, du, elmap, nodalloads, dirichletbc, lambda, mu, N, dNdchi, ipoints, iweights)
.. function:: solve_elasticity_increment!(X, u, du, elmap, nodalloads, dirichletbc, lambda, mu, N, dNdchi, ipoints, iweights)
source:
-------
[JuliaFEM/src/elasticity_solver.jl:278](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/elasticity_solver.jl#L278)
Solve one increment of elasticity problem
docs/api/JuliaFEM.interfaces.rst
+6 -6
View File
@@ -4,13 +4,13 @@ JuliaFEM.interfaces
Exported
--------
.. function:: solve_elasticity_interface!()
.. jl:function:: solve_elasticity_interface!()
This is generic interface that reads data from data model, solves elasticity
problem and updates model.
:param model : to be defined
source:
-------
[JuliaFEM/src/interfaces.jl:22](https://github.com/JuliaFEM/JuliaFEM.jl/tree/08077a783a6459336cf07b02e323df0e9977c74f/src/interfaces.jl#L22)
Parameters
----------
model : to be defined
docs/build_api.jl
+1 -8
View File
@@ -41,7 +41,6 @@ const modules = append_list
# run(`cp $main_folder/$file_ $this_folder`)
cd(dirname(@__FILE__)) do
<<<<<<< HEAD
# Run the doctests *before* we start to generate *any* documentation.
# for m in modules
# failures = failed(doctest(m))
@@ -51,15 +50,12 @@ cd(dirname(@__FILE__)) do
# exit(1) # Bail when doctests fail.
# end
# end
=======
>>>>>>> origin/master
# Generate and save the contents of docstrings as markdown files.
index = Index()
for mod in modules
Lexicon.update!(index, save(joinpath(api_directory, "$(mod).rst"), mod))
end
<<<<<<< HEAD
#save(joinpath(api_directory, "index.rst"))
save(joinpath(api_directory, "index.rst"), index)
# Add a reminder not to edit the generated files.
# open(joinpath(api_directory, "README.md"), "w") do f
@@ -68,10 +64,7 @@ cd(dirname(@__FILE__)) do
# all changes to the originating docstrings/files rather than these ones.
# """)
# end
=======
# save(joinpath(api_directory, "index.rst"), index; md_subheader = :category)
>>>>>>> origin/master
# info("Adding all documentation changes in $(api_directory) to this commit.")
# success(`git add $(api_directory)`) || exit(1)