dyna在计算sph模拟的流固耦合问题时,总是中止计算,是什么原因?
Input file: G:\sph\1.k
The native file format : 32-bit small endian
Memory size from command line: 100000000
Binary file : 32 bit IEEE selected
expanding memory to 100000000 starting ls-dyna
******** notice ******** notice ******** notice ********
* *
* This is the LS-DYNA Finite Element code. *
* *
* Neither LSTC nor the authors assume any responsibility for *
* the validity, accuracy, or applicability of any results *
* obtained from this system. The user must verify his own *
* results. *
* *
* LSTC endeavors to make the LS-DYNA code as complete, *
* accurate and easy to use as possible. *
* Suggestions and comments are welcomed. Please report any *
* errors encountered in either the documentation or results *
* immediately to LSTC through your site focus. *
* *
* Copyright (C) 1990-2009 *
* by Livermore Software Technology Corp. (LSTC) *
* All rights reserved *
* *
******** notice ******** notice ******** notice ********
contracting memory to 1 beginning of keyword reader
expanding memory to 21002
expanding memory to 126706
expanding memory to 171706
expanding memory to 201706
expanding memory to 406402
expanding memory to 428084
expanding memory to 480316
expanding memory to 532548
expanding memory to 553301
*** Warning extra card found in *CONTACT read
The following card will be ignored.
Each Contact definition must start with a
*CONTACT_... card
*** Warning extra card found in *CONTACT read
The following card will be ignored.
Each Contact definition must start with a
*CONTACT_... card
*** Warning extra card found in *CONTACT read
The following card will be ignored.
Each Contact definition must start with a
*CONTACT_... card
expanding memory to 574401
expanding memory to 594401
expanding memory to 614401
expanding memory to 634701
expanding memory to 654701
expanding memory to 681219
expanding memory to 879435
expanding memory to 903551
expanding memory to 956243
expanding memory to 978741
expanding memory to 1000741
expanding memory to 1046418
Memory required to process keyword : 1046418
contracting memory to 1 end of keyword reader
LS-DYNA will perform a structural only analysis
ls-dyna ls971s.50638 R4.2 date 02/04/2009
expanding memory to 2833 material, node and element data
expanding memory to 2920 material properties
expanding memory to 91801 element data
expanding memory to 249626 gradient of damage data
expanding memory to 707968 user control for contact interfaces
expanding memory to 707973 for boundary cond. vectorization
expanding memory to 720006 arbitrary element numbering input
*** Warning warpage angle of element 2136129 is 2.9010E+01 degrees
*** Warning warpage angle of element 2136130 is 2.9004E+01 degrees
*** Warning warpage angle of element 2136131 is 2.9010E+01 degrees
*** Warning warpage angle of element 2136132 is 2.9031E+01 degrees
*** Warning warpage angle of element 2136148 is 2.9010E+01 degrees
*** Warning warpage angle of element 2136149 is 2.9010E+01 degrees
*** Warning warpage angle of element 2136150 is 2.9010E+01 degrees
*** Warning warpage angle of element 2136151 is 2.9010E+01 degrees
*** Warning warpage angle of element 2140129 is 2.8989E+01 degrees
*** Warning warpage angle of element 2140130 is 2.9004E+01 degrees
*** Warning warpage angle of element 2140131 is 2.9007E+01 degrees
*** Warning warpage angle of element 2140132 is 2.9009E+01 degrees
*** Warning warpage angle of element 2140148 is 2.9010E+01 degrees
*** Warning warpage angle of element 2140149 is 2.9010E+01 degrees
*** Warning warpage angle of element 2140150 is 2.9010E+01 degrees
*** Warning warpage angle of element 2140151 is 2.9010E+01 degrees
*** Warning Number of elements exceed warpage angle = 16
min warpage angle = 0.2898871E+02
max warpage angle = 0.2903063E+02
expanding memory to 787246 spc cards
expanding memory to 1181451 process element material data
expanding memory to 1183457 load curve data
contracting memory to 1182077 regularize time history data
contracting memory to 1182077 rigidwall definitions
expanding memory to 1215114 nodal constraint sets
expanding memory to 1216333 contact interfaces
expanding memory to 1250369 implicit contact 1
expanding memory to 1250371 implicit contact 2
expanding memory to 1525043 contact segment data 1
expanding memory to 1642991 contact segment data 2
expanding memory to 1763113 contact segment data 3
expanding memory to 2502086 contact definition
expanding memory to 2552084 rigid body
expanding memory to 2931438 contact friction
expanding memory to 2945388 rigid body b.c.
expanding memory to 3323694 joints
contracting memory to 3323694 shell solid interfaces
expanding memory to 3835753 SPH elements
expanding memory to 3847638 discrete element
contracting memory to 3847638 accelerometers
expanding memory to 3847738 rigid bodies for disp. termination
contracting memory to 3847738 transformation to move nodes
expanding memory to 3847746 contact data
expanding memory to 3850008 contact data
expanding memory to 3850012 thermal expansion 0
expanding memory to 3850016 flexible body 2
contracting memory to 3850016 implicit joints
expanding memory to 3898248 initialization phase
slave surface of interface = 2000153
type = a 3
surface timestep = 6.573E-05
controlling slave node ID = 6081274
part ID of slave node = 2000804
current minimum time step = 6.573E-05
master surface of interface = 2000153
type = a 3
surface timestep = 2.838E-05
controlling master node ID = 4320934
part ID of master node = 2000805
current minimum time step = 2.838E-05
slave surface of interface = 2000145
type = a 3
surface timestep = 1.000E+16
controlling slave node ID = 12024
part ID of slave node = 2327507
current minimum time step = 2.838E-05
master surface of interface = 2000145
type = a 3
surface timestep = 1.000E+16
controlling master node ID = 12024
part ID of master node = 2327507
current minimum time step = 2.838E-05
slave surface of interface = 2000154
type = a 5
surface timestep = 1.251E+02
controlling slave node ID = 6364249
part ID of slave node = 2327507
current minimum time step = 2.838E-05
master surface of interface = 2000154
type = a 5
surface timestep = 2.838E-05
controlling master node ID = 4320934
part ID of master node = 2000805
current minimum time step = 2.838E-05
slave surface of interface = 2000139
type = 13
surface timestep = 2.838E-05
controlling slave node ID = 4320934
part ID of slave node = 2000805
current minimum time step = 2.838E-05
The LS-DYNA time step size should not exceed 2.838E-05
to avoid contact instabilities. If the step size is
bigger then scale the penalty of the offending surface.
s t o r a g e a l l o c a t i o n
Memory required to begin solution : 3898248
Additional dynamically allocated memory: 1081800
Total: 4980048