Dislotwin_TWIP-TRIP Parameters, Martensite and Twinning VFs

[yigitnd]

   I have several questions:
   
   1) In many papers, p_sl and q_sl values are 1.15 and 1.0 respectively. However, DAMASK does not allow me to do so because these values are out of the range of those parameters. If the model works with p=1.15 and q_sl=1.0, why is there such limitation and how can I make DAMASK work with p_sl=1.15 and q=1.0?
   
   2) Dislotwin model parameters used in my simulation at 323K is shown below. Model creates twin VF and martensite VF around 2-3% at 35-40% strain, which is too low compared with the experiments and simulations in the papers. I tried many different parameter sets (including changing interaction matrices as well) but unable to solve this problem. What may be the cause?
   
    N_sl: [12]
    f_edge: [1.0]
    b_sl: [2.56e-10]
    Q_sl: [3.5e-19]
    p_sl: [1]
    q_sl: [1.15]
    B: [0.001]
    i_sl: [30.0]
    v_0: [0.0001]
    tau_0: [123000000.0]
    D_a: 2.0
    Q_cl: 3.0e-19
    rho_mob_0: [50000000000.0]
    rho_dip_0: [50000000000.0]
    h_sl-sl: [0.122, 0.122, 0.625, 0.07, 0.137, 0.137, 0.122]
    N_tw: [12]
    b_tw: [1.2e-10]
    L_tw: 9e-08
    i_tw: 2.0
    t_tw: [5.0e-08]
    p_tw: [2]
    h_tw-tw: [1.0, 1.0]
    h_sl-tw: [1.0, 1.0, 1.0]
    N_tr: [12]
    b_tr: [1.2e-10]
    L_tr: 9e-08
    i_tr: 2.0
    t_tr: [1.0e-07]
    p_tr: [15]
    V_mol: 7.09e-06
    c/a_hP: 1.633
    Delta_G: 177.69
    h_tr-tr: [1.0, 1.0]
    h_sl-tr: [1.5, 1.5, 1.5]
    T_ref: 323
    Gamma_sf: 0.02645
    x_c: 1.0e-09
    V_cs: 1.67e-29
    D: 3.0e-06
    
    
    3) The parameter set above is used at different temperatures by changing temperature dependent parameters like Gamma_sf, Delta_G, T_ref and  tau_0 with the calculated values at those temperatures in the papers. Martensite VF is expected to decrease with increasing temperature (due to high SFE), but the opposite happens. I am sharing the results. Can you make comments on this situation?
    
    Thank you in advance and appreciate your time.