非叠加焊接d3plot例子基础是用固态元素焊接模拟例中非相联分析表示温度字段先用热分析计算,然后用单机结构分析绘制节点以计算变形热分析binout文件引用为温度字段源//www.cannycloudtech.com/thermal/welding-new/uncoupled-d3plot//www.cannycloudtech.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
非叠加焊接d3plot
例子基础是用固态元素焊接模拟例中非相联分析表示温度字段先用热分析计算,然后用单机结构分析绘制节点以计算变形热分析binout文件引用为温度字段源
*BOUNDARY_CONVECTION_SET *BOUNDARY_RADIATION_SET *BOUNDARY_SPC_SET *BOUNDARY_THERMAL_WELD_TRAJECTORY *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_THERMAL_ID *CONTACT_TIED_SURFACE_TO_SURFACE_ID *CONTACT_TIED_SURFACE_TO_SURFACE_THERMAL_ID *CONTROL_ACCURACY *CONTROL_IMPLICIT_AUTO *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_SOLUTION *CONTROL_SOLUTION *CONTROL_TERMINATION *CONTROL_THERMAL_SOLVER *CONTROL_THERMAL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_GLSTAT *DATABASE_MATSUM *DATABASE_RCFORC *DATABASE_SPCFORC *DEFINE_CURVE *DEFINE_CURVE_TITLE *ELEMENT_SOLID *END *INCLUDE *INITIAL_TEMPERATURE_SET *KEYWORD_JOBID *LOAD_BODY_Z *LOAD_THERMAL_D3PLOT *MAT_CWM_TITLE *MAT_ELASTIC_PLASTIC_THERMAL_TITLE *MAT_THERMAL_CWM_TITLE *MAT_THERMAL_ISOTROPIC *NODE *PART *SECTION_SOLID *SET_NODE_LIST_TITLE *SET_SEGMENT_GENERAL *SET_SEGMENT_TITLE *TITLE
开工热分析
*KEYWORD_JOBID thermalstep $ $=============================CONTROL DEFINITIONS ============================== $ *TITLE welding steel parts - uncoupled thermal step *CONTROL_SOLUTION $ soln 1 *CONTROL_TERMINATION $ endtim 99 *CONTROL_ACCURACY $ osu inn pidosu iacc 1 *CONTROL_THERMAL_SOLVER $ atype ptype solver 1 1 11 $ msglvl maxitr 0 50 *CONTROL_THERMAL_TIMESTEP $ ts tip its tmin tmax dtemp tscp lcts 1 1.00 0.2 -902 -901 5000 0.50 *DEFINE_CURVE_TITLE minimal thermal time step size $ lcid sidr sfa sfo offa offo dattyp 902 0 1.0 1.0 0.000 0.000 0 $ a1 o1 0.0 0.1 7.0 0.1 10.0 2.0 100.0 4.0 *DEFINE_CURVE_TITLE maximal thermal step size $ lcid sidr sfa sfo offa offo dattyp 901 0 1.0 1.0 0.000 0.000 0 $ a1 o1 0.0 0.25 7.0 0.25 10.0 4.0 100.0 10.0 $================================ OUTPUTDATA =================================== *DATABASE_BINARY_D3PLOT $ dt 0.01 *DATABASE_GLSTAT $ dt 0.01 *DATABASE_MATSUM $ dt 0.01 *DATABASE_RCFORC $ dt 0.01 *DATABASE_SPCFORC $ dt 0.01 $ $============================== PART DEFINITIONS =============================== $ *PART $ title profile $ pid secid mid eosid hgid grav adpopt tmid 1 1 1 *PART $ title plate $ pid secid mid eosid hgid grav adpopt tmid 2 1 1 *PART $ title weldingseam 3 1 2 $ $============================= SECTION PROPERTIES ============================== $ *SECTION_SOLID $ secid elform 1 -2 $ $================== MECHANICAL AND THERMAL MATERIAL PROPERTIES ================= $ *MAT_THERMAL_ISOTROPIC $ tmid tro tgrlc tgmult tlat hlat 1 7.85e-9 $ hc tc 4.60e8 40.0 *MAT_THERMAL_CWM_TITLE weld metal $ tmid tro tgrlc tgmult hdeath tdeath 2 7.85e-9 0 0.0 511e6 0.001 $ lchc lctc tlstart tlend tistart tiend hghost tghost 106 107 -1001 -1000 -2.0 -1.0 511e6 0.001 *DEFINE_CURVE_TITLE lchc temp dependent specific heat $ lcid sidr sfa sfo offa offo dattyp 106 0 1.0 1e6 0.000 0.000 0 $ a1 o1 273.0 500 473.0 550 10000 750 *DEFINE_CURVE_TITLE lctc temp dependent thermal conductivity $ lcid sidr sfa sfo offa offo dattyp 107 0 1.0 1.0 0.000 0.000 0 $ a1 o1 273.0 15.0 473.0 20.0 10000.0 30.0 $ $=======================NODE/ELEMENT/SET/SEGMENT DEFINTIONS ==================== $ *INCLUDE 05_welding_solid_parts.k $ $ *SET_NODE_LIST_TITLE nodes_weldingseam_path sid 1 $ *SET_NODE_LIST_TITLE nodes_plate_edge_side1 sid 2 $ *SET_NODE_LIST_TITLE nodes_plate_edge_side2 sid 3 $ $ *SET_SEGMENT_TITLE segments_weldingseam_path sid 1 $ *SET_SEGMENT_TITLE segments_rad_and_conv sid 2 $ *SET_SEGMENT_TITLE segments_weldingseam_contact sid 3 $ *SET_SEGMENT_GENERAL $ sid 4 $ option e1 e2 part 1 2 $ $=========================== CONTACT DEFINITIONS =============================== $ *CONTACT_TIED_SURFACE_TO_SURFACE_THERMAL_ID 1,contact_tied $ ssid msid sstyp mstyp sboxid mboxid spr mpr 3 4 0 0 $ fs fd dc vc vdc penchk bt dt $ sfs sfm sst mst sfst sfmt fsf vsf -0.1 -0.1 $ k frad ho lmin lmax ftoslv bc_flg algo 0.026 0.0 50.0 0.5 2.0 0.5 1 0 $ card a $ card b $ igap 2 *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_THERMAL_ID 2,contact_automatic $ ssid msid sstyp mstyp sboxid mboxid spr mpr 1 2 3 3 $ fs fd dc vc vdc penchk bt dt $ sfs sfm sst mst sfst sfmt fsf vsf $ k frad ho lmin lmax ftoslv bc_flg algo 0.026 0.0 1.0 0.5 2 0.5 1 0 $ $================= THERMAL BOUNDARY AND INITIAL CONDITIONS ===================== $ *INITIAL_TEMPERATURE_SET $ nsid temp loc 0 293.15 *BOUNDARY_THERMAL_WELD_TRAJECTORY $ pid ptyp nsid1 spd1 nsid2 spd2 ncyc relvel 3 1 1 -301 -1 -301 5 $ iform lcid q lcrot lcmov lclat disc 2 302 5.0e6 $ p1 p2 p3 p4 p5 p6 p7 p8 3.0 2.5 2.0 8.0 0.5 1.0 *DEFINE_CURVE $ lcid sidr sfa sfo offa offo dattyp 301 0 1.0 1.0 0.000 0.000 0 $ time,temp 0.0,0.0 0.1,6.0 6.9,6.0 7.0,0.0 101.0,0.0 *DEFINE_CURVE $ lcid sidr sfa sfo offa offo dattyp 302 0 1.0 1.0 0.000 0.000 0 $ time,temp 0.0,0.0 0.1,1.0 6.9,1.0 7.0,0.0 101.0,0.0 *BOUNDARY_CONVECTION_SET $ ssid 2 $ hlcid hmult tlcid tmult loc 0 0.01 0 293.15 *BOUNDARY_RADIATION_SET $ ssid type 2 1 $ flcid fmult tlcid tmult loc 03.4022E-11 0 293.15 $ *END
二叉结构分析
*KEYWORD_ID structuralstep $ $=============================CONTROL DEFINITIONS ============================== $ *TITLE welding steel parts - solids *CONTROL_SOLUTION $ soln 0 *CONTROL_TERMINATION $ endtim 98 *CONTROL_IMPLICIT_GENERAL $ imflag dt0 imform nsbs igs cnstn form zero_v 1 0.2 2 1 2 0 0 1 *CONTROL_IMPLICIT_SOLUTION $ nsolv 12 *CONTROL_IMPLICIT_AUTO $ iauto iteopt itewin dtmin dtmax 1 11 5 0.01 -901 *CONTROL_ACCURACY $ osu inn pidosu iacc 1 *DEFINE_CURVE_TITLE maximal structural step size $ lcid sidr sfa sfo offa offo dattyp 901 0 1.0 1.0 0.000 0.000 0 $ a1 o1 0.0 0.25 7.0 0.25 10.0 4.0 100.0 10.0 $ $================================ OUTPUTDATA =================================== $ *DATABASE_BINARY_D3PLOT $ dt 0.01 *DATABASE_GLSTAT $ dt 0.01 *DATABASE_MATSUM $ dt 0.01 *DATABASE_RCFORC $ dt 0.01 *DATABASE_SPCFORC $ dt 0.01 $ $============================== PART DEFINITIONS =============================== $ *PART $ title profile $ pid secid mid eosid hgid grav adpopt tmid 1 1 1 *PART $ title plate $ pid secid mid eosid hgid grav adpopt tmid 2 1 1 *PART $ title weldingseam 3 1 2 $ $============================= SECTION PROPERTIES ============================== $ *SECTION_SOLID $ secid elform 1 -2 $ $================== MECHANICAL AND THERMAL MATERIAL PROPERTIES ================= $ *MAT_ELASTIC_PLASTIC_THERMAL_TITLE steel $ mid ro 1 7.85e-9 $ t1 t2 t3 t4 273.0 493.0 1273.0 10000.0 $ e1 e2 e3 e4 210000 210000 75000 1000 $ pr1 pr2 pr3 pr4 0.285 0.285 0.3 0.45 $ alpha1 alpha2 alpha3 alpha4 12e-6 12e-6 14e-6 0.0 $ sigy1 sigy2 sigy3 sigy4 435.0 100.0 20 1 $ etan1 etan2 etan3 etan4 1000.0 1000.0 250.0 1 *MAT_CWM_TITLE weld metal $ mid ro lcem lcpr lcsy lchr lcat beta 2 7.85e-09 101 102 103 104 105 1.0 $ tastart taebd tlstart tlend eghost pghost aghost 1300 1400 1200 1400 10000.0 0.49 0.0 $ t2phase t1phase 800 500 *DEFINE_CURVE_TITLE lecm temp dependent youngs modulus $ lcid sidr sfa sfo offa offo dattyp 101 0 1.0 1.0 0.000 0.000 0 $ a1 o1 273.0 210000 493.0 210000 10000.0 1000 *DEFINE_CURVE_TITLE lcpr temp dependent poissons ratio $ lcid sidr sfa sfo offa offo dattyp 102 0 1.0 1.0 0.000 0.000 0 $ a1 o1 273.0 0.30 473.0 0.30 10000.0 0.49 *DEFINE_CURVE_TITLE lcsy temp dependent yield stress $ lcid sidr sfa sfo offa offo dattyp 103 0 1.0 1.0 0.000 0.000 0 $ a1 o1 273.0 240 473.0 240 10000.0 5 *DEFINE_CURVE_TITLE lchr temp dependent hardening modulus $ lcid sidr sfa sfo offa offo dattyp 104 0 1.0 1.0 0.000 0.000 0 $ a1 o1 273.0 700 473.0 700 10000.0 5 *DEFINE_CURVE_TITLE lcat temp dependent thermal expansion coefficient $ lcid sidr sfa sfo offa offo dattyp 105 0 1.0 1.0e-06 0.000 0.000 0 $ a1 o1 273.0 17.00 473.0 17.00 1000.0 22.00 10000.0 0.00 $ $=======================NODE/ELEMENT/SET/SEGMENT DEFINTIONS ==================== $ *INCLUDE 05_welding_solid_parts.k $ $ *SET_NODE_LIST_TITLE nodes_weldingseam_path sid 1 $ *SET_NODE_LIST_TITLE nodes_plate_edge_side1 sid 2 $ *SET_NODE_LIST_TITLE nodes_plate_edge_side2 sid 3 $ $ *SET_SEGMENT_TITLE segments_weldingseam_path sid 1 $ *SET_SEGMENT_TITLE segments_rad_and_conv sid 2 $ *SET_SEGMENT_TITLE segments_weldingseam_contact sid 3 $ *SET_SEGMENT_GENERAL $ sid 4 $ option e1 e2 part 1 2 $ $=========================== CONTACT DEFINITIONS =============================== $ *CONTACT_TIED_SURFACE_TO_SURFACE_ID 1,contact_tied $ ssid msid sstyp mstyp sboxid mboxid spr mpr 3 4 0 0 $ fs fd dc vc vdc penchk bt dt $ sfs sfm sst mst sfst sfmt fsf vsf -0.1 -0.1 $ card a $ card b $ igap 2 *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_ID 2,contact_automatic $ ssid msid sstyp mstyp sboxid mboxid spr mpr 1 2 3 3 $ fs fd dc vc vdc penchk bt dt 0.3 $ sfs sfm sst mst sfst sfmt fsf vsf $ $====================== MECHANICAL BOUNDARY CONDITIONS ========================= $ *BOUNDARY_SPC_SET $ nid/nsid cid dofx dofy dofz dofrx dofry dofrz 2 1 1 1 *BOUNDARY_SPC_SET $ nid/nsid cid dofx dofy dofz dofrx dofry dofrz 3 1 1 1 *LOAD_BODY_Z $ lcid sf lciddr xc yc zc cid 201 9.81e3 *DEFINE_CURVE $ lcid 201 $ time,temp 0,1.0 101.0,1.0 *LOAD_THERMAL_D3PLOT $ load the previously calculated temperature fields from the file which is $ referenced in the execution expression with "T=thermalstep.d3plot" $ *END
开工热分析 |
二叉结构分析 |
|
|
开工热分析 |
二叉结构分析 |
|
|
