减输入
*节点卡删除
*KEYWORD $ $=============================CONTROL DEFINITIONS ============================== $ *TITLE Thermal contact between heating plate and parts *CONTROL_SOLUTION $ thermal analysis $ soln nlq isnan lcint lcacc ncdcf 1 *CONTROL_THERMAL_SOLVER $ transient thermal analysis $ atype ptype solver cgtol gpt eqheat fwork sbc 1 1 11 $ msglvl maxitr abstol reltol omega tsf *CONTROL_THERMAL_TIMESTEP $ ts tip its tmin tmax dtemp tscp lcts 1 0.1 *CONTROL_TERMINATION $ endtim endcyc dtmin endeng endmas nosol 10 $ $================================ OUTPUTDATA =================================== $ *DATABASE_BINARY_D3PLOT $ dt 0.1 *DATABASE_GLSTAT $ dt 0.1 *DATABASE_MATSUM $ dt 0.1 *DATABASE_TPRINT $ dt 0.1 $ $============================== PART DEFINITIONS =============================== $ *PART $ title heating_plate $ pid secid mid eosid hgid grav adpopt tmid 1 1 0 0 0 0 0 2 *PART $ title cube_contact $ pid secid mid eosid hgid grav adpopt tmid 2 1 0 0 0 0 0 1 *PART $ title cube_close $ pid secid mid eosid hgid grav adpopt tmid 3 1 0 0 0 0 0 1 *PART $ title cube_far $ pid secid mid eosid hgid grav adpopt tmid 4 1 0 0 0 0 0 1 $ $============================= SECTION PROPERTIES ============================== $ *SECTION_SOLID $ secid elform 1 -1 $ $======================== THERMAL MATERIAL PROPERTIES ========================== $ *MAT_THERMAL_ISOTROPIC $ zinc $ tmid tro tgrlc tgmult tlat hlat 1 7.13e-9 $ hc tc 3.88e8 110.0 *MAT_THERMAL_ISOTROPIC $ steel $ tmid tro tgrlc tgmult tlat hlat 2 7.85e-9 $ hc tc 4.60e8 40.0 $ $=======================NODE/ELEMENT/SET/SEGMENT DEFINTIONS ==================== $ *INCLUDE 06_heating_plate_parts.k $ $*SET_SEGMENT_TITLE segments_plate_bottom sid 1 $ *SET_PART_LIST $ sid 5 $ pid1 pid2 pid3 2 3 4 $ $=========================== CONTACT DEFINITIONS =============================== $ *CONTACT_SURFACE_TO_SURFACE_THERMAL_ID 1,contact $ ssid msid sstyp mstyp sboxid mboxid spr mpr 5 1 2 3 $ fs fd dc vc vdc penchk bt dt $ sfs sfm sst mst sfst sfmt fsf vsf $ k frad h0 lmin lmax ftoslv 0.1 6.3E-12 50.0 0.01 2.0 $ frad based on emissivity e1=e2=0.2 $==================== THERMAL BOUNDARY and INITIAL CONDITIONS ================== $ *INITIAL_TEMPERATURE_SET $ nid/sid temp loc 0 293.15 *BOUNDARY_FLUX_SET $ ssid 1 $ lcid mlc1 mlc2 mlc3 mlc4 loc nhisv 101 -1.0 -1.0 -1.0 -1.0 *DEFINE_CURVE $ lcid sidr sfa sfo offa offo dattyp 101 0 1 1 0 0 0 $ time , heat flux 0.0,0.0 1.0,200.0 11.0,200.0 $ *END