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. *

* *

* by Livermore Software Technology Corp. (LSTC) *

* *

******** 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