  The detection of GW170817/GRB170817A implied the strong association between
short gamma-ray bursts (SGRBs) and binary neutron star (BNS) mergers which
produce gravitational waves (GWs). More evidence is needed to confirm the
association and reveal the physical processes of BNS mergers. The upcoming High
Altitude Detection of Astronomical Radiation (HADAR) experiment, excelling in a
wide field of view (FOV) and a large effective area above tens of GeV, is a
hope for the prompt detection of very-high-energy (VHE; > 10 GeV) SGRBs. The
aim of this paper is to simulate and analyse GW/SGRB joint detections by future
GW detector networks in synergy with HADAR, including the second generation
LIGO, Virgo and KAGRA and the third generation ET and CE. We provide a brief
introduction of the HADAR experiment for SGRB simulations and its expected SGRB
detections. For GW simulations, we adopt a phenomenological model to describe
GWs produced by BNS mergers and introduce the signal-noise ratios (SNRs) as
detector responses. Following a theoretical analysis we compute the
redshift-dependent efficiency functions of GW detector networks. We then
construct the simulation of GW detection by Monte Carlo sampling. We compare
the simulated results of LIGO-Virgo O2 and O3 runs with their actual detections
as a check. The combination of GW and SGRB models is then discussed for joint
detection, including parameter correlations, triggered SNRs and efficiency
skymaps. The estimated joint detection rates are 0.09-2.52 per year for LHVK
network with HADAR under different possible configurations, and approximately
0.27-7.89 per year for ET+CE network with HADAR.
