Commit Graph

296 Commits

Author SHA1 Message Date
Jukka Aho e13d6482ce Use PkgTestSuite for CI (#131)
This standardizes the CI process between JuliaFEM packages
2017-07-21 16:28:03 +03:00
Jukka Aho a666bb4bd8 Use package AbaqusReader.jl (#127)
Source code related to read and parse ABAQUS .inp files is now living in
it's own repository `AbaqusReader.jl` and in this commit we cleanup the
same files from this repository.

- add AbaqusReader to .travis.yml because it's not registered package yet
- initialize Mesh from AbaqusReader.jl dict
- remove ABAQUS tests and files moved to AbaqusReader.jl
- remove references to old module Abaqus
- move ABAQUS code to preprocess.jl (what is left)
- close issue #122
- close issue #55
2017-07-21 00:40:52 +03:00
Jukka Aho 9ad66be08b Make code 0.6 compatible (#128)
* running v0.6 conversion code proposed by @ovainola in #108.
* change travis so that build is done using 0.6
* documentation is build from 0.6
* fix most of deprecation warnings
* fix test to pass 0.6
2017-07-20 20:46:57 +03:00
Jukka Aho 49b9642e6a Measure testing time for each test file separately (#126)
Uses TimerOutputs.
2017-07-20 17:11:37 +03:00
Jukka Aho 5e3a0c644e add CheckHeader and CheckTabs (#125)
* add `CheckHeader.jl` and `CheckTabs.jl` to `travis.yml`
* fix tab -> 4 spaces
2017-07-20 17:10:15 +03:00
Jukka Aho 436bf119b2 Set up documentation + lint (#121)
- remove some automatically generated stuff not should even be in
repository
- set up lint + Documents.jl in same way it is defined in freshly started projects
- add lint + doctest to after_success so that build pass, these needs to be fixed later
- build is failing on nightly (0.7) but it's ok for release (0.5.2)
- Documenter.jl supports doctests, so this closes least #23
- build system is now on Travis-CI completely, so this closes also #68
2017-07-19 07:34:38 +03:00
Jukka Aho c0bf19687a Reading Code Aster med mesh properly
- MED file actually supports multiple element sets or node sets for
  single element/node. This is now implemented.
- Minor cleanup of code, etc..
- Fixes issue #111.
2017-05-30 06:26:12 +03:00
Andreas Dörfler d50f5fadf0 Add Pyr5 elements, including med support and tests 2017-05-28 00:45:19 +02:00
Jukka Aho 5327459ba4 Use TimerOUtputs.jl to measure time and memory performance
- Added TimerOutputs.jl to REQUIRE file
- JuliaFEM has now TimerOutput handle `to` and function `print_statistics()` to get summary of time and memory usage
- TimerOutputs is used in LinearSolver at the moment and tested in file `test_elasticity_2d_linear_with_surface_load.jl`
2017-04-11 10:39:18 +03:00
Jukka Aho b3f0531746 Add postprocessing features (#100)
* refactored code for solvers.

* Added elementary tests for least-squares fitting of strain and stress fields

* A more realistic postprocess + Xdmf writing test

* removed debug keyword argument from test

* Rewrite update_xdmf!

New function to update Xdmf file no longer takes Solver object but
xdmf, problem, time and fields to write, for example

julia> update_xdmf!(xdmf, problem, 0.0, ["displacement", "temperature"])

All problems are written separately and put together into one
SpatialCollection, allowing to have more structured Xdmf and making
it easier to write complicated field configurations. Support for Xdmf
API 3.0 added.

* Support for Tensor6 field writing

* moved update_xdmf! to io.jl

* Removed some empty files

* Not use old Postprocessor, obsolete code.

* Not use old XDMF (obsolete code). Fixed test.

* removed some postprocessing to pass test, maybe we should drop abaqus.jl from code as obsolete

* add function get_temporal_collection back, it's used by update_xdmf of modal solver

* postprocess of boundary problems also

* added test for contact pressure. dl+quad test output was written in wrong file, fixed.

* postprocess for contact pressure

* contact pressure postprocess

* with boundary problems always store also the primary unknown field

* Change "reaction force" -> "lambda"

* testing postprocess of reaction force also

* sign convention
2017-03-21 08:36:18 +02:00
Jukka Aho a0d18568ce Contact algorithms testing & develpoment (#95)
* contact 3d patch test, standard lagrange, small sliding, linear tet4 elements

* tet4 dual basis contact patch test pass

* contact 3d patch test, standard lagrange, small sliding, linear tet4 elements

* tet4 dual basis contact patch test pass

* Patch test for linear elements standard lagrange / dual lagrange pass now

* Patch test for quadratic contact surfaces for standard + dual basis pass

* refactoring

* renamed files

* Improvements to preprocess scripts

* convert several elements to node sets in one command

* possibility to find particular node from mesh filtered by node set

* 2d small sliding contact patch test, linear elements

* Added backward compatibility

* 2d contact algorithms pass patch tests

* test data for 2d contacts

* no common models in different tests. testing generalized alpha stabilization

* Preprocess tests

* moved tests from test_preprocess_aster_reader.jl to test_preprocess.jl

* generalized-alpha time integration, alpha=0.0 by default

* Improvements to logging

* JuliaFEM.jl: can set environment variable to one of logging levels: OFF, CRITICAL, ERROR, WARNING, INFO, DEBUG

* problems_contact_2d_autodiff.jl: do not loop over nodes if logging level != DEBUG
2017-03-02 08:43:58 +02:00
Jukka Aho 93b520fa2c Added backward compatibility 2017-02-27 08:27:08 +02:00
Jukka Aho 3fafcf243e Improvements to preprocess scripts
* convert several elements to node sets in one command

* possibility to find particular node from mesh filtered by node set
2017-02-27 07:34:04 +02:00
Jukka Aho f275ce3767 Bug/mortar discretization (#88)
* new mortar segmentation tests which are failing

* test_problems_mortar_3d.jl: first test (Tet4) pass

* solvers.jl: diagonal of A is now properly filled, if that option is used. Another option is to remove zero rows from matrix system, which is on by default

* problems_mortar.jl: added new function diagnose_interface to calculate quantities from interface hopefully revealing bugs in calculation

* problems_mortar_3d.jl: added docstring for check_orientation! and removed flooding debug messages not helping to debug anything

* solvers.jl: Another way to solve Ax = b

* Refactored code to make implementation of Tri6 assemble! easier

* Patch test with linear Tet4 elements and quadratic Tet10 elements pass

When using quadratic elements, in polygon clipping algorithm element is divided to linear sub-elements as proposed in [Puso2008]. Interpolation of Lagrange multiplier space is done using quadratic shape functions.

References
----------

[Puso2008] Puso, Michael A., T. A. Laursen, and Jerome Solberg. "A segment-to-segment mortar contact method for quadratic elements and large deformations." Computer Methods in Applied Mechanics and Engineering 197.6 (2008): 555-566.

* increased coverage by adding diagnose_interface

* test using dual basis, failing for unknown reason

* Fixed dual basis construction for Mortar/Tet4

The coefficient matrix Ae for one particular slave element e is the result performing numerical integration on *all* integration cells associated with this element [Popp2013]. Ae cannot be calculated "cell-wise" like it was done before. Now patch test will pass also using `interface.properties.dual_basis = true` option. Partially integrated slave elements are supported as well.

References
----------

[Popp2013] Popp, Alexander, et al. "Improved robustness and consistency of 3D contact algorithms based on a dual mortar approach." Computer Methods in Applied Mechanics and Engineering 264 (2013): 67-80.

* Minor modifications to preprocess.jl

- removed two functions which are unimplemented (but maybe planned in future)
- added function create_node_set_from_element_set!, which can be used, like name suggests, to create a node set from nodes belonging to some set of elements.

* solvers.jl: now prints a list of overconstrained nodes which can be easily copy-pasted to problem.assembly.removed_dofs list to solver overconstrained situation manually

* Increase code coverage

Added a new test which tests dual basis 3d mortar + adjust option when using Tet4 in elasticity problem.

* Tet10 + Dual basis still failing, others are working

* mortar 3d low level tests

* linear surface element projection tests pass

* Introduced basis transform constant alpha

Tet10 + dual basis patch test still failing, but single element low level routine tests gives expected results with alpha=0.2

* added new integration rule FPG12 for triangular elements

* added drop_tolerance option to remove very small values from constraint matrices

* Introduced a basis transform matrix T

Constructing bi-orthogonal basis for quadratic surfaces is ill-conditioned. By doing a basis transform N' = N*T for slave side displacement vector it's possible to construct a bi-orthogonal basis in a same way than with linear elements. Setting alpha=0.2 ensures that quadratic basis functions are strictly positive in practical cases.

* fix 3d clipping test routine, accepts only 3d vertices

* dropped number of integration poitns from 12 to 7 in quadratic mortar surfaces intrestingly gives more accurate results, maybe something numerical error in FPG12 integration rule..?

* added two displacement patch tests + output writing for all cases

* %s/Int64/Int/g

* Changed test data location

* Fine tuning of logging levels
2017-02-25 18:40:14 +02:00
Jukka Aho 0ca7efd631 Testing/code coverage (#83)
Change the code coverage to green. 

* removed duplicate code

* Removed unused code

* removed unmaintained code

* DCTI + DVTI refactored

* discrete fields refactored and tested

* fields are now tested quite well.

* Removed obsolete code not used anywhere

* Element descriptions to common dictionary

* size in global const dictionary also

* Added coverage to sparse tools and removed couple unused functions

* get nonzero rows from SparseMatrixCSC

* bugfix: extending element basis now working and tested

* Removed two unused functions from elements.jl

* removed useless function

* Useless conversion

* remove elasticity assembly using ForwardDiff because it's not used anywhere'

* Added basic testing for NURBS. Fixed bug in NSolid interpolation.

* removed unused functions

* Removed some debug stuff

* renamed file

* removed field assembly posthook, i think not good idea at all

* test for nnz(K) == 0 and automatic determination of dofs

* Testing that solver is throwing error if having problems with boundary assembly

* Removed some unused options. Refactoring.

* Moved solver non-related code to elements.jl

* Removed custom exception (no need)

* unneeded postprocess code

* More tests for NURBS elements.

* Removed unfinished .mail parser

* proper use of Logging package

* also read results

* renamed test file

* create_surface_elements accepts surface name in String now

* bugfix: remove zero rows from constraint matrix after manually removing dofs from some boundary assemblies.

* New test, displacement 3d patch test

* skip displacement field in surface element splitting if not defined

* test element splitting and linear surface elements, fails.

* Bugfix: Xdmf, not XDMF

* removed nonworking tests, requires bugfix

* abaqus_read_results is not working -> bug
2017-01-30 12:28:33 +02:00
Jukka Aho 40ef9388c3 Fix FactCheck dependencies (issue #80)
- Remove some obsolete tests related to Field.
- Fixed a test checking is header information found from source files
  (was using FactCheck and commented out).
- Add header information for all files (was missing from
  materials_plasticity.jl)
2017-01-27 10:15:52 +02:00
Jukka Aho 500fbbcf9f Fixed multiple dispatch issues on io.jl
Related to issue #74, io.jl is no more modifying functionality of
LightXML
2017-01-21 12:32:53 +02:00
Jukka Aho b2be50757b Inceased coverage (#77)
- added Xdmf() to solver for two tests to test also creating Xdmf after solution.
2017-01-21 09:57:52 +02:00
Jukka Aho 28c31b4c5a fixed test 2017-01-09 13:21:07 +02:00
Jukka Aho 3eba94a52d removed obsolete test files 2017-01-09 11:26:26 +02:00
Jukka Aho cf762c4f25 fixed test 2017-01-09 09:32:15 +02:00
Jukka Aho 0913b23d01 possibility to manually remove dofs from problems by adding dof number to problem.assembly.removed_dofs vector 2017-01-09 08:56:31 +02:00
Jukka Aho d0023c5dd4 fixed test 2017-01-06 07:24:58 +02:00
Jukka Aho 3f4bcfbf7e fixed test 2017-01-06 07:22:09 +02:00
Jukka Aho ea94f3a942 io tests 2017-01-05 11:29:29 +02:00
Jukka Aho 1630dd1ec0 dict fields 2017-01-05 11:03:19 +02:00
Jukka Aho cde0f15dcd no id in constructor 2017-01-05 09:32:43 +02:00
Jukka Aho 3ba1046e0b fixed element initialization 2017-01-05 09:30:41 +02:00
Jukka Aho 1d4a150022 fixed 2d finite sliding contact for julia 0.5 2017-01-05 09:09:11 +02:00
Jukka Aho 5601009bd7 SparseVectorCOO 2017-01-05 09:08:32 +02:00
Jukka Aho 78da7a9ecb connectivity information needs now to be given every element 2017-01-05 07:22:15 +02:00
Jukka Aho 7d46d33b6a removed obsolete code 2017-01-05 07:20:17 +02:00
Jukka Aho 1c8bc9bc64 removed null space related code as obsolete 2017-01-05 07:14:54 +02:00
Jukka Aho 89645be3bd added license string. this is not a test (tests missing) 2017-01-05 05:32:13 +02:00
Frondelius 46ccba4cf1 Improving test coverage: test_lagrange.jl fixing copy-paste 2017-01-03 22:38:32 +02:00
Frondelius 36a7e17bf2 Improving test coverage: test_lagrange.jl 2017-01-03 22:11:27 +02:00
Frondelius c4d5185b27 Improving test coverage: test_preprocess.jl 2017-01-03 18:30:17 +02:00
Frondelius 0fefd2791c Improving test coverage: test_preprocess_abaqus_reader.jl 2017-01-03 17:11:10 +02:00
Jukka Aho 2914d39330 added test for zero eigenvalue model (no dirichlet boundaries defined) 2016-12-12 07:40:18 +02:00
Jukka Aho f383bd7140 test stress tensor at integration points compared to Ca solution 2016-12-02 10:09:27 +02:00
Jukka Aho 319a0174c6 test stress extrapolation to nodes 2016-12-01 11:14:32 +02:00
Jukka Aho 05bd562cb4 first attempts to verify stress calculation to nodes 2016-11-30 16:19:40 +02:00
Jukka Aho 4c4130e2f9 fixed whitespace bug in reading *NODE block 2016-11-28 10:09:48 +02:00
Tero Frondelius 53a7d665c5 try catch to ignore *ELSET creation from existing sets 2016-11-17 21:35:40 +02:00
Tero Frondelius b7f3900c03 missing testdata 2016-11-17 20:40:02 +02:00
Tero Frondelius 88a72f7bde Fixed abaqus_read_mesh() 1.0E+00 problem 2016-11-17 19:43:20 +02:00
Tero Frondelius 2ef93e54a6 fixing the tests after last update 2016-11-17 18:54:08 +02:00
Jukka Aho e64ed7f6a3 v0.5 fixes 2016-11-15 22:01:00 +02:00
Jukka Aho 95200a81e1 fixed deprecation warnings for 0.5 2016-11-13 13:24:08 +02:00
Olli Väinölä c3bc708f01 Ideal plasticity converged, both 2D and 3D 2016-10-28 15:46:04 +03:00