contact constraint example

About the development of the FEM module/workbench.

Moderator: bernd

User avatar
bernd
Posts: 8485
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

contact constraint example

Postby bernd » Fri Jan 27, 2017 12:30 pm

I downloaded the contact constraint example from here: viewtopic.php?t=16749#p132166 applied a new material Calculix-Steel and a new Solver calculix. But ccx breaks with some error. May be I use the wrong iteration parameter ... The FreeCAD file is attached.

Code: Select all

FemTools.setup_working_dir()  -->  self.working_dir = /tmp/FEMWB/
FemInputWriterCcx --> self.dir_name  -->  /tmp/FEMWB/
FemInputWriterCcx --> self.main_file_name  -->  Fusion_Mesh.inp
FemInputWriterCcx --> self.file_name  -->  /tmp/FEMWB//Fusion_Mesh.inp
('  ReferenceShape : ', 'Face', ', ', 'Fusion', ', ', u'Fusion', ' --> ', 'Face5')
('  ReferenceShape : ', 'Face', ', ', 'Fusion', ', ', u'Fusion', ' --> ', 'Face8')
('  ReferenceShape : ', 'Face', ', ', 'Fusion', ', ', u'Fusion', ' --> ', 'Face6')
('  ReferenceShape : ', 'Face', ', ', 'Fusion', ', ', u'Fusion', ' --> ', 'Face10')
len femnodes_ele_table:2983
('  ReferenceShape : ', 'Face', ', ', 'Fusion', ', ', u'Fusion', ' --> ', 'Face9')
len femnodes_ele_table:2983
len node_set: 577
len bit_pattern_dict:1316
('found Faces: ', 252)
('faces: ', [[1322, 2], [1516, 1], [1701, 1], [1727, 1], [1728, 1], [1729, 1], [1733, 3], [1745, 2], [1746, 3], [1751, 3], [1754, 2], [1756, 1], [1757, 4], [1759, 4], [1761, 1], [1762, 1], [1763, 2], [1765, 1], [1773, 2], [1775, 1], [1778, 2], [1779, 2], [1781, 1], [1782, 4], [1786, 3], [1803, 2], [1804, 2], [1805, 2], [1806, 4], [1807, 2], [1809, 3], [1810, 4], [1812, 4], [1813, 3], [1815, 1], [1817, 2], [1821, 4], [1824, 3], [1825, 3], [1826, 2], [1828, 4], [1830, 2], [1834, 1], [1836, 2], [1839, 1], [1840, 3], [1842, 1], [1846, 2], [1847, 4], [1848, 2], [1852, 4], [1854, 1], [1866, 4], [1869, 4], [1875, 3], [1881, 2], [1882, 1], [1883, 2], [1885, 3], [1886, 3], [1888, 3], [1895, 1], [1896, 2], [1897, 2], [1902, 2], [1905, 1], [1914, 2], [1921, 4], [1933, 2], [1935, 2], [1940, 3], [1941, 2], [1946, 1], [1947, 3], [1953, 1], [1958, 3], [1961, 3], [1966, 2], [1976, 3], [1983, 2], [1989, 2], [1993, 1], [1994, 4], [1996, 4], [1998, 2], [2001, 2], [2008, 1], [2010, 3], [2018, 4], [2024, 1], [2030, 4], [2031, 4], [2042, 3], [2044, 1], [2046, 4], [2047, 2], [2056, 4], [2057, 2], [2071, 1], [2075, 1], [2076, 4], [2079, 3], [2080, 1], [2081, 1], [2082, 2], [2083, 3], [2088, 1], [2089, 2], [2090, 2], [2091, 2], [2093, 2], [2107, 3], [2108, 2], [2116, 4], [2119, 4], [2122, 1], [2125, 3], [2128, 4], [2135, 1], [2136, 3], [2140, 3], [2141, 4], [2146, 1], [2147, 2], [2151, 3], [2152, 4], [2154, 4], [2155, 3], [2156, 3], [2160, 3], [2164, 3], [2166, 4], [2167, 1], [2168, 2], [2169, 3], [2171, 2], [2172, 1], [2174, 4], [2176, 2], [2179, 4], [2183, 1], [2189, 3], [2191, 1], [2192, 4], [2197, 3], [2198, 3], [2200, 1], [2204, 2], [2206, 2], [2211, 3], [2216, 2], [2218, 3], [2220, 2], [2222, 4], [2228, 4], [2229, 1], [2231, 1], [2235, 2], [2238, 4], [2241, 2], [2246, 1], [2247, 2], [2250, 1], [2251, 2], [2259, 2], [2262, 2], [2264, 3], [2266, 4], [2268, 4], [2269, 2], [2272, 4], [2273, 4], [2276, 3], [2278, 4], [2279, 1], [2286, 1], [2290, 3], [2292, 3], [2296, 3], [2307, 3], [2310, 4], [2312, 3], [2314, 1], [2318, 1], [2321, 3], [2323, 2], [2325, 4], [2327, 1], [2336, 4], [2344, 2], [2346, 4], [2347, 3], [2349, 2], [2355, 2], [2358, 3], [2361, 2], [2363, 4], [2369, 3], [2372, 2], [2380, 4], [2390, 1], [2392, 4], [2396, 4], [2397, 4], [2400, 3], [2404, 3], [2406, 2], [2408, 3], [2425, 1], [2430, 2], [2436, 2], [2439, 2], [2442, 3], [2451, 1], [2452, 3], [2462, 1], [2464, 1], [2472, 1], [2477, 4], [2478, 1], [2482, 3], [2488, 3], [2491, 1], [2493, 1], [2496, 2], [2499, 1], [2505, 3], [2512, 2], [2514, 3], [2521, 1], [2526, 4], [2530, 1], [2534, 2], [2536, 3], [2537, 4], [2542, 1], [2546, 2], [2548, 4], [2550, 3], [2552, 4], [2556, 3], [2557, 2], [2559, 2], [2560, 2], [2570, 1], [2580, 2], [2584, 2], [2585, 1], [2586, 1], [2589, 1], [2592, 3], [2596, 2]])
Writing time input file: 0.657094 

