"especially for transformers")
kappa = ith or ip
net["_options"] = {}
_add_ppc_options(net,
calculate_voltage_angles=False,
trafo_model="pi",
check_connectivity=check_connectivity,
mode="sc",
switch_rx_ratio=2,
init_vm_pu="flat",
init_va_degree="flat",
enforce_q_lims=False,
recycle=None)
_add_sc_options(net,
fault=fault,
case=case,
lv_tol_percent=lv_tol_percent,
tk_s=tk_s,
topology=topology,
r_fault_ohm=r_fault_ohm,
kappa_method=kappa_method,
x_fault_ohm=x_fault_ohm,
kappa=kappa,
ip=ip,
ith=ith,
branch_results=branch_results,
return_all_currents=return_all_currents)
init_results(net, "sc")
if fault == "3ph":
_calc_sc(net)
logger.warning("Branch results are in beta mode and might not always be reliable, "
"especially for transformers")
# Convert bus to numpy array
if bus is None:
bus = net.bus.index.values
else:
bus = np.array([bus]).ravel()
kappa = ith or ip
net["_options"] = {}
_add_ppc_options(net, calculate_voltage_angles=False, trafo_model="pi",
check_connectivity=check_connectivity, mode="sc", switch_rx_ratio=2,
init_vm_pu="flat", init_va_degree="flat", enforce_q_lims=False,
recycle=None)
_add_sc_options(net, fault=fault, case=case, lv_tol_percent=lv_tol_percent, tk_s=tk_s,
topology=topology, r_fault_ohm=r_fault_ohm, kappa_method=kappa_method,
x_fault_ohm=x_fault_ohm, kappa=kappa, ip=ip, ith=ith,
branch_results=branch_results, return_all_currents=return_all_currents,
inverse_y=inverse_y)
init_results(net, "sc")
if fault in ("2ph", "3ph"):
_calc_sc(net, bus)
elif fault == "1ph":
_calc_sc_1ph(net, bus)
else:
raise ValueError("Invalid fault %s" % fault)
def _calc_current(net, ppci_orig, bus):
