如题,导出K文件并进行求解的时候出现了如上所示的错误,但是在之前用的另一个模型上使用相似的设置进行求解并没有问题发生。定义的curve为x从0 到0.5,y为恒定值3.08的一段直线。
部分K文件内容如下:
*CONTROL_CONTACT $# slsfac rwpnal islchk shlthk penopt thkchg orien enmass 2.000000 0.000 1 0 0 0 1 0 $# usrstr usrfrc nsbcs interm xpene ssthk ecdt tiedprj 0 0 0 0 4.000000 0 0 0 $# sfric dfric edc vfc th th_sf pen_sf 0.000 0.000 0.000 0.000 0.000 0.000 0.000 $# ignore frceng skiprwg outseg spotstp spotdel spothin 0 0 0 0 0 0 0.000 $# isym nserod rwgaps rwgdth rwksf icov swradf ithoff 0 0 0 0.000 1.000000 0 0.000 0 $# shledg pstiff ithcnt tdcnof ftall unused shltrw 0 0 0 0 0 0 0.000 *CONTROL_ENERGY $# hgen rwen slnten rylen 2 2 2 2 *CONTROL_SHELL $# wrpang esort irnxx istupd theory bwc miter proj 20.000000 1 -1 1 2 2 1 0 $# rotascl intgrd lamsht cstyp6 tshell nfail1 nfail4 psnfail 1.000000 0 0 1 0 0 0 0 $# psstupd irquad cntco itsflg irquad 0 0 0 0 2 *CONTROL_TERMINATION $# endtim endcyc dtmin endeng endmas 0.500000 0 0.000 0.000 0.000 *CONTROL_TIMESTEP $# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st 0.000 0.900000 0 0.000 0.000 0 0 0 $# dt2msf dt2mslc imscl 0.000 0 0 *DATABASE_BINARY_D3PLOT $# dt lcdt beam npltc psetid 0.005000 0 0 0 0 $# ioopt 0 *DATABASE_BINARY_D3THDT $# dt lcdt beam npltc psetid 5.0000E-4 0 0 0 0 *DATABASE_FORMAT $# iform ibinary 2 0 *DATABASE_RCFORC_MOMENT $# cid nodes nodem 1 0 0 *DATABASE_EXTENT_BINARY $# neiph neips maxint strflg sigflg epsflg rltflg engflg 0 0 3 1 1 1 1 1 $# cmpflg ieverp beamip dcomp shge stssz n3thdt ialemat 0 0 4 1 1 1 2 0 $# nintsld pkp_sen sclp unused msscl therm intout nodout 0 0 1.000000 0 0 0STRESS STRESS *BOUNDARY_PRESCRIBED_MOTION_RIGID $# pid dof vad lcid sf vid death birth 2 1 0 2 1.000000 01.0000E+28 0.000 *CONTACT_ERODING_SURFACE_TO_SURFACE_ID $# cid title 1 $# ssid msid sstyp mstyp sboxid mboxid spr mpr 1 2 3 3 0 0 0 0 $# fs fd dc vc vdc penchk bt dt 0.300000 0.250000 0.000 0.000 0.000 0 0.000 1.0000E+7 $# sfs sfm sst mst sfst sfmt fsf vsf 1.000000 1.000000 0.000 0.000 1.000000 1.000000 1.000000 1.000000 $# isym erosop iadj 0 1 0 *PART $# title Part 1 for Mat 2 and Elem Type 1 $# pid secid mid eosid hgid grav adpopt tmid 1 1 2 0 0 0 0 0 *SECTION_SOLID $# secid elform aet 1 1 0 *MAT_PLASTIC_KINEMATIC $# mid ro e pr sigy etan beta 2 0.007800 2.0700E+5 0.250000 0.835000 80900.000 1.000000 $# src srp fs vp 0.000 0.000 0.800000 0.000 *PART $# title Part 2 for Mat 1 and Elem Type 2 $# pid secid mid eosid hgid grav adpopt tmid 2 2 1 0 0 0 0 0 *SECTION_SOLID $# secid elform aet 2 16 0 *MAT_RIGID $# mid ro e pr n couple m alias 1 0.003500 9.6000E+5 0.200000 0.000 0.000 0.000 $# cmo con1 con2 1.000000 5 7 $# lco or a1 a2 a3 v1 v2 v3 0.000 0.000 0.000 0.000 0.000 0.000 *INITIAL_VELOCITY_RIGID_BODY $# pid vx vy vz vxr vyr vzr icid 2 3.080000 0.000 0.000 0.000 0.000 0.000 0 *DEFINE_CURVE $# lcid sidr sfa sfo offa offo dattyp 1 0 1.000000 1.000000 0.000 0.000 0 $# a1 o1 0.000 3.080000 0.500000 3.080000
|