Analysis failed with exit code 201
--------start of stderr-------

--------end of stderr---------
--------start of stdout-------

************************************************************

CalculiX Version 2.11, Copyright(C) 1998-2015 Guido Dhondt
CalculiX comes with ABSOLUTELY NO WARRANTY. This is free
software, and you are welcome to redistribute it under
certain conditions, see gpl.htm

************************************************************

You are using an executable made on So 31. Jul 13:26:31 CEST 2016

  The numbers below are estimated upper bounds

  number of:

   nodes:         2983
   elements:         2596
   one-dimensional elements:            0
   two-dimensional elements:            0
   integration points per element:            4
   degrees of freedom per node:            3
   layers per element:            1

   distributed facial loads:          252
   distributed volumetric loads:            0
   concentrated loads:            0
   single point constraints:          312
   multiple point constraints:            1
   terms in all multiple point constraints:            1
   tie constraints:            1
   dependent nodes tied by cyclic constraints:            0
   dependent nodes in pre-tension constraints:            0

   sets:            6
   terms in all sets:         9304

   materials:            2
   constants per material and temperature:            4
   temperature points per material:            1
   plastic data points per material:            0

   orientations:            0
   amplitudes:            2
   data points in all amplitudes:            2
   print requests:            2
   transformations:            0
   property cards:            0

 *INFO: time control parameters set to:
        i0 =            4
        ir =            8
        ip =            9
        ic =          200
        il =           10
        ig =          400
        is =            0
        ia =          200
        ij =            0
        it =            0
        df =   0.25000000000000000     
        dc =   0.50000000000000000     
        db =   0.75000000000000000     
        da =   0.84999999999999998     
        ds =    0.0000000000000000     
        dh =    0.0000000000000000     
        dd =    1.5000000000000000     
        wg =    0.0000000000000000     

 STEP            1

 Static analysis was selected

 Newton-Raphson iterative procedure is active

 Decascading the MPC's

 Determining the structure of the matrix:
 number of equations
 8637
 number of nonzero lower triangular matrix elements
 269628

 increment 1 attempt 1 
 increment size= 1.000000e-01
 sum of previous increments=0.000000e+00
 actual step time=1.000000e-01
 actual total time=1.000000e-01

 iteration 1

 Number of contact spring elements=0

 Determining the structure of the matrix:
 number of equations
 8637
 number of nonzero lower triangular matrix elements
 269628

 Using up to 2 cpu(s) for the stress calculation.

 Using up to 2 cpu(s) for the symmetric stiffness/mass contributions.

 Factoring the system of equations using the symmetric spooles solver
 Using up to 2 cpu(s) for spooles.

 Using up to 2 cpu(s) for the stress calculation.

 average force= 2095.935420
 time avg. forc= 2095.935420
 largest residual force= 0.000000 in node 2899 and dof 2
 largest increment of disp= 5.586546e+00
 largest correction to disp= 5.586546e+00 in node 1738 and dof 2

 no convergence

 iteration 2

 Number of contact spring elements=0

 Determining the structure of the matrix:
 number of equations
 8637
 number of nonzero lower triangular matrix elements
 269628

 Using up to 2 cpu(s) for the stress calculation.

 Using up to 2 cpu(s) for the symmetric stiffness/mass contributions.

 Factoring the system of equations using the symmetric spooles solver
 Using up to 2 cpu(s) for spooles.

 Using up to 2 cpu(s) for the stress calculation.

 average force= 2095.935420
 time avg. forc= 2095.935420
 largest residual force= 0.000000 in node 2903 and dof 2
 largest increment of disp= 5.586546e+00
 largest correction to disp= 3.813563e-12 in node 709 and dof 2

 convergence

 Using up t
--------end of stdout---------
CalculiX failed ccx finished with error 201
OS: Debian GNU/Linux 8.7 (jessie)
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.17.9803 (Git)
Build type: Unknown
Branch: master
Hash: dd63fe92107a87f6003f33518e387743a2090a1c
Python version: 2.7.9
Qt version: 4.8.6
Coin version: 4.0.0a
OCC version: 7.0.0
contact-example-udated.fcstd
(356.43 KiB) Downloaded 62 times
Last edited by bernd on Sun Jan 29, 2017 7:34 pm, edited 1 time in total.
fandaL
Posts: 343
Joined: Thu Jul 24, 2014 8:29 am

Re: contact constraint example

Postby fandaL » Sat Jan 28, 2017 10:35 am

bernd wrote:
contact-example-udated.fcstd
With your file I get results. I didn't change any parameter.

OS: Windows 7
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.17.9822 (Git)
Build type: Release
Branch: master
Hash: 0df4af4ba5e6edce825af6ecf6d024c704292633
Python version: 2.7.8
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 7.0.0
User avatar
bernd
Posts: 8485
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: contact constraint example

Postby bernd » Sat Jan 28, 2017 11:07 am

fandaL wrote:
bernd wrote:contact-example-udated.fcstd
With your file I get results. I didn't change any parameter
I get resuts too, but they are different from the ones included in the file. These are the original results.
fandaL
Posts: 343
Joined: Thu Jul 24, 2014 8:29 am

Re: contact constraint example

Postby fandaL » Sat Jan 28, 2017 11:31 am

Original absolute displacement:
29,7740 mm
After recomputing (in the last result object "CalculiX_static_time_0_0_results009"):
29.7562 mm
bernd wrote:...new material Calculix-Steel...
There is "Steel-Generic" with E = 200 000 MPa, I didn't change it.
User avatar
bernd
Posts: 8485
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: contact constraint example

Postby bernd » Sat Jan 28, 2017 12:17 pm

fandaL wrote:Original absolute displacement:
29,7740 mm
After recomputing (in the last result object "CalculiX_static_time_0_0_results009"):
29.7562 mm
bernd wrote:...new material Calculix-Steel...
There is "Steel-Generic" with E = 200 000 MPa, I didn't change it.
which version of ccx do you use?
User avatar
bernd
Posts: 8485
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: contact constraint example

Postby bernd » Sat Jan 28, 2017 12:27 pm

would you post your *.inp file?
fandaL
Posts: 343
Joined: Thu Jul 24, 2014 8:29 am

Re: contact constraint example

Postby fandaL » Sat Jan 28, 2017 3:08 pm

same results with ccx 2.11 and ccx 2.8p2
Fusion_Mesh.inp.txt
(178.35 KiB) Downloaded 62 times
User avatar
bernd
Posts: 8485
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: contact constraint example

Postby bernd » Sat Jan 28, 2017 8:18 pm

fandaL wrote:same results with ccx 2.11 and ccx 2.8p2
Fusion_Mesh.inp.txt
It's crazy, if I run the file on linux with 2.11 it breaks after second time step, if I run the file on windows with 2.10 I get the same results as you. I get 10 results and the last one has max deformation 29.76 mm, really strange.

could one of other linux guys test if the file runs ? I'd be curious if it only breaks for me or if it breaks for all linux binaries.

cheers bernd
User avatar
bernd
Posts: 8485
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: contact constraint example

Postby bernd » Sat Jan 28, 2017 8:31 pm

fandaL wrote:
bernd wrote:contact-example-udated.fcstd
With your file I get results. I didn't change any parameter.

OS: Windows 7
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.17.9822 (Git)
Build type: Release
Branch: master
Hash: 0df4af4ba5e6edce825af6ecf6d024c704292633
Python version: 2.7.8
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 7.0.0
same for me on Windows ... and the old ccx 2.8.2 which is included in FreeCAD

OS: Windows 7
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.17.9822 (Git)
Build type: Release
Branch: master
Hash: 0df4af4ba5e6edce825af6ecf6d024c704292633
Python version: 2.7.8
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 7.0.0

Thanks fandaL for testing.
User avatar
makkemal
Posts: 395
Joined: Wed Apr 29, 2015 12:41 pm
Location: South Africa
Contact:

Re: contact constraint example

Postby makkemal » Sun Jan 29, 2017 8:52 am

I don't have a problem this model runs.
ccx 2.10

Code: Select all

OS: Ubuntu 16.04.1 LTS
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.17.9820 +11 (Git)
Build type: Unknown
Branch: Beam3Doutput
Hash: fe6786daa60732912946ee09d347f435a9505ca8
Python version: 2.7.12
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 6.8.0.oce-0.17
contact_example.fcstd
(292.37 KiB) Downloaded 62 times
contact.png
contact.png (329.39 KiB) Viewed 1644 